Acetylcholinesterase: Role of the enzyme's charge distribution in steering charged ligands toward the active site

The electrostatic steering of charged ligands toward the active site of Torpedo californica acetylcholinesterase is investigated by Brownian dynamics simulations of wild type enzyme and several mutated forms, in which some normally charged residues are neutralized. The simulations reveal that the total ligand influx through a surface of 42 Å radius centered in the enzyme monomer and separated from the protein surface by I-14 Å is not significantly influenced by electrostatic interactions. Electrostatic effects are visible for encounters with a surface of 32 Å radius, which is partially hidden inside the protein, but mostly within the solvent. A clear accumulation of encounter events for that sphere is observed in the area directly above the entrance to the active site gorge. In this area, the encounter events are increased by 40% compared to the case of a neutral ligand. However, the differences among the encounter rates for the various mutants considered here are not pronounced, all rate constants being within ±10% of the average value. The enzyme charge distribution becomes more important as the charged ligand moves toward the bottom of the gorge, where the active site is located. We show that neither the enzyme's total charge, nor its dipole moment, fully account for the electrostatic steering of ligand to the active site. Higher moments of the enzyme's charge distribution are also important. However, for a series of mutations for which the direction of the enzyme dipole moment is constant within a few degrees, one observes a gradual decrease in the diffusional encounter rate constant with the number of neutralized residues. On the other hand, for other mutants that change the direction of the dipole moment from that of the wild type, the calculated encounter rate constants can be very close to that of the wild type. The present work yields two new insights to the kinetics of acetylcholinesterase. First, evolution appears to have built a redundant electrostatic steering capability into this important enzyme through the overall distribution of its thousands of partially charged atoms. And second, roughly half of the rate enhancement due to electrostatics arises from steering of the substrate outside the enzyme; the other half of the rate enhancement arises from improved trapping of the substrate after it has entered the gorge. The computational results reproduce qualitatively, and help to rationalize, many surprising experimental results obtained recently for human acetylcholinesterase. © 1996 John Wiley & Sons, Inc.     

Tumor necrosis factor-alpha-induced reactive oxygen species formation is mediated by JNK1-dependent ferritin degradation and elevation of labile iron pool

Tumor necrosis factor alpha induces increased reactive oxygen species (ROS) generation in different experimental models. However, the nature of this phenomenon is still unknown. We hypothesized that TNF-induced ROS formation is due to JNK-regulated ferritin degradation and an increase in labile iron pool (LIP). We used as a model human prostate cancer cells, DU145. TNF treatment induced ROS formation, which was reduced to the control level in cells pretreated with desferrioxamine, an iron chelator. TNF induced a drop in light chain of the ferritin level, as judged by immunoblotting and an increase in LIP, evaluated by calcein fluorescence. Moreover, we observed that the JNK inhibitor SP600125 abolished TNF-induced changes in LIP, which suggests that JNK kinases are involved in this process. To explore which one of the JNK kinases is responsible for these effects, DU145 cells were transiently transfected with plasmids encoding inactive mutants of JNK1 or JNK2. The cells expressing inactive JNK1 mutant, but not cells expressing JNK2 mutant or possessing an empty vector, were completely resistant to TNF-induced ROS generation, ferritin degradation, and an increase in LIP. These data suggest that TNF-induced ROS formation is mediated by JNK1, which regulates ferritin degradation and thus the level of highly reactive iron.     

Evolution of the 28S ribosomal RNA gene in anurans: regions of variability and their phylogenetic implications

Fifteen restriction sites were mapped to the 28S ribosomal RNA gene of individuals representing 54 species of frogs, two species of salamanders, a caecilian, and a lungfish. Eight of these sites were present in all species examined, and two were found in all but one species. Alignment of these conserved restriction sites revealed, among anuran 28S rRNA genes, five regions of major length variation that correspond to four of 12 previously identified divergent domains of this gene. One of the divergent domains (DD8) consists of two regions of length variation separated by a short segment that is conserved at least throughout tetrapods. Most of the insertions, deletions, and restriction-site variations identified in the 28S gene will require sequence-level analysis for a detailed reconstruction of their history. However, an insertion in DD9 that is coextensive with frogs in the suborder Neobatrachia, a BstEII site that is limited to representatives of two leptodactylid subfamilies, and a deletion in DD10 that is found only in three ranoid genera are probably synapomorphies.      

Prediction of titration properties of structures of a protein derived from molecular dynamics trajectories.

