|
Originally published In Press as doi:10.1074/jbc.M006437200 on October 3, 2000
J. Biol. Chem., Vol. 275, Issue 52, 40992-40999, December 29, 2000
Structure and Orientation of Two
Voltage-dependent Anion-selective Channel
Isoforms
AN ATTENUATED TOTAL REFLECTION FOURIER-TRANSFORM INFRARED
SPECTROSCOPY STUDY*
Helge
Abrecht ,
Erik
Goormaghtigh§¶,
Jean-Marie
Ruysschaert§, and
Fabrice
Homblé
From the Laboratoire de Physiologie
Végétale, CP 206/2, Faculté des Sciences,
Université Libre de Bruxelles, Bld du Triomphe, B-1050 Brussels,
Belgium and the § Laboratoire de Chimie Physique des
Macromolécules aux Interfaces, CP 206/2, Faculté des
Sciences, Université Libre de Bruxelles, Bld du Triomphe,
B-1050 Brussels, Belgium
Received for publication, July 19, 2000, and in revised form, September 28, 2000
 |
ABSTRACT |
Two VDAC (voltage-dependent
anion-selective channel) isoforms were purified from seed cotyledons of
Phaseolus vulgaris by chromatofocusing chromatography.
Attenuated total reflection Fourier-transform infrared (ATR-FTIR)
spectroscopy was used to study the structural properties of the two
isoforms reconstituted in a mixture of asolectin and 5% stigmasterol.
The IR spectra of the two VDAC isoforms were highly similar indicating
50 to 53% anti-parallel -sheet. The orientation of the -strands
relative to the barrel axis was calculated from the experimentally
obtained dichroic ratios of the amide I -sheet component and the
amide II band. Comparing the IR spectra of the reconstituted VDAC
isoforms with the IR spectra of the bacterial porin OmpF, for which a
high resolution structure is available, provided evidence for a general
structural organization of the VDAC isoforms similar to that of
bacterial porins. Hydrogen-deuterium exchange measurements indicated
that the exchange of the amide protons occurs to a higher extent in the
two VDAC isoforms than in the OmpF porin.
 |
INTRODUCTION |
Mitochondria are surrounded by two distinct membranes. The
mitochondrial outer membrane is freely permeable for small hydrophilic solutes up to approximately 6 kDa. This property has been attributed to
the presence of a pore forming protein, the
voltage-dependent anion-selective channel
(VDAC)1 or mitochondrial
porin (1). Biophysical properties of VDAC channels inserted in planar
lipid bilayers have been studied in a large variety of organisms. VDAC
channels switch from a high conducting fully open state (3.6 to 4.5 nS
in 1 M KCl) at low voltages to several low-conducting
substates upon application of voltages higher than ±20 mV. VDAC is
slightly anion selective in the fully open state but switches to cation
selectivity in the low-conducting substates (2). In these substates,
the flow of negatively charged metabolites such as succinate, citrate, phosphate (3), and adenine nucleotides (4, 5) through the channel is
strongly reduced compared with the fully open state. Therefore, VDAC is
thought to regulate the mitochondrial metabolism by regulating the flux
of metabolites across the mitochondrial outer membrane.
VDAC channels are involved in different cellular events like the
induced release of cytochrome c, which constitutes an early step in apoptosis (6). VDAC is the binding site for cytoplasmic enzymes
like glycerol kinase and hexokinase (7) and for the cytoskeleton (8).
The different VDAC isoforms that exist in yeast, plants, and mammals
could be involved in specific functions. Two human VDAC isoforms have
different binding affinities for hexokinase (9). Plant VDAC isoforms
have slightly different electrophysiological properties (10), and
mitochondria isolated from yeast vdac minus mutants
show differences in their outer membrane permeability to NADH depending
on which of the three mouse VDAC isoforms was expressed (11).
Therefore, it seems reasonable to assume that the multitude of
functional properties that are attributed to VDAC could arise from the
presence of different VDAC isoforms in mitochondria.
VDAC proteins from different far related organisms generally exhibit
poor sequence homologies (~30%) but VDAC from plants seem to have
more conserved amino acid sequences (12). Despite poor sequence
homologies, it is assumed that all VDACs have a rather conserved
-barrel structure similar to that of bacterial porins. This
assumption arose from secondary structure predictions (13-15) that
suggested anti-parallel -strands as the main structural feature. CD
measurements of purified yeast VDACs (16, 17) indicated a high
-sheet content. Electron microscopic studies of two-dimensional
crystals of fungal VDAC and low-resolution images of human VDAC
crystals suggested that VDAC proteins form a pore (18, 19). Although
the molecular structure of several bacterial porins has been resolved
at the atomic level (20, 21), no high-resolution VDAC crystals are available.
Here, we report on the secondary structure, orientation, and
accessibility to the water phase of two VDAC isoforms purified from
seeds of P. vulgaris. A possible structural organization of
the VDAC isoforms is proposed and is compared with that of the
Escherichia coli porin OmpF.
 |
