A novel small molecule inhibitor of p32 mitochondrial protein overexpressed in glioma

BACKGROUND

The mitochondrial protein p32 is a validated therapeutic target of cancer overexpressed in glioma. Therapeutic targeting of p32 with monoclonal antibody or p32-binding LyP-1 tumor-homing peptide can limit tumor growth. However, these agents do not specifically target mitochondrial-localized p32 and would not readily cross the blood-brain barrier to target p32-overexpressing gliomas. Identifying small molecule inhibitors of p32 overexpressed in cancer is a more rational therapeutic strategy. Thus, in this study we employed a pharmacophore modeling strategy to identify small molecules that could bind and inhibit mitochondrial p32.

METHODS

A pharmacophore model of C1q and LyP-1 peptide association with p32 was used to screen a virtual compound library. A primary screening assay for inhibitors of p32 was developed to identify compounds that could rescue p32-dependent glutamine-addicted glioma cells from glutamine withdrawal. Inhibitors from this screen were analyzed for direct binding to p32 by fluorescence polarization assay and protein thermal shift. Affect of the p32 inhibitor on glioma cell proliferation was assessed by Alamar Blue assay, and affect on metabolism was examined by measuring lactate secretion.

RESULTS

Identification of a hit compound (M36) validates the pharmacophore model. M36 binds directly to p32 and inhibits LyP-1 tumor homing peptide association with p32 in vitro. M36 effectively inhibits the growth of p32 overexpressing glioma cells, and sensitizes the cells to glucose depletion.

CONCLUSIONS

This study demonstrates a novel screening strategy to identify potential inhibitors of mitochondrial p32 protein overexpressed in glioma. High throughput screening employing this strategy has potential to identify highly selective, potent, brain-penetrant small molecules amenable for further drug development.

Recognition of dextran-superparamagnetic iron oxide nanoparticle conjugates (Feridex) via macrophage scavenger receptor charged domains

Dextran-coated superparamagnetic iron oxide nanoparticles (dextran-SPIO conjugates) offer the attractive possibility of enhancing MRI imaging sensitivity so that small or diffuse lesions can be detected. However, systemically injected SPIOs are rapidly removed by macrophages. We engineered embryonic cells (HEK293T) to express major macrophage scavenger receptor (SR) subtypes including SR-AI, MARCO, and endothelial receptor collectin-12. These SRs possess a positively charged collagen-like (CL) domain and they promoted SPIO uptake, while the charge neutral lipoprotein receptor SR-BI did not. In silico modeling indicated a positive net charge on the CL domain and a net negative charge on the cysteine-rich (CR) domain of MARCO and SR-AI. In vitro experiments revealed that CR domain deletion in SR-AI boosted uptake of SPIO 3-fold, while deletion of MARCO's CR domain abolished this uptake. These data suggest that future studies might productively focus on the validation and further exploration of SR charge fields in SPIO recognition.

Synthesis of 5-(6-hydroxy-7H-purine-8-ylthio)- 2-(N-hydroxyformamido)pentanoic acid

We have developed a synthetic route for the preparation of a hybrid bisubstrate small molecule based on a nucleoside. A prototype compound was designed and docked into the catalytic domain of the AdSS enzyme bridging the region between the magnesium center of the protein to the nucleoside region. The synthesis involves coupling a brominated peptide fragment capable of complexing magnesium to a thiolated nucleoside to obtain the hybrid model compound.

Coherent expression chromosome cluster analysis reveals differential regulatory functions of amino-terminal and distal parathyroid hormone-related protein domains in prostate carcinoma

Parathyroid hormone-related protein (PTHrP) has a number of cancer-related actions. While best known for causing hypercalcemia of malignancy, it also has effects on cancer cell growth, apoptosis, and angiogenesis. Studying the actions of PTHrP in human cancer is complicated because there are three isoforms and many derived peptides. Several peptides are biologically active at known or presumed cell surface receptors; in addition, the PTHrP-derived molecules can exert effects at the cell nucleus. To address this complexity, we studied gene expression in a DU 145 prostate cancer cell line that was stably transfected with control vector, PTHrP 1-173 and PTHrP 33-173. With this model, regulatory effects of the amino-terminal portion of PTHrP would result only from transduction with the full-length molecule, while effects pertaining to distal sequences would be evident with either construct. Analysis of the expression profiles by microarrays demonstrated nonoverlapping groups of differentially expressed genes. Amino-terminal PTHrP affected groups of genes involved in apoptosis, prostaglandin and sex steroid metabolism, cell-matrix interactions, and cell differentiation, while PTHrP 33-173 caused substantial increases in MHC class I antigen expression. This work demonstrates the distinct biological actions of the amino-terminus compared to distal mid-molecule or carboxy-terminal sequences of PTHrP in prostate carcinoma cells and provides targets for further study of the malignant process.

Role of synucleins in Alzheimer's disease

Alzheimer's disease (AD) and Parkinson's disease (PD) are the most common causes of dementia and movement disorders in the elderly. While progressive accumulation of oligomeric amyloid-beta protein (Abeta) has been identified as one of the central toxic events in AD leading to synaptic dysfunction, accumulation of alpha-synuclein (alpha-syn) resulting in the formation of oligomers has been linked to PD. Most of the studies in AD have been focused on investigating the role of Abeta and Tau; however, recent studies suggest that alpha-syn might also play a role in the pathogenesis of AD. For example, fragments of alpha-syn can associate with amyloid plaques and Abeta promotes the aggregation of alpha-syn in vivo and worsens the deficits in alpha-syn tg mice. Moreover, alpha-syn has also been shown to accumulate in limbic regions in AD, Down's syndrome, and familial AD cases. Abeta and alpha-syn might directly interact under pathological conditions leading to the formation of toxic oligomers and nanopores that increase intracellular calcium. The interactions between Abeta and alpha-syn might also result in oxidative stress, lysosomal leakage, and mitochondrial dysfunction. Thus, better understanding the steps involved in the process of Abeta and alpha-syn aggregation is important in order to develop intervention strategies that might prevent or reverse the accumulation of toxic proteins in AD.

Interaction of organic cations with organic anion transporters

Studies of the organic anion transporters (Oats) have focused mainly on their interactions with organic anionic substrates. However, as suggested when Oat1 was originally identified as NKT (Lopez-Nieto, C. E., You, G., Bush, K. T., Barros, E. J., Beier, D. R., and Nigam, S. K. (1997) J. Biol. Chem. 272, 6471-6478), since the Oats share close homology with organic cation transporters (Octs), it is possible that Oats interact with cations as well. We now show that mouse Oat1 (mOat1) and mOat3 and, to a lesser degree, mOat6 bind a number of "prototypical" Oct substrates, including 1-methyl-4-phenylpyridinium. In addition to oocyte expression assays, we have tested binding of organic cations to Oat1 and Oat3 in ex vivo assays by analyzing interactions in kidney organ cultures deficient in Oat1 and Oat3. We also demonstrate that mOat3 transports organic cations such as 1-methyl-4-phenylpyridinium and cimetidine. A pharmacophore based on the binding affinities of the tested organic cations for Oat3 was generated. Using this pharmacophore, we screened a chemical library and were able to identify novel cationic compounds that bound to Oat1 and Oat3. These compounds bound Oat3 with an affinity higher than the highest affinity compounds in the original set of prototypical Oct substrates. Thus, whereas Oat1, Oat3, and Oat6 appear to function largely in organic anion transport, they also bind and transport some organic cations. These findings could be of clinical significance, since drugs and metabolites that under normal physiological conditions do not bind to the Oats may undergo changes in charge and become Oat substrates during pathologic conditions wherein significant variations in body fluid pH occur.