This paper explores the dependence of the molecular dynamics (MD) trajectory of a protein molecule on the titration state assigned to the molecule. Four 100-ps MD trajectories of bovine pancreatic trypsin inhibitor (BPTI) were generated, starting from two different structures, each of which was held in two different charge states. The two starting structures were the X-ray crystal structure and one of the solution structures determined by NMR, and the charge states differed only in the ionization state of N terminus. Although it is evident that the MD simulations were too short to sample fully the equilibrium distribution of structures in each case, standard Poisson-Boltzmann titration state analysis of the resulting configurations shows general agreement between the overall titration behavior of the protein and the charge state assumed during MD simulation: at pH 7, the total net charge of the protein resulting from the titration analysis is consistently lower for the protein with the N terminus assumed to be neutral than for the protein with the N terminus assumed to be charged. For most of the ionizable residues, the differences in the calculated pKaS among the four trajectories are statistically negligible and remain in good agreement with the data obtained by crystal structure titration and by experiment. The exceptions include the N terminus, which responds directly to the change of its imposed charge; the C terminus, which in the NMR structure interacts strongly with the former; and a few other residues (Arg 1, Glu 7, Tyr 35, and Arg 42) whose pKaS reflect the initial structure and the limited trajectory lengths. This study illustrates the importance of the careful assignment of protonation states at the start of MD simulations and points to the need for simulation methods that allow for the variation of the protonation state in the calculation of equilibrium properties.     

Acetylcholinesterase: diffusional encounter rate constants for dumbbell models of ligand.

For some enzymes, virtually every substrate molecule that encounters the entrance to the active site proceeds to reaction, at low substrate concentrations. Such diffusion-limited enzymes display high apparent bimolecular rate constants ((kcat/KM)), which depend strongly upon solvent viscosity. Some experimental studies provide evidence that acetylcholinesterase falls into this category. Interestingly, the asymmetric charge distribution of acetylcholinesterase, apparent from the crystallographic structure, suggests that its electrostatic field accelerates the encounter of its cationic substrate, acetylcholine, with the entrance to the active site. Here we report simulations of the diffusion of substrate in the electrostatic field of acetylcholinesterase. We find that the field indeed guides the substrate to the mouth of the active site. The computed encounter rate constants depend upon the particular relative geometries of substrate and enzyme that are considered to represent successful encounters. With loose reaction criteria, the computed rates exceed those measured experimentally, but the rate constants vary appropriately with ionic strength. Although more restrictive reaction criteria lower the computed rates, they also lead to unrealistic variation of the rate constants with ionic strength. That these simulations do not agree well with experiment suggests that the simple diffusion model is incomplete. Structural fluctuations in the enzyme or events after the encounter may well contribute to rate limitation.     

Arabidopsis TCH4, regulated by hormones and the environment, encodes a xyloglucan endotransglycosylase.

Adaptation of plants to environmental conditions requires that sensing of external stimuli be linked to mechanisms of morphogenesis. The Arabidopsis TCH (for touch) genes are rapidly upregulated in expression in response to environmental stimuli, but a connection between this molecular response and developmental alterations has not been established. We identified TCH4 as a xyloglucan endotransglycosylase by sequence similarity and enzyme activity. Xyloglucan endotransglycosylases most likely modify cell walls, a fundamental determinant of plant form. We determined that TCH4 expression is regulated by auxin and brassinosteroids, by environmental stimuli, and during development, by a 1-kb region. Expression was restricted to expanding tissues and organs that undergo cell wall modification. Regulation of genes encoding cell wall-modifying enzymes, such as TCH4, may underlie plant morphogenetic responses to the environment.     

Characterization of cross-bridge elasticity and kinetics of cross-bridge cycling during force development in single smooth muscle cells

Force development in smooth muscle, as in skeletal muscle, is believed to reflect recruitment of force-generating myosin cross-bridges. However, little is known about the events underlying cross-bridge recruitment as the muscle cell approaches peak isometric force and then enters a period of tension maintenance. In the present studies on single smooth muscle cells isolated from the toad (Bufo marinus) stomach muscularis, active muscle stiffness, calculated from the force response to small sinusoidal length changes (0.5% cell length, 250 Hz), was utilized to estimate the relative number of attached cross-bridges. By comparing stiffness during initial force development to stiffness during force redevelopment immediately after a quick release imposed at peak force, we propose that the instantaneous active stiffness of the cell reflects both a linearly elastic cross-bridge element having 1.5 times the compliance of the cross-bridge in frog skeletal muscle and a series elastic component having an exponential length-force relationship. At the onset of force development, the ratio of stiffness to force was 2.5 times greater than at peak isometric force. These data suggest that, upon activation, cross-bridges attach in at least two states (i.e., low-force-producing and high-force-producing) and redistribute to a steady state distribution at peak isometric force. The possibility that the cross-bridge cycling rate was modulated with time was also investigated by analyzing the time course of tension recovery to small, rapid step length changes (0.5% cell length in 2.5 ms) imposed during initial force development, at peak force, and after 15 s of tension maintenance. The rate of tension recovery slowed continuously throughout force development following activation and slowed further as force was maintained. Our results suggest that the kinetics of force production in smooth muscle may involve a redistribution of cross-bridge populations between two attached states and that the average cycling rate of these cross-bridges becomes slower with time during contraction.     