EXPERIMENTAL PROCEDURES |
Materials--
Seeds from P. vulgaris var.
"Streamline" were purchased from a local store. Percoll,
chromatofocusing media PBE 94, and Polybuffer PB 96 were from
Amersham Pharmacia Biotech. Asolectin (soybean phosphatidylcholine) and
stigmasterol were obtained from Sigma. Genapol-X-080 was purchased from
Fluka and D2O from Merck. Bio-Beads were from Bio-Rad. The
bacterial porin, OmpF from E. coli was a kind gift from Dr.
T. Schirmer (University of Basel, Basel, Switzerland).
VDAC Purification--
The two VDAC isoforms were purified from
seeds of P. vulgaris var. "Streamline" following the
procedure described previously with minor modifications. The two VDAC
proteins were purified and separated in a single chromatofocusing
chromatography step, i.e. the HTP batch described in the
previous protocol (22) was omitted. Briefly, seeds from P. vulgaris were soaked in tap water for 24 h and mitochondria
were isolated from the cotyledons by differential centrifugation steps
and further purified on a 28% Percoll gradient (23). Then, the final
mitochondrial pellet (300-400 mg of protein, 20-25 mg/ml) was mixed
with an equal volume of 2-fold concentrated chromatofocusing start
buffer (see below) supplemented with 4% (v/v) Genapol-X-080 to
solubilize the membrane proteins. This suspension was aliquoted (2 ml)
and frozen at 20 °C until use. A 10-ml (1.6 × 5 cm)
chromatofocusing column was prepared with polybuffer exchanger media
(PBE 94, Amersham Pharmacia Biotech) and equilibrated with the start
buffer containing 25 mM ethanolamine, 4 M urea,
0.1% (v/v) Genapol-X-080, pH 9.4. A 2-ml mitochondria suspension in
start buffer (containing 2% (v/v) Genapol-X-080) was thawed at room
temperature and incubated for 30 min at 4 °C. Insoluble material was
removed by centrifugation at 10,000 × g for 5 min
(Beckman, Allegra 21R). The supernatant, containing approximately 20 mg
of protein, was loaded onto the chromatofocusing column and proteins
were eluted with 10% (v/v) Polybuffer 96 (PB 96, Amersham Pharmacia
Biotech), 4 M urea, 0.1% (v/v) Genapol-X-080 at pH 7.0. The fractions containing the individual VDAC isoforms from 3 chromatofocusing columns (about 25 ml) were pooled and concentrated
about 20-fold (Centriprep-30, Amicon). The elution buffer was exchanged
against 1 mM Hepes, 0.1% (v/v) Genapol-X-080, pH 7.2, by
gel filtration (Sephadex G-75, Amersham Pharmacia Biotech) at a flow
rate of 7.5 cm/h. The proteins eluted at a concentration of about 0.07 to 0.12 mg/ml.
Reconstitution of the Purified VDAC Isoforms--
950 µg of
purified asolectin (24) and 50 µg of stigmasterol were dissolved in
chloroform and the solvent was evaporated under a stream of
N2 to have a thin lipid film which was dried further
overnight under vacuum. The dry lipids were resuspended in 2 ml of
protein solution (~0.1 mg of protein/ml). The
lipid/protein/detergent mixture was incubated 30 min at room
temperature under constant agitation followed by 3 freeze/thawing
cycles. The removal of detergent was achieved by 4 subsequent additions
of 40 mg of Bio-Beads (previously washed with methanol and water).
Incubation times were 45, 45, 30, and 30 min under constant agitation
at 4 °C. After removal of the final Bio-Beads, the proteoliposome
suspension was diluted with 1 mM Hepes, pH 7.2, to a total
volume of 4.5 ml. The proteins associated with the lipid vesicles were
pelleted at 37,000 rpm for 2 h (Beckman L7, SW-60 rotor). This
washing step was repeated once. Finally, the proteoliposomes were
suspended in 20-30 µl of 1 mM Hepes, pH 7.2. The
proteoliposomes had a lipid to protein ratio of about 10:1 (w/w).
Protein recovery was in the range of 25-35%.
Bacterial porin was diluted in 2 ml of 1 mM Hepes, 0.1%
(v/v) Genapol, pH 7.2, to a similar concentration as were the VDAC proteins and reconstitution was then carried out as above described for
VDAC. Liposome blanks were prepared with 2 ml of buffer in the absence
of protein.
In order to demonstrate the association of VDAC proteins with the
lipids, the proteoliposomes were centrifuged on a 32 to 5% (w/v)
sucrose gradient with a 40% (w/v) sucrose cushion for 16 h at
35,000 rpm (Beckman L7, SW-60 rotor). The gradient was then
fractionated from the bottom of the tube. Turbidity of each fraction
was measured at 405 nm and protein was detected with the BCA assay (Pierce).
IR Spectroscopy--
Attenuated total reflection Fourier
transform infrared (ATR-FTIR) spectra were obtained on a PerkinElmer
1720 FTIR spectrophotometer equipped with a liquid nitrogen-cooled
mercury cadmium telluride detector as described elsewhere (25).
The spectra were recorded at a nominal resolution of 4 cm 1 and 128 scans were averaged for each measurement. The
internal reflection element was a germanium ATR plate (50 × 20 × 2 mm, Harrick EJ2121) with an aperture angle of 45°,
yielding 25 internal reflections. The spectrophotometer was
continuously purged with dry air. Spectra with polarized incident light
and kinetics of hydrogen/deuteration exchange were measured on a Bruker
IFS-55 spectrophotometer equipped with liquid nitrogen-cooled mercury cadmium telluride detector and a polarizer mount assembly (KRS-5 polarizer).
Sample Preparation--
20-30 µl of reconstituted protein
were spread on one side of the ATR plate. The solvent was slowly
evaporated under a continuous stream of nitrogen. So, 15-30 µg of
protein were deposited on the ATR plate.
Analysis of the Secondary Structure--
The determination of
protein secondary structure was based on the multivariate statistical
analysis method for band shape recognition similar as previously
reported (26). This analysis used a commercially available partial
least-square package to determine fractional secondary structure
composition (PLSPlus 2.1 from Spectra Calc, Galactic Industries, Salem
NH). A reference set (RaSP50) was generated by collecting ATR-FTIR
spectra of 50 proteins with known x-ray structures selected to
represent a wide range of -helix and -sheet compositions as well
as 60 different protein domain folds. Before analysis, all spectra were
normalized so that the PLS algorithm was forced to use band shapes of
the amide I and amide II bands (between 1720 and 1500 cm 1) rather than absolute intensities. Cross-validation
procedures using the RaSP50 protein spectra have shown RMS structure
determination errors of ±4.5% for -helix and ±6.3% for
-sheet.
Orientation of the Secondary Structure--
ATR-FTIR
spectroscopy provides information on the orientation of the peptide
secondary structure embedded in the hydrated lipid film where the lipid
acyl chains are oriented perpendicular to the surface of the ATR plate
(27). The protein spectra were recorded with parallel and perpendicular
polarized light with respect to the plane of incidence, averaging 512 scans at a resolution of 2 cm 1. To determine the dichroic
ratio of the different secondary structure components, the amide I band
(1700 to 1600 cm 1) was decomposed by a least-square
curve-fitting procedure using a Cauchy (Lorenzian/Gaussian) function as
described previously (25). The integrated areas corresponding to the
-sheet components (near 1630 cm 1) were determined, and
the ratio of the integrated areas from the two polarized spectra then
yielded the dichroic ratio of the -sheet,
R ATR. The same procedure was
employed to determine the dichroic ratio, R ATR, of the -helix
component (at ~1655 cm 1). The dichroic ratio of the
amide II band was determined from the ratio of the integrated areas
between 1590 and 1505 cm 1 from the two polarized spectra.
The dichroic ratio of the lipid (C=O), called below
RATRiso, was obtained from the parallel and perpendicular
polarized spectra. RATRiso was used to compute the film
thickness and the values of the electric field components
EX, EY, and EZ at 1631 cm 1 as described previously (28). For these
calculations refractive indexes of 4.0 and 1.44 were used for the
Ge-plate and the sample film, respectively.
In a -sheet the amide I transition dipole moment is oriented
parallel to the C=O bond which is oriented perpendicular to the
-strand axis. The N-H bending contributes mainly to the amide II
transition dipole, which is oriented parallel to the strand axis.
Because two symmetry axes are needed to describe the sheet orientation,
the combination of the dichroic ratios from both amide I and amide II
bands are required to determine the orientation of the anti-parallel
-strands. In a barrel structure, axial symmetry around the barrel
symmetry axis can be partially obtained (29-31). The orientation of
the -strands was calculated according to Marsh (30, 31).
Kinetics of the Hydrogen/Deuterium Exchange--
Deuteration
kinetics were carried out as described previously (32, 33). A computer
program controlled the kinetic measurements. Prior to each experiment
10 spectra were recorded to verify the stability and the
reproducibility of the system. At time 0 the sample was flushed with
D2O-saturated N2 at a flow rate of 4 liters/min. For each spectrum 24 scans were accumulated at a resolution
of 2 cm 1. Background and water vapor spectra were
subtracted from the kinetic spectra as reported (32).
The areas of the amide I and amide II bands were calculated by
integration of the spectra between 1700 and 1600 cm 1 and
1590 and 1505 cm 1, respectively. The area of the amide II
band decreases as deuteration of the protein proceeds. To take into
account variations of the overall spectral intensities related to film
swelling upon hydration in the first minutes of the measurement, the
amide II/amide I ratio was determined for each spectrum. Thus, the H/D
exchange was monitored as the evolution of the ratio of amide II/amide I areas, expressed as percent of non-exchanged amide protons. The amide
II/amide I ratios of the non-deuterated spectra were taken as 100%,
whereas the 0% value corresponds to a zero absorption in the amide II
region (34).
We observed in our experiments a considerable absorbance of the lipids
in the amide I and amide II region that overlapped the protein spectra.
To obtain the pure protein spectra, the spectra of H/D exchange kinetic
of a liposome blank were recorded under identical conditions. A
computer program subtracted the spectra of a liposome blank from each
corresponding sample spectrum. To account for varying absorbance
intensities in sample and lipid spectra, a subtraction coefficient was
used to cancel the lipid (C=O) area integrated between 1705 and 1770 cm 1.
 |