Analysis of metagene portraits reveals distinct transitions during kidney organogenesis

Organogenesis is a multistage process, but it has been difficult, by conventional analysis, to separate stages and identify points of transition in developmentally complex organs or define genetic pathways that regulate pattern formation. We performed a detailed time-series examination of global gene expression during kidney development and then represented the resulting data as self-organizing maps (SOMs), which reduced more than 30,000 genes to 650 metagenes. Further clustering of these maps identified potential stages of development and suggested points of stability and transition during kidney organogenesis that are not obvious from either standard morphological analyses or conventional microarray clustering algorithms. We also performed entropy calculations of SOMs generated for each day of development and found correlations with morphometric parameters and expression of candidate genes that may help in orchestrating the transitions between stages of kidney development, as well as macro- and micropatterning of the organ.

Thermal denaturation of wild type and mutant recombinant acetylcholinesterase from amphioxus: effects of the temperature of in vitro expression and of reversible inhibitors

We have studied the thermal inactivation at 37 degrees C of wild type and mutant ChE2 (C310A, F312I, C466A, C310A/F312I, and C310A/C466A) from amphioxus (Branchiostoma floridae) expressed in vitro in COS-7 monkey cells under three sets of conditions: 30 degrees C for 48 h, 30 degrees C for 24 h and 37 degrees C for 24 h, and 37 degrees C for 48 h. We found biphasic denaturation curves for all enzymes and conditions, except wild type and C310A ChE2 expressed at 30 degrees C for 48 h. Generally, single mutants are more unstable than wild type, and the double mutants are even more unstable. We propose a model involving stable and unstable conformations of the enzymes to explain these results, and we discuss the implications of the model. We also found a correlation between the melting temperature of the ChEs and the rates at which they denature at 37 degrees C, with the denaturation of the unstable conformation dominating the relationship. Reversible cholinergic inhibitors protect the ChEs from thermal denaturation, and in some cases produce monophasic denaturation curves; we also propose a model to explain this stabilization.

Fibrinogen signal transduction as a mediator and therapeutic target in inflammation: lessons from multiple sclerosis

The blood protein fibrinogen as a ligand for integrin and non-integrin receptors functions as the molecular nexus of coagulation, inflammation and immunity. Studies in animal models and in human disease have demonstrated that extravascular fibrinogen that is deposited in tissues upon vascular rupture is not merely a marker, but a mediator of diseases with an inflammatory component, such as rheumatoid arthritis, multiple sclerosis, sepsis, myocardial infarction and bacterial infection. The present article focuses on the recent discoveries of specific cellular targets and receptors for fibrinogen within tissues that have extended the role of fibrinogen from a coagulation factor to a regulator of inflammation and immunity. Fibrinogen has the potential for selective drug targeting that would target its proinflammatory properties without affecting its beneficial effects in hemostasis, since it interacts with different receptors to mediate blood coagulation and inflammation. Strategies to target receptors for fibrinogen and fibrin within the tissue microenvironment could reveal selective and disease-specific agents for therapeutic intervention in a variety of human diseases associated with fibrin deposition.

p75 neurotrophin receptor regulates tissue fibrosis through inhibition of plasminogen activation via a PDE4/cAMP/PKA pathway

Clearance of fibrin through proteolytic degradation is a critical step of matrix remodeling that contributes to tissue repair in a variety of pathological conditions, such as stroke, atherosclerosis, and pulmonary disease. However, the molecular mechanisms that regulate fibrin deposition are not known. Here, we report that the p75 neurotrophin receptor (p75^NTR), a TNF receptor superfamily member up-regulated after tissue injury, blocks fibrinolysis by down-regulating the serine protease, tissue plasminogen activator (tPA), and up-regulating plasminogen activator inhibitor-1 (PAI-1). We have discovered a new mechanism in which phosphodiesterase PDE4A4/5 interacts with p75^NTR to enhance cAMP degradation. The p75^NTR-dependent down-regulation of cAMP results in a decrease in extracellular proteolytic activity. This mechanism is supported in vivo in p75^NTR-deficient mice, which show increased proteolysis after sciatic nerve injury and lung fibrosis. Our results reveal a novel pathogenic mechanism by which p75^NTR regulates degradation of cAMP and perpetuates scar formation after injury.

Cystic fibrosis transmembrane conductance regulator (CFTR) functionality is dependent on coatomer protein I (COPI)

BACKGROUND INFORMATION

Cystic fibrosis results from mutations in the ABC transporter CFTR (cystic fibrosis transmembrane conductance regulator), which functions as a cAMP-regulated anion channel. The most prevalent mutation in CFTR, the Phe(508) deletion, results in the generation of a trafficking and functionally deficient channel. The cellular machineries involved in modulating CFTR trafficking and function have not been fully characterized. In the present study, we identified a role for the COPI (coatomer protein I) cellular trafficking machinery in the development of the CFTR polypeptide into a functional chloride channel. To examine the role of COPI in CFTR biosynthesis, we employed the cell line ldlF, which harbours a temperature-sensitive mutation in epsilon-COP, a COPI subunit, to inhibit COPI function and then determined whether the CFTR polypeptide produced from the transfected gene developed into a cAMP-regulated chloride channel.

RESULTS

When COPI was inactivated in the ldlF cells by an elevated temperature pulse (39 degrees C), the CFTR polypeptide was detected on the cell surface by immunofluorescence microscopy and cell-surface biotinylation. Therefore, CFTR proceeded upstream in the secretory pathway in the absence of COPI function, a result demonstrated previously by others. In contrast, electrophysiological measurements indicated an absence of cAMP-stimulated chloride efflux, suggesting that channel function was impaired. In comparison, expression of CFTR at the same elevated temperature (39 degrees C) in an epsilon-COP-rescued cell line [ldlF(ldlF)], which has an introduced wild-type epsilon-COP gene in addition to the mutant epsilon-COP gene, showed restoration of cAMP-stimulated channel activity, confirming the requirement of COPI for channel function.

CONCLUSIONS

These results therefore suggest that generation of the folded-functional conformation of CFTR requires COPI.

Synaptic arrangement of the neuroligin/beta-neurexin complex revealed by X-ray and neutron scattering

Neuroligins are postsynaptic cell-adhesion proteins that associate with their presynaptic partners, the neurexins. Using small-angle X-ray scattering, we determined the shapes of the extracellular region of several neuroligin isoforms in solution. We conclude that the neuroligins dimerize via the characteristic four-helix bundle observed in cholinesterases, and that the connecting sequence between the globular lobes of the dimer and the cell membrane is elongated, projecting away from the dimer interface. X-ray scattering and neutron contrast variation data show that two neurexin monomers, separated by 107 A, bind at symmetric locations on opposite sides of the long axis of the neuroligin dimer. Using these data, we developed structural models that delineate the spatial arrangements of different neuroligin domains and their partnering molecules. As mutations of neurexin and neuroligin genes appear to be linked to autism, these models provide a structural framework for understanding altered recognition by these proteins in neurodevelopmental disorders.

A tryptophan in the bottleneck of the catalytic gorge of an invertebrate acetylcholinesterase confers relative resistance to carbamate and organophosphate inhibitors

