Critical amino acids responsible for conferring calcium channel characteristics are located on the surface and aroundβ-turn potentials of channel proteins

Calcium ion is thought to be one of the initial signals in the process of synaptic modification. Various reports have described that the critical amino acids responsible for determining calcium permeability of ion channels are glutamic acid, glutamine, arginine, and asparagine. By using a computational method (MacPROT) distinguishing transmembrane, globular, and surface sequences of proteins, the present work predicts that the critical amino acids exist within surface regions of the proteins. Furthermore, occurrence ofβ-turn probabilities can be predicted around these critical residues by the protein conformational prediction method of Chou and Fasman. The results suggest that the critical amino acids exist at hydrophilic spaces or canals of membranous channel proteins and that the redirection potential of the protein chain induced by the turn structures provides the conformational change requisite for the ion selectivity and gating (opening/closing) of the channels.     

Structural effects of amino acid substitutions on the P21 proteins: Evidence for a malignant conformation

The conformational effects of different amino acid substitutions for Gly at position 12 in theras-oncogene-encoded P21 proteins have been investigated using conformational energy calculations. Mutations that cause amino acid substitutions for Gly 12 result in a protein that produces malignant transformation of cells. It had previously been shown that substitution of Val, Lys, or Ser for Gly at position 12 results in a major conformational change, and that the preferred lowest energy structure for each of the substituted peptides is identical. It is now found that substitution for Gly 12 of other amino acids that have widely disparate helix-nucleating potentials and completely different side chains (Asp, Asn, Cys, Phe, Tle, Leu, and Ala) all produce this identical lowest energy conformation. This finding is consistent with the recent results of site-specific mutagenesis experiments showing that P21 proteins containing these amino acids at position 12 all promote malignant transformation of cells and suggests the existence of a malignancy-causing conformation for the P21 proteins.     

Comparison of the low energy conformations of an oncogenic and a non-oncogenic p21 protein, neither of which binds GTP or GDP

Oncogenic p21 protein, encoded by theras-oncogene, that causes malignant transformation of normal cells and many human tumors, is almost identical in sequence to its normal protooncogene-encoded counterpart protein, except for the substitution of arbitrary amino acids for the normally occurring amino acids at critical positions such as Gly 12 and Gin 61. Since p21 is normally activated by the binding of GTP in place of GDP, it has been postulated that oncogenic forms must retain bound GTP for prolonged time periods. However, two multiply substituted p21 proteins have been cloned, neither of which binds GDP or GTP. One of these mutant proteins with Val for Gly 10, Arg for Gly 12, and Thr for Ala 59 causes cell transformation, while the other, similar protein with Gly 10, Arg 12, Val for Gly 13 and Thr 59 does not transform cells. To define the critical conformational changes that occur in the p21 protein that cause it to become oncogenic, we have calculated the low energy conformations of the two multiply substituted mutant p21 proteins using a new adaptation of the electrostatically driven Monte Carlo (EDMC) technique, based on the program ECEPP. We have used this method to explore the conformational space available to both proteins and to compute the average structures for both using statistical mechanical averaging. Comparison of the average structures allows us to detect the major differences in conformation between the two proteins. Starting structures for each protein were calculated using the recently deposited x-ray crystal coordinates for the p21 protein, that was energy-refined using ECEPP, and then perturbed using the EDMC method to compute its average structure. The specific amino acid substitutions for both proteins were then generated into the lowest energy structure generated by this procedure, subjected to energy minimization and then to full EDMC perturbations. We find that both mutant proteins exhibit major differences in conformation in specific regions, viz., residues 35–47, 55–78, 81–93, 96–110, 115–126, and 123–134, compared with the EDMC-refined x-ray structure of the wild-type protein. These regions have been found to be the most flexible in the p21 protein bound to GDP from prior molecular dynamics calculations (Dykeset al., 1993). Comparison of the EDMC-average structure of the transforming mutant with that of the nontransforming mutant reveals major structural differences at residues 10–16, 32–40, and 60–68. These structural differences appear to be the ones that are critical in activation of the p21 protein. Analysis of the correlated motions of the different regions of the two mutant proteins reveals that changes in the conformation of regions in the carboxyl half of the protein are caused by changes in conformation around residues 10–16 and are transmitted by means of residues around Gln 61. The latter region therefore constitutes a molecular switch unit, in agreement with conclusions from prior work.On leave from the Department of Chemistry, University of Gdask, ul. Sobieskiego 18, 80-952 Gdask, Poland.     

Use of Chou-Fasman amino acid conformational parameters to analyze the organization of the genetic code and to construct protein genealogies