RESULTS |
VDAC Purification--
As a modification to the previously
reported protocol (22), purification of the two VDAC isoforms was
achieved in a single purification step. Mitochondria were isolated from
the cotyledons of 24-h imbibed kidney beans (P. vulgaris)
seeds and solubilized with 2% (v/v) Genapol X-080 and 4 M
urea. The solubilized proteins were loaded on and eluted from a
chromatofocusing column using a pH gradient from pH 9.4 to 7.0. Fig.
1 shows a Coomassie Blue-stained SDS gel
with the two individual VDAC proteins which were separated and purified
by chromatofocusing chromatography (Fig. 1, lanes 1 and
3). The two isoforms are called VDAC 31 and VDAC 32 according to their respective apparent molecular mass. SDS gels loaded
with high amounts of protein (10-20 µg) confirm the high purity of the protein preparation (Fig. 1, lanes 2 and 4).
The staining intensities reflect the abundance of each isoform in the
seed mitochondria (22).

View larger version (22K):
[in this window]
[in a new window]
|
Fig. 1.
Purification of two VDAC isoforms. VDAC
31, lane 1, and VDAC 32, lane 3, were purified
from mitochondria of P. vulgaris cotyledons using a 10 ml of
PBE column and a pH gradient from pH 9.4 to 7.0 in the presence of 4 M urea. Fractions from 3 columns were concentrated.
Lanes 2 and 4 are VDAC 31 and VDAC 32, respectively, recovered from Ge plates after IR spectra have been
recorded to demonstrate that no degradation occurred during the
reconstitution process and IR measurements. Note at the bottom of
lane 1, where likely the Coomassie Blue-stained ampholytes
appeared as a large spot that disappeared after the gel filtration
(compare with lane 2). The 12% SDS gel was Coomassie Blue
stained.
|
|
Reconstitution--
The purified VDAC isoforms were reconstituted
in asolectin supplemented with 5% (w/w) stigmasterol. The detergent
was removed from the mixed lipid-detergent-protein micelles by
adsorption on Bio-Beads to promote the formation of proteoliposomes.
Formation of proteoliposomes was demonstrated by the co-migration of
lipids and proteins during a sucrose gradient centrifugation (Fig.
2). No protein aggregates were detected
at the bottom of the tube. Fig. 2 shows the distribution pattern for
the VDAC 31, VDAC 32, and bacterial OmpF porin. One population of
proteoliposomes was obtained with each protein, indicating a quite
homogenous distribution of protein within the liposomes. The
functional channel properties of the reconstituted VDAC isoforms were
quite similar to those reported for other VDACs (22) (data not
shown).

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 2.
Sucrose density gradient centrifugation of
reconstituted VDAC 31 (A), VDAC 32 (B), and OmpF (C). The proteins
were reconstituted into liposomes as described. The proteoliposomes
were centrifuged on a discontinuous 40, 32-5% (w/v) sucrose gradient
at 120,000 × g for 16 h at 4 °C. Fractions
were collected from the bottom of the tube. Turbidity
(circles) of each fraction was determined spectrometrically
at 405 nm. Protein concentration (squares) was determined
with the BCA assay.
|
|
ATR-FTIR Spectroscopy--
We compared the IR spectra of the
reconstituted VDAC isoforms with the IR spectrum of the bacterial porin
OmpF for which a high-resolution crystal structure is available (Fig.
3A). Fourier deconvolution of
the spectra allows a more detailed analysis of the spectral composition
(Fig. 3B). Spectra of VDAC 31 and VDAC 32 are almost
identical in the region of the amide I band between 1700 and 1600 cm 1. The amide I bands are characterized by a major
component located around 1631 cm 1 which can be assigned
to -sheet. The position of the major component at 1631 cm 1 and a weak component at about 1694 cm 1
in the amide I region, together with the major amide II band at 1535 cm 1, indicate that the -strands are in an
anti-parallel configuration (35). The shape of the amide I bands of the
two VDAC isoforms are very similar to that of the OmpF porin,
suggesting a similar secondary structure composition. The IR spectrum
of the OmpF is virtually identical to that previously reported (36,
37). The secondary structure content of the VDAC isoforms and of the OmpF is reported in Table I. The
secondary structure of all three proteins is very similar. It must be
noted here that the variations observed are in the range of the error
made in evaluating the secondary structure. The secondary structure
determined here for OmpF closely matches the crystal structure (Table
I) (21).

