Jingyu Zhan
Di Xia- Laboratory of Cell Biology, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, MD, United States
The mitochondrial respiratory chain has long been a primary target for the development of fungicides for its indispensable role in various cellular functions including energy metabolism. Over the years, a wide range of natural and synthetic fungicides and pesticides targeting the respiratory chain complexes have been discovered or developed and used in agriculture and in medicine, which brought considerable economic gains but was also accompanied by the emergence of resistance to these compounds. To delay and overcome the onset of resistance, novel targets for fungicides development are actively being pursued. Mitochondrial AAA protein Bcs1 is necessary for the biogenesis of respiratory chain Complex III, also known as cyt bc1 complex, by delivering the last essential iron-sulfur protein subunit in its folded form to the cyt bc1 precomplex. Although no report on the phenotypes of knock-out Bcs1 has been reported in animals, pathogenic Bcs1 mutations cause Complex III deficiency and respiratory growth defects, which makes it a promising new target for the development of fungicides. Recent Cryo-EM and X-ray structures of mouse and yeast Bcs1 revealed the basic oligomeric states of Bcs1, shed light on the translocation mechanism of its substrate ISP, and provided the basis for structure-based drug design. This review summarizes the recent progress made on understanding the structure and function of Bcs1, proposes the use of Bcs1 as an antifungal target, and provides novel prospects for fungicides design by targeting Bcs1.
1 Introduction
1.1 Mitochondrial respiratory chain components as effective drug targets
Mitochondrial respiratory chain components are validated targets for disease treatment in medicine, and fungal and pest controls in agriculture by natural or synthetic compounds (Fisher et al., 2020). The respiratory Complex III, for example, has been targeted for treatment of malarial parasites and Pneumocystis jiroveci [Pneumocystis carinii] pneumonia with atovaquone (Baggish and Hill, 2002). Complex III is also widely targeted by antifungals in agriculture with synthetic azoxystrobin and famoxadone (Gao et al., 2002; Esser et al., 2004; Esser et al., 2019). The effective disease control by targeting the respiratory chain components reflects their fundamental and indispensable role in cellular function. However, like all chemical therapies, resistance to these treatments is inevitable, rendering drugs useless. Studies of resistance mechanisms indicated that pathogenic organisms are equipped with a plethora of tools to evade attacks by drugs. These tools include mutations to target drug binding site, deploying a large number of ABC transporters that have broad substrate specificity on the cell surface to fend off invading drugs, reducing the intracellular concentration of the drug by degradation, chemical modification, and storage into intracellular vesicles, and going into hibernation (Prasad et al., 2016). To circumvent or delay the onset of resistance, various approaches have been used such as designing derivatives of existing compounds, using a combination of different drugs against several targets simultaneously, and identifying new targets.
1.2 New targets are needed to overcome drug resistance
Despite the effort, drug resistance remains a major issue in treating microbial infections. New targets are actively being sought after. Identification of new targets is by no means easy because the targets must be verified extensively for its lethality to the pathogens and its safety to the hosts. One approach to efficiently identify new targets is to find proteins that are connected and essential to the validated targets, which includes protein assembly factors, translocations factors, and quality control factors. Inactivation of these factors adversely affect the structural and functional integrity of validated targets, thus achieving the goal of disease control. In this article, we review a mitochondrial protein factor Bcs1 that is essential to the assembly of respiratory Complex III and its potential as a drug target.
1.3 Cellular protein translocation pathways that are required for Complex III assembly
Nearly all human proteins are encoded by the nuclear genes and synthesized in the cytosol, except for 13 proteins that are encoded by mitochondrial genome and produced by the mitochondrial protein synthesis machinery (Boengler et al., 2011). A significant portion of nuclear-encoded proteins has to be either translocated across or inserted into the membrane through complex and energy-dependent translocation machineries. Most proteins are exported in unfolded states through the universal Sec pathway to their designated locations where they subsequently become folded to acquire their functional states (Park and Rapoport, 2012). As the main route for protein translocation, components of the Sec apparatus are evolutionarily conserved, as found in Sec61 in eukaryotes (Deshaies and Schekman, 1987) and SecYEG complex in bacteria (Bieker et al., 1990) and chloroplasts (Laidler et al., 1995). Nuclear encoded mitochondrial proteins are imported as unfolded precursors via the TOM (translocase of the outer membrane) and TIM (translocase of the inner membrane) complexes (Wasilewski et al., 2017). However, some proteins need to be translocated in their folded states, which in turn requires specialized translocons. This can happen when specific cofactors (Brandt et al., 1993; Lill, 2009) or metal ions (Tottey et al., 2008) have to be inserted into the protein active sites before translocation, or when a multi-subunit protein complex is directed to the translocation system by the signal sequence only present in one subunit of the complex (Sauve et al., 2007). In bacteria, this highly unusual task of transporting folded protein is carried out by the twin arginine translocation (Tat) pathway (Palmer and Berks, 2012). For example, some bacterial periplasmic proteins contain non-covalently bound redox cofactors that can only be inserted in the cytosol prior to translocation, and they are exported to the periplasm by the Tat system in the folded form with the cofactor installed.
Unlike Sec pathway, the Tat system is not ubiquitous, found in most, but not all bacteria and archaea (Palmer and Berks, 2012). The Tat system is also found in the thylakoid membrane of chloroplasts in plants and is essential for membrane translocation of many substrates including the assembly of the cytochrome cyt b6f complex in the photosynthetic pathway by translocating the folded Rieske iron-sulfur protein (ISP) with the iron-sulfur cluster (Fe2S2) installed (Kurisu et al., 2003). The Tat apparatus is broadly conserved between prokaryotes and chloroplasts, and requires minimally three integral membrane proteins to function: TatA, TatB and TatC (Rollauer et al., 2012; Ramasamy et al., 2013; Rodriguez et al., 2013; Zhang et al., 2014; Palmer and Stansfeld, 2020). The Tat pathway substrates share a consensus signal sequence at the N-terminus featuring a critical pair of consecutive Arginine residues, which is recognized by TatC, the core component the Tat system (Ramasamy et al., 2013). Substrate translocation by the Tat apparatus is fueled by the cross-membrane proton motive force (PMF). Although the structure for the TatC pore (Rollauer et al., 2012; Ramasamy et al., 2013) and the helical regions of TatA (Rodriguez et al., 2013) and TatB (Zhang et al., 2014) have been experimentally determined individually, the stoichiometry of these components in the active Tat complex is not established and the structure of the substrate-activated Tat complex is not known. The mechanism by which folded substrates with various sizes get translocated across biomembranes without compromising membrane integrity remains elusive. By contrast, the Tat pathway is no longer preserved in fungal and animal mitochondria, although the need to translocate the folded ISP subunit for the assembly of the respiratory chain cyt bc1 complex (also known as Complex III), a homolog to the cyt b6f in mitochondria, is still present. It became clear a decade ago that an AAA (ATPases Associated with diverse cellular Activities) protein Bcs1 was identified to perform this role in mitochondria in substitution for the Tat pathway (Fernandez-Vizarra et al., 2007; Wagener et al., 2011).
2 Function of mitochondrial AAA protein Bcs1
2.1 Bcs1 translocates folded ISP across the mitochondrial inner membrane for the assembly of cyt bc1 complex
The cyt bc1 complex is the central segment of the respiratory chain that performs reaction of electron transfer from NADH to molecular oxygen, and couples this reaction to proton pumping across the mitochondrial inner membrane (IM), creating a Proton Motive Force (PMF) across the membrane that fuels ATP synthesis and is used by many other cellular activities. Depending on the biological complexity, cyt bc1 complex can have as few as the only three essential subunits in bacteria (Yang and Trumpower, 1986), or up to 11 subunits in vertebrates (Xia et al., 1997). The three essential subunits are cytochrome c1 (cyt c1, or Cyt1 in yeast), cytochrome b (cyt b), and the Rieske iron-sulfur protein (ISP or Rip1 in yeast), all equipped with redox cofactors and involved in electron transfer. For eukaryotic cyt bc1 complex, cyt b is the only subunit that is coded by the mitochondrial DNA and plays a nucleating role for the cyt bc1 complex biogenesis (OBrien et al., 2009; Needs et al., 2021). Upon translation in mitochondrial matrix, the nascent cyt b peptide is co-translationally inserted into the IM, and the redox cofactor is subsequently installed from the matrix side (Smith et al., 2012). Nuclear-encoded cyt c1 is imported into the matrix via the TOM-TIM23 complex as an unfolded precursor (Neupert, 2015), which matures through IM insertion and heme c1 co-factor incorporation at intermembrane space (IMS) (Smith et al., 2012). Cyt b, cyt c1 and other non-catalytic subunits, also called supernumerary subunits, form the core of the bc1 precomplex, which is still inactive due to the lack of the essential subunit ISP (Figure 1).
FIGURE 1. Bcs1 is involved in the late stage of bc1 complex biogenesis. The biogenesis of yeast bc1 complex (3CX5.pdb) is shown here to represent the overall process for eukaryotic bc1 complex. It starts with the translation of mitochondrial DNA encoded cyt b. Via stepwise addition of several non-catalytic subunits and catalytic cyt c1, the biogenesis of bc1 proceeds from early precomplex to late precomplex. ISP is the last catalytic subunit to be assembled into the bc1 complex, which is delivered in its folded form to bc1 precomplex by the translocase Bcs1.
ISP is the third catalytic subunit of cyt bc1 complex, with a N-terminal unstructured random coil integrated into the core region of the cyt bc1 complex, a middle transmembrane helix segment anchored to the IM, and a globular C-terminal domain (ISP-ED, extrinsic domain) remaining in the IMS (Xia et al., 1997). The iron-sulfur cluster (Fe2S2) is embedded in the ISP-ED and mediates electron transfer between ubiquinol and cyt c1 (Xia et al., 1997). Due to the matrix localization of the Fe2S2 cluster biosynthesis machinery (Lill, 2009), nuclear encoded ISP is first imported into the mitochondrial matrix as a precursor peptide via the TOM-TIM23 complex (Neupert, 2015), which subsequently becomes folded upon installation of the Fe2S2 cluster in the matrix. Then, the fully folded ISP-ED of the ISP precursor, stabilized additionally by a matrix protein Mzm1 (Cui et al., 2012), is translocated across the IM to IMS for the assembly of the cyt bc1 complex (Figure 1). In bacteria and plant chloroplasts where the Tat system is preserved, folded ISP-ED is exported to the periplasm and the intrathylakoid space, respectively, by the Tat system (Aldridge et al., 2008). In mitochondria, the AAA protein Bcs1 is shown to replace the Tat system to perform this role.
As the central component of the respiratory chain, the cyt bc1 complex has long been a popular target for the development of antibiotics (Fisher et al., 2020). Natural compounds such as strobilurin (Gewehr and Sauter, 2019) and synthetic compounds such as azoxystrobin (Esser et al., 2004; Esser et al., 2019) have been used world-wide as agricultural fungicides, which brought enormous economic gains. However, there are major concerns about these fungicides over emerging drug resistance (Lucas et al., 2015), as well as the toxicity towards the non-target terrestrial animals and aquatic lives. The sequence and structure of the catalytic subunits in cyt bc1 complex are highly conserved throughout evolution (Table 1), which makes it very difficult to design inhibitors with high selectivity. Furthermore, increasing number of fungal pathogens are developing resistance to the current inhibitors by evolving mutations in the drug target, which gets selected and amplified upon fungicide usage (Lucas et al., 2015). To delay and overcome the onset of resistance, novel targets with potentially better inhibitor specificity for fungicides development are needed.
TABLE 1. Sequence conservation to human homologue.
2.2 Bcs1 mutations are linked to cyt bc1 complex deficiencies
Bcs1 was first identified as an integral membrane protein in yeast mitochondria, for its indispensable role in incorporating ISP subunit into cyt bc1 complex (Nobrega et al., 1992). Bcs1 mutant in yeast led to a selective loss of ISP subunit in cyt bc1 complex and consequent enzymatic defect, which could be rescued/reverted by overexpressing Bcs1 but not ISP (Nobrega et al., 1992; Wagener et al., 2011). Further investigations into the ISP biogenesis pathway found that Bcs1 acts as a translocase that could recognize and translocate folded ISP-ED across the mitochondrial IM for it to assemble into the cyt bc1 complex (Wagener et al., 2011). Orthologs of Bcs1 were identified in nearly all eukaryotes, and similar cyt bc1 complex deficiencies were found to be linked with Bcs1 ortholog mutants (De Meirleir et al., 2003).
Human ortholog to Bcs1 is encoded by BCS1L (Bcs1-like) gene, and pathogenic BCS1L variants cause autosomal recessive inheritance with different prognoses (Hinson et al., 2007; Olahova et al., 2019). Infants born with the GRACILE syndrome show characteristic symptoms including growth retardation, aminoaciduria, cholestasis, iron overload, lactic acidosis, and early death. Because of its severeness, newborns develop liver and kidney failure quickly and could not survive for over a few months (Akduman et al., 2016; Guo et al., 2022). In milder cases such as the Björnstad syndrome, reported phenotypes include twisted hair shafts (pili torti) and hearing loss (Hinson et al., 2007). Clinical manifestations caused by pathogenic BCS1L variants can land anywhere between the two extremes.
Among various mitochondria diseases, BCS1L mutations remain to be the most frequent causes of mitochondrial respiratory chain deficiency (Fernandez-Vizarra and Zeviani, 2015), which in turn validates the indispensable role Bcs1 plays in the assembly of the cyt bc1 complex. If the Bcs1 activity of a eukaryotic pathogen could be selectively inhibited, the growth of the pathogen would be severely impaired due to the loss of cyt bc1 function, making Bcs1 a potential target for eukaryotic pathogens including fungus. The recent structural studies of mouse and yeast Bcs1 (Kater et al., 2020; Tang et al., 2020; Xia, 2021) revealed differences that laid the foundation to identify vulnerabilities in the enzyme suitable for structure-based design of selective drugs.
3 Structure of Bcs1
Bcs1 belongs to the AAA-ATPase protein family that features a conserved AAA-ATPase domain at the C-terminus and a unique N-terminal region for its specific function (Figures 2A, B). The N-terminal region could be further divided into a single transmembrane helix (TM) that anchors Bcs1 on the mitochondrial IM, and a Bcs1-specific domain that resides in the matrix with the ATPase domain. The C-terminal ATPase domain bears sequence similarity to AAA family proteins, with a typical nucleotide binding pocket formed at the interface between the large subdomain (AAA-LD) and small subdomain (AAA-SD). The AAA-ATPase domain harness the energy released from ATP hydrolysis to power the specific function of the protein (Gates and Martin, 2020).
FIGURE 2. Overview of heptameric Bcs1 in Apo state. (A) Schematic diagram showing the domain organization of Bcs1. EM density map of yeast Bcs1 shown in side view (B) and top view (C) from the matrix side. (D) Cut side view of yBcs1 in Apo 1 state with a conical cavity in the mitochondrial inner membrane and a cylindric cavity in the matrix. The two cavities are separated by the middle seal loops (EMD-10193). (E) Cut side view of yBcs1 in Apo 2 state where the middle seal loops were pointing towards the matrix thus connecting the two cavities (EMD-10194). The cylindric cavity in the matrix can accommodate a folded ISP-ED (shown in yellow), as is shown in cut side view (F) and top view (G).
Recently, structures of mouse and yeast Bcs1 (mBcs1, yBcs1) have been determined using Cryo-EM and X-ray crystallography in different nucleotide bound states (Kater et al., 2020; Tang et al., 2020), and both form a homo-heptameric ring with a central pore regardless of the nucleotide states (Figure 2; Figure 3). The rare heptameric assembly of Bcs1 is quite unique compared to the canonical AAA proteins which mostly form homo-hexameric rings. The yeast and mouse Bcs1 overall structures closely resemble each other with differences in certain aspects, thus the features discussed below will apply to both unless otherwise specified.
FIGURE 3. Nucleotide-dependent conformational changes of Bcs1. Ribbon diagram of mBcs1 structures determined in Apo, ATPγS and ADP states respectively, from left to right. (A) Side view of the heptameric Bcs1 ring. (B) Top view from the matrix side of the heptameric association of Bcs1-specific region and AAA region. (C) Isolated Bcs1 protomer models in different nucleotide states highlighting the domain organization. Bcs1 can be divided into TM region, Bcs1-specific region and the AAA region, and the AAA region can be further divided into a large subdomain (AAA-LD) and a small subdomain (AAA-SD). The nucleotide binding pocket lies in the interface between AAA-LD and AAA-SD.
3.1 Structures of Apo Bcs1
In Apo state, the heptameric assembly of Bcs1 has an overall shape of a mushroom when viewed from the matrix side along the central 7-fold axis. The molecule can be divided into three regions along the symmetry axis: the TM region, the Bcs1-specific region and the AAA region; they are similar in height and each measured about 30–40 Å. Within the Bcs1-specific region that connects the AAA region, there is a less-dense part, about 15 Å in height, which is defined as interstitial gap (Figure 2B) (Tang et al., 2020; Xia, 2021). At the mushroom cap, there is a large cylindrical cavity formed by the ATPase domain accessible from the matrix side (Figures 2C–E). This cylindrical cavity has an average diameter of ∼35 Å, a depth of 40 Å, and an estimated size of 28,000 Å3 (Tang et al., 2020), which is large enough to accommodate a fully folded ISP-ED (Figures 2F, G). This is consistent with the observation that Bcs1 can interact with ISP-ED in the absence of nucleotides (Wagener et al., 2011; Wagener and Neupert, 2012). On the opposite side of the matrix cavity, there is a smaller conical cavity encircled by the TM helices in the IM. In yBcs1, the TM helices extend through the entire IM and intersect at the intermembrane space side (IMS), thus seal the conical cavity from the IMS side. While in Apo mBcs1 (PDB:6UKP), only half of the TM helices were modelled (Figure 3). The two cavities are separated by the middle seal loops in the Bcs1-specific domain, reaching to the center of the pore.
It is worth noting that there are two Apo structures of yBcs1 obtained from 3D classification and refinement of one Cryo-EM dataset, named Apo1 and Apo2 respectively; their major differences lie in the arrangement of the seal loop in the Bcs1-specific region. Similar to mBcs1, yBcs1 in Apo1 state has the 7 seal loops tightly packed together on the membrane plane that separates the two cavities (Figure 2D), while yBcs1 in Apo2 state has the seal loops pointing away from membrane plane towards the matrix, which consequently opens up the middle seal pore and connects the two cavities (Figure 2E) (Kater et al., 2020).
3.2 Structure of Bcs1 with bound ADP
It has been shown that Bcs1 migrates differently in different nucleotide states on a Blue-Native PAGE (Hinson et al., 2007; Wagener et al., 2011; Tang et al., 2020), with the Apo state and the uniform ADP-bound state migrating slower compared to the uniform ATPγS-bound state, suggesting nucleotide-dependent conformational changes that alter its apparent molecular size. The distinct conformations were revealed in detail when mBcs1 structures were determined with uniform ATPγS bound (PDB:6UKS), and ADP bound (PDB:6UKO) (Tang et al., 2020). The yBcs1 structure with uniform ADP bound was also determined (PDB:6SH3) (Kater et al., 2020).
For both mBcs1 and yBcs1, their structures for the ADP state are very similar to those of Bcs1 in the Apo state (Apo1 state for yBcs1), which features a large substrate-binding cavity in the AAA region opening towards the matrix, compatible with ISP-ED binding; and an interstitial gap between the Bcs1-specific region and the AAA region. The substrate-binding cavity is tightly sealed from the IM side by the middle seal loops in the Bcs1-specific domain. The TM helices in mBcs1 are partially built, while in yBcs1 they are complete. The close resemblance between Bcs1 in Apo state and in ADP state also explains their similar migration patterns in the Blue-native PAGE and their ability to bind ISP (Hinson et al., 2007; Wagener et al., 2011; Tang et al., 2020).
3.3 ATPγS binding induced conformational changes in Bcs1
The ATPγS-bound state is only available for mBcs1 (Tang et al., 2020) (Figure 3). Upon interaction with ATPγS, mBcs1 underwent conformational changes to a more compact bullet shape compared to the Apo/ADP state. While the TM helices could not be determined and were not modelled, the overall height of the matrix region (Bcs1-specific domain and AAA domain) reduced from 70 Å in the Apo/ADP state to 63 Å in the ATPγS state (Figure 3A). Especially, the interstitial gap between Bcs1-specific region and the AAA region disappeared. The cylindrical cavity formed by the ATPase domain shrank by two thirds in size to a diameter of 20 Å and a depth of 25 Å (Figure 3B), which could no longer accommodate the folded ISP-ED (Tang et al., 2020). Nonetheless, it has been shown that in the presence of AMP-PNP, Bcs1 could bind the N-terminal segment of the ISP presumably after the translocation of the ISP-ED (Wagener et al., 2011). The fact that in the ATPγS-bound state, the TM helices could not be modelled due to poor density suggests a lack of interaction in the TM region, which may open up the seal to the IMS (Figures 3A, C) and allow ISP-ED to pass through. Thus, the ATPγS state is deemed as a post-translocation state (Tang et al., 2020).
4 The underlying translocation mechanism
The translocation of ISP-ED is accomplished through ATP hydrolysis-powered conformational changes of Bcs1. The biochemical and structural data of Bcs1 in different nucleotide states provided valuable insights into its translocation mechanism. Unlike the canonical AAA proteins with hexameric association, the heptameric assembly affords Bcs1 a large enough cavity for substrate entry and translocation. Based on biochemical and structural analyses, there are two potential binding sites on Bcs1 for ISP, one is in the AAA domain that is selectively occupied by ISP-ED from the matrix side when Bcs1 is in Apo/ADP state, the other one is present inside the IM accessible from the IMS side that interacts with the N-terminus of ISP after the translocation of ISP-ED, when Bcs1 is in the ATPγS state (Wagener et al., 2011; Kater et al., 2020; Tang et al., 2020).
Therefore, a complete translocation cycle of ISP-ED by Bcs1 could be depicted as follows: (1) In Apo/ADP state, fully folded ISP-ED in mitochondria matrix enters the large substrate-binding cavity in Bcs1 AAA domain. (2) Upon interaction with ATP, Bcs1 protomers go through conformational changes to ATPγS state and push ISP-ED to the IMS, while the N-terminal region of ISP remains in the Bcs1 central pore. (3) Subsequent ATP hydrolysis resets the Bcs1 back to the ADP/Apo state, which also abolishes its interaction with the N-terminal region of ISP and releases the entire ISP for the assembly of the cyt bc1 complex (Figure 4).
FIGURE 4. Model for translocation of ISP by Bcs1. The TM region, Bcs1-specific domain, AAA-LD and AAA-SD of Bcs1 are colored in cyan, green, purple and blue, respectively. In a translocation cycle, Bcs1 populates two major conformations: Apo/ADP state and ATP state. (1) In Apo/ADP state, fully folded ISP-ED in mitochondria matrix enters the substrate-binding cavity in Bcs1 AAA domain. (2) Upon interaction with ATP, Bcs1 undergoes conformational change to the more compact ATP state. The substrate-binding cavity contracts and pushes ISP-ED to the IMS, while the N-terminal region of ISP remains in the Bcs1 central pore. (3) Subsequent ATP hydrolysis resets Bcs1 back to the ADP/Apo state, which also abolishes the interaction between Bcs1 central pore and the N-terminus of ISP. The entire ISP is released from Bcs1 for the assembly of bc1 complex.
Despite all the progress made on understanding the biogenesis of the mitochondrial cyt bc1 complex and the role Bcs1 plays in it, there are still important questions unanswered. First of all, it remains unclear how the seven subunits of Bcs1 coordinate with one another during the translocation process. For the substrate-binding cavity in the AAA region to accommodate the folded ISP-ED, Bcs1 has to be in either uniform ADP state or Apo state. However, in subsequent steps where ATP binding and hydrolysis cause major conformational changes of Bcs1, it is unclear if the transition occurs simultaneously among the seven subunits or sequentially around the ring. Canonical hexameric AAA proteins such as p97 (Cooney et al., 2019; Pan et al., 2021) and ClpAP (Lopez et al., 2020) follow the sequential mechanism where subunits within the ring hydrolyze ATP one by one, causing co-existing nucleotide states and a spiral conformation (Gates and Martin, 2020). Considering Bcs1 is a non-canonical AAA protein with heptameric assembly and a folded protein substrate, it remains to be seen if it is an exception. Secondly, how the entire ISP gets released from the Bcs1 at the end of the translocation cycle remains unknown. Post translocation of ISP-ED, the N-terminal region of ISP remains in the Bcs1 central pore. The N-terminal region of ISP is composed of a TM helix which is about 35 residues long, and a matrix region which could have 30 residues in yeast or over 100 residues in human. It has to either exit through the central pore of Bcs1 to the IMS and re-insert into the IM, or diffuse along the IM after the Bcs1 ring laterally opens up. Thirdly, the proton motive force across the IM generated by mitochondrial respiration must be preserved at all times during the translocation process. The detailed mechanism of ISP translocation and release by Bcs1 without disrupting the membrane permeability remains to be uncovered.
5 Bcs1 as a potential drug target for microbial infections
Most targets for antibiotics against bacterial infections take advantage of the evolutionary long distance between bacteria and animals, making it one of the most successful treatment options in modern medicine. Microbial infections caused by fungi or oomycetes, both in agriculture and in animals, are difficult to treat as the pathogens are much closer to their hosts evolutionarily and share more common biological pathways. This awkward situation is reflected by the fact that only a handful classes of antifungals are available for use in medicine, which is in sharp contrast to the numerous antibiotics against bacterial infections. Therefore, a successful drug candidate should meet many criteria including the following two: (1) its inhibition leads to irreversible harm to the pathogen, and (2) ideally, the target should be unique to the pathogen to avoid off-target effect to its hosts.
5.1 Bcs1 is connected to the function of the cellular respiratory chain
Although a complete inactivation of Bcs1 in yeast was accomplished, allowing its role in Complex III assembly to be uncovered, Bcs1 gene knock-out in animals has not been reported. Based on the yeast studies, full inactivation of Bcs1L in animals is likely to be lethal, making it challenging to study in higher organisms. Nevertheless, pathogenic Bcs1L mutations identified in clinical cases (De Meirleir et al., 2003; Hinson et al., 2007; Kasapkara et al., 2014; Fernandez-Vizarra and Zeviani, 2015; Olahova et al., 2019; Guo et al., 2022) point to the vulnerable regions on Bcs1 that are important to maintain its proper function. If these regions are perturbed, there might be severe consequences resembling the clinical symptoms of the mutations, which also makes them potential targets for drug intervention. By mapping the pathogenic mutations on mBcs1 structure, which shares 94% sequence identity with human BCS1L (Table 1), it is revealed that the pathogenic mutations are clustered in Bcs1-specific region and the AAA-ATPase region (Figures 5A, B). More specifically, the severe GRACILE syndrome is frequently caused by mutations in Bcs1-specific region, while the milder Björnstad syndrome is linked to mutations in the AAA region. Combined with the structural information on Bcs1 and available phenotypes of pathogenic mutations, it has been proposed that the structural integrity of Bcs1-specific region is critical for persevering the proton motive force across the IM, and the fetal GRACILE syndrome could be the result of proton leakage caused by defective Bcs1-specific domain on top of the failed cyt bc1 complex assembly (Tang et al., 2020). Thus, using Bcs1 as a target could potentially interfere with both the activity of Complex III in the respiratory chain and the function of mitochondrial proton motive force across IM.
FIGURE 5. Pathogenic mutations on Bcs1 and potential target for fungicide design. (A, B) Pathogenic mutations of BCS1L are mapped onto the ribbon diagram of mBcs1 as red spheres. The TM region, Bcs1-specific domain, AAA-LD and AAA-SD of Bcs1 are colored in cyan, green, purple and blue, respectively. Cut side view is shown in (A) and Bcs1 protomer is shown in (B). (C–E) Ribbon diagram of mBcs1 colored by sequence conservation, with most conserved residues colored in red and the least conserved residues in blue. The cut side view and bottom view of Bcs1 ring are shown in (C) and (E) respectively, and Bcs1 protomer is shown in (D). (F) Zoom in view of an unconserved pocket in AAA-SD in the close vicinity of nucleotide binding pocket of Bcs1. The bound ADP is shown as yellow sticks. The surface of Bcs1 is also colored by sequence conservation.
5.2 Identifications of unique features that distinguish Bcs1 of pathogens from hosts
Although it is considered unconventional or even dangerous to target mitochondrial respiratory chain because it is a shared and conserved pathway between pathogens and hosts (Table 1), exceptions to this rule exist. For example, the FDA approved anti-malaria drug atovaquone specifically targets the Complex III of the pathogen’s respiratory chain (Baggish and Hill, 2002). Atovaquone is shown to be safe for prophylactic use in humans. Structural study of famoxadone, a Complex III specific antifungal widely used in crop disease control, revealed organism-specific binding to the active site with a IC50 value 400 times smaller against bacterial cyt bc1 when compared to that against bovine mitochondrial complex (Gao et al., 2002). Thus, it is possible to identify a unique locus in the target, which gives a significant difference between pathogen and host even the target is part of the shared pathway.
To find such locations on Bcs1-specific and AAA regions suitable for drug targeting, a multisequence alignment of Bcs1 from difference species is provided (Figure 6) and mapped onto yBcs1 structure (Figures 5C–E). In the close vicinity of the nucleotide-binding site, there is a pocket in the AAA-SD that appears to be highly divergent in the sequence alignment (Figure 5F). Potentially, a molecule that specifically targets the ATPase site while having a tail stretching to explore this pocket in pathogenic fungal Bcs1 could block the ATP binding and hydrolysis, thus serve as a fungicide with sufficient specificity to the pathogen. Alternatively, molecules that bind to the interstitial gap in the Bcs1-specific region may be able to block the overall conformational change of Bcs1 upon ATP binding, thus have the Bcs1 stuck in Apo/ADP conformation and unable to finish its translocation cycle.
FIGURE 6. Multiple sequence alignment of Bcs1. Residues with pathogenic mutations implicated in diseases such as GRACILE syndrome, Björnstad syndrome and bc1 complex deficiencies are marked with a green dot. Protein sequences were retrieved from the Uniprot with the accession numbers shown in parentheses: Human (Q9Y276), Bovine (Q5E9H5), Mouse (Q9CZP5), African clawed frog (Q7ZTL7), Zebrafish (Q7ZV60), House spider (A0A2L2YIU0), African Malaria Mosquito (A0A182RX69), Parasitic wasp (A0A7M7PZM7), Fruit fly (Q9VL22), Baker’s yeast (P32839), Malaria parasite P. vivax (A0A564ZXW0).
Since Bcs1 is found in nearly all eukaryotes, it could also be targeted in other eukaryotic pathogens. For instance, both malaria parasite and its transmission vector African malaria mosquito have Bcs1 orthologues which share a sequence identity of 39% and 60% with human BCS1L, respectively (Table 1; Figure 6). Small molecules that are tailored to specifically target the pathogenic Bcs1 could be designed to control these eukaryotic pathogens.
5.3 Methods that can be employed to identify Bcs1 inhibitors
Currently, there are two in vitro assays of Bcs1 that can be employed for screening inhibitor libraries. One is the ATPase assay that measures ATP hydrolysis by isolated Bcs1. ATPase assay has been widely used to screen inhibitors for other AAA ATPases such as the human AAA protein p97 (Chou et al., 2014), which was quite successful to bring at least one anti-cancer drug to clinical trial (Tang et al., 2019; Xia and Tang, 2020). Another assay that can potentially be used is based on the conformational change induced by binding of ATP or its analogs. The conformational change of Bcs1 can be detected by the Blue-Native PAGE method as illustrated in the previous publications (Hinson et al., 2007; Wagener et al., 2011; Tang et al., 2020), allowing detection of inhibitors that block this conformational transition. These assays can potentially be automated for high throughput screening (HTS).
With the available high-resolution Bcs1 structures of both mBcs1 and yBcs1 in different conformations, we are also in a position to carry out structure-based in silico screen of compound database to select promising candidates. The HTS approach realizes sensitive detection and automatic operation, but its positive rate is low, thus the cost of HTS is relatively high. However, the virtual screening has a relatively high hit rate and inexpensive. By combining the results of virtual screening with experimental screening, the latter can be more targeted, has a higher hit rate, thus shortening the research cycle and reducing the cost of drug development.
Author contributions
All authors listed have made a substantial, direct, and intellectual contribution to the work and approved it for publication.
Funding
This research was supported by the Intramural Research Program (Z01 Number: ZIA BC 010600) of the Center for Cancer Research, National Cancer Institute, National Institutes of Health.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
AAA, ATPases Associated with diverse cellular Activities; Cyt1, yeast cyt c1; cyt c1, cytochrome c1; cyt b, cytochrome b; Fe2S2, iron-sulfur cluster; HTS, high throughput screening; ISP, Rieske iron-sulfur protein; ISP-ED, extrinsic domain of the Rieske iron-sulfur protein; IM, mitochondrial inner membrane; IMS, mitochondrial intermembrane space; OM, mitochondrial outer membrane; PMF, proton motive force; TAT, twin arginine translocation; TIM, translocase of the inner mitochondrial membrane; TM, transmembrane; TOM, translocase of the outer mitochondrial membrane; mBcs1, mouse Bcs1; yBcs1, yeast Bcs1.
References
Akduman, H., Eminoglu, T., Okulu, E., Erdeve, O., Atasay, B., and Arsan, S. (2016). A neonate presenting with gracile syndrome and nornstad phenotype associated with Bcs1l mutation. Genet. Couns. 27 (4), 509–512.
Aldridge, C., Spence, E., Kirkilionis, M. A., Frigerio, L., and Robinson, C. (2008). Tat-dependent targeting of Rieske iron-sulphur proteins to both the plasma and thylakoid membranes in the cyanobacterium Synechocystis PCC6803. Mol. Microbiol. 70 (1), 140–150. doi:10.1111/j.1365-2958.2008.06401.x
Baggish, A. L., and Hill, D. R. (2002). Antiparasitic agent atovaquone. Antimicrob. Agents Ch 46 (5), 1163–1173. doi:10.1128/aac.46.5.1163-1173.2002
Bieker, K. L., Phillips, G. J., and Silhavy, T. J. (1990). The Sec and prl genes of escherichia-coli. J. Bioenerg. Biomembr. 22 (3), 291–310. doi:10.1007/bf00763169
Boengler, K., Heusch, G., and Schulz, R. (2011). Nuclear-encoded mitochondrial proteins and their role in cardioprotection. Bba-Mol Cell Res. 1813 (7), 1286–1294. doi:10.1016/j.bbamcr.2011.01.009
Brandt, U., Yu, L., Yu, C. A., and Trumpower, B. L. (1993). The mitochondrial targeting presequence of the Rieske iron-sulfur protein is processed in a single step after insertion into the cytochrome-bc1 complex in mammals and retained as a subunit in the complex. J. Biol. Chem. 268 (12), 8387–8390. doi:10.1016/s0021-9258(18)52883-0
Chou, T. F., Bulfer, S. L., Weihi, C. C., Li, K. L., Lis, L. G., Walters, M. A., et al. (2014). Specific inhibition of p97/VCP ATPase and kinetic analysis demonstrate interaction between D1 and D2 ATPase domains. J. Mol. Biol. 426 (15), 2886–2899. doi:10.1016/j.jmb.2014.05.022
Cooney, I., Han, H., Stewart, M. G., Carson, R. H., Hansen, D. T., Iwasa, J. H., et al. (2019). Structure of the Cdc48 segregase in the act of unfolding an authentic substrate. Science 365 (6452), 502–505. doi:10.1126/science.aax0486
Cui, T. Z., Smith, P. M., Fox, J. L., Khalimonchuk, O., and Winge, D. R. (2012). Late-stage maturation of the Rieske Fe/S protein: Mzm1 stabilizes Rip1 but does not facilitate its translocation by the AAA ATPase Bcs1. Mol. Cell Biol. 32 (21), 4400–4409. doi:10.1128/mcb.00441-12
De Meirleir, L., Seneca, S., Damis, E., Sepulchre, B., Hoorens, A., Gerlo, E., et al. (2003). Clinical and diagnostic characteristics of complex III deficiency due to mutations in the BCS1L gene. Am. J. Med. Genet. A 121A (2), 126–131. doi:10.1002/ajmg.a.20171
Deshaies, R. J., and Schekman, R. (1987). A yeast mutant defective at an early stage in import of secretory protein precursors into the endoplasmic-reticulum. J. Cell Biol. 105 (2), 633–645. doi:10.1083/jcb.105.2.633
Esser, L., Quinn, B., Li, Y. F., Zhang, M. Q., Elberry, M., Yu, L., et al. (2004). Crystallographic studies of quinol oxidation site inhibitors: A modified classification of inhibitors for the cytochrome bc(1) complex. J. Mol. Biol. 341 (1), 281–302. doi:10.1016/j.jmb.2004.05.065
Esser, L., Zhou, F., Yu, C. A., and Xia, D. (2019). Crystal structure of bacterial cytochrome bc(1) in complex with azoxystrobin reveals a conformational switch of the Rieske iron-sulfur protein subunit. J. Biol. Chem. 294 (32), 12007–12019. doi:10.1074/jbc.ra119.008381
Fernandez-Vizarra, E., Bugiani, M., Goffrini, P., Carrara, F., Farina, L., Procopio, E., et al. (2007). Impaired complex III assembly associated with BCS1L gene mutations in isolated mitochondrial encephalopathy. Hum. Mol. Genet. 16 (10), 1241–1252. doi:10.1093/hmg/ddm072
Fernandez-Vizarra, E., and Zeviani, M. (2015). Nuclear gene mutations as the cause of mitochondrial complex III deficiency. Front. Genet. 6, 134. doi:10.3389/fgene.2015.00134
Fisher, N., Meunier, B., and Biagini, G. A. (2020). The cytochrome bc(1) complex as an antipathogenic target. Febs Lett. 594 (18), 2935–2952. doi:10.1002/1873-3468.13868
Gao, X. G., Wen, X. L., Yu, C. A., Esser, L., Tsao, S., Quinn, B., et al. (2002). The crystal structure of mitochondrial cytochrome bc1 in complex with famoxadone: The role of aromatic-aromatic interaction in inhibition. Biochemistry-Us 41 (39), 11692–11702. doi:10.1021/bi026252p
Gates, S. N., and Martin, A. (2020). Stairway to translocation: AAA+ motor structures reveal the mechanisms of ATP-dependent substrate translocation. Protein Sci. 29 (2), 407–419. doi:10.1002/pro.3743
Gewehr, M., and Sauter, H. (2019). “Strobilurins and other complex III inhibitors,” in Modern crop protection compounds (United States: Wiley), 634–681.
Guo, W. C., Shao, Y. F., Lang, Y. H., Wang, H., Lin, Y., Liu, X. Y., et al. (2022). Identification of two novel variants of <i>BCS1L</i> gene in a patient with classical GRACILE syndrome. Nephrology 27 (10), 810–814. doi:10.1111/nep.14086
Hinson, J. T., Fantin, V. R., Schonberger, J., Breivik, N., Siem, G., McDonough, B., et al. (2007). Missense mutations in the BCS1L gene as a cause of the Bjornstad syndrome. New Engl. J. Med. 356 (8), 809–819. doi:10.1056/nejmoa055262
Kasapkara, C. S., Tumer, L., Ezgu, F. S., Kucukcongar, A., and Hasanoglu, A. (2014). BCS1L gene mutation causing GRACILE syndrome: Case report. Ren. Fail. 36 (6), 953–954. doi:10.3109/0886022x.2014.900422
Kater, L., Wagener, N., Berninghausen, O., Becker, T., Neupert, W., and Beckmann, R. (2020). Structure of the Bcs1 AAA-ATPase suggests an airlock-like translocation mechanism for folded proteins. Nat. Struct. Mol. Biol. 27 (2), 142–149. doi:10.1038/s41594-019-0364-1
Kurisu, G., Zhang, H., Smith, J. L., and Cramer, W. A. (2003). Structure of the cytochrome b6f complex of oxygenic photosynthesis: Tuning the cavity. Science 302 (5647), 1009–1014. doi:10.1126/science.1090165
Laidler, V., Chaddock, A. M., Knott, T. G., Walker, D., and Robinson, C. (1995). A SecY homolog in Arabidopsis thaliana. J. Biol. Chem. 270 (30), 17664–17667. doi:10.1074/jbc.270.30.17664
Lill, R. (2009). Function and biogenesis of iron-sulphur proteins. Nature 460 (7257), 831–838. doi:10.1038/nature08301
Lopez, K. E., Rizo, A. N., Tse, E., Lin, J., Scull, N. W., Thwin, A. C., et al. (2020). Conformational plasticity of the ClpAP AAA+ protease couples protein unfolding and proteolysis. Nat. Struct. Mol. Biol. 27 (5), 406–416. doi:10.1038/s41594-020-0409-5
Lucas, J. A., Hawkins, N. J., and Fraaije, B. A. (2015). The evolution of fungicide resistance. Adv. Appl. Microbiol. 90, 29–92. doi:10.1016/bs.aambs.2014.09.001
Needs, H. I., Protasoni, M., Henley, J. M., Prudent, J., Collinson, I., and Pereira, G. C. (2021). Interplay between mitochondrial protein import and respiratory complexes assembly in neuronal Health and degeneration. Life-Basel. 11 (5), 432. doi:10.3390/life11050432
Neupert, W. (2015). A perspective on transport of proteins into mitochondria: A myriad of open questions. J. Mol. Biol. 427 (6), 1135–1158. doi:10.1016/j.jmb.2015.02.001
Nobrega, F. G., Nobrega, M. P., and Tzagoloff, A. (1992). BCS1, a novel gene required for the expression of functional Rieske iron-sulfur protein in Saccharomyces cerevisiae. EMBO J. 11 (11), 3821–3829. doi:10.1002/j.1460-2075.1992.tb05474.x
Obrien, E. A., Zhang, Y., Wang, E., Marie, V., Badejoko, W., Lang, B. F., et al. (2009). Gobase: An organelle genome database. Nucleic Acids Res. 37, D946–D950. doi:10.1093/nar/gkn819
Olahova, M., Berti, C. C., Collier, J. J., Alston, C. L., Jameson, E., Jones, S. A., et al. (2019). Molecular genetic investigations identify new clinical phenotypes associated with BCS1L-related mitochondrial disease. Hum. Mol. Genet. 28 (22), 3766–3776. doi:10.1093/hmg/ddz202
Palmer, T., and Berks, B. C. (2012). The twin-arginine translocation (Tat) protein export pathway. Nat. Rev. Microbiol. 10 (7), 483–496. doi:10.1038/nrmicro2814
Palmer, T., and Stansfeld, P. J. (2020). Targeting of proteins to the twin-arginine translocation pathway. Mol. Microbiol. 113 (5), 861–871. doi:10.1111/mmi.14461
Pan, M., Yu, Y., Ai, H., Zheng, Q., Xie, Y., Liu, L., et al. (2021). Mechanistic insight into substrate processing and allosteric inhibition of human p97. Nat. Struct. Mol. Biol. 28 (7), 614–625. doi:10.1038/s41594-021-00617-2
Park, E., and Rapoport, T. A. (2012). Mechanisms of sec61/SecY-mediated protein translocation across membranes. Annu Rev. Biophys. 41, 21–40. doi:10.1146/annurev-biophys-050511-102312
Prasad, R., Shah, A. H., and Rawal, M. K. (2016). Antifungals: Mechanism of action and drug resistance. Adv. Exp. Med. Biol. 892, 327–349. doi:10.1007/978-3-319-25304-6_14
Ramasamy, S., Abrol, R., Suloway, C. J. M., and Clemons, W. M. (2013). The glove-like structure of the conserved membrane protein TatC provides insight into signal sequence recognition in twin-arginine translocation. Structure 21 (5), 777–788. doi:10.1016/j.str.2013.03.004
Rodriguez, F., Rouse, S. L., Tait, C. E., Harmer, J., De Riso, A., Timmel, C. R., et al. (2013). Structural model for the protein-translocating element of the twin-arginine transport system. P Natl. Acad. Sci. USA 110 (12), E1092–E1101. doi:10.1073/pnas.1219486110
Rollauer, S. E., Tarry, M. J., Graham, J. E., Jaaskelainen, M., Jager, F., Johnson, S., et al. (2012). Structure of the TatC core of the twin-arginine protein transport system. Nature 492 (7428), 210–214. doi:10.1038/nature11683
Sauve, V., Bruno, S., Berks, B. C., and Hemmings, A. M. (2007). The SoxYZ complex carries sulfur cycle intermediates on a peptide swinging arm. J. Biol. Chem. 282 (32), 23194–23204. doi:10.1074/jbc.m701602200
Smith, P. M., Fox, J. L., and Winge, D. R. (2012). Biogenesis of the cytochrome bc(1) complex and role of assembly factors. Biochim. Biophys. Acta 1817 (2), 276–286. doi:10.1016/j.bbabio.2011.11.009
Tang, W. K., Borgnia, M. J., Hsu, A. L., Esser, L., Fox, T., de Val, N., et al. (2020). Structures of AAA protein translocase Bcs1 suggest translocation mechanism of a folded protein. Nat. Struct. Mol. Biol. 27 (2), 202–209. doi:10.1038/s41594-020-0373-0
Tang, W. K., Odzorig, T., Jin, W., and Xia, D. (2019). Structural basis of p97 inhibition by the site-selective anticancer compound CB-5083. Mol. Pharmacol. 95 (3), 286–293. doi:10.1124/mol.118.114256
Tottey, S., Waldron, K. J., Firbank, S. J., Reale, B., Bessant, C., Sato, K., et al. (2008). Protein-folding location can regulate manganese-binding versus copper- or zinc-binding. Nature 455 (7216), 1138–1142. doi:10.1038/nature07340
Wagener, N., Ackermann, M., Funes, S., and Neupert, W. (2011). A pathway of protein translocation in mitochondria mediated by the AAA-ATPase Bcs1. Mol. Cell 44 (2), 191–202. doi:10.1016/j.molcel.2011.07.036
Wagener, N., and Neupert, W. (2012). Bcs1, a AAA protein of the mitochondria with a role in the biogenesis of the respiratory chain. J. Struct. Biol. 179 (2), 121–125. doi:10.1016/j.jsb.2012.04.019
Wasilewski, M., Chojnacka, K., and Chacinska, A. (2017). Protein trafficking at the crossroads to mitochondria. Biochim. Biophys. Acta Mol. Cell Res. 1864, 125–137. doi:10.1016/j.bbamcr.2016.10.019
Xia, D. (2021). Structural snapshots of the cellular folded protein translocation machinery Bcs1. FEBS J. 288 (9), 2870–2883. doi:10.1111/febs.15576
Xia, D., and Tang, W. K. (2020). Structural basis of p97 inhibition by the anti-cancer compound CB-5083. Biophys. J. 118 (3), 49a–a. doi:10.1016/j.bpj.2019.11.446
Xia, D., Yu, C. A., Kim, H., Xia, J. Z., Kachurin, A. M., Zhang, L., et al. (1997). Crystal structure of the cytochrome bc1 complex from bovine heart mitochondria. Science 277 (5322), 60–66. doi:10.1126/science.277.5322.60
Yang, X. H., and Trumpower, B. L. (1986). Purification of a three-subunit ubiquinol-cytochrome c oxidoreductase complex from Paracoccus denitrificans. J. Biol. Chem. 261 (26), 12282–12289. doi:10.1016/s0021-9258(18)67236-9
Keywords: Bcs1, Complex III, cytochrome bc1 complex, iron-sulfur protein, respiratory chain, AAA protein
Citation: Zhan J and Xia D (2023) Bcs1, a novel target for fungicide. Front. Chem. 11:1146753. doi: 10.3389/fchem.2023.1146753
Received: 17 January 2023; Accepted: 03 March 2023;
Published: 13 March 2023.
Edited by:
Peta Harvey, The University of Queensland, AustraliaReviewed by:
Ashutosh Kumar Pandey, The State University of New Jersey, United StatesRashmi Panigrahi, University of Alberta, Canada
Copyright © 2023 Zhan and Xia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Di Xia, xiad@mail.nih.gov