Amphioxus, an invertebrate chordate, has two acetylcholinesterases (AChEs): cholinesterase 1 (ChE1) and cholinesterase 2 (ChE2). ChE1 is up to 329-fold more resistant to a variety of carbamate and organophosphate inhibitors, including a number of insecticides, when compared with ChE2. One difference between the two enzymes is at the position homologous to Phe331 in Torpedo AChE. In Torpedo AChE, this residue is a component of the hydrophobic subsite and defines one side of the bottleneck in the catalytic gorge of the enzyme. In ChE1, the homologous residue is Trp353; in ChE2, it is Phe353. We used site-directed mutagenesis to investigate the proposal that the resistance of ChE1 to inhibition by carbamates and organophosphates was due to this difference, creating a ChE1 W353F mutant to widen the bottleneck. The mutation virtually abolishes the difference in sensitivity to the inhibitors. The ChE1 W353F mutant is only 2- to 3-fold more resistant than ChE2 to carbamates and is actually 2.5- to 10-fold more sensitive to inhibition by organophosphates. The differences in resistance are due to different affinities of the enzymes for the inhibitors, not different reactivities. Molecular modeling supports the proposal that the difference in inhibition is due to the width of the bottleneck of the gorge. Our results have implications for insecticide resistance in insects, in particular mosquitoes and aphids.

The crystal structure of beryllofluoride Spo0F in complex with the phosphotransferase Spo0B represents a phosphotransfer pretransition state

A number of regulatory circuits in biological systems function through the exchange of phosphoryl groups from one protein to another. Spo0F and Spo0B are components of a phosphorelay that control sporulation in the bacterium Bacillus subtilis through the exchange of a phosphoryl group. Using beryllofluoride as a mimic for phosphorylation, we trapped the interaction of the phosphorylated Spo0F with Spo0B in the crystal lattice. The transition state of phosphoryl transfer continues to be a highly debated issue, as to whether it is associative or dissociative in nature. The geometry of Spo0F binding to Spo0B favors an associative mechanism for phosphoryl transfer. In order to visualize the autophosphorylation of the histidine kinase, KinA, and the subsequent phosphoryl transfer to Spo0F, we generated in silico models representing these reaction steps.

Alpha-conotoxin OmIA is a potent ligand for the acetylcholine-binding protein as well as alpha3beta2 and alpha7 nicotinic acetylcholine receptors

The molluskan acetylcholine-binding protein (AChBP) is a homolog of the extracellular binding domain of the pentameric ligand-gated ion channel family. AChBP most closely resembles the alpha-subunit of nicotinic acetylcholine receptors and in particular the homomeric alpha7 nicotinic receptor. We report the isolation and characterization of an alpha-conotoxin that has the highest known affinity for the Lymnaea AChBP and also potently blocks the alpha7 nAChR subtype when expressed in Xenopus oocytes. Remarkably, the peptide also has high affinity for the alpha3beta2 nAChR indicating that alpha-conotoxin OmIA in combination with the AChBP may serve as a model system for understanding the binding determinants of alpha3beta2 nAChRs. alpha-Conotoxin OmIA was purified from the venom of Conus omaria. It is a 17-amino-acid, two-disulfide bridge peptide. The ligand is the first alpha-conotoxin with higher affinity for the closely related receptor subtypes, alpha3beta2 versus alpha6beta2, and selectively blocks these two subtypes when compared with alpha2beta2, alpha4beta2, and alpha1beta1deltaepsilon nAChRs.

Inactivation of an invertebrate acetylcholinesterase by sulfhydryl reagents: the roles of two cysteines in the catalytic gorge of the enzyme

We have used site-directed mutagenesis and molecular modeling to investigate the inactivation of an invertebrate acetylcholinesterase (AChE), ChE2 from amphioxus, by the sulfhydryl reagents 5,5'-dithiobis(2-nitrobenzoic acid) (DTNB) and N-ethylmaleimide (NEM), creating various mutants, including C310A and C466A, and the double mutants C310A/C466A and C310A/F312I, to assess the relative roles of the two cysteines and a proposal that the increased rate of inactivation in the F312I mutant is due to increased access to Cys310. Our results suggest that both cysteines may be involved in inactivation by sulfhydryl reagents, but that the cysteine in the vicinity of the acyl pocket is more accessible. We speculate that the inactivation of aphid AChEs by sulfhydryl reagents is due to the presence of a cysteine homologous to Cys310. We also investigated the effects of various reversible cholinergic ligands, which bind to different subsites of the active site of the enzyme, on the rate of inactivation by DTNB of wild type ChE2 and ChE2 F312I. For the most part the inhibitors protect the enzymes from inactivation by DTNB. However, a notable exception is the peripheral site ligand propidium, which accelerates inactivation in the wild type ChE2, but retards inactivation in the F312I mutant. We propose that these opposing effects are the result of an altered allosteric signal transduction mechanism in the F312I mutant compared to the wild type ChE2.

Dissection of Synapse Induction by Neuroligins

To study synapse formation by neuroligins, we co-cultured hippocampal neurons with COS cells expressing wild type and mutant neuroligins. The large size of COS cells makes it possible to test the effect of neuroligins presented over an extended surface area. We found that a uniform lawn of wild type neuroligins displayed on the cell surface triggers the formation of hundreds of uniformly sized, individual synaptic contacts that are labeled with neurexin antibodies. Electron microscopy revealed that these artificial synapses contain a presynaptic active zone with docked vesicles and often feature a postsynaptic density. Neuroligins 1, 2, and 3 were active in this assay. Mutations in two surface loops of neuroligin 1 abolished neuroligin binding to neurexin 1beta, a presumptive presynaptic binding partner for postsynaptic neuroligins, and blocked synapse formation. An analysis of mutant neuroligins with an amino acid substitution that corresponds to a mutation described in patients with an autistic syndrome confirmed previous reports that these mutant neuroligins have a compromised capacity to be transported to the cell surface. Nevertheless, the small percentage of mutant neuroligins that reached the cell surface still induced synapse formation. Viewed together, our data suggest that neuroligins generally promote artificial synapse formation in a manner that is associated with beta-neurexin binding and results in morphologically well differentiated synapses and that a neuroligin mutation found in autism spectrum disorders impairs cell-surface transport but does not completely abolish synapse formation activity.

Identification of Molecular Determinants That Modulate Trafficking of ΔF508 CFTR, the Mutant ABC Transporter Associated With Cystic Fibrosis

Cystic fibrosis is a life-shortening inherited disorder associated primarily with a three-base in frame deletion that eliminates Phe508 in the ABC transporter, cystic fibrosis transmembrane conductance regulator (CFTR). Mutant CFTR, designated deltaF508 CFTR, is misprocessed and retained intracellularly. It is unclear what causes the trafficking impairment despite extensive investigative effort and the disease's prevalence. We hypothesize that the trafficking impairment is mediated by "receptors" of the cellular trafficking machinery that at three sequential "trafficking checkpoints" govern (1) exit from the endoplasmic reticulum (ER), (2) Golgi to the ER retrieval, and (3) targeting from post-Golgi compartments to lysosomes. We propose that, because of the Phe508 deletion and polypeptide misfolding: (1) a forward-directing signal recognized by the sec24 component of the COPII complex that mediates ER exit is eliminated; (2) a basic amino acid signal recognized by the COPI machinery involved in Golgi to ER retrieval becomes activated; and (3) a tyrosine-based sorting signal that targets to the lysosomes likewise becomes activated. We employed recently reported crystal structures of CFTR nucleotide binding domain 1 and sec24 in computational docking models to identify the most plausible CFTR-sec24 recognition domain. Site-directed mutagenesis and heterologous expression were also used to identify amino acid sequences that operate in Golgi to ER and post-Golgi to lysosome targeting. The importance of considering a multiple checkpoint model for trafficking is that rationale design of pharmaceutical interventions would require abrogation of all major checkpoints to deliver deltaF508 CFTR to the cell surface.

Conformational preferences and activities of peptides from the catecholamine release-inhibitory (catestatin) region of chromogranin A