View larger version (23K):
[in this window]
[in a new window]
|
Fig. 3.
A, ATR-FTIR spectra of VDAC 31 (a), VDAC 32 (b), and OmpF (c)
reconstituted in asolectin and 5% (w/w) stigmasterol, representing the
amide I and the amide II region. From each spectrum the spectrum a
liposome blank has been subtracted to obtain the pure protein spectrum
(a subtraction coefficient was chosen to cancel the (C=O) band
between 1770 and 1705 cm 1). B, Fourier
self-deconvolution of the original spectra in A with an
enhancement factor of K = 1.8. The ordinate scales in A
and B correspond to the spectrum of VDAC 31 (a). The spectra
of VDAC 32 (b) and OmpF (c) have been rescaled
and are off-set for better comparison and clarity.
|
|
View this table:
[in this window]
[in a new window]
|
Table I
Secondary structure of the reconstituted VDAC isoforms and bacterial
porin determined by ATR-FTIR
Quantitative determination of the secondary structure content was based
on spectral band shape recognition using the FTIR spectra of 50 known
proteins as the reference set. The values in parentheses are deduced
from the OmpF crystal structure.
|
|
Orientation of the Secondary Structure--
The ATR-FTIR
spectroscopy technique provides information about the orientation of
proteins inserted in a lipid bilayer. Protein spectra were recorded
with parallel and perpendicular polarized incident light (Fig.
4A). The dichroic spectra
(i.e. parallel minus perpendicular polarized spectrum) of
the two VDAC isoforms and of the OmpF are illustrated in Fig.
4B. The three spectra show a positive dichroism signal for
the main component of the amide I band at ~1630 cm 1 and
a positive deviation of the amide II band indicating that the -sheet
structure has a preferred orientation. A similar dichroic spectrum of
porin has been reported previously (37). To determine the orientation
of the -strands, the amide I band was decomposed into its different
components by a least-square curve-fitting procedure to obtain the
dichroic ratio R ATR of the
-sheet component. For a quantitative determination of the -strand
orientation, the dichroic ratios of both the amide I band (here
R ATR) and the amide II band
have to be determined because the transition dipole moment of
-strands does not exhibit uni-axial symmetry (29). The orientation
I of the amide I transition dipole moment with respect
to normal to the membrane is related to the strand tilt, , by
I = ( /2) . For the amide II transition
dipole moment, which is oriented parallel to the strand axis,
II = . For a regular -barrel, the transition
dipole moments are distributed with axial symmetry around the barrel
axis and the orientation of the transition dipole moments is related to
the dichroic ratios RATR by,
|
(Eq. 1)
|
where P2(x) = (1/2)(3x2 1) is a second order Legendre
polynomial and is the orientation of the barrel axis relative to the membrane normal (30, 31). This equation can be written for either
the amide I band or the amide II band, so that two equations are
obtained to calculate the order parameters. Ex, Ey, and Ez, are the amplitudes of
the electric field in a coordinate system, where z is
perpendicular to the ATR plate surface, x and y
are in the ATR-plate plane with x in the incidence plane and
y perpendicular to the incidence plane. The values of
Ex, Ey, and Ez were computed considering the
real film thickness (see "Experimental Procedures"). The tilt
angle, , was calculated from the experimental dichroic ratios of
both the amide I (R ATR) and
the amide II bands using Equation 1. The -strands in all three
proteins are tilted at a very similar angle of ~45° to 47° relative to the barrel axis (Table II). A
value of 45° for the strand tilt in OmpF is in good agreement with
the mean strand orientation deduced from the crystal structure (21,
31). The order parameter corresponding to the angle is similar for
the two VDAC isoforms ( P2(cos ) = 0.36 and 0.45 for
VDAC 31 and VDAC 32, respectively) but is significantly smaller than
for OmpF P2(cos ) = 0.66). This point will be
discussed later. The value for OmpF agrees well with data reported by
Marsh (31). The contribution of the -helical component of the amide
I band was too low (~10%) to have enough accuracy to calculate the
orientation of the -helix.

