Biomedical subjects

Steen E Pedersen

Publications and source records attributed to Steen E Pedersen.

8 recordsLinked to original sources

PubMed2006

Kinetic modeling and fitting software for interconnected reaction schemes: VisKin.

Reaction kinetics for complex, highly interconnected kinetic schemes are modeled using analytical solutions to a system of ordinary differential equations. The algorithm employs standard linear algebra methods that are implemented using MatLab functions in a Visual Basic interface. A graphical user interface for simple entry of reaction schemes facilitates comparison of a variety of reaction schemes. To ensure microscopic balance, graph theory algorithms are used to determine violations of thermodynamic cycle constraints. Analytical solutions based on linear differential equations result in fast comparisons of first order kinetic rates and amplitudes as a function of changing ligand concentrations. For analysis of higher order kinetics, we also implemented a solution using numerical integration. To determine rate constants from experimental data, fitting algorithms that adjust rate constants to fit the model to imported data were implemented using the Levenberg-Marquardt algorithm or using Broyden-Fletcher-Goldfarb-Shanno methods. We have included the ability to carry out global fitting of data sets obtained at varying ligand concentrations. These tools are combined in a single package, which we have dubbed VisKin, to guide and analyze kinetic experiments. The software is available online for use on PCs.

Computer Simulation↗

PubMed2006

Nicotinic acetylcholine receptor channel electrostatics determined by diffusion-enhanced luminescence energy transfer.

The electrostatic potentials within the pore of the nicotinic acetylcholine receptor (nAChR) were determined using lanthanide-based diffusion-enhanced fluorescence energy transfer experiments. Freely diffusing Tb3+ -chelates of varying charge constituted a set of energy transfer donors to the acceptor, crystal violet, a noncompetitive antagonist of the nAChR. Energy transfer from a neutral Tb3+ -chelate to nAChR-bound crystal violet was reduced 95% relative to the energy transfer to free crystal violet. This result indicated that crystal violet was strongly shielded from solvent when bound to the nAChR. Comparison of energy transfer from positively and negatively charged chelates indicate negative electrostatic potentials of -25 mV in the channel, measured in low ionic strength, and -10 mV measured in physiological ionic strength. Debye-Hückel analyses of potentials determined at various ionic strengths were consistent with 1-2 negative charges within 8 A of the crystal violet binding site. To complement the energy transfer experiments, the influence of pH and ionic strength on the binding of [3H]phencyclidine were determined. The ionic strength dependence of binding affinity was consistent with -3.3 charges within 8 A of the binding site, according to Debye-Hückel analysis. The pH dependence of binding had an apparent pKa of 7.2, a value indicative of a potential near -170 mV if the titratable residues are constituted of aspartates and glutamates. It is concluded that long-range potentials are small and likely contribute little to selectivity or conductance whereas close interactions are more likely to contribute to electrostatic stabilization of ions and binding of noncompetitive antagonists within the channel.

Binding Sites↗

PubMed2006

Computed pore potentials of the nicotinic acetylcholine receptor.

Electrostatic surface potentials in the vestibule of the nicotinic acetylcholine receptor (nAChR) were computed from structural models using the University of Houston Brownian Dynamics program to determine their effect on ion conduction and ionic selectivity. To further determine whether computed potentials accurately reflect the electrostatic environment of the channel, the potentials were used to predict the rate constants for diffusion-enhanced fluorescence energy transfer; the calculated energy transfer rates are directly comparable with those determined experimentally (see companion article by Meltzer et al. in this issue). To include any effects on the local potentials by the bound acceptor fluorophore crystal violet, its binding site was first localized within the pore by fluorescence energy transfer measurements from dansyl-C6-choline bound to the agonist sites and also by simulations of binding using Autodock. To compare the computed potentials with those determined experimentally, we used the predicted energy transfer rates from Tb3+ chelates of varying charge to calculate an expected potential using the Boltzmann relationship. This expected potential (from -20 to -40 mV) overestimates the values determined experimentally (from -10 to -25 mV) by two- to fourfold at similar conditions of ionic strength. Although the results indicate a basic discrepancy between experimental and computed surface potentials, both methods demonstrate that the vestibular potential has a relatively small effect on conduction and selectivity.