Previous modeling (PDB 1cfk) of the catecholamine release-inhibitory "catestatin" region of chromogranin A (CgA) suggested a beta-strand/loop/beta-strand active conformation, displaying an electropositive Arg-rich loop (R(351)AR(353)GYGFR(358)). To explore this possibility, we studied NMR structures of linear and cyclic synthetic catestatin, bovine (bCgA(344-364)) or human (hCgA(352-372)). By 2-D (1)H-NMR, the structure of linear catestatin (hCgA(352-372)) exhibited the NOE pattern of a coiled loop (PDB 1lv4). We then constrained the structure, cyclizing the putative Arg-rich loop connecting the beta-strands: cyclic bCgA(350-362) ([C(0)]F(350)RARGYGFRGPGL(362)[C(+14)]). Favored conformations of cyclic bCgA(350-362) were determined by (1)H-NMR and (13)C-NMR spectroscopy. Cyclic bCgA(350-362) conformers (PDB 1n2y) adopted a "twisted-loop" conformation. Alignment between the homology model and the cyclic NMR structure showed that, while portions of the NMR structure's mid-molecule and carboxy-terminus were congruent with the homology model (RMSD, 1.61-1.91 A), the amino-terminal "twisted loop" coiled inward and away from the model (RMSD, 3.36 A). Constrained cyclic bCgA(350-362) did not exert nicotinic cholinergic antagonist activity (IC(50)>10 microM), when compared to full-length linear (IC(50) approximately 0.42-0.56 microM), or cyclic (IC(50) approximately 0.74 microM) catestatin. Thus, loss of activity in the small, constrained peptide did not result from either [Cys]-extension or cyclization, per se. While linear catestatin displays coiled character, a small cyclic derivative lost biological activity, perhaps because its amino-terminal domain deviated sharply from the predicted active conformation. These results refine the relationship between structure and function in catestatin, and suggest goals in future peptidomimetic syntheses, in particular attempts to constrain the correct amino-terminal shape for biological activity.

Amino acids defining the acyl pocket of an invertebrate cholinesterase

Amphioxus (Branchiostoma floridae) cholinesterase 2 (ChE2) hydrolyzes acetylthiocholine (AsCh) almost exclusively. We constructed a homology model of ChE2 on the basis of Torpedo californica acetylcholinesterase (AChE) and found that the acyl pocket of the enzyme resembles that of Drosophila melanogaster AChE, which is proposed to be comprised of Phe330 (Phe290 in T. californica AChE) and Phe440 (Val400), rather than Leu328 (Phe288) and Phe330 (Phe290), as in vertebrate AChE. In ChE2, the homologous amino acids are Phe312 (Phe290) and Phe422 (Val400). To determine if these amino acids define the acyl pocket of ChE2 and its substrate specificity, and to obtain information about the hydrophobic subsite, partially comprised of Tyr352 (Phe330) and Phe353 (Phe331), we performed site-directed mutagenesis and in vitro expression. The aliphatic substitution mutant F312I ChE2 hydrolyzes AsCh preferentially but also butyrylthiocholine (BsCh), and the change in substrate specificity is due primarily to an increase in k(cat) for BsCh; K(m) and K(ss) are also altered. F422L and F422V produce enzymes that hydrolyze BsCh and AsCh equally due to an increase in k(cat) for BsCh and a decrease in k(cat) for AsCh. Our data suggest that Phe312 and Phe422 define the acyl pocket. We also screened mutants for changes in sensitivity to various inhibitors. Y352A increases the sensitivity of ChE2 to the bulky inhibitor ethopropazine. Y352A decreases inhibition by BW284c51, consistent with its role as part of the choline-binding site. Aliphatic replacement mutations produce enzymes that are more sensitive to inhibition by iso-OMPA, presumably by increasing access to the active site serine. Y352A, F353A and F353V make ChE2 considerably more resistant to inhibition by eserine and neostigmine, suggesting that binding of these aromatic inhibitors is mediated by pi-pi or cation-pi interactions at the hydrophobic site. Our results also provide information about the aromatic trapping of the active site histidine and the inactivation of ChE2 by sulfhydryl reagents.

The Arg451Cys-neuroligin-3 mutation associated with autism reveals a defect in protein processing

The neuroligins are a family of postsynaptic transmembrane proteins that associate with presynaptic partners, the beta-neurexins. Neurexins and neuroligins play a critical role in initiating formation and differentiation of synaptic junctions. A recent study reported that a mutation of neuroligin-3 (NL3), an X-linked gene, was found in siblings with autistic spectrum disorder in which two affected brothers had a point mutation that substituted a Cys for Arg451. To characterize the mutation at the biochemical level, we analyzed expression and activity of the mutated protein. Mass spectrometry comparison of the disulfide bonding pattern between the native and the mutated proteins indicates the absence of aberrant disulfide bonding, suggesting that the secondary structure of the mutated protein is conserved. However, the mutation separately affects protein expression and activity. The Cys mutation causes defective neuroligin trafficking, leading to retention of the protein in the endoplasmic reticulum. This, in turn, decreases the delivery of NL3 to the cell surface. Also, the small fraction of protein that reaches the cell membrane lacks or has markedly diminished beta-neurexin-1 (NX1beta) binding activity. Other substitutions for Arg451 allow for normal cellular expression but diminished affinity for NX1beta. Our findings reveal a cellular phenotype and loss of function for a congenital mutation associated with autistic spectrum disorders.

Structural Characterization of Recombinant Soluble Rat Neuroligin 1: Mapping of Secondary Structure and Glycosylation by Mass Spectrometry

Neuroligins (NLs) are a family of transmembrane proteins that function in synapse formation and/or remodeling by interacting with beta-neurexins (beta-NXs) to form heterophilic cell adhesions. The large N-terminal extracellular domain of NLs, required for beta-NX interactions, has sequence homology to the alpha/beta hydrolase fold superfamily of proteins. By peptide mapping and mass spectrometric analysis of a soluble recombinant form of NL1, several structural features of the extracellular domain have been established. Of the nine cysteine residues in NL1, eight are shown to form intramolecular disulfide bonds. Disulfide pairings of Cys 117 to Cys 153 and Cys 342 to Cys 353 are consistent with disulfide linkages that are conserved among the family of alpha/beta hydrolase proteins. The disulfide bond between Cys 172 and Cys 181 occurs within a region of the protein encoded by an alternatively spliced exon. The disulfide pairing of Cys 512 and Cys 546 in NL1 yields a structural motif unique to the NLs, since these residues are highly conserved. The potential N-glycosylation sequons in NL1 at Asn 109, Asn 303, Asn 343, and Asn 547 are shown occupied by carbohydrate. An additional consensus sequence for N-glycosylation at Asn 662 is likely occupied. Analysis of N-linked oligosaccharide content by mass matching paradigms reveals significant microheterogeneous populations of complex glycosyl moieties. In addition, O-linked glycosylation is observed in the predicted stalk region of NL1, prior to the transmembrane spanning domain. From predictions based on sequence homology of NL1 to acetylcholinesterase and the molecular features of NL1 established from mass spectrometric analysis, a novel topology model for NL three-dimensional structure has been constructed.

Identification of the protein kinase A regulatory RIalpha-catalytic subunit interface by amide H/2H exchange and protein docking

An important goal after structural genomics is to build up the structures of higher-order protein-protein complexes from structures of the individual subunits. Often structures of higher order complexes are difficult to obtain by crystallography. We have used an alternative approach in which the structures of the individual catalytic (C) subunit and RIalpha regulatory (R) subunit of PKA were first subjected to computational docking, and the top 100,000 solutions were subsequently filtered based on amide hydrogen/deuterium (H/2H) exchange interface protection data. The resulting set of filtered solutions forms an ensemble of structures in which, besides the inhibitor peptide binding site, a flat interface between the C-terminal lobe of the C-subunit and the A- and B-helices of RIalpha is uniquely identified. This holoenzyme structure satisfies all previous experimental data on the complex and allows prediction of new contacts between the two subunits.