View larger version (27K):
[in this window]
[in a new window]
|
Fig. 4.
A, polarized ATR-FTIR spectra of
reconstituted VDAC 31 recorded with parallel polarized incident light
(dashed line) and recorded with perpendicular polarized
incident light (line). RATRiso, the
film thickness and the electric field amplitudes were determined as
described under "Experimental Procedures" from the two polarized
spectra. For VDAC 31 RATRiso was determined to
be 1.58 (± 0.01), the film thickness, d, was 0.33 (±0.01) µm and
the electric field components EX2,
EY2, and
EZ2 were 1.96 (±0.002), 2.23 (±0.007),
and 1.55 (±0.01), respectively, for three independent determinations.
For VDAC 32 RATRiso = 1.62 (±0.04), d = 0.4 (±0.05) µm, EX2 = 1.97 (±0.002),
EY2 = 2.23 (±0.007),
EZ2 = 1.62 (±0.004) (n = 3). For OmpF RATRiso = 1.45 (±0.05), d = 0.22 (±0.04) µm, EX2 = 1.98 (±0.003), EY2 = 2.2 (±0.01),
EZ2 = 1.21 (±0.15), n = 2). B, dichroic spectra of reconstituted VDAC 31 (line), VDAC 32 (dashed line), and OmpF
(dashed-dotted line). The dichroic spectra were obtained by
subtracting the perpendicular spectra from the parallel spectra.
Before subtraction, the perpendicular spectra were multiplied by the
corresponding RATRiso for scaling. For clarity,
the spectra of OmpF and VDAC 31 are off-set by +2 and 2
milliabsorption units, respectively.
|
|
View this table:
[in this window]
[in a new window]
|
Table II
Orientation of the -strands in VDAC and OmpF
The dichroic ratio of the -sheet component of the amide I band and
the dichroic ratio of the amide II band were determined from parallel
and perpendicular polarized ATR-FTIR spectra of the different proteins
in films of asolectin + 5% stigmasterol. The tilt angle, , of
the -strands with respect to the barrel axis were derived from
Equation 1. The data are presented as the mean of three experiments
(for OmpF, n = 2), ± SE.
|
|
Kinetics of the Hydrogen/Deuterium Exchange--
The kinetic of
H/D exchange of the amide protons is monitored by evaluating the amide
II/amide I area ratio. As amide II arises predominantly from the
peptide N-H bending, the rate of H/D exchange is related to both the
solvent accessibility of the NH amide groups of the peptide bond (which
is related to the tertiary structure) and to the stability of secondary
structure elements. The kinetics of H/D exchange was followed by
monitoring the decrease of the area of the amide II band as a function
of time to D2O exposure. A series of ATR-FTIR spectra
recorded upon H/D exchange of reconstituted VDAC 31 appears in Fig.
5. Clearly, the amide II area (near 1530 cm 1) decreases as deuteration proceeds due to a shift of
the amide II band upon deuteration to near 1445 cm 1
(called amide II'). The deuteration kinetics of the two VDAC isoforms
and of the OmpF porin are compared in Fig.
6. The H/D exchange occurs to different
extents in VDAC and OmpF. For both VDAC isoform more than 60% of the
amide hydrogens are exchanged after 10 h. On the contrary, under
the same experimental conditions, only about 35% of the OmpF amide
protons are exchanged which is consistent with previously reported
values (36). Although the decrease of the amide II intensity measures
the global H/D exchange, the exchange kinetics can quantitatively be
analyzed as follows. The exchange curves follow a multiexponential
decay. The H/D exchange is a first-order reaction involving different
groups, i, of amide protons characterized by a common period
Ti,

View larger version (83K):
[in this window]
[in a new window]
|
Fig. 5.
Deuteration of VDAC 31. The spectra were
recorded between 1800 and 1400 cm 1 as a function of time
exposure to D2O saturated N2 flow. Integration
of the area of the amide I and amide II bands was performed using the
base lines (dotted lines) drawn between the intersection of
the spectra (vertical lines). Spectra of a liposome blank
recorded under identical conditions have been subtracted from each
corresponding sample spectrum. For subtraction, a coefficient was
computed to cancel the (C=O) band between 1770 and 1705 cm 1 (see "Experimental Procedures").
|
|

View larger version (21K):
[in this window]
[in a new window]
|
Fig. 6.
Deuterium/hydrogen exchange kinetics of VDAC
31 (squares), VDAC 32 (circles), and
OmpF (triangles) reported as percent non-exchanged
hydrogen. Deuteration percentage was evaluated from the evolution
of the area amide II/area amide I ratio upon exposure to
D2O. The vertical bars represent the standard
deviation of the mean of two (VDAC 31) or three (VDAC 32) independent
experiments.
|
|
|
(Eq. 2)
|
where ai is the number of amino acids of the
group i. Within the experimental reproducibility limits, a
nonlinear fitting of the experimental curves (Fig. 6) could be
performed with three exponentials indicating that the exchange can be
described as the exchange of three groups of amide protons with the
proportions a1, a2, and a3
characterized respectively by the periods T1,
T2, and T3. Each time constant (T1,
T2, and T3) was very similar for the two VDAC
isoforms. To compare the proportions of each group of amide protons,
the Ti values of the two VDAC kinetics were replaced by
their mean values (T1 = 1.1, T2 = 10, and
T3 = 2300 min) and a second fitting was performed using
Ti constrained to their mean values (32). The number of
amide protons (in % of total amide protons) in each group are reported
in Table III. The three groups of the
OmpF amide protons had time constants of T1 = 1.3, T2 = 28, and T3 = 5637 min (Table III). Only
17% of the OmpF amide protons but 30% of the VDAC amide protons
exhibit a fast exchange rate with comparable time constants of
T1 = 1.3 and 1.1 min, respectively. On the contrary,
two-thirds of the OmpF amide protons exhibit a very slow exchange rate
(T3 = 5637 min). In comparison, only 48% of the VDAC amide
protons exhibit a slow exchange rate but still faster than in OmpF
(Table III).
View this table:
[in this window]
[in a new window]
|
Table III
Analysis of the hydrogen-deuterium exchange kinetics
a1, a2, and a3 represent the proportions of
amide protons in percent of the total amide protons of the two VDAC
isoforms and OmpF. The groups of the VDACs amide protons are
characterized by time constants of T1 = 1.1, T2 = 10, and T3 = 2300 min. Those of the OmpF porin have time constants
of T1 = 1.3, T2 = 28, T3 = 5637 min. The
proportions were obtained from the analysis of the H/D exchange curves
recorded for the two VDAC isoforms or the bacterial porin OmpF as
described in the text.
|
|
 |
DISCUSSION |
Two VDAC isoforms were purified from kidney bean cotyledons in a
single purification step using chromatofocusing chromatography as
previously reported (22). Highly purified proteins were reconstituted in a mixture of asolectin and stigmasterol, a major sterol of plant
membranes. Sterols have been found to be associated with VDAC isolated
from Neurospora crassa and bovine heart mitochondria (38,
39). It has been suggested recently that sterols are essential for the
proper folding of VDAC (40).
The IR spectra of the reconstituted VDAC isoforms are strikingly
similar displaying characteristics typical for an anti-parallel -pleated sheet. The two spectra could indeed almost be superimposed, strongly suggesting that the two isoforms have an almost identical secondary structure. Quantitative analysis indicated that about 50 to
53% of the amino acid residues are in -sheet conformation. Moreover, the percentages of the different secondary structure elements
calculated from the amide I and amide II bands associated to OmpF and
VDAC were almost identical. The agreement of the secondary structure of
OmpF determined by IR analysis and from the crystal structure (Table I)
strongly supports the validity of the structures proposed for VDAC.
Despite the remarkable similarity between the VDAC and the OmpF IR
spectra a few differences can be observed. The main peak of the VDACs
amide I bands (1631 cm 1) is slightly shifted toward a
higher wavenumber compared with the peak of the OmpF spectrum (1629 cm 1). Kleffel et al. (36) correlated the
position of the main -sheet peak with the length (i.e.
the number of amino acids) of the -strands. It has been found that
the peak at around 1630 cm 1 shifted toward higher
wavenumbers for shorter -strands. The slight shift in the VDAC
spectra could therefore be an indication for shorter -strands in
VDAC than in OmpF (~12 residues on the average). The -helix
content in VDAC (~9 to 13%) seems slightly higher than in OmpF,
visible as a shoulder in the VDAC IR spectra at about 1660 cm 1 (Fig. 3). This observation could correlate with the
assumption that the N terminus of VDAC forms an amphipathic -helix.
It has been suggested that this -helical N terminus might play a
particular role in protein import into the mitochondrial outer membrane
and in the voltage-dependent gating of the channel
(41-43).
The possibility to obtain information on the orientation of the
different secondary structures is certainly one of the greatest advantages of ATR-FTIR spectroscopy. Information on the orientation of
the dipole can be obtained by measuring the IR spectra with different
polarized incident light (polarized parallel and perpendicular with
respect to the plane of incidence). The accuracy of the results obtained with this method has been recently demonstrated with bacterial
porins, the atomic structures of which are known (31). The same method
was applied here to get information on the structure of two membrane
proteins for which no high-resolution structures were available. For a
quantitative analysis of the -strand tilt, the dichroic ratio of the
lipid (C=O) band was chosen as representative of a dichroic ratio
that would correspond either to an unordered dipole or to a dipole
oriented at the magic angle with respect to the Ge-plate normal. It is
called below RATRiso. With the value of
RATRiso the film thickness and the electric
field amplitudes were computed. This procedure assumes that the
orientation of the lipid (C=O) transition dipole moment is closed to
the magic angle (54.7°). In fact, NMR and x-ray studies carried out
on several phospholipids have suggested that the average orientation of
the two fatty acyl chain (C=O) transition dipole moments are close
to the magic angle (44). The validity to determine
RATRiso from the lipid (C=O) band of
polarized ATR-FTIR spectra has been demonstrated recently (45).
In this work, we calculated the orientation of the main secondary
structure, i.e. anti-parallel -strands, from the
experimentally obtained dichroic ratio of the -sheet component of
the amide I band and the dichroic ratio of the amide II band. Our
results show that the -strands of VDAC and the OmpF porin are tilted at a very similar angle with respect to the barrel axis ( = 45 to 47°, Table II). This result was obtained from expressions derived for regular -barrels, which have an axial distribution of transition dipole moments around the barrel symmetry axis. It can be compared with
the topology of planar -sheets which display non-axial symmetry and
for which the orientation of the -strands can be calculated according to Marsh (29). Applying this approach to the dichroic ratios
given in Table II we obtain a tilt angle of the -strands ~ 45° for both the two VDAC isoforms and the OmpF porin. Regarding OmpF, these results are in very good agreement with the value derived
from the crystal structure (~45°) where 11 strands are tilted at
about 35° and 5 strands exhibit a more oblique tilt (21).
The structural geometries of a regular -barrel are mainly determined
by the number, n, of -strands and the tilt angle of
the -strands (46). The radius, R, of a cylindrical barrel is related to the strand tilt by,
|
(Eq. 3)
|
where d = 0.472 nm is the distance between two
residues in a -strand. From the knowledge of the barrel diameter and
the mean orientation of the -strands determined here by polarized ATR-FTIR, the number of -strands in VDAC can be calculated. A diameter of VDAC of 3.7 nm has been determined by electron microscopy studies of yeast VDAC two-dimensional crystals and from low-resolution crystals of human VDAC (18, 19). Applying Equation 3 to this diameter
and a mean -strand tilt of 46° to 47° results in a barrel that
is formed from 17 -strands. Because -barrels have an even number
of -strands VDAC could have either 16 or 18 -strands (a tilt of
48° would give 16 -strands and 44° 18 -strands). Secondary
structure predictions have suggested both, 16 and 18 -strands (14,
47). The 16 (or 18) -strands would have a mean length of
approximately 9 (or 8) residues per strand, when 50 to 53% of the
total residues (~280) are in -sheet conformation. From the strand
tilt determined for OmpF (45°) and a mean diameter of 3.4 nm 16 -strands are calculated which corresponds to that observed in the
crystal structure.
Apparently, the VDAC -barrels are somewhat more inclined in the
lipid membrane than the OmpF barrels because the order parameters determined for the VDAC isoforms are significantly smaller than the one
measured for OmpF. Assuming a similar ordering of the membrane in both
systems, it can be calculated that the VDAC barrels are tilted at a
higher angle than the OmpF barrel.
Hydrogen/deuterium exchange has long been used for the analysis of
protein structure and dynamics. The H/D exchange data contain information regarding the strength of H-bonding, the solvent
accessibility, and dynamic structure of proteins. Because the exchange
of only the amide protons is monitored by IR spectroscopy, the exchange data are directly proportional to the number of amino acid residues in
the protein. The exchange data suggest that the amide protons of VDAC
and OmpF can be grouped in three populations with different time
kinetics corresponding to fast, intermediate and slow exchange rates.
Analysis of the data provide evidence that about 69% of the total OmpF
amino acids exchange very slowly, characterized by T3 = 5637 min. Only about 48% of the VDAC amide protons exhibit a slow
exchange rate, albeit 2.4 times faster than the OmpF amide protons. In
contrast, the proportion of the amide protons that has a fast exchange
rate is larger in VDAC (30%) than in OmpF (17%). The exchange
proceeded to a larger extent in VDAC than in OmpF. More than 60% of
the VDAC amide protons were exchanged after 10 h compared to only
35% of the OmpF amide protons. Since both VDAC and OmpF form water
filled pores, solvent accessibility should not be a limiting factor for
the H/D exchange. In fact, it has been shown that high exchange rates
correlates with the occurrence of water filled pores in a
K+-channel and the aquaporin CHIP28 (48, 49). Apparently,
strong H bonding causes the low exchange in OmpF whereas the exchange data indicate that less strong H bonding may occur in VDAC. When compared with the bacterial porin, a less stable sheet structure can be
correlated with shorter strands in agreement with the first, the
determination of the number and mean length of the strands in the
barrel (see above) and second, the higher frequency of the sheet
contribution to the amide I band as discussed above.
In conclusion, we show in this report that the two VDAC isoforms from
P. vulgaris have a very similar secondary structure and
orientation of the -sheet and furthermore, their H/D exchange is
undistinguishable. Taken together, this suggests that the two isoforms
have a very similar overall structure. Furthermore, our data provide an
experimental support to a general structural organization of both VDAC
isoforms similar to the bacterial porins. VDAC could form a somewhat
tilted -barrel (Fig. 7) composed of 16 (or 18) -strands with a mean length of 9 (or 8) residues per strand. The strands are tilted at about 46 to 47° relative to the barrel axis
to form a barrel with a diameter of 3.7 nm.

View larger version (29K):
[in this window]
[in a new window]
|
Fig. 7.
Model of a VDAC
-barrel. About 50-53% of the protein's
amino acids are involved in the barrel formation. The -barrel is
formed from 16 (or 18) anti-parallel -strands which have a mean
length of 9 (or 8) amino acids. The triangles and
circles represent alternating hydrophobic and hydrophilic
residues facing either the membrane lipids or the barrel lumen. H
bonding between neighboring -strands is indicated. The strands are
tilted at ~ 46° relative to the barrel axis to form a
barrel (backbone) diameter of 3.7 nm. The -barrel is tilted in the
lipid bilayer at the angle .
|
|
 |
ACKNOWLEDGEMENTS |
We thank Dr. T. Schirmer from the University
of Basel, Basel, Switzerland, for the OmpF porin gift and G. Vandenbussche for scientific discussion and advice.
 |
FOOTNOTES |
*
This work was supported by a grant of the 178 Communnauté Française de Belgique-Actions de Recherches
Concertées 178 and TMR Marie Curie Research Training Grant
ERBFMBICT971977 of the European Commission (to H. A.).The costs of publication of this article were defrayed in part by the
payment of page charges. The article
must therefore be hereby marked
"advertisement" in accordance with 18 U.S.C. Section
1734 solely to indicate this fact.
¶
Research Director from the National Fund for Scientific
Research (Belgium).
Research Director from the National Fund for Scientific
Research (Belgium). To whom correspondence should be addressed:
Université Libre de Bruxelles, Laboratoire de Physiologie
Végétale, CP 206/2, Bld du Triomphe, B-1050 Brussels,
Belgium. Tel.: 32-2-650-5383; Fax: 32-2-650-5382; E-mail:
fhomble@ulb.ac.be.
Published, JBC Papers in Press, October 3, 2000, DOI 10.1074/jbc.M006437200
 |
ABBREVIATIONS |
The abbreviations used are:
VDAC, voltage-dependent anion-selective channel;
ATR-FTIR, attenuated total reflection Fourier-transform infrared.
 |
REFERENCES |
| 1.
|
Benz, R.
(1994)
Biochim. Biophys. Acta
1197,
167-196
|
| 2.
|
Benz, R.,
Kottke, M.,
and Brdiczka, D.
(1990)
Biochim. Biophys. Acta
1022,
311-318
|
| 3.
|
Hodge, T.,
and Colombini, M.
(1997)
J. Membr. Biol.
157,
271-279
|
| 4.
|
Benz, R.,
Wojtczak, L.,
Bosch, W.,
and Brdiczka, D.
(1988)
FEBS Lett.
231,
75-80
|
| 5.
|
Rostovtseva, T.,
and Colombini, M.
(1997)
Biophys. J.
72,
1954-1962
|
| 6.
|
Shimizu, S.,
Narita, M.,
and Tsujimoto, Y.
(1999)
Nature
399,
483-487
|
| 7.
|
Adams, V.,
Griffin, L.,
Towbin, J.,
Gelb, B.,
Worley, K.,
and McCabe, E. R.
(1991)
Biochem. Med. Metab. Biol.
45,
271-291
|
| 8.
|
Linden, M.,
and Karlsson, G.
(1996)
Biochem. Cell Biol. Commun.
218,
833-836
|
| 9.
|
Blachly-Dyson, E.,
Zambronicz, E. B., Yu, W. H.,
Adams, V.,
McCabe, E. R.,
Adelman, J.,
Colombini, M.,
and Forte, M.
(1993)
J. Biol. Chem.
268,
1835-1841
|
| 10.
|
Elkeles, A.,
Breiman, A.,
and Zizi, M.
(1997)
J. Biol. Chem.
272,
6252-6260
|
| 11.
|
Xu, X.,
Decker, W.,
Sampson, M. J.,
Craigen, W. J.,
and Colombini, M.
(1999)
J. Membr. Biol.
170,
89-102
|
| 12.
|
Roosens, N.,
Al Bitar, F.,
Jacobs, M.,
and Homblé, F.
(2000)
Biochim. Biophys. Acta
1463,
470-476
|
| 13.
|
Forte, M.,
Guy, H. R.,
and Mannella, C. A.
(1987)
J. Bioenerg. Biomembr.
19,
341-350
|
| 14.
|
Rauch, G.,
and Moran, O.
(1994)
Biochem. Biophys. Res. Commun.
200,
908-915
|
| 15.
|
Mannella, C. A.,
Neuwald, A. F.,
and Lawrence, C. E.
(1996)
J. Bioenerg. Biomembr.
28,
163-169
|
| 16.
|
Shao, L.,
Kinnally, K. W.,
and Mannella, C. A.
(1996)
Biophys. J.
71,
778-786
|
| 17.
|
Koppel, D. A.,
Kinnally, K. W.,
Masters, P.,
Forte, M.,
Blachly, D. E.,
and Mannella, C. A.
(1998)
J. Biol. Chem.
273,
13794-13800
|
| 18.
|
Guo, X. W.,
Smith, P. R.,
Cognon, B.,
D'Arcangelis, D.,
Dolginova, E.,
and Mannella, C. A.
(1995)
J. Struct. Biol.
114,
41-59
|
| 19.
|
Dolder, M.,
Zeth, K.,
Tittmann, P.,
Gross, H.,
Welte, W.,
and Wallimann, T.
(1999)
J. Struct. Biol.
127,
64-71
|
| 20.
|
Weiss, M. S.,
Kreusch, A.,
Schiltz, E.,
Nestel, U.,
Welte, W.,
Weckesser, J.,
and Schulz, G. E.
(1991)
FEBS Lett.
280,
379-382
|
| 21.
|
Cowan, S. W.,
Schirmer, T.,
Rummel, G.,
Steiert, M.,
Ghosh, R.,
Pauptit, R. A.,
Jansonius, J. N.,
and Rosenbusch, J. P.
(1992)
Nature
358,
727-733
|
| 22.
|
Abrecht, H.,
Wattiez, R.,
Ruysschaert, J. M.,
and Homblé, F.
(2000)
Plant Physiol.
124,
1181-1190
|
| 23.
|
Douce, R.,
Bourguignon, J.,
Brouquisse, R.,
and Neuburger, M.
(1987)
Methods Enzymol.
148,
403-415
|
| 24.
|
Kagawa, Y.,
and Racker, E.
(1971)
J. Biol. Chem.
246,
5477-5487
|
| 25.
|
Goormaghtigh, E.,
Cabiaux, V.,
and Ruysschaert, J. M.
(1990)
Eur. J. Biochem.
193,
409-420
|
| 26.
|
Cabiaux, V.,
Oberg, K. A.,
Pancoska, P.,
Walz, T.,
Agre, P.,
and Engel, A.
(1997)
Biophys. J.
73,
406-417
|
| 27.
|
Fringeli, U. P.,
and Gunthard, H. H.
(1981)
Mol. Biol. Biochem. Biophys.
31,
270-332
|
| 28.
|
Goormaghtigh, E.,
Raussens, V.,
and Ruysschaert, J. M.
(1999)
Biochim. Biophys. Acta
1422,
105-185
|
| 29.
|
Marsh, D.
(1997)
Biophys. J.
72,
2710-2718
|
| 30.
|
Marsh, D.
(1998)
Biophys. J.
75,
354-358
|
| 31.
|
Marsh, D.
(2000)
J. Mol. Biol.
297,
803-808
|
| 32.
|
Goormaghtigh, E.,
Vigneron, L.,
Scarborough, G. A.,
and Ruysschaert, J. M.
(1994)
J. Biol. Chem.
269,
27409-27413
|
| 33.
|
Sonveaux, N.,
Shapiro, A. B.,
Goormaghtigh, E.,
Ling, V.,
and Ruysschaert, J. M.
(1996)
J. Biol. Chem.
271,
24617-24624
|
| 34.
|
de Jongh, H. H.,
Goormaghtigh, E.,
and Ruysschaert, J. M.
(1995)
Biochemistry
34,
172-179
|
| 35.
|
Goormaghtigh, E.,
Cabiaux, V.,
and Ruysschaert, J. M.
(1994)
Subcell. Biochem.
23,
329-362
|
| 36.
|
Kleffel, B.,
Garavito, R. M.,
Baumeister, W.,
and Rosenbusch, J. P.
(1985)
EMBO J.
4,
1589-1592
|
| 37.
|
Nabedryk, E.,
Garavito, R. M.,
and Breton, J.
(1988)
Biophys. J.
53,
671-676
|
| 38.
|
Freitag, H.,
Genchi, G.,
Benz, R.,
Palmieri, F.,
and Neupert, W.
(1982)
FEBS Lett.
145,
72-76
|
| 39.
|
De Pinto, V.,
Benz, R.,
and Palmieri, F.
(1989)
Eur. J. Biochem.
183,
179-187
|
| 40.
|
Popp, B.,
Schmid, A.,
and Benz, R.
(1995)
Biochemistry
34,
3352-3361
|
| 41.
|
Kleene, R.,
Pfanner, N.,
Pfaller, R.,
Link, T. A.,
Sebald, W.,
Neupert, W.,
and Tropschug, M.
(1987)
EMBO J.
6,
2627-2633
|
| 42.
|
Pfaller, R.,
Kleene, R.,
and Neupert, W.
(1990)
Experientia
46,
153-161
|
| 43.
|
Thomas, L.,
Blachly, D. E.,
Colombini, M.,
and Forte, M.
(1993)
Proc. Natl. Acad. Sci. U. S. A.
90,
5446-5449
|
| 44.
|
Hauser, H.,
Pascher, I.,
and Sundell, S.
(1988)
Biochemistry
27,
9166-9174
|
| 45.
|
Bechinger, B.,
Ruysschaert, J. M.,
and Goormaghtigh, E.
(1999)
Biophys. J.
76,
552-563
|
| 46.
|
Murzin, A. G.,
Lesk, A. M.,
and Chothia, C.
(1994)
J. Mol. Biol.
236,
1369-1381
|
| 47.
|
Mannella, C. A.
(1997)
J. Bioenerg. Biomembr.
29,
525-531
|
| 48.
|
Tatulian, S. A.,
Cortes, D. M.,
and Perozo, E.
(1998)
FEBS Lett.
423,
205-212
|
| 49.
|
Harris, A. L.,
Walter, A.,
and Zimmerberg, J.
(1989)
J. Membr. Biol.
109,
243-250
|
Copyright © 2000 by The American Society for Biochemistry and Molecular Biology, Inc.

CiteULike Complore Connotea Del.icio.us Digg Reddit Technorati What's this?
Copyright © 2000 by the American Society for Biochemistry and Molecular Biology.
|
Advertisement
Advertisement
|