Cell Membrane↗

PubMed2006

Electrostatic steering at acetylcholine binding sites.

The electrostatic environments near the acetylcholine binding sites on the nicotinic acetylcholine receptor (nAChR) and acetylcholinesterase were measured by diffusion-enhanced fluorescence energy transfer (DEFET) to determine the influence of long-range electrostatic interactions on ligand binding kinetics and net binding energy. Changes in DEFET from variously charged Tb3+ -chelates revealed net potentials of -20 mV at the nAChR agonist sites and -14 mV at the entrance to the AChE active site, in physiological ionic strength conditions. The potential at the alphadelta-binding site of the nAChR was determined independently in the presence of d-tubocurarine to be -14 mV; the calculated potential at the alphagamma-site was approximately threefold stronger than at the alphadelta-site. By determining the local potential in increasing ionic strength, Debye-Hückel theory predicted that the potentials near the nAChR agonist binding sites are constituted by one to three charges in close proximity to the binding site. Examination of the binding kinetics of the fluorescent acetylcholine analog dansyl-C6-choline at ionic strengths from 12.5 to 400 mM revealed a twofold decrease in association rate. Debye-Hückel analysis of the kinetics revealed a similar charge distribution as seen by changes in the potentials. To determine whether the experimentally determined potentials are reflected by continuum electrostatics calculations, solutions to the nonlinear Poisson-Boltzmann equation were used to compute the potentials expected from DEFET measurements from high-resolution models of the nAChR and AChE. These calculations are in good agreement with the DEFET measurements for AChE and for the alphagamma-site of the nAChR. We conclude that long-range electrostatic interactions contribute -0.3 and -1 kcal/mol to the binding energy at the nAChR alphadelta- and alphagamma-sites due to an increase in association rates.

Acetylcholine↗

PubMed2006

Site specificity of agonist-induced opening and desensitization of the Torpedo californica nicotinic acetylcholine receptor.

Agonist-binding kinetics to the nicotinic acetylcholine receptor (AChR) from Torpedo californica were measured using sequential-mixing stopped-flow fluorescence methods to determine the contribution of each individual site to agonist-induced opening and desensitization. Timed dansyl-C6-choline (DC6C) binding followed by its dissociation upon mixing with high, competing agonist concentrations revealed four kinetic components: an initial, fast fluorescence decay, followed by a transient increase, and then two characteristic decays that reflect dissociation from the desensitized agonist sites. The transient increase resulted from DC6C binding to the open-channel based on its prevention by proadifen, a noncompetitive antagonist. Further characterization of DC6C channel binding by the inhibition of [3H]phencyclidine binding and by equilibrium measurements of DC6C fluorescence yielded KD values of 2-4 microM for the desensitized AChR and approximately 600 microM for the closed state. At this site, DC6C displayed a strongly blue-shifted emission spectrum, higher intrinsic fluorescence, and weaker energy transfer from tryptophans than when bound to either agonist site. The initial, fast fluorescence decay was assigned to DC6C dissociation from the alphadelta site of the AChR in its closed conformation, on the basis of inhibition with the site-selective antagonists d-tubocurarine and alpha-conotoxin MI. Fast decay amplitude data indicated an apparent affinity of 0.9 microM for the closed-state alphadelta site; the closed-state alphagamma-site affinity is inferred to be near 100 microM. These values and the known affinities for the desensitized conformation show that the alphagamma site drives AChR desensitization to a approximately 40-fold greater extent than the alphadelta site, undergoes energetically larger conformational changes, and is the primary determinant of agonist potency.

Animals↗

PubMed2004

Calmodulin interactions with IQ peptides from voltage-dependent calcium channels.