Finding needles in haystacks: Reranking DOT results by using shape complementarity, cluster analysis, and biological information

We present an evaluation of our results for the first Critical Assessment of PRedicted Interaction (CAPRI). The methods used include the molecular docking program DOT, shape analysis tool FADE, cluster analysis and filtering based on biological data. Good results were obtained for most of the seven CAPRI targets, and for two systems, submissions having the highest number of correctly predicted contacts were produced.

Tryptophan fluorescence reveals conformational changes in the acetylcholine binding protein

The recent characterization of an acetylcholine binding protein (AChBP) from the fresh water snail, Lymnaea stagnalis, shows it to be a structural homolog of the extracellular domain of the nicotinic acetylcholine receptor (nAChR). To ascertain whether the AChBP exhibits the recognition properties and functional states of the nAChR, we have expressed the protein in milligram quantities from a synthetic cDNA transfected into human embryonic kidney (HEK) cells. The protein secreted into the medium shows a pentameric rosette structure with ligand stoichiometry approximating five sites per pentamer. Surprisingly, binding of acetylcholine, selective agonists, and antagonists ranging from small alkaloids to larger peptides results in substantial quenching of the intrinsic tryptophan fluorescence. Using stopped-flow techniques, we demonstrate rapid rates of association and dissociation of agonists and slow rates for the alpha-neurotoxins. Since agonist binding occurs in millisecond time frames, and the alpha-neurotoxins may induce a distinct conformational state for the AChBP-toxin complex, the snail protein shows many of the properties expected for receptor recognition of interacting ligands. Thus, the marked tryptophan quenching not only documents the importance of aromatic residues in ligand recognition, but establishes that the AChBP will be a useful functional as well as structural surrogate of the nicotinic receptor.

Residues in the epsilon subunit of the nicotinic acetylcholine receptor interact to confer selectivity of waglerin-1 for the alpha-epsilon subunit interface site

Waglerin-1 (Wtx-1) is a 22-amino acid peptide that competitively antagonizes muscle nicotinic acetylcholine receptors (nAChRs). Previous work demonstrated that Wtx-1 binds to mouse nAChRs with higher affinity than receptors from rats or humans, and distinguished residues in alpha and epsilon subunits that govern the species selectivity. These studies also showed that Wtx-1 binds selectively to the alpha-epsilon binding site with significantly higher affinity than to the alpha-delta binding site. Here we identify residues at equivalent positions in the epsilon, gamma, and delta subunits that govern Wtx-1 selectivity for one of the two binding sites on the nAChR pentamer. Using a series of chimeric and point mutant subunits, we show that residues Gly-57, Asp-59, Tyr-111, Tyr-115, and Asp-173 of the epsilon subunit account predominantly for the 3700-fold higher affinity of the alpha-epsilon site relative to that of the alpha-gamma site. Similarly, we find that residues Lys-34, Gly-57, Asp-59, and Asp-173 account predominantly for the high affinity of the alpha-epsilon site relative to that of the alpha-delta site. Analysis of combinations of point mutations reveals that Asp-173 in the epsilon subunit is required together with the remaining determinants in the epsilon subunit to achieve Wtx-1 selectivity. In particular, Lys-34 interacts with Asp-173 to confer high affinity, resulting in a DeltaDeltaG(INT) of -2.3 kcal/mol in the epsilon subunit and a DeltaDeltaG(INT) of -1.3 kcal/mol in the delta subunit. Asp-173 is part of a nonhomologous insertion not found in the acetylcholine binding protein structure. The key role of this insertion in Wtx-1 selectivity indicates that it is proximal to the ligand binding site. We use the binding and interaction energies for Wtx-1 to generate structural models of the alpha-epsilon, alpha-gamma, and alpha-delta binding sites containing the nonhomologous insertion.

Hidden Markov models-based system (HMMSPECTR) for detecting structural homologies on the basis of sequential information

HMMSPECTR is a tool for finding putative structural homologs for proteins with known primary sequences. HMMSPECTR contains four major components: a data warehouse with the hidden Markov models (HMM) and alignment libraries; a search program which compares the initial protein sequences with the libraries of HMMs; a secondary structure prediction and comparison program; and a dominant protein selection program that prepares the set of 10-15 "best" proteins from the chosen HMMs. The data warehouse contains four libraries of HMMs. The first two libraries were constructed using different HHM preparation options of the HAMMER program. The third library contains parts ("partial HMM") of initial alignments. The fourth library contains trained HMMs. We tested our program against all of the protein targets proposed in the CASP4 competition. The data warehouse included libraries of structural alignments and HMMs constructed on the basis of proteins publicly available in the Protein Data Bank before the CASP4 meeting. The newest fully automated versions of HMMSPECTR 1.02 and 1.02ss produced better results than the best result reported at CASP4 either by r.m.s.d. or by length (or both) in 64% (HMMSPECTR 1.02) and 79% (HMMSPECTR 1.02ss) of the cases. The improvement is most notable for the targets with complexity 4 (difficult fold recognition cases).

Electrostatic environment surrounding the activation loop phosphotyrosine in the oncoprotein v-Fps

Autophosphorylation of Tyr-1073 in the activation loop of the oncoprotein v-Fps enhances the phosphoryl transfer reaction without influencing substrate, ATP, or metal ion binding affinities [Saylor, P., et al. (1998) Biochemistry 37, 17875-17881]. A structural model of v-Fps, generated from the insulin receptor, indicates that pTyr-1073 chelates two arginines. Mutation of these residues to alanine (R1042A and R1066A) results in weakly phosphorylated enzymes, indicating that one electropositive center is insufficient for attaining maximum loop phosphorylation and concomitant high catalytic activity. While the turnover rate for R1066A is similar to that for a mutant lacking a phosphorylatable residue in the activation loop, the rate for R1042A is 50-fold slower. While solvent perturbation studies suggest that the former is due to a slow phosphoryl transfer step, the latter effect results from a slow conformational change in the mutant, potentially linked to motions in the catalytic loop. Binding of a stoichiometric quantity of Mg(2+) is essential for ATP binding and catalysis, while binding of an additional Mg(2+) ion activates further the wild-type enzyme. The affinity of the R1066A enzyme for the second Mg(2+) ion is 23-fold higher than that of the phosphorylated or unphosphorylated form of wild-type v-Fps, with substrate binding unaffected. Conversely, the affinity of R1066A for a substrate mimic lacking a phosphorylation site is 12-fold higher than that for the phosphorylated or unphosphorylated form of wild-type v-Fps, with binding of the second Mg(2+) ion unaffected. A comparison of these enzyme-independent parameters indicates that Arg-1042 and Arg-1066 induce strain in the active site in the repressed form of the enzyme. While this strain is not relieved in the phosphorylated form, the improvements in catalysis in activated v-Fps compensate for reduced metal and substrate binding affinities.

Subunit interface selective toxins as probes of nicotinic acetylcholine receptor structure

