|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
J. Biol. Chem., Vol. 279, Issue 24, 25364-25373, June 11, 2004
Functional Characterization of a Second Porin Isoform in Drosophila melanogasterDMPORIN2 FORMS VOLTAGE-INDEPENDENT CATION-SELECTIVE PORES*![]() ![]() ![]() ![]() ![]() ![]() ![]()
From the
Received for publication, September 24, 2003 , and in revised form, March 29, 2004.
Mitochondrial porins or voltage-dependent anion-selective channels are channel-forming proteins mainly found in the mitochondrial outer membrane. Genome sequencing of the fruit fly Drosophila melanogaster revealed the presence of three additional porin-like genes. No functional information was available for the different gene products. In this work we have studied the function of the gene product closest to the known Porin gene (CG17137 coding for DmPorin2). Its coding sequence was expressed in Escherichia coli. The recombinant DmPorin2 protein is able to form channels similar to those formed by DmPorin1 reconstituted in artificial membranes. Furthermore, DmPorin2 is clearly voltage-independent and cation-selective, whereas its counterpart isoform 1 is voltage-dependent and anion-selective. Sequence comparison of the two porin isoforms indicates the exchange of four lysines in DmPorin1 for four glutamic acids in DmPorin2. We have mutated two of them (Glu-66 and Glu-163) to lysines to investigate their role in the functional features of the pore. The mutants E163K and E66K/E163K are endowed with an almost full inversion of the ion selectivity. Both single mutations partially restore the voltage dependence of the pore. We found that an additional effect with the double mutant E66K/E163K was the restoration of voltage dependence. Protein structure predictions highlight a 16 -strand pattern, typical for porins. In a three-dimensional model of DmPorin2, Glu-66 and Glu-163 are close to the rim of the channel, on two opposite sides. DmPorin2 is expressed in all the fly tissues and in all the developmental stages tested. Our main conclusions are as follows. 1) The CG17137 gene may express a porin with a functional role in D. melanogaster. 2) We have identified two amino acids of major relevance for the voltage dependence of the porin pore.
Mitochondrial porins or voltage-dependent anion-selective channels (VDACs)1 are pore-forming proteins mainly found in the mitochondrial outer membrane, where they are able to form large hydrophilic channels (13). Although many features of this protein family have been known for more than 20 years, the atomic structure is still lacking, and only models have been computed according to secondary structure predictions (47). All models indicate that the transmembrane segments typical for these channels are formed by antiparallel, amphipathic -strands, similar to the situation in bacterial porins (8). Porin or VDAC is the major pore-forming protein of the mitochondrial outer membrane and has a unique role in energetic metabolism since all metabolites, including ADP and ATP, have to pass through the porin channels (9, 10). It is the receptor where kinases bind to the mitochondrial surface, which is also important for the metabolic energy channeling in the cell (11). The role of mitochondrial porin isoform 1 or VDAC1 in apoptosis was also thoroughly investigated in the last few years (1216). Its participation in the permeability transition pore has been reported (1719). A recent report (20) focused on the VDAC2 isoform as an important player in the cell death pathway.
In the last years it has been demonstrated that the chromosomes of eukaryotic cells contain several genes coding for porins. In Saccharomyces cerevisiae the deletion of the porin gene is not lethal for the cells unless they were forced to live under strictly mitochondrially driven, aerobic conditions (21). This observation led to the discovery of a second gene, coding for a porin analogue and termed VDAC2 by sequence homology (22). Although the yeast VDAC2 is unable to form channels in vitro, its overexpression complements the inability to grow on non-fermentable carbon sources, the cellular defect induced by the VDAC1 deletion (23). The genomes of humans, rats, and mice contain three different porin genes coding for distinctly expressed isoforms (2325). Human and mouse VDAC1 and VDAC2 isoforms show channel forming activity in vitro and can complement the yeast VDAC1 deficiency, whereas VDAC3 can only partially complement this defect (25). An interesting picture of the porin genes in a higher eukaryotic organism has been revealed in the fruit fly Drosophila melanogaster. The porin gene, which codes for DmPorin1, was cloned, sequenced, and mapped at 32B3-4 on the second chromosome. It has an identity of about 60% compared with those of the porin isoforms in mammals (26, 27). The organization of the region adjacent to the porin locus was elucidated after completion of the fruit fly genome. Most surprising, it revealed the close proximity of three additional genes with significant homology to the porin gene called CG17137, CG17139, and CG17140 (28). These three genes were predicted in silico by alignment of the genomic sequences with the EST sequences in the Berkeley Drosophila Genome Project collection. They are on the same strand as the porin gene, and they have a similar organization with two introns interrupting the coding sequence (28). The 5'-untranslated region is present only in the porin gene. The nucleotide distances separating the four porin-like genes are extremely reduced. EST sequences have been reported for CG17137 and for CG17139 and CG17140, albeit considerably less frequent than for the porin gene. In a most recent update of the Flybase Data Bank (FLYBASE, FBgn0051722; LocusID, 34497), the putative genes CG17139 and CG17140 have been considered as a single transcription unit termed CG31722, since no 3'-EST is available for the gene CG17139. The expression of CG17137 and CG31722 is established at the transcriptional level, since the corresponding ESTs have been identified. Alignment of the deduced protein sequences shows sequence conservation in pairs. The protein sequence deduced from CG17137 is 42% identical to D. melanogaster porin isoform 1, but it is only 26% identical to the deduced amino acid sequence of CG17140 or CG31722 (28). In this paper we have cloned the coding sequence corresponding to CG17137 and expressed it in Escherichia coli. The expressed protein was purified and refolded. It showed pore forming activity with features resembling but not overlapping those of the isoform 1. In addition, the functional characterization revealed new features that could be of general interest to the understanding of the function of this protein. Site-directed mutagenesis of amino acids that are different in both porin sequences allowed interesting insights into the role of these amino acids in voltage dependence of eukaryotic porin.
PCR Amplification of D. melanogaster Porin-coding SequencesThe CG17137 (FLYBASE, FBgn0032308; LocusID, 34499) coding sequence was amplified from a D. melanogaster male adult cDNA library in gt10 using the primers DmPorin2for and DmPorin2rev (Table I). The primer DmPorin2for introduces a BamHI restriction site and the ATG start codon into the PCR product, whereas the primer DmPorin2rev introduces a HindIII restriction site and the TAA stop codon. 1.5 µl of phage lysate from the library were added to 28.5 µl of H2O, warmed at 70 °C for 5 min, and incubated on ice. The PCR was performed according to standard protocol in a volume of 50 µl, and the final concentration of dNTPs was 0.2 mM each; the MgCl2 final concentration was 1.5 mM, and the final concentration of the primers was 0.5 µM each. The PCR product was cloned into the linear TOPO-TA vector (Invitrogen) and sequenced (GenBankTM accession number AJ580916
[GenBank]
).
The control coding sequence of the DmPorin1 (GenBankTM accession number X92408 [GenBank] ) was amplified from a cDNA clone isolated in our laboratory (26). The upstream primer DmPorin1for introduced the BamHI restriction site and the ATG start codon into the PCR product; similarly, the downstream primer DmPorin1rev introduced a HindIII site together with the TAA stop codon. 10 ng of plasmid were used as template, and the reaction was performed in a volume of 50 µl. All the other conditions were the same as described above. Cloning was in the TOPO-TA vector (Invitrogen). PCR Amplification of CG17137 from D. melanogaster Developmental Stage cDNA LibrariesThe same PCR was performed to check whether the CG17137 gene was transcribed in other developmental stages of the fly. Eight different cDNA libraries were used as templates: embryo 03 h; embryo 1224 h; embryo III instar larval stages; early III instar larval stage; late III instar larval stage; pupal instar stage; and late pupal instar stage. The protocol for PCR and the amplification parameters were the same as described above. PCR was also designed on Rp49, a ribosomal protein, and was used to amplify the underlying sequence in the various libraries to standardize the template amount (not shown). Cloning in the pQE30 VectorThe fragments released from the TOPO TA vector after digestion with BamHI and HindIII were ligated in the expression vector pQE30 (Qiagen) to obtain the plasmid pQE30DmPorin2. pQE30DmPorin2 also codes for an additional six histidines (His tag) upstream of the polyclonal site containing the cDNA coding for DmPorin2. XL10 gold ultracompetent cells (Stratagene) were transformed using 110 ng of pQE30DmPorin2, following the protocol provided by the manufacturer. The cells were spread on an ampicillin- and kanamycin-containing agar plate. Positive clones were picked, and the plasmid DNA was extracted from a single colony, digested, and checked on a 1% agarose gel for the right fragment size and was finally sequenced. Site-directed MutagenesisThe QuickChangeTM site-directed mutagenesis kit (Stratagene) was used to change two glutamic acids (Glu-66 and Glu-163) into two lysines in the DmPorin2 sequence by PCR. The mutant plasmid pQE30DmPorin2e66k was obtained by using the DmPorin2 cDNA cloned in the plasmid pQE30 as template and the primers E66Kfor and E66Krev (Table I). Similarly, the mutant plasmid pQE30DmPorin2e163k was obtained by mutagenesis of pQE30DmPorin2 with the primers E163Kfor and E163Krev (Table I). A double mutant, harboring both the E66K and E163K mutations, has been created using pQE30DmPorin2e66k as template in a further PCR with the primers E163Kfor and E163Krev. The PCR was performed according to the protocol of the manufacturer with Pfu Turbo DNA polymerase (2.5 units) and 125 ng of the primers. The amount of template used was 10 ng. The PCR product was digested using DpnI (10 units) at 37 °C for 1h.1 µl of the digested DNA was transformed in competent E. coli M15 cells. The DNA purified from single colonies was sequenced to check the products.
Heterologous Expression of DmPorin GenesThe expression of the different genes was performed in E. coli strains devoid of the major porin genes (genotype: BL21(DE3), Renaturation of the ProteinThe proteins purified under denaturing conditions (8 M) have lost their structure under these conditions and needed to be refolded to get channel-forming activity in lipid bilayer experiments (29). The eluted proteins were dialyzed at 4 °C overnight against a solution containing 1% Genapol X-080 (Fluka) supplemented with 0.5% (w/v) cholesterol. The proteins were kept at 4 °C before use in the black lipid bilayer measurements. Lipid Bilayer ExperimentsThe method used for the black lipid bilayer experiments has been described previously (30). Membranes were formed from a 1% (w/v) solution of diphytanoylphosphatidylcholine (Avanti Polar Lipids, Alabaster AL) in n-decane by painting onto a circular hole (surface area about 0.4 mm2) separating the two compartments of a Teflon cell. For standard single channel conductance experiments, the Teflon chamber was filled with an unbuffered salt solution. The voltage across the membrane was applied through silver/silver chloride electrodes (with salt bridges) inserted into the aqueous compartments on both sides of the membrane. The membrane current was measured with a current amplifier (Keithley 427 current amplifier). The amplified signal was monitored with a storage oscilloscope and recorded on a strip chart recorder. Zero current potential measurements were performed using a high impedance electrometer (Keithley 617 current amplifier) as described previously (31). Voltage Dependence StudiesThe refolded recombinant protein was added at a concentration of 500 ng/ml to both sides of a black diphytanoylphosphatidylcholine/n-decane membrane. After about 20 min, the reconstitution of porin channels into the membrane reached equilibrium. Then different potentials were applied to both sides of the membrane starting with ±10 mV (data not shown). These experiments were repeated with ±20 to ±80 mV in steps of ±10 mV. The experiments were analyzed in the following way. The membrane conductance (G) as a function of voltage, Vm, was measured when the opening and closing of channels reached an equilibrium, i.e. after the exponential decay of the membrane current following the voltage step Vm. G was divided by the initial value of the conductance (G0, which was a linear function of the voltage) obtained immediately after the onset of the voltage. Ion Selectivity DeterminationDiphytanoylphosphatidylcholine membranes were formed in 50 mM KCl, and concentrated porin solutions were added to the aqueous phase when the membranes were in the black state. After incorporation of 10100 channels into a membrane, salt gradients were established by addition of small amounts of concentrated KCl solution to one side of the membrane. Analysis of the zero current membrane potential was performed using the Goldman-Hodgkin-Katz equation (31). Three-dimensional Modeling of the DmPorin2 Amino Acid SequenceThe three-dimensional model of DmPorin2 was computed by using a procedure similar to that adopted previously for porin (isoform 1 or VDAC1) (7). In summary, the protein topography starting from the sequence was predicted with a neural network-based method and a hidden Markov model-based method (32, 33).
Both predictors suggested 16 transmembrane strands along the protein sequence, with the same reliability index. An The target was then modeled on templates taken from the Protein Data Bank (codes 2OMF [PDB] , 1PRN [PDB] , and 2POR [PDB] ) with the same number of strands in the barrel. MODELLER 6.2 (34) was used to perform essentially a building by homology procedure. The N-terminal portion was not considered, since the expert-driven alignment between the target and the templates is strictly based on the transmembrane regions, constraining gaps in the loop regions, external to the membrane phase. Small differences in the alignment led to slightly different models; the model with the best PROCHECK (35) score (minimum percentage of residues in disallowed conformation) was considered as the most reliable model.
Expression and Purification of DmPorin1 and DmPorin2 Fig. 1A shows an agarose gel of the PCR amplification from the D. melanogaster male adult cDNA library using the primers DmPorin2for and reverse specific for CG17137 (DmPorin2) and the primers DmPorin1for and reverse specific for porin (DmPorin1). Both cDNAs were cloned in a TA vector and then in the pQE30 expression vector to yield the plasmids pQE30DmPorin1 and pQE30DmPorin2. Both plasmids were transferred into an E. coli BL21 strain, which lacks all major outer membrane proteins. The expressed proteins were purified from inclusion bodies. They had the expected apparent molecular masses (Fig. 1B). The expression level in bacteria was around 30 µg of protein/ml of bacterial culture for both proteins, and the purified protein was free of bacterial contaminants as judged by SDS-PAGE.
PCR Amplification of the Gene for DmPorin2 in Different RNA-derived cDNA LibrariesA similar PCR amplification was performed on a set of D. melanogaster libraries derived by poly(A)+ RNA purified at different developmental stages of the fly (Fig. 2). The coding sequence delimited by the PCR primers (DmPorin2for and rev) was specifically amplified in all the libraries tested and had a similar level. This result demonstrates that the gene of DmPorin2 was transcribed in all developmental stages of D. melanogaster.
Channel Forming Activity and Voltage Dependence of the Recombinant PorinsThe electrophysiological properties of the purified and refolded DmPorin1 were studied in black lipid bilayer membranes and compared with those of the native porin purified from D. melanogaster mitochondria, which have been reported previously (36). Single channel experiments performed with small amounts of refolded recombinant protein (110 ng/ml) in 1 M KCl showed current increases in discrete steps with a predominant single channel conductance of about 4.5 nS (Fig. 3), similar to native porin isolated from D. melanogaster mitochondria. This result also indicated that the His tag attached to the recombinant protein did not interfere with the channel conductance.
In the next step the voltage dependence of DmPorin1 was analyzed. Voltage dependence is a typical feature of mitochondrial porins (3). Starting with application of about 2030 mV, the membrane current decreased in an exponential fashion for both positive and negative potentials. This result indicated symmetrical response to the applied voltage. The data of Fig. 3D show the plot of the relative conductance G/G0 as a function of the applied voltage, i.e. it corresponds to the voltage dependence of DmPorin1. It is noteworthy that the voltage dependence of the recombinant porin was identical to that measured previously with the native mitochondrial porin of D. melanogaster (36). These results indicate that recombinant DmPorin1 was successfully re-natured and that its properties were undistinguishable from those of native porin. Similar experiments were performed with the refolded recombinant DmPorin2. It also formed channels in lipid bilayer membranes with a long lifetime at low voltage. DmPorin2 also formed preferentially channels with a single channel conductance of 4.5 nS in 1 M KCl (Fig. 3A). However, the channel forming activity was much smaller than that of DmPorin1 under otherwise identical conditions. At least a 10 times higher concentration of DmPorin2 was needed to obtain the same number of reconstituted channels as compared with DmPorin1. It is unknown whether this is a normal feature of DmPorin2 or is because of an incomplete refolding process. The results of single channel experiments suggest that both porins form channels with very similar properties (Fig. 3). The voltage dependence of DmPorin2 was also investigated and compared with that of DmPorin1. Most surprising, DmPorin2 did not display any voltage dependence as clearly evidenced by Fig. 3, B and D. The single channel conductance of DmPorin1 and DmPorin2 in different salt solutions was also studied. The recombinant DmPorin1 and DmPorin2 proteins formed channels that had the same conductance only in 1 M KCl. However, in other salt solutions with different mobility of the ions the conductance differed. When the salt was composed of the highly mobile Cl- anion and the less mobile Li+ cation, the conductance of DmPorin1 was about twice that recorded with DmPorin2. The opposite was the case when a less mobile anion (CH3COO-) was combined with a mobile cation (K+). This influence on the conductance could be explained by the effect of charges influencing ion transport through the channel, i.e. DmPorin1 and DmPorin2 could exhibit different ion selectivity (Fig. 3C). This possibility was studied by zero current membrane potential measurements, which allow the calculation of the permeability ratio Pcation divided by Panion (Pc/Pa) for channels in multichannel experiments. In all experiments with DmPorin1, the more diluted side of the membrane became negative, which indicated preferential movement of anions through the channel. The zero current membrane potential for a 4-fold gradient of KCl was about -2.2 mV (average of three experiments) on the more diluted side (Table II). The results were in good agreement with data derived previously from native porin isolated from mitochondria (36). Similar experiments were also performed with recombinant DmPorin2. This channel was found to be cation-selective because the more dilute side had a positive potential of about 18 mV for an 8-fold KCl gradient, and the permeability coefficient Pcation divided by Panion was 3.8. This result indicated that the two isoforms of D. melanogaster porins exhibit different ion selectivity.
Electrophysiological Properties of Three DmPorin2 MutantsThe data from the reconstitution experiments clearly demonstrated a different functional behavior of the two porins. Analysis of the amino acid composition and sequence comparison of the two proteins supported this view. DmPorin1 has a smaller number of acidic amino acids (in particular, 11 glutamic acids versus 22 in DmPorin2) and a similar number of basic amino acids (25 lysines and 4 arginines against 21 and 7, respectively, in DmPorin2). The alignment of DmPorin1 and DmPorin2 showed that 10 lysines are in a conserved position, but four lysines in DmPorin1 are exchanged by glutamic acids in DmPorin2 (Fig. 4). In particular lysines 135 and 176 (sequence position of DmPorin2) are not well conserved among porins from different sources (see Fig. 4) and were excluded from our study. We thus selected glutamic acid 66 and 163 in DmPorin2, corresponding to Lys-65 and Lys-162 in DmPorin1, for site-directed mutagenesis to study the role of these amino acids in voltage dependence and ion selectivity.
To check whether Glu-66 and Glu-163 may influence the ion selectivity and the voltage dependence of DmPorin2, they were replaced by lysines through site-directed mutagenesis. The single mutants E66K and E163K were constructed together with the double mutant E66K/E163K. The electrophysiological properties of the three mutant proteins were tested in lipid bilayer experiments. The mutations resulted in dramatic changes of the functional behavior of the protein. The mutant E66K had a higher channel forming activity in low salt solutions possibly because the ionic strength interfered with the protein stability; thus, decreasing the ionic strength from 1 M to 150 mM or 50 mM KCl resulted in well defined channels (Fig. 5). In 1 M KCl it behaved in a manner similar to DmPorin2. The most frequent single channel conductance was 4 nS; nevertheless, the conductance was very disturbed, and the incorporation events were scarce.
Channel forming activity of E163K was lower than that of DmPorin2, and the single channels showed higher current noise even in low ionic strength buffer (data not shown). The most frequent single channel conductance in 1 M KCl was around 2.53 nS (Fig. 5). The smaller channel forming activity of the mutant as compared with DmPorin2 may be caused by the exchange of Glu-163 by lysine, which seems to interfere with the stability of the protein. Similar results were obtained for the double mutant, E66K/E163K, which also showed high current noise in the reconstitution experiments. The most frequent single channel conductance was 2.53 nS in 1 M KCl (Fig. 5). The behavior of this protein could be explained if the role of both glutamic acids is taken into consideration. Glutamic acid 66 mutation did not alter channel forming activity, but the substitution of glutamate 163 by lysine did. The double mutant showed the characteristic of the two single amino acid replacements and in particular that of the E163K mutant. The Replacement of Glutamic Acids by Lysines Confers Voltage Dependence to rDmPorin2DmPorin1 has a net positive charge and is voltage-dependent, whereas DmPorin2 has a net negative charge and is voltage-independent. Since it has been proposed that the voltage sensor of VDAC is formed by a cluster of positively charged groups (37), it is possible that the amino acids Lys-66 and Lys-163 of DmPorin1 may be involved in the voltage sensor. This means that the presence of glutamic acid residues at the same position in DmPorin2 can explain, at least in part, the absence of any voltage dependence in this protein. To check this possibility, the voltage dependence of DmPorin2 mutants was investigated. The ratio G/G0 calculated for DmPorin2 E66K was 0.89 at a voltage of ± 50 mV and 0.8 at a voltage of ± 80 mV, which suggests that this mutant exhibited some voltage dependence. Similar results were obtained with the other mutants, which were also slightly voltage-dependent (see Table II and Fig. 5). It is noteworthy that the double mutant E66K/E163K exhibited the highest voltage dependence (G/G0 was 0.61 at ± 50 mV and 0.58 at ± 80 mV), modifying the voltage independence shown for DmPorin2. Besides the voltage dependence we also studied the ion selectivity of the mutants. This was done to check whether the amino acids Glu-66 and Glu-163 were also involved in channel selectivity. The results are also included in Table II. The ratio Pc/Pa derived for DmPorin2 E66K was similar (indeed even higher) than that of DmPorin2, which means that this protein also formed cation-selective channels. The ratio Pc/Pa derived for the E163K and the E66K/E163K mutants in KCl was about 1.2 in both cases, i.e. the channels formed by these proteins were still slightly cation-selective although much less than DmPorin2 (Pc/Pa = 3.8). These results suggest that the mutation E66K has no effect on ion selectivity. The mutation E163K had a higher impact on the selectivity of the channels, but they were still cation-selective, which means that other charged amino acids also have an influence on the ion selectivity of these channels. Predicting the Structure of the DmPorin2 PoreThe DmPorin2 coding sequence (GenBankTM accession number AJ580916 [GenBank] ) obtained by amplification of an adult cDNA library was identical to the sequence deduced from the genome. In comparison with the other fruit fly porins, DmPorin2 contains an additional 10 amino acids (EQQGGDVVVN) toward the C terminus (Fig. 4).
The sequence of DmPorin2 was modeled by a refined neural network tool to predict a secondary structure and, in turn, a three-dimensional model of the protein (Fig. 6). The neural network software indicates the presence of 16 amphipathic
Porin/VDAC is a multifunctional protein involved in several cellular processes, but its biological significance has not been fully characterized (3842). In particular, the functions of the various VDACs or eukaryotic porins and their relationships are still open questions. In recent years we used the D. melanogaster to elucidate some of these questions in a relatively simple eukaryotic model organism (2628, 36). More porin genes have been found in the fly after the completion of the fly genome project (28). No information is available about the possible expression of the encoded proteins or the specific function of the gene products. Very recently Komarov et al. (43) have reported on the physiological properties of D. melanogaster porin-like proteins.
In the present work we have shown for the first time that the transcript corresponding to the gene termed CG17137 (28) is widely present in the various developmental stages of the fly. We have thus amplified its coding sequence from a cDNA library and expressed it in E. coli. After refolding in the presence of sterol, lipids, and detergents (29, 44), this protein was used for reconstitution experiments in black lipid membrane and showed pore forming activity. The formation of channels by DmPorin2 exhibited features typical for eukaryotic porins. Secondary structure predictions suggest the presence of amphipathic A detailed analysis of the channel formation of DmPorin2 revealed interesting differences between the homologous DmPorin1 protein. DmPorin2 is cation-selective instead of anion-selective protein, and it is voltage-independent. The charged amino acid residues in or near the channel wall are likely responsible for the weak ion selectivity of the channels formed by the eukaryotic porins (46). It is clear from the amino acid composition that DmPorin2 is a more acidic protein than DmPorin1. We assumed that exchanges of oppositely charged amino acids could cause the change of electrophysiological properties of the protein. We selected Glu-66 and Glu-163 in DmPorin2, corresponding to Lys-65 and Lys-162 in DmPorin1, for a site-directed mutagenesis study by assuming these amino acids were relevant for voltage dependence and ion selectivity. The E66K mutation did not change the preference for cations of the protein, whereas E163K decreased the cation selectivity. Neither E163K nor the double mutant E66K/E163K were able to completely reverse the ion selectivity to anion-selective as in DmPorin1. It clearly has to be concluded that only the position corresponding to Glu-163 is relevant in the ion selectivity. Voltage dependence of the mutants was studied next. As described above, DmPorin2 did not show any voltage dependence, which is a typical feature of these proteins, consequently also called VDAC. Most surprising, both single mutations E66K and E163K added a slight voltage dependence to DmPorin2. The double mutant E66K/E163K showed a higher voltage dependence, comparable in size to the voltage dependence measured for recombinant DmPorin1. It is tempting to speculate that the effect of the charge changes at positions 66 and 163 was additive.
For additional hints regarding the reason for such evident influence exerted by these mutations on the electrophysiological properties of the DmPorin2, we used predictive modeling (7) to speculate about the position of residues Glu-66 and Glu-163 in the protein. Most interesting, in our model both residues 66 and 163 are located close to the predicted rim of the pore, on opposite sides. In particular, residue 66 is located in a tight turn between two
We examined the literature to compare these results with other predictive models reported by other groups. In a 12 Voltage dependence is a typical feature of eukaryotic porins; it has the effect of partially closing the pore. Even though it has been detected only in reconstituted systems, it is thought to be of metabolic significance. Metabolites such as succinate, malate, and ATP, normally able to diffuse through the pore, should be delayed in crossing the outer membrane after the voltage-dependent switching of porin to closed states (9, 10).
Thomas et al. (37) proposed that a positively charged domain in the pore-forming protein, called voltage-sensor, distributed over a large region of the protein should be responsible for sensing voltage changes across the membrane and thus for triggering the partial closure of the pore. If this is correct, increasing the positive charge in the voltage sensor should increase the voltage dependence as it has indeed been demonstrated by the same authors with the substitution of lysine for glutamate at the end of a putative transmembrane
In a very recent report Komarov et al. (43) expressed in In conclusion, DmPorin2 shares with DmPorin1 the ability to form channels with a relatively large diameter. From electrophysiological studies we can conclude that they should have a similar function in the cell. This is confirmed by the pattern of transcription and by the structural predictions. The different although slight ionic selectivity might be related to a subtle preference for metabolites. More important could be the difference in voltage dependence. A porin not responding to changes in the voltage across the membrane could give rise to a non-regulated channel, always open to hydrophilic permeation. A similar channel-forming protein might thus be useful in emergency situations to rescue damages or impairment to the main porin DmPorin1. Ultimately, we cannot exclude the possibility that the second isoform might be involved in other different functions not yet known. For example, a recent report (20) suggests a primary function of VDAC2 isoform in mammals as a specific inhibitor of BAK-dependent mitochondrial apoptosis. We did not characterize the other predicted genes CG17139 and CG17140 from D. melanogaster because there is not a general consensus about their relevance. The compilators of the Berkeley Drosophila Genome Project have recently joined these two putative genes in a single one (CG31722). In comparison with the other predicted porin-like genes, DmPorin2 (or CG17137) has a higher probability to be a true gene, since it also contains a putative promoter and other untranslated regions. In addition, there is a large number of reported EST sequences related to it in the data bases, and its sequence does not show the unexpected 60-amino acid long N-terminal addition predicted for the putative products of genes CG17139 and CG17140 (28). This means that the CG17137 has functional features making them a biologically relevant new member of the Porin/VDAC gene family. The position of the porin protein in the eukaryotic cell and its involvement in the mitochondrial energy metabolism and in other cellular pathways such as programmed cell death may have raised the need for additional isoforms able to compensate for one another. This could have originated by modern gene duplications in D. melanogaster. Today porin isoforms may have similar but not fully coincident functions.
* This work was supported in part by University of Catania Grant COFIN2003058409_004 (to V. D. P.) and Deutsche Forschungsgemeinschaft Project Be 865/10 (to R. B.). The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked "advertisement" in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.
|| Recipient of a fellowship for the Ministry for Instruction, University, Research project "Hydrolases from Thermophiles: Structure, Function, and Homologous and Heterologous Expression" and supported by a project on Functional Genomics (Consiglio Nazionale delle Ricerche).
** Recipient of grants "Development and Implementation of Algorithms for Predicting Protein Structure" (Ministry for University, Research in Science and Technology), Molecular Genetics (Consiglio Nazionale delle Ricerche), and Postgenomics PNR 2001-2003, FIRB Article 8.
1 The abbreviations used are: VDAC, voltage-dependent anion-selective channel; EST, expressed sequence tag; for, forward; rev, reverse; nS, nanosiemens; Dm, D. melanogaster.
The nine D. melanogaster cDNA libraries used in this work were the kind gifts of Dr. Ruggiero Caizzi, Department of Genetics, Bari, Italy. We thank Prof. George Wolf, Dept. of Nutritional Sciences, University of California, Berkeley, for revision of the manuscript.
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
|
Advertisement | |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||