Crystallization and overview of the overall fold of the Long Chain Acyl-CoA Dehydrogenase
Structures were obtained for wildtype LCAD complexed with and without acetoacetyl-CoA to 2.8 Å and 2.5Å resolution, respectively. In addition, the structure of the E291Q mutant complexed with lauric acid has been determined to 2.0 Å resolution (See Table 1). In accord with UniProKB, the residue numbering system in this report is based on the 430 amino acid preprotein that includes the 30-residue mitochondrial transport peptide. However, the protein used for our crystallization was the 400 amino acid mature protein without the transport peptide.
Table 1
Data Collection and Refinement Statistics
Data Collection
|
Wt-LCAD
|
Wt-LCAD
|
LCAD (E291Q)
|
PDB code
|
8W0T
|
8W0U
|
8W0Z
|
Bound ligand
|
None
|
Acetoacetyl-CoA
|
Lauric Acid
|
No. of measured reflections
|
252,360
|
309,146
|
851,736
|
Unique Refs.
|
52,001(2,227)
|
45,139 (4,365)
|
245,928 (11,015)
|
Resolution (Å)
|
50 − 2.5 (2.59–2.5)
|
50 − 2.8 (2.59–2.80)
|
50 − 2.0 (2.03-2.0)
|
Space Group
|
P212121
|
P21
|
P1
|
Unit Cell a (Å)
|
101.1
|
86.0
|
86.3
|
b
|
102.1
|
95.3
|
94.8
|
c
|
173.6
|
118.8
|
119.2
|
α (º)
|
90
|
90
|
89.24
|
β
|
90
|
106.52
|
74.88
|
γ
|
90
|
90
|
88.35
|
Rsym (last shell)
|
0.071(0.265)
|
0.078(0.474)
|
0.080(0.458)
|
Completeness (%)
|
82.8 (47.9)
|
98.6 (95.8)
|
96.7(87.1)
|
I/σI
|
15.5
|
11.7
|
14
|
Refinement
|
Resolution range
|
30 − 2.5
|
50 − 3.0
|
50 − 2.0
|
R factor (%)
|
20.0
|
22.5
|
21.4
|
Rfree (%)
|
26.2
|
29.7
|
25.6
|
No of water molecules
|
244
|
91
|
1859
|
Average B factor, Å2 (proteins)
|
49.3
|
63.2
|
38.0
|
Average B factor, Å2 (solvent)
|
36.0
|
35.9
|
40.4
|
Stereochemical ideality (rms deviation)
|
Bonds (Å)
|
0.007
|
0.006
|
0.007
|
Angles (°)
|
1.3
|
1.4
|
1.3
|
Ramachandran analysis
|
Residues in preferred regions (%)
|
90.7
|
87.6
|
95.4
|
Residues in allowed regions (%)
|
9.0
|
12.1
|
4.3
|
Residues in disallowed regions (%)
|
0.3
|
0.3
|
0.3
|
Values in parentheses are for the highest resolution shell. |
The overall structure of LCAD is similar to that of the other ACADs whose three-dimensional structures have been determined previously21–25. LCAD is a homotetramer and the tetrameric arrangement of the monomers is a dimer of dimers (Fig. 2a). Similar to the other known structures of acyl-CoA dehydrogenases, each LCAD monomer consists of three domains. The ribbon diagram of an LCAD monomer, showing the N-terminal α-helical domain (helices A-F), a β-sheet domain (strands 1–7) and a C-terminal α-helical domain (helices G-K), is presented in Fig. 2b. The first few residues at the N-terminus were not resolved in the wild type and E291Q structures. The chain starts at Glu34 in all four monomers of the mutant LCAD structure. In the wild type structure, the chain starts at Glu34 (monomer A) or Leu36 (monomers B, C, and D). The C-terminal residue Lys430 is observed in all the monomers in the 2.8 Å structure of the wild type protein complexed with acetoacetyl-CoA, but not in the 2.5 Å wild type ligand-free and 2.0 Å (E291Q mutant) structures.
Although no substrates or inhibitors could be co-crystallized with wildtype crystals, the inhibitor acetoacetyl-CoA could be soaked into these crystals. Interestingly, only two monomers (Molecules A and B) of the tetramer have acetoacetyl-CoA bound, presumably due to the crystal packing. Binding of acetoacetyl-CoA to Molecules C and D would result in steric clashes between the adenine ring of FAD in Mol C and loop residues located between Glu151 and His156 of symmetry related Mol A, and between the FAD adenine of Mol D and the loop comprised of Glu105-His106-Leu107 of symmetry related Mol B.
The crystal packing of the LCAD E291Q P1 form is very similar to the crystal packing of the wild type LCAD P21 crystal, except that, in the P1 space group, the adenine ring binding site now has moved slightly away from the symmetry related A and B molecules to avoid steric clashes, transforming the P21 space group of wild type LCAD to the P1 space group of the E291Q mutant crystal. Although we obtained crystals of the LCAD E291Q mutant co-crystallized with substrate lauroyl-CoA (C12-CoA), density for the CoA moiety is not observed, and only weak electron density of the alkyl chain portion is observed, likely due to hydrolysis of C12-CoA during the crystallization process.
As shown in the overlay of LCAD with MCAD and other ACAD structures (Supplementary Figure S1), the overall fold and topology of LCAD are similar to the other known ACAD structures. The root mean square deviations (RMSDs) between the main-chain atoms of human LCAD and rat SCAD, pig MCAD, and human VLCAD are 1.63Å, 1.78Å, and 1.83Å, respectively.
Catalytic residue and stereo-specificity of proton abstraction
The structure of LCAD confirms that Glu291 functions as the catalytic base, consistent with previous molecular modeling and site directed mutagenesis studies26. The structure of the catalytic mutant E291Q with substrate lauroyl-CoA (C12-CoA) bound in the active site shows that Gln291, extending out from helix G, is located at the active site close to the Cα-Cβ bond of the substrate (Fig. 3a). The distance between the carboxylate of the Glu291 and Cα of the CoA is 2.8 Å, and the distance between the N5 of the FAD and Cβ of the CoA is 3.8 Å. This arrangement of the isoalloxazine ring of the FAD, the Cα and Cβ of the thioester substrate, and the carboxylate of Glu291 is ideally suited for abstraction of the pro-R hydrogen as a proton and transfer of the Cβ pro-R hydrogen to the N-5 atom of the FAD as a hydride ion27, consistent with α-proton abstraction by Glu291. Crystallographic, mutagenesis, and chemical modification studies established Glu405 of MCAD (Glu376 of the mature protein) as the catalytic residue (see Thorpe and Kim review28). Sequence comparisons among all the members of the acyl-CoA dehydrogenase family show that the location of the catalytic base (Glu405 in MCAD) in the loop between helices J and K is conserved in all but two of the known acyl-CoA dehydrogenases, LCAD and IVD (Supplementary Figure S2). The crystal structure of human IVD complexed with CoA persulfide confirmed its catalytic residue as Glu283 (Glu254 of the mature protein)23. Although the catalytic base of MCAD, Glu405, is separated by more than 100 residues from those of LCAD and IVD, these residues are topologically conserved in the 3D structures of these proteins, and the proR-proR stereospecificity of the α,β- dehydrogenation reaction is maintained.
As in the structures of other acyl-CoA dehydrogenases, the carbonyl oxygen of the fatty acyl-CoA moiety of the substrate makes hydrogen bonds to the 2’-hydroxyl of the FAD ribityl chain (3.1Å) and to the peptide backbone of (Gly412-N, 2.6Å) (Fig. 3a). The corresponding distances in human MCAD are 2.8 and 2.9 Å. These interactions are not only responsible for the orientation of the substrate, but they are also crucial for polarization of the substrate and the lowering of the pKa of the Cα proton for abstraction by the catalytic base, Glu-29129,30.
Comparison between the Ligand-bound and Unbound subunits in wild-type LCAD and between wild-type and mutant subunits:
Among the known ACAD crystal structures, MCAD21, IBD25, and VLCAD31 are the only ACADs whose structures reported have been determined both with and without bound substrate or inhibitor. These ligand-bound and unbound ACAD structures have shown that, although subtle changes in the conformations of side chains lining the substrate binding cavity are observed, there are no large conformational changes accompanying ligand binding. Likewise, the overall structure of ligand-bound E291Q LCAD is similar to those of both the ligand bound and unliganded wild-type structures.
In all LCAD structures lacking bound CoA derivatives (all eight monomers of mutant LCAD, all four monomers of ligand-free wild type LCAD, and the C and D monomers of the ligand-bound wild type LCAD), residues Ser179 through Gln182 in the vicinity of the CoA moiety of the substrate are relatively disordered in the absence of ligand. His228 and nearby residues Asp180 and Leu181 are highly mobile, judging from their weak, diffused electron densities. The side chain of Tyr282 is rotated away from the conformation of the bound structure by approximately 125 degrees. New hydrogen bonding interactions contributing to significant CoA binding affinity are observed upon ligand binding (Fig. 3b). Arg424, situated at the entrance of the binding pocket, reorients from a position ~ 5.0Å from the CoA phosphate into the active site, forming a hydrogen bond (3.3Å) with the phosphate of the CoA. The side chain of His228, which lies within a hydrogen bonding distance (3.0Å) of the CoA phosphate, moves up and rotates away (Fig. 3b and Supplementary Figure S3). His228-Nε makes a H-bond with the phosphate of CoA. Having rotated approximately 125 degrees from its position in the unbound structure, the hydroxyl group of Tyr282 makes an H-bond with His228-ND1, so that the phenyl ring of Tyr282 stacks with the adenine ring of CoA.
Structural basis for substrate chain length and branched chain specificity:
LCAD is active toward substrates with chain length ranging from C10 to C18-CoA with C12- and C14-CoA as the optimum substrates32. Figure 4a shows the largely hydrophobic amino acid residues that line the fatty acid binding pocket of human LCAD. The LCAD binding cavity is deeper and wider compared to those of other mammalian ACADs, accounting for the lack of specific interactions between the fatty acyl chain of C12-CoA (bound substrate) and cavity-forming residues. In all eight monomers of E291Q LCAD, electron densities for the substrate fatty-acyl portion are weak and interactions with cavity-forming hydrophobic residues are not seen. This lack of specific interactions may account for the decreased activity of LCAD, in general, compared with MCAD and SCAD2.
Structural basis for metabolism of THC-CoA by LCAD
Consistent with previous reports of metabolism of THC-CoA by LCAD2, inspection of the active site cavity reveals that it is considerably larger than required for C12-CoA (Fig. 4a) and large enough to accommodate THC-CoA. Figure 4b shows a model of THC-CoA bound to the LCAD active site pocket. Hydrogen bonding interactions between THC-CoA and residues Ala125, Asn128, Gln408, and Tyr411 are shown. Furthermore, the structural basis for the unique substrate binding cavity can be seen in the long E helix, where Pro132 is located in a position to cause unwinding of about two turns of the E helix, thus creating a large cavity. Figure 5a shows overlays of the LCAD substrate binding cavity and E-helix with those of SCAD, MCAD, and VLCAD, showing that unwinding of the helix due to Pro132 allows expansion of the cavity to accommodate longer-chain fatty acids and the steroid moiety of THC-CoA. In addition, small residues in the immediate vicinity of Pro132 (Cys129, Ser130, Gly131, and Gly133) contribute to widening of the cavity (Fig. 5b).
Stabilization of LCAD tetramer
As with other ACADs, FAD is bound at the interface of the LCAD dimer and binding is dependent upon stability of the tetrameric structure. The K333Q polymorphism has been reported to be associated with decreased enzyme activity, protein stability, and FAD content33. This residue is located in a chain of salt bridges between two monomers that can stabilize the tetramer conformation (Fig. 6). Disruption of this salt bridge may destabilize the tetramer and account for both the decreased FAD content and enzyme activity observed in the K333Q form of human LCAD.
Interaction between ETF and Long-chain acyl-CoA dehydrogenase
In the mammalian mitochondrial matrix, electron transfer flavoprotein (ETF) links the activity of at least nine different acyl-CoA dehydrogenases to the respiratory chain by accepting and subsequently transferring electrons to the membrane bound ETF-ubiquinone oxidoreductase (ETF-QO) (Fig. 1).
An overlay of the LCAD:ETF (not shown) complex structure based on the previously solved structure of MCAD:ETF34 showed that the side chain of βLeu195 of ETF nestles into a hydrophobic pocket formed by LCAD residues Phe60, Gly97, Leu98, Val101, Ile110, Gly112 and Val120, corresponding to the hydrophobic pocket formed by residues Phe 52, Gly89, Leu90, Thr93, Leu102, Leu104 and Ile112, in the MCAD:ETF complex (Figure S2)34. Although the side chain of Phe60 of LCAD is oriented ~ 180° away from that seen in MCAD:ETF complex, it forms close contacts with βLeu195 of ETF. The ionic interaction between Glu212 of MCAD and αArg249 of ETF plays a key role to stabilize electron transfer competent state34. Leu248 is found at the corresponding position in LCAD; however, the nearby residue Glu260 is in a position to form the salt bridge with ETF αArg249 (Supplementary Fig. S2).
Circadian regulation of LCAD and MCAD has been demonstrated and associated with SIRT3-mediated deacetylation35. LCAD Lys318 and Lys322 have been identified as targets of SIRT3 and are conserved in MCAD and LCAD (Fig. S2) as SIRT3 targets36. The location of these residues, on the surface of the protein and in a position to interact with ETF34, provides a structural basis for circadian regulation of LCAD and MCAD via SIRT3-mediated deacetylation.
Expanded binding site cavities in mammalian and bacterial ACADs.
Identification of the role of Pro132 in partial unwinding of the E-helix prompted a search for other ACADs that might share this expanded substrate-binding cavity. Alignment of LCAD with ACAD10, ACAD11, and ACAD12 reveals the presence of a proline at residues 782 and 463, respectively, of ACAD10 and ACAD11 (Supplementary Figure S2). Examination of ACAD11 structure (PDB code:2wbi) reveals unwinding of the helix at Pro463, permitting binding of a bulky substrate such as cholesterol (Fig. 7a). In addition, the ACAD11 binding site is further expanded to accommodate longer-chain substrates, consistent with the reported ability of ACAD11 to oxidize substrates as long as C244. No structure of ACAD10 has been determined, but sequence similarity to LCAD and ACAD11 suggests the presence of a large substrate binding cavity. Although metabolism of 2-methyl C15-CoA by ACAD10 has been demonstrated, metabolism of THC-CoA has not been tested4. ACAD12 does not contain a proline residue at this position, suggesting that it does not share this large substrate binding cavity.
Bacteria, notably Mycobacterium tuberculosis, are known to utilize steroids as carbon sources20,37. Thomas et al have identified a unique α2β2 heterotetrameric ACAD, encoded by separate genes and distinct from the homotetrameric eucaryotic ACADs38, that carry out β-oxidation of the cholesterol side chain. The crystal structures of two of these α2β2 ACADs, ChsE4-ChsE5 (FadE26-FadE27) and ChsE1-ChsE2 (FadE28-FadE29)20, have been determined, showing that the α subunit binds the FAD isoalloazine, contains the catalytic Glu residue, and substrate binding site; while the β subunit interacts with the ADP-ribose moiety of FAD. Comparison of the structure of ChsE4 with that of LCAD demonstrates the large active site cavity, showing helix unwinding at Pro 91 of ChsE4, corresponding to Pro132 of LCAD (Fig. 7b).
Evolutionary origins of LCAD
LCAD is historically classified as a member of the Acyl- CoA dehydrogenase family responsible for mitochondrial β-oxidation of fatty acyl-CoA thioesters, that includes SCAD, MCAD, LCAD, VLCAD, and ACAD9. Phylogenetic analysis has shown that individual eucaryotic ACADs bear greater sequence similarity to bacterial ACADs than to other eucaryotic ACADs, suggesting that clades containing these five eucaryotic ACADs diverged well before the common ancestor of Archaea, Bacteria, and Eucarya39. Although LCAD shares limited sequence identity (24%-34%) with other eucaryotic ACADs, BLAST searching of bacterial and archaeal genomes revealed that an ACAD from Sneathiella chungangensis (WP_161340043.1) displayed the greatest homology to human LCAD, with 59% amino acid sequence identity, while M. tuberculosis contains ACADs with as high as 52% sequence identity. Examination of the aligned sequences of human and M. tuberculosis ACADs (Supplementary Figure S2) reveals a complex evolutionary pattern with preservation of functional residues. The catalytic LCAD Glu291 is conserved in IVD and the M. tuberculosis ACADs. Glu405 of MCAD is conserved in MCAD, VLCAD, SCAD, and IBD, while ACADs 10, 11, and 12 and Mtb4HR3 contain aspartate at the position corresponding to MCAD Glu405. Residues forming the hydrophobic ETF binding pocket are largely conserved in all the ACADs. The salt bridge with ETF αArg249 is preserved in human ACADs, with conservation of LCAD Glu260 in LCAD, IVD, and ACADs 10–12 and conservation of MCAD Glu241 in MCAD, SCAD, and IBD. Interestingly, Glu is found at both positions in VLCAD and ACAD9.
Pro132, responsible for helix unwinding and expansion of the substrate binding cavity, is conserved in LCAD and the four M. tuberculosis ACADs, including MtbChsE4 and MtbChsE2 that are known to metabolize cholesterol (Fig. 7b and Supplementary Figure S2, and Yang et al.20). Although substrate specificities of ACAD10, ACAD11, Mtb4HR3_A, and MtbCNF74574 are unknown, the presence of the proline residue in the E helix suggests that they may also metabolize sterol substrates.
Physiological function of LCAD
All reported human cases of long chain fatty acid dehydrogenase deficiency have been linked to VLCAD mutations and no human case of LCAD deficiency has been reported. However, although LCAD substrate specificity overlaps with those of VLCAD and MCAD, the minor allele frequency of the E291K polymorphism, resulting in elimination of the catalytic base, is very low (1.10e-5, https://gnomad.broadinstitute.org/variant/2-211068168-C-T), suggesting an essential function for LCAD that cannot be replaced by VLCAD or MCAD. A noncoding polymorphism in the ACADL gene has been associated with alterations in serum 2,6-dimethylheptanoic acid levels40,41. The ability of LCAD to accommodate bulky substrates such as branched chain fatty acids and sterol substrates is not shared by SCAD, MCAD, VLCAD, or ACAD9, suggesting that metabolism of one of these compounds is an essential function of LCAD.
Fasting-induced hypoketotic hypoglycemia and cardiac hypertrophy seen in both LCAD- and VLCAD-null mice are similar to symptoms seen in human VLCAD deficiency. LCAD knockout mice also exhibit fasting-induced hepatosteatosis and cardiac dysfunction in association with altered branched-chain amino acid metabolism, decreased anaplerosis and activation of the integrated stress response12,15. Although VLCAD mRNA levels are 85 times that of LCAD in human heart and skeletal muscle2, suggesting that VLCAD is primarily responsible for long chain fatty acid metabolism in these tissues, LCAD has been identified as a modulator of cardiac remodeling in a mouse model of stress-induced hypertrophy42. Circadian regulation of LCAD via Sirt3-mediated deacetylation suggests a role in energy production in response to feeding and fasting cycles35,36,43 and in vitro studies have implicated LCAD as a tumor suppressor in human cancer17–19,44,45.
In human lung, LCAD is localized to ATII cells responsible for synthesis and secretion of pulmonary surfactant, a mixture of phospholipids, cholesterol, and surfactant proteins, and LCAD knockout mice exhibit surfactant dysfunction, decreased lung compliance and increased susceptibility to influenza infection7. Although regulation of cholesterol content is essential for surfactant structure, mechanisms of cholesterol synthesis, uptake and clearance are not well understood46. Single cell transcriptomic analysis of pulmonary fibrosis has identified a cholesterol metabolic process and downregulation of LCAD47, raising the possibility that sterol metabolism may be a function of LCAD in alveolar ATII cells.
In conclusion, we present structural evidence showing that LCAD is a prototype of a distinct class of eucaryotic ACADs that has evolved an enlarged substrate-binding cavity suitable for β oxidation of bulky substrates including branched chain fatty acyl-CoAs and sterol derivatives. This structure, where the presence of proline in the E helix leads to formation of a large active-site cavity, has also been found in steroid-metabolizing ACADs in bacteria and ACAD11 in eukaryotes. Sequence similarity suggests that substrate binding site of ACAD10 may share this structure. The expanded substrate specificity of LCAD raises the possibility of multiple functions for this enzyme in normal physiology, as well as a mechanism for metabolic switching during pathological states such as cardiomyopathy and tumorigenesis.