The pentametric assembly of the nicotinic acetylcholine receptor with two of the five subunit interfaces serving as a ligand binding sites offers an opportunity to distinguish features on the surfaces of the subunits, and their ligand specificity characteristics. The receptor from mammalian muscle, with its circular order of homologous subunits (alphagamma alphadelta beta), assembles in a unique arrangement. The residues governing assembly can be ascertained through mutagenesis. Selectivity of certain natural toxins is sufficient to distinguish between sites at the alphagamma and alphadelta subunit interfaces. By interchanging residues on the gamma and delta subunits through mutagenesis, and ascertaining how they interact with the alpha subunit, determinants forming the binding sites can be delineated. The alpha-conotoxins show a 10,000-fold preference for the alphadelta over alphagamma subunit interface with alphaepsilon falling in between. The waglerins show a 2,000-fold preference for alphaepsilon over the alphagamma and alphadelta interfaces. Finally, the alpha-neurotoxin from N. mossambica mossambica shows a 10,000-fold preference for the alphagamma and alphadelta interfaces over alphaepsilon. Identification of interactive residues through mutagenesis, when coupled with homology modeling of domains and site-directed residue modification, has revealed important elements of receptor structure.

Protein docking using continuum electrostatics and geometric fit

The computer program DOT quickly finds low-energy docked structures for two proteins by performing a systematic search over six degrees of freedom. A novel feature of DOT is its energy function, which is the sum of both a Poisson-Boltzmann electrostatic energy and a van der Waals energy, each represented as a grid-based correlation function. DOT evaluates the energy of interaction for many orientations of the moving molecule and maintains separate lists scored by either the electrostatic energy, the van der Waals energy or the composite sum of both. The free energy is obtained by summing the Boltzmann factor over all rotations at each grid point. Three important findings are presented. First, for a wide variety of protein-protein interactions, the composite-energy function is shown to produce larger clusters of correct answers than found by scoring with either van der Waals energy (geometric fit) or electrostatic energy alone. Second, free-energy clusters are demonstrated to be indicators of binding sites. Third, the contributions of electrostatic and attractive van der Waals energies to the total energy term appropriately reflect the nature of the various types of protein-protein interactions studied.

Catalytic assessment of the glycine-rich loop of the v-Fps oncoprotein using site-directed mutagenesis

The three glycine residues in the glycine-rich loop of the oncoprotein, v-Fps, were mutated to determine the function of these highly conserved residues in catalysis. The kinase domains of six mutants (G928A,S, G930A,S, and G933A,S) and the wild-type enzyme were expressed and purified as fusion proteins of glutathione-S-transferase in Escherichia coli, and their catalytic properties were assessed using steady-state kinetic, inhibition, viscosity and autophosphorylation studies. Although both G928A and G930A had no detectable activity toward the substrate peptide (EAEIYEAIE), the other mutants had apparent, but varying activities. G930S lowered the rate of phosphoryl transfer by 130-fold while G928S and G933S had smaller (6-9-fold) reductions in this step. These effects on catalytic function parallel the reductions in turnover and autophosphorylation but, for G933S and G933A, net product release is still rate limiting at saturating substrate and ATP concentrations. On the basis of K(I) measurements, the effects on turnover for these mutants may be due to improved ADP affinity. While ADP affinity is reduced 2- and 3-fold for G928S and G930S, the affinity of this product is increased by 22- and 7-fold for G933S and G933A. In contrast, ATP affinity is enhanced by 5-fold for G928S and G933S and is reduced by less than 2-fold for G930S. These complex, differential effects on nucleotide binding indicate that the glycines influence the relative affinities of ADP and ATP. On the basis of the results of serine replacements, Gly-928 and Gly-930 enhance ADP affinity by 9- and 2-fold compared to ATP affinity whereas Gly-933 diminishes ADP affinity by approximately 4-fold compared to ATP affinity. These findings demonstrate that the functions of the loop lie not only in modulating the rate of the phosphoryl transfer step but also in balancing the relative affinities of ATP and ADP. These effects on nucleotide specificity may be a contributing element for the stabilization of the phosphoryl transition state and may also facilitate quick release of bound products.

Common EF-hand motifs in cholinesterases and neuroligins suggest a role for Ca2+ binding in cell surface associations

Comparisons of protein sequence via cyclic training of Hidden Markov Models (HMMs) in conjunction with alignments of three-dimensional structure, using the Combinatorial Extension (CE) algorithm, reveal two putative EF-hand metal binding domains in acetylcholinesterase. Based on sequence similarity, putative EF-hands are also predicted for the neuroligin family of cell surface proteins. These predictions are supported by experimental evidence. In the acetylcholinesterase crystal structure from Torpedo californica, the first putative EF-hand region binds the Zn2+ found in the heavy metal replacement structure. Further, the interaction of neuroligin 1 with its cognate receptor neurexin depends on Ca2+. Thus, members of the alpha,beta hydrolase fold family of proteins contain potential Ca2+ binding sites, which in some family members may be critical for heterologous cell associations.

600 ps molecular dynamics reveals stable substructures and flexible hinge points in cAMP dependent protein kinase

Molecular dynamics simulations of the catalytic subunit of cAMP dependent protein kinase (cAPK) have been performed in an aqueous environment. The relations among the protein hydrogen-bonding network, secondary structural elements, and the internal motions of rigid domains were examined. The values of fluctuations of protein dihedral angles during dynamics show quite distinct maxima in the regions of loops and minima in the regions of alpha-helices and beta-strands. Analyses of conformation snapshots throughout the run show stable subdomains and indicate that these rigid domains are constrained during the dynamics by a stable network of hydrogen bonds. The most stable subdomain during the dynamics was in the small lobe including part of the carboxy-terminal tail. The most significant flexible region was the highly conserved glycine-rich loop between beta strands 1 and 2 in the small lobe. Many of the main chain dihedral angle changes measured in a comparison of the crystallographic structures of "open" and "closed" conformations of cAPK correspond to the highly flexible residues found during dynamics.

Analysis of cholinesterase inactivation and reactivation by systematic structural modification and enantiomeric selectivity

We show here with a congeneric series of Rp- and Sp-alkoxymethyl phosphonothiolates of known absolute stereochemistry that chiral selectivity in their reaction with acetylcholinesterase can be described in terms of discrete orientational and steric requirements. Stereoselectivity depends on acyl pocket dimensions, which govern leaving group orientation and a productive association of the phosphonyl oxygen in the oxyanion hole. Overall geometry is consistent with a pentavalent intermediate where the attacking serine and leaving group are at apical positions. Oxime reactivation of the phosphonylated enzyme occurs through a similar associative intermediate presumably forming an oxime phosphonate. The oximes of differing structure show distinct angles of attacking the phosphate where the attack angles and access to the phosphorus are constrained in the sterically impacted gorge. Hence, efficacy of oxime reactivation is dependent on both oxime and conjugated phosphonate structures.

Delineation of selective cyclic GMP-dependent protein kinase Ialpha substrate and inhibitor peptides based on combinatorial peptide libraries on paper

Peptide libraries on cellulose paper have proven to be valuable tools for the a priori determination of substrate specificities of cyclic AMP- and cyclic GMP-dependent protein kinases (cAMP-kinase and cGMP-kinase) on the basis of octa-peptide sequences. Here, we report the extension of our peptide library screens to 12-mer and 14-mer peptide sequences, resulting in highly cGMP-kinase Ialpha selective peptides. The sequences TQAKRKKSLAMA-amide and TQAKRKKSLAMFLR-amide, with Km values for cGMP-kinase Ialpha of 0.7 and 0.26 microM and Vmax values of 11.5 and 10.9 micromol/min/mg, respectively, display a high specificity for this enzyme. Furthermore, replacing the phosphate acceptor residue serine with alanine in TQAKRKKSLAMA-amide resulted in the highly cGMP-kinase Ialpha selective inhibitor peptide TQAKRKKALAMA-amide, with inhibitor constants for cGMP-kinase Ialpha and cAMP-kinase of 7.5 microM and 750 microM, respectively. Selective cGMP-kinase inhibitors have the potential to play an important role in the elucidation of the distinct cellular functions of cGMP-kinase separate from those activated by cAMP-kinases, and, therefore, may play an important role as pharmaceutical targets. Molecular docking experiments of the most cGMP-kinase selective sequences on a molecular model of the catalytic domain of cGMP-kinase Ialpha suggest that they adopt unique conformations, which differ significantly from those observed for the cAMP-kinase-specific inhibitor PKI(5-24). Our results suggest that despite their structural similarities, cAMP-kinase and cGMP-kinase use distinct peptide substrate and inhibitor conformations, which could account for their unique substrate specificities. These findings are further supported by cAMP- and cGMP-kinase-selective inhibitor analogs with (D)-Ala residues at the inhibitory positions.