Summary Chou-Fasman parameters, measuring preferences of each amino acid for different conformational regions in proteins, were used to obtain an amino acid difference index of conformational parameter distance (CPD) values. CPD values were found to be significantly lower for amino acid exchanges representing in the genetic code transitions of purines, GA than for exchanges representing either transitions of pyrimidines, CU, or transversions of purines and pyrimidines. Inasmuch as the distribution of CPD values in these non GA exchanges resembles that obtained for amino acid pairs with double or triple base differences in their underlying codons, we conclude that the genetic code was not particularly designed to minimize effects of mutation on protein conformation. That natural selection minimizes these changes, however, was shown by tabulating results obtained by the maximum parsimony method for eight protein genealogies with a total occurrence of 4574 base substitutions. At the beginning position of the codons GA transitions were in very great excess over other base substitutions, and, conversely, CU transitions were deficient. At the middle position of the codons only fast evolving proteins showed an excess of GA transitions, as though selection mainly preserved conformation in these proteins while weeding out mutations affecting chemical properties of functional sites in slow evolving proteins. In both fast and slow evolving proteins the net direction of transitions and transversions was found to be from G beginning codons to non-G beginning codons resulting in more commonly occurring amino acids, especially alanine with its generalized conformational properties, being replaced at suitable sites by amino acids with more specialized conformational and chemical properties. Historical circumstances pertaining to the origin of the genetic code and the nature of primordial proteins could account for such directional changes leading to increases in the functional density of proteins.In order to further explore the course of protein evolution, a modified parsimony algorithm was developed for constructing protein genealogies on the basis of minimum CPD length. The algorithm's ability to judge with finer discrimination that in protein evolution certain pathways of amino acid substitution should occur more readily than others was considered a potential advantage over strict maximum parsimony. In developing this CPD algorithm, the path of minimum CPD length through intermediate amino acids allowed by the genetic code for each pair of amino acids was determined. It was found that amino acid exchanges representing two base changes have a considerably lower average CPD value per base substitution than the amino acid exchanges representing single base changes. Amino acid exchanges representing three base changes have yet a further marked reduction in CPD per base change. This shows how extreme constraining effects of stabilizing selection can be circumvented, for by way of intermediate amino acids almost any amino acid can ultimately be substituted for another without damage to an evolving protein's conformation during the process.     

Measuring Plasma Membrane Protein Endocytic Rates by Reversible Biotinylation

Plasma membrane proteins are a large, diverse group of proteins comprised of receptors, ion channels, transporters and pumps. Activity of these proteins is responsible for a variety of key cellular events, including nutrient delivery, cellular excitability, and chemical signaling. Many plasma membrane proteins are dynamically regulated by endocytic trafficking, which modulates protein function by altering protein surface expression. The mechanisms that facilitate protein endocytosis are complex and are not fully understood for many membrane proteins. In order to fully understand the mechanisms that control the endocytic trafficking of a given protein, it is critical that the protein s endocytic rate be precisely measured. For many receptors, direct endocytic rate measurements are frequently achieved utilizing labeled receptor ligands. However, for many classes of membrane proteins, such as transporters, pumps and ion channels, there is no convenient ligand that can be used to measure the endocytic rate. In the present report, we describe a reversible biotinylation method that we employ to measure the dopamine transporter (DAT) endocytic rate. This method provides a straightforward approach to measuring internalization rates, and can be easily employed for trafficking studies of most membrane proteins.Download video file.^{(127M, mp4)}     

Presynaptic Calmodulin targets: lessons from structural proteomics

Introduction: Calmodulin (CaM) is a highly conserved Ca²⁺-binding protein that is exceptionally abundant in the brain. In the presynaptic compartment of neurons, CaM transduces changes in Ca²⁺ concentration into the regulation of synaptic transmission dynamics.

Areas covered: We review selected literature including published CaM interactor screens and outline established and candidate presynaptic CaM targets. We present a workflow of biochemical and structural proteomic methods that were used to identify and characterize the interactions between CaM and Munc13 proteins. Finally, we outline the potential of ion mobility-mass spectrometry (IM-MS) for conformational screening and of protein-protein cross-linking for the structural characterization of CaM complexes.

Expert commentary: Cross-linking/MS and native MS can be applied with considerable throughput to protein mixtures under near-physiological conditions, and thus effectively complement high-resolution structural biology techniques. Experimental distance constraints are applicable best when obtained by combining different cross-linking strategies, i.e. by using cross-linkers with different spacer length and reactivity, and by using the incorporation of unnatural photo-reactive amino acids. Insights from structural proteomics can be used to generate CaM-insensitive mutants of CaM targets for functional studies in vitro or ideally in vivo.     