Calmodulin (CaM) functions as a Ca(2+) sensor for inactivation and, in some cases, facilitation of a variety of voltage-dependent Ca(2+) channels. A crucial determinant for CaM binding to these channels is the IQ motif in the COOH-terminal tail of the channel-forming subunit. The binding of CaM to IQ peptides from Lc-, P/Q-, and R-type, but not N-type, voltage-dependent Ca(2+) channels increases the Ca(2+) affinity of both lobes of CaM, producing similar N- and C-lobe Ca(2+) affinities. Ca(2+) associates with and dissociates from the N-lobe much more rapidly than the C-lobe when CaM is bound to the IQ peptides. Compared with the other IQ peptides, CaM-bound Lc-IQ has the highest Ca(2+) affinity and the most rapid rates of Ca(2+) association at both lobes, which is likely to make Ca(2+) binding to CaM, bound to this channel, less sensitive than other channels to intracellular Ca(2+) buffers. These kinetic differences in Ca(2+) binding to the lobes of CaM when bound to the different IQ motifs may explain both the ability of CaM to perform multiple functions in these channels and the differences in CaM regulation of the different voltage-dependent Ca(2+) channels.

Amino Acid Sequence↗

PubMed2003

Site-selective agonist binding to the nicotinic acetylcholine receptor from Torpedo californica.

Fluorescent energy transfer measurements of dansyl-C6-choline binding to the nicotinic acetylcholine receptor (AChR) from Torpedo californica were used to determine binding characteristics of the alpha gamma and alpha delta binding sites. Equilibrium binding measurements show that the alpha gamma site has a lower fluorescence than the alpha delta site; the emission difference is due to differences in the intrinsic fluorescence of the bound fluorophores rather than differences in energy transfer at the two sites. Stopped-flow fluorescence kinetics showed that dissociation of dansyl-C6-choline from the AChR in the desensitized conformation occurs 5-10-fold faster from the alpha gamma site than from the alpha delta site. The dissociation rates are robust for distinct protein preparations, in the presence of noncompetitive antagonists, and over a broad range of ionic strengths. Equilibrium fluorescent binding measurements show that dansyl-C6-choline binds with higher affinity to the alpha delta site (K = 3 nM) than to the alpha gamma site (K = 9 nM) when the AChR is desensitized. Similar affinity differences were observed for acetylcholine itself. The distinct dissociation rates permit the extent of desensitization to be measured at each site during the time course of binding. This sequential mixing method of measuring the desensitized state population at each agonist site can be applied to study the mechanism of AChR activation and subsequent desensitization in detail.

Animals↗

PubMed2002

Orientation of d-tubocurarine in the muscle nicotinic acetylcholine receptor-binding site.

Ligand modification and receptor site-directed mutagenesis were used to examine binding of the competitive antagonist, d-tubocurarine (dTC), to the muscle-type nicotinic acetylcholine receptor (AChR). By using various dTC analogs, we measured the interactions of specific dTC functional groups with amino acid positions in the AChR gamma-subunit. Because data for mutations at residue gammaTyr(117) were the most consistent with direct interaction with dTC, we focused on that residue. Double mutant thermodynamic cycle analysis showed apparent interactions of gammaTyr(117) with both the 2-N and the 13'-positions of dTC. Examination of a dTC analog with a negative charge at the 13'-position failed to reveal electrostatic interaction with charged side-chain substitutions at gamma117, but the effects of side-chain substitutions remained consistent with proximity of Tyr(117) to the cationic 2-N of dTC. The apparent interaction of gammaTyr(117) with the 13'-position of dTC was likely mediated by allosteric changes in either dTC or the receptor. The data also show that cation-pi electron stabilization of the 2-N position is not required for high affinity binding. Molecular modeling of dTC within the binding pocket of the acetylcholine-binding protein places the 2-N in proximity to the residue homologous to gammaTyr(117). This model provides a plausible structural basis for binding of dTC within the acetylcholine-binding site of the AChR family that appears consistent with findings from photoaffinity labeling studies and with site-directed mutagenesis studies of the AChR.

Animals↗