Mechanism of action of chromogranin A on catecholamine release: molecular modeling of the catestatin region reveals a beta-strand/loop/beta-strand structure secured by hydrophobic interactions and predictive of activity

A novel fragment of chromogranin A, known as 'catestatin' (bovine chromogranin A344-364), inhibits catecholamine release from chromaffin cells and noradrenergic neurons by acting as a non-competitive nicotinic cholinergic antagonist, and may therefore constitute an endogenous autocrine feedback regulator of sympathoadrenal activity. To characterize how this activity depends on the peptide's structure, we searched for common 3-dimensional motifs for this primary structure or its homologs. Catestatin's primary structure bore significant (29-35.5% identity, general alignment score 44-57) sequence homology to fragment sequences within three homologs of known 3-dimensional structures, based on solved X-ray crystals: 8FAB, IPKM, and 2IG2. Each of these sequences exists in nature as a beta-strand/loop/beta-strand structure, stabilized by hydrophobic interactions between the beta-strands. The catestatin structure was stable during molecular dynamics simulations. The catestatin loop contains three Arg residues, whose electropositive side chains form the terminus of the structure, and give rise to substantial uncompensated charge asymmetry in the molecule. A hydrophobic moment plot revealed that catestatin is the only segment of chromogranin A predicted to contain amphiphilic beta-strand. Circular dichroism in the far ultraviolet showed substantial (63%) beta-sheet structure, especially in a hydrophobic environment. Alanine-substitution mutants of catestatin established a crucial role for the three central arginine residues in the loop (Arg351, Arg353, and Arg358), though not for two arginine residues in the strand region toward the amino-terminus. [125I]Catestatin bound to Torpedo membranes at a site other than the nicotinic agonist binding site. When the catestatin structure was 'docked' with the extracellular domain of the Torpedo nicotinic cholinergic receptor, it interacted principally with the beta and delta subunits, in a relatively hydrophobic region of the cation pore extracellular orifice, and the complex of ligand and receptor largely occluded the cation pore, providing a structural basis for the non-competitive nicotinic cholinergic antagonist properties of the peptide. We conclude that a homology model of catestatin correctly predicts actual features of the peptide, both physical and biological. The model suggests particular spatial and charge features of the peptide which may serve as starting points in the development of non-peptide mimetics of this endogenous nicotinic cholinergic antagonist.

Determining ligand orientation and transphosphonylation mechanisms on acetylcholinesterase by Rp, Sp enantiomer selectivity and site-specific mutagenesis

Acetylcholinesterase, an enzyme of the serine hydrolase family, catalyzes the rapid hydrolysis of certain carboxyl esters. Other acyl esters efficiently transacylate the enzyme with a subsequent, slow deacylation step. Of these, the phosphoryl and phosphonyl esters are perhaps of greatest mechanistic interest since individual enantiomers of known absolute stereochemistry can be isolated and their interactions with the dissymmetric enzyme active site examined. We describe here studies of a series of enantiomeric Rp- and Sp-alkylphosphonates interacting with mouse acetylcholinesterase. Since the acetylcholinesterase is generated by recombinant DNA methods, mutant enzymes can be made with specific replacements of individual amino acid side chains. Individual amino acid replacements in the acyl pocket, the choline subsite and at the active center gorge entry have been generated, and the reaction kinetics of the mutant enzymes analyzed. These studies have shown that substitution of aliphatic amino acids for phenylalanines 295 and 297 in the acyl pocket diminishes, and in some cases, actually inverts chiral preferences. The combined structure-activity approach, where both ligand and enzyme are modified systematically, has enabled us to show that the restricted dimensions of the acyl pocket in the active center dictate enantiomeric selectivity. Moreover, the reactions of compounds of known absolute stereochemistry show three requirements for efficient transphosphonylation: (a) apposition of the phosphate with the gamma-oxygen on Ser 203 to form a pentavalent, presumed trigonal bipyramidal intermediate, (b) polarization of the phosphonyl oxygen bond by its positioning in the oxyanion hole, and (c) positioning the leaving group towards the gorge exit.

Kinetic analyses of mutations in the glycine-rich loop of cAMP-dependent protein kinase

The conserved glycines in the glycine-rich loop (Leu-Gly50-Thr-Gly52-Ser-Phe-Gly55-Arg-Val) of the catalytic (C) subunit of cAMP-dependent protein kinase were each mutated to Ser (G50S, G52S, and G55S). The effects of these mutations were assessed here using both steady-state and pre-steady-state kinetic methods. While G50S and G52S reduced the apparent affinity for ATP by approximately 10-fold, substitution at Gly55 had no effect on nucleotide binding. In contrast to ATP, only mutation at position 50 interfered with ADP binding. These three mutations lowered the rate of phosphoryl transfer by 7-300-fold. The combined data indicate that G50 and G52 are the most critical residues in the loop for catalysis, with replacement at position 52 being the most extreme owing to a larger decrease in the rate of phosphoryl transfer (29 vs 1.6 s-1 in contrast to 500 s-1 for wild-type C). Surprisingly, all three mutations lowered the affinity for Kemptide by approximately 10-fold, although none of the loop glycines makes direct contact with the substrate. The inability to correlate the rate constant for net product release with the dissociation constant for ADP implies that other steps may limit the decomposition of the ternary product complex. The observations that G52S (a) selectively affects ATP binding and (b) significantly lowers the rate of phosphoryl transfer without making direct contact with either the nucleotide or the peptide imply that this residue serves a structural role in the loop, most likely by positioning the backbone amide of Ser53 for contacting the gamma-phosphate of ATP. Energy-minimized models of the mutant proteins are consistent with the observed kinetic consequences of each mutation. The models predict that only mutation of Gly52 will interfere with the observed hydrogen bonding between the backbone and ATP.

Identification of pairwise interactions in the alpha-neurotoxin-nicotinic acetylcholine receptor complex through double mutant cycles

alpha-Neurotoxins are potent inhibitors of the nicotinic acetylcholine receptor (nAChR), binding with high affinity to the two agonist sites located on the extracellular domain. Previous site-directed mutagenesis had identified three residues on the alpha-neurotoxin from Naja mossambica mossambica (Lys27, Arg33, and Lys47) and four residues on the mouse muscle nAChR alpha-subunit (Val188, Tyr190, Pro197, and Asp200) as contributing to binding. In this study, thermodynamic mutant cycle analysis was applied to these sets of residues to identify specific pairwise interactions. Amino acid variants of alpha-neurotoxin from N. mossambica mossambica at position 33 and of the nAChR at position 188 showed strong energetic couplings of 2-3 kcal/mol at both binding sites. Consistently smaller yet significant linkages of 1.6-2.1 kcal/mol were also observed between variants at position 27 on the toxin and position 188 on the receptor. Additionally, toxin residue 27 coupled to the receptor residues 190, 197, and 200 at the alphadelta binding site with observed coupling energies of 1.5-1.9 kcal/mol. No linkages were found between toxin residue Lys47 and the receptor residues studied here. These results provide direct evidence that the two conserved cationic residues Arg33 and Lys27, located on loop II of the toxin structure, are binding in close proximity to the alpha-subunit region between residues 188-200. The toxin residue Arg33 is closer to Val188, where it is likely stabilized by adjacent negative or aromatic residues on the receptor structure. Lys27 is positioned closer to Tyr190, Pro197, and Asp200, where it is likely stabilized through electrostatic interaction with Asp200 and/or cation/pi interactions with Tyr190.