Charge-charge interactions are key determinants of the pK values of ionizable groups in ribonuclease Sa (pI=3.5) and a basic variant (pI=10.2)

The pK values of the titratable groups in ribonuclease Sa (RNase Sa) (pI=3.5), and a charge-reversed variant with five carboxyl to lysine substitutions, 5K RNase Sa (pI=10.2), have been determined by NMR at 20 degrees C in 0.1M NaCl. In RNase Sa, 18 pK values and in 5K, 11 pK values were measured. The carboxyl group of Asp33, which is buried and forms three intramolecular hydrogen bonds in RNase Sa, has the lowest pK (2.4), whereas Asp79, which is also buried but does not form hydrogen bonds, has the most elevated pK (7.4). These results highlight the importance of desolvation and charge-dipole interactions in perturbing pK values of buried groups. Alkaline titration revealed that the terminal amine of RNase Sa and all eight tyrosine residues have significantly increased pK values relative to model compounds.A primary objective in this study was to investigate the influence of charge-charge interactions on the pK values by comparing results from RNase Sa with those from the 5K variant. The solution structures of the two proteins are very similar as revealed by NMR and other spectroscopic data, with only small changes at the N terminus and in the alpha-helix. Consequently, the ionizable groups will have similar environments in the two variants and desolvation and charge-dipole interactions will have comparable effects on the pK values of both. Their pK differences, therefore, are expected to be chiefly due to the different charge-charge interactions. As anticipated from its higher net charge, all measured pK values in 5K RNase are lowered relative to wild-type RNase Sa, with the largest decrease being 2.2 pH units for Glu14. The pK differences (pK(Sa)-pK(5K)) calculated using a simple model based on Coulomb's Law and a dielectric constant of 45 agree well with the experimental values. This demonstrates that the pK differences between wild-type and 5K RNase Sa are mainly due to changes in the electrostatic interactions between the ionizable groups. pK values calculated using Coulomb's Law also showed a good correlation (R=0.83) with experimental values. The more complex model based on a finite-difference solution to the Poisson-Boltzmann equation, which considers desolvation and charge-dipole interactions in addition to charge-charge interactions, was also used to calculate pK values. Surprisingly, these values are more poorly correlated (R=0.65) with the values from experiment. Taken together, the results are evidence that charge-charge interactions are the chief perturbant of the pK values of ionizable groups on the protein surface, which is where the majority of the ionizable groups are positioned in proteins.     

Protonation free energy levels in complex molecular systems

All proteins, nucleic acids, and other biomolecules contain residues capable of exchanging protons with their environment. These proton transfer phenomena lead to pH sensitivity of many molecular processes underlying biological phenomena. In the course of biological evolution, Nature has invented some mechanisms to use pH gradients to regulate biomolecular processes inside cells or in interstitial fluids. Therefore, an ability to model protonation equilibria in molecular systems accurately would be of enormous value for our understanding of biological processes and for possible rational influence on them, like in developing pH dependent drugs to treat particular diseases. This work presents a derivation, by thermodynamic and statistical mechanical methods, of an expression for the free energy of a complex molecular system at arbitrary ionization state of its titratable residues. This constitutes one of the elements of modeling protonation equilibria. Starting from a consideration of a simple acid-base equilibrium of a model compound with a single tritratable group, we arrive at an expression which is of general validity for complex systems. The only approximation used in this derivation is the postulating that the interaction energy between any pair of titratable sites does not depend on the protonation states of all the remaining ionizable groups.     

Methyl-beta-cyclodextrin induces mitochondrial cholesterol depletion and alters the mitochondrial structure and bioenergetics

There is growing evidence of mitochondrial membrane raft-like microdomains that are involved in the apoptotic pathway. The aim of this study was to investigate the effect of methyl-beta-cyclodextrin (MβCD), being a well-known lipid microdomain disrupting agent and cholesterol chelator, on the structure and bioenergetics of rat liver mitochondria (RLM). We observed that MβCD decreases the function of RLM, induces changes in the mitochondrial configuration state and decreases the calcium chloride-induced swelling. These data suggest that disruption of mitochondrial raft-like microdomains by cholesterol efflux on one hand impairs mitochondrial bioenergetics, but on the other hand it protects the mitochondria from swelling.