Lytic activity localized to membrane-spanning region of ?X174 E protein

Summary Lytic activity of the X174E (lysis) protein had previously been localized to the amino terminal 51 amino acids (a.a.) of the molecule (Blasi and Lubitz 1985). This E gene lytic activity has here been further localized to the amino terminal 29 a.a., a region of the protein which is thought to just span the cell membrane (Young and Young 1982). X174 E gene fusions to both the lacZ gene and the chloramphenicol acetyl transferase (CAT) gene resulted in fusion proteins with lytic activity. Fusion to a third protein, trpE, did not result in lytic activity. These results support a model of oligomerization of the X174 E protein for lytic activity since both -galactosidase and CAT exist as tetramers in their native state. A difference in the composition of the charged amino acids at the cytoplasmic boundary between the various fusion proteins could also account for these results, since these amino acids may play a role in proper anchoring of the E protein in the cell membrane. In a spontaneous E gene mutant, which introduces a proline residue at position 9 of the E protein, lytic activity of the E protein was decreased, but not abolished. The presence of the helix-breaking proline at this position may interfere with insertion of the lysis protein into the cell membrane, leading to the decreased functional activity of the protein.     

The regularities of the changes of amino acid physico-chemical properties within the genetic code

Summary All the codons of the genetic code can be arranged into the closed one-step mutation ring, containing three periods of the same sequence of mutations (2,3,3,3,1,3,3,3,1,3,3,3,1,3,3,3,2,3,3,3). The codons of Gly play a role of the connecting element between the end of the third, and the beginning of the first period of the genetic code. The reactivity of amino acids, expressed by the reaction rates of aminolysis reaction of N-hydroxysuccinimide esters of protected amino acids with p-anisidine, changes periodically with the respect to the mutation periods of the genetic code. Chou-Fasman P as well as P conformational parameters of amino acids, and also the compositional frequencies of amino acids in proteins, demonstrate the pseudosymmetry pattern with respect to the center of one-step mutation ring, which is situated between Thr ACY and ACR codons.     

A conservation and rigidity based method for detecting critical protein residues

Background

Certain amino acids in proteins play a critical role in determining their structural stability and function. Examples include flexible regions such as hinges which allow domain motion, and highly conserved residues on functional interfaces which allow interactions with other proteins. Detecting these regions can aid in the analysis and simulation of protein rigidity and conformational changes, and helps characterizing protein binding and docking. We present an analysis of critical residues in proteins using a combination of two complementary techniques. One method performs in-silico mutations and analyzes the protein's rigidity to infer the role of a point substitution to Glycine or Alanine. The other method uses evolutionary conservation to find functional interfaces in proteins.

Results

We applied the two methods to a dataset of proteins, including biomolecules with experimentally known critical residues as determined by the free energy of unfolding. Our results show that the combination of the two methods can detect the vast majority of critical residues in tested proteins.

Conclusions

Our results show that the combination of the two methods has the potential to detect more information than each method separately. Future work will provide a confidence level for the criticalness of a residue to improve the accuracy of our method and eliminate false positives. Once the combined methods are integrated into one scoring function, it can be applied to other domains such as estimating functional interfaces.

     

Ion channels formed by amphipathic helical peptides

Channel forming peptides (CFPs) are amphipathic peptides, of length ca. 20 residues, which adopt an -helical conformation in the presence of lipid bilayers and form ion channels with electrophysiological properties comparable to those of ion channel proteins. We have modelled CFP channels as bundles of parallel trans-bilayer helices surrounding a central ion-permeable pore. Ion-channel interactions have been explored via accessible surface area calculations, and via evaluation of changes in van der Waals and electrostatic energies as a K⁺ ion is translated along the length of the pore. Two CFPs have been modelled: (a) zervamicin-A1-16, a synthetic apolar peptaibol related to alamethicin, and (b) -toxin from Staphylococcus aureus. Both of these CFPs have previously been shown to form ion channels in planar lipid bilayers, and have been shown to have predominantly helical conformations. Zervamicin-A1-16 channels were modelled as bundles of 4 to 8 parallel helices. Two related helix bundle geometries were explored. K⁺channel interactions have been shown to involve exposed backbone carbonyl oxygen atoms. -Toxin channels were modelled as bundles of 6 parallel helices. Residues Q3, D11 and D18 generate favourable K⁺-channel interactions. Rotation of W15 about its C-C bond has been shown to be capable of occluding the central pore, and is discussed as a possible model for sidechain conformational changes in relation to ion channel gating.