Toxins selective for subunit interfaces as probes of nicotinic acetylcholine receptor structure

The pentameric structure of the nicotinic acetylcholine receptor with two of the five subunit interfaces serving as ligand binding sites offers an opportunity to distinguish features on the surfaces of the subunits and their ligand specificity characteristics. This problem has been approached through the study of assembly of subunits and binding characteristics of selective peptide toxins. The receptor, with its circular order of homologous subunits (alpha gamma alpha delta beta), assembles in only one arrangement, and through mutagenesis, the residues governing assembly can be ascertained. Selectivity of certain toxins is sufficient to readily distinguish between sites at the alpha gamma and alpha delta interfaces. By interchanging residues on the gamma and delta subunits, and ascertaining how they interact with the alpha-subunit, determinants forming the binding sites can be delineated. The alpha-conotoxins, which contain two disulfide loops and 12-14 amino acids, show a 10,000-fold preference for the alpha delta over the alpha gamma subunit interface with alpha epsilon falling between the two. The waglerins, as 22-24 amino acid peptides with a single core disulfide loop, show a 2000-fold preference for alpha epsilon over the alpha gamma and alpha delta interfaces. Finally, the 6700 Da short alpha-neurotoxin from N. mossambica mossambica shows a 10,000-fold preference for the alpha gamma and alpha delta interfaces over alpha epsilon. Selective mutagenesis enables one to also distinguish alpha-neurotoxin binding at the alpha gamma and alpha delta subunits. This information, when coupled with homology modeling of domains and site-directed residue modification, reveals important elements of receptor structure and conformation.

Role of the glycine triad in the ATP-binding site of cAMP-dependent protein kinase

A glycine-rich loop in the ATP-binding site is one of the most highly conserved sequence motifs in protein kinases. Each conserved glycine (Gly-50, Gly-52, and Gly-55) in the catalytic (C) subunit of cAMP-dependent protein kinase (cAPK) was replaced with Ser and/or Ala. Active mutant proteins were expressed in Escherichia coli, purified to apparent homogeneity, separated into phosphoisoforms, and characterized. Replacing Gly-55 had minimal effects on steady-state kinetic parameters, whereas replacement of either Gly-50 or Gly-52 had major effects on both Km and kcat values consistent with the prediction of the importance of the tip of the glycine-rich loop for catalysis. Substitution of Gly-50 caused a 5-8-fold reduction in Km (ATP), an 8-12-fold increase in Km (peptide), and a 3-5-fold drop in kcat. The Km (ATP) and Km (peptide) values of C(G52S) were increased 8- and 18-fold, respectively, and the kcat was decreased 6-fold. In contrast to catalytic efficiency, the ATPase rates of C(G50S) and C(G52S) were increased by more than an order of magnitude. The thermostability of each mutant was slightly increased. Unphosphorylated C(G52S) was characterized as well as several isoforms phosphorylated at a single site, Ser-338. All of these phosphorylation-defective mutants displayed a substantial decrease in both enzymatic activity and thermal stability that correlated with the missing phosphate at Thr-197. These results are correlated with the crystal structure, models of the respective mutant proteins, and conservation of the Glys within the protein kinase family.

A model of the nicotinic receptor extracellular domain based on sequence identity and residue location

We have modeled the extracellular domains of individual subunits (amino acids 31-200) in the nicotinic acetylcholine receptor using sequence homology with copper binding proteins of known crystal structure, plastocyanin and pseudoazurin, and data from recent site-specific mutagenesis, antibody mapping, and site-directed labelling studies. These data formed an initial model that was refined using molecular dynamics and mechanics as well as electrostatic and solvation energy calculations. The sequences between residues 31 and 164 in the alpha 1-subunit and corresponding residues in homologous receptor subunits show similarity with the core sequence of the cation binding site in plastocyanin and pseudoazurin, a region in the template proteins characterized by multiple hairpin loops. In addition to defining the subunit interfaces that comprise the site for agonist and competitive antagonist binding in more detail, the findings show that negatively charged residues cluster in domains arranged to diminish electrostatic free energy of the complex. Electrostatic factors also appear to distinguish the ligand binding interfaces, alpha gamma and alpha delta, from the other three interfaces on the pentameric receptor.

Solution structure of synthetic peptide inhibitor and substrate of cAMP-dependent protein kinase. A study by 2D H NMR and molecular dynamics

Peptides derived from the inhibitor of cAMP-dependent protein kinase. PKI, have been studied by 2D 1H NMR techniques. These include the inhibitor PKI(6-22), the substrate [Ala20-Ser21]PKI(5-24), and a phosphorylated form of the latter [Ala20-Ser21P]PKI(5-24). A homologous fold was found in the three peptides which consisted of an N-terminal segment in helical conformation to residue 13 and a C-terminal segment poorly defined conformationally. A parallel study was carried out by molecular dynamics (MD) for the inhibitor peptide PKI(5-24). The N-terminal helix, as observed in the crystal structure of the catalytic subunit-PKI(5-24) complex, was conserved in the MD simulations with the enzyme-free inhibitor. Similarly the Gly14-Gly17 turn was apparent in all MD structures, whereas the C-terminal region, residues 18-24, was directed towards the N-terminal helix in contrast to the extended conformation of this segment pointing away from the N-terminal helix in the crystal structure. This is primarily due to ionic interaction between Asp9 and Arg15. Indeed, a detailed analysis of the NOE contacts by NOESY at low temperature (2 degrees C) shows the occurrence of pH-dependent contacts with Phe10. We conclude that the binding of short inhibitors, such as PKI(5-24), to the enzyme involves a conformational rearrangement of the C-terminal region. The substrate [Ala20-Ser21]PKI(5-24) and the product [Ala20-Ser21P]PKI(5-24), give very similar structures with local rearrangements involving some of the side chains.

Synergistic binding of nucleotides and inhibitors to cAMP-dependent protein kinase examined by acrylodan fluorescence spectroscopy

We have engineered an acrylodan-modified derivative of the catalytic subunit of cyclic AMP-dependent protein kinase (cAPK) whose fluorescence emission signal has allowed the synergistic binding between nucleotides and physiological inhibitors of cAPK to be examined (Whitehouse, S., and Walsh, D. A. (1983) J. Biol. Chem. 258, 3682-3692). In the presence of the regulatory subunit, RI, the affinity of cAPK for adenosine, ADP, AMPPNP (adenosine 5'-(beta, gamma-imino)triphosphate), or ATP was 5-, 50-, 120-, and 15,000-fold enhanced, while in the presence of the heat-stable inhibitor protein of cAPK (PKI), there was a 3-, 20-, 33-, and 2000-fold enhancement in the binding of these nucleotides, respectively. A short inhibitor peptide, PKI-(14-22), enhanced the binding of ADP to the same degree as did full-length PKI (20-fold) but, in contrast, did not significantly enhance the binding of ATP or AMPPNP. The full binding synergism between PKI and either ATP (2000-fold) or AMPPNP (33-fold) to cAPK could, however, be mimicked by a longer peptide, PKI-(5-24), suggesting that the PKI NH2 terminus (residues 5-13) is most likely critical. Since this region is remote from the ATP gamma-phosphate, the binding synergism must arise through an extended network communication mechanism between the PKI NH2 terminus and the ATP binding site.