Ins and outs of AlphaFold2 transmembrane protein structure predictions

Transmembrane (TM) proteins are major drug targets, but their structure determination, a prerequisite for rational drug design, remains challenging. Recently, the DeepMind’s AlphaFold2 machine learning method greatly expanded the structural coverage of sequences with high accuracy. Since the employed algorithm did not take specific properties of TM proteins into account, the reliability of the generated TM structures should be assessed. Therefore, we quantitatively investigated the quality of structures at genome scales, at the level of ABC protein superfamily folds and for specific membrane proteins (e.g. dimer modeling and stability in molecular dynamics simulations). We tested template-free structure prediction with a challenging TM CASP14 target and several TM protein structures published after AlphaFold2 training. Our results suggest that AlphaFold2 performs well in the case of TM proteins and its neural network is not overfitted. We conclude that cautious applications of AlphaFold2 structural models will advance TM protein-associated studies at an unexpected level.

Conformational Changes in a Macrolide Antibiotic Binding Protein From Mycobacterium smegmatis Upon ADP Binding

Rv3197 (MABP-1), a non-canonical ABC protein in Mycobacterium tuberculosis, has ATPase activity and confers inducible resistance to the macrolide family of antibiotics. Here we have shown that MSMEG_1954, the homolog of Rv3197 in M. smegmatis, has a similar function of conferring macrolide resistance. Crystal structures of apo-MSMEG_1954 (form1 and form 2) and MSMEG_1954 in complex with ADP have been determined. These three structures show that MSMEG_1954 has at least two different conformations we identify as closed state (MSMEG_1954-form 1) and open state (MSMEG_1954-form 2 and MSMEG_1954-ADP). Structural superimposition shows that the MSMEG_1954-form 2 and MSMEG_1954-ADP complex have similar conformation to that observed for MABP-1 and MABP-1-erythromicin complex structure. However, the antibiotic binding pocket in MSMEG_1954-form 1 is completely blocked by the N-terminal accessory domain. When bound by ADP, the N-terminal accessory domain undergoes conformational change, which results in the open of the antibiotic binding pocket. Because of the degradation of N terminal accessory domain in MSMSG_1954-form 2, it is likely to represent a transitional state between MSMEG_1954-form 1 and MSMEG_1954-ADP complex structure.

Increase in cardioprotective SUR2A does not alter heart rate and heart rate regulation by physical activity and diurnal rhythm

SUR2A is an ABC protein serving as a regulatory subunit of ATP-sensitive (KATP) channels. An increase in SUR2A levels is cardioprotective and it is a potential therapeutic strategy against ischaemic heart disease, heart failure and other diseases. However, whether overexpression of this protein has any adverse effects is yet to be fully understood. Objectives Here, we examined the heart rate and the heart rate diurnal variation in mice overexpressing SUR2A (SUR2A+) and their littermate controls (WT) using ECG telemetry that was continuously recorded for 14 days (days 8–23 post-radiotransmitter implantation). Methods Using SigmaPlot 14.0 and Microsoft Excel, Area Under the Curve (AUC) for each parameter was calculated and plotted in a graph. Results Both WT and SUR2A+ mice were more physically active during nights and there were no significant differences between two phenotypes. Physical activity was associated with increased heart rate in both phenotypes, but there were no differences in heart rate between phenotypes irrespective of physical activity or time of the day. A diurnal heart rate variation was preserved in the SUR2A+ mice. As area under the curve (AUC) analysis has the potential to reveal differences that are invisible with other statistical methods, we compared AUC of heart rate in SUR2A+ and WT mice. This analysis did not yield anything different from traditional analysis. Conclusions We conclude that increased SUR2A levels are not associated with changes in physical activity, heart rate and/or circadian rhythm influence on the heart rate. This lack of adverse effects supports a notion that manipulation with SUR2A levels is a promising cardioprotective strategy.

AlphaFold2 transmembrane protein structure prediction shines

Transmembrane (TM) proteins are major drug targets, indicated by the high percentage of prescription drugs acting on them. For a rational drug design and an understanding of mutational effects on protein function, structural data at atomic resolution are required. However, hydrophobic TM proteins often resist experimental structure determination and in spite of the increasing number of cryo-EM structures, the available TM folds are still limited in the Protein Data Bank. Recently, the DeepMind's AlphaFold2 machine learning method greatly expanded the structural coverage of sequences, with high accuracy. Since the employed algorithm did not take specific properties of TM proteins into account, the validity of the generated TM structures should be assessed. Therefore, we investigated the quality of structures at genome scales, at the level of ABC protein superfamily folds, and also in specific individual cases. We tested template-free structure prediction also with a new TM fold, dimer modeling, and stability in molecular dynamics simulations. Our results strongly suggest that AlphaFold2 performs astoundingly well in the case of TM proteins and that its neural network is not overfitted. We conclude that a careful application of its structural models will advance TM protein associated studies at an unexpected level.

Translation of ABCE1 Is Tightly Regulated by Upstream Open Reading Frames in Human Colorectal Cells

ATP-binding cassette subfamily E member 1 (ABCE1) belongs to the ABC protein family of transporters; however, it does not behave as a drug transporter. Instead, ABCE1 actively participates in different stages of translation and is also associated with oncogenic functions. Ribosome profiling analysis in colorectal cancer cells has revealed a high ribosome occupancy in the human ABCE1 mRNA 5′-leader sequence, indicating the presence of translatable upstream open reading frames (uORFs). These cis-acting translational regulatory elements usually act as repressors of translation of the main coding sequence. In the present study, we dissect the regulatory function of the five AUG and five non-AUG uORFs identified in the human ABCE1 mRNA 5′-leader sequence. We show that the expression of the main coding sequence is tightly regulated by the ABCE1 AUG uORFs in colorectal cells. Our results are consistent with a model wherein uORF1 is efficiently translated, behaving as a barrier to downstream uORF translation. The few ribosomes that can bypass uORF1 (and/or uORF2) must probably initiate at the inhibitory uORF3 or uORF5 that efficiently repress translation of the main ORF. This inhibitory property is slightly overcome in conditions of endoplasmic reticulum stress. In addition, we observed that these potent translation-inhibitory AUG uORFs function equally in cancer and in non-tumorigenic colorectal cells, which is consistent with a lack of oncogenic function. In conclusion, we establish human ABCE1 as an additional example of uORF-mediated translational regulation and that this tight regulation contributes to control ABCE1 protein levels in different cell environments.

The Association of Estrogen Receptor Activity, Interferon Signaling, and MHC Class I Expression in Breast Cancer

Purpose: The expression of major histocompatibility complex class I (MHC I) has previously been reported to be negatively associated with estrogen receptor (ER) expression. Furthermore, the expression of MHC I, level of tumor-infiltrating lymphocytes (TILs), and expression of interferon (IFN) mediator MxA are positively associated with one another in human breast cancers. This study aimed to investigate the mechanisms of association of MHC I with estrogen and IFN signaling. Methods: The human leukocyte antigen (HLA)-ABC protein expression was analyzed in breast cancer cell lines. The expressions of HLA-A and MxA mRNAs were analyzed in MCF-7 cells in Gene Expression Omnibus (GEO) data. ER and HLA-ABC expressions and TIL levels in tumor tissue were also analyzed in ER+/ human epidermal growth factor receptor 2 (HER2)- breast cancer patients who randomly received either neoadjuvant chemotherapy or estrogen modulator treatment following surgical resection. Results: HLA-ABC protein expression was decreased after β-estradiol treatment or hESR-GFP transfection and increased after fulvestrant or IFN-γ treatment in breast cancer cell lines. In GEO data, HLA-A and MxA expression was increased after ESR1 shRNA transfection. When comparing the two patient groups, ER Allred score was significantly lower and the HLA-ABC expression and TIL levels were significantly higher in the estrogen modulator treated group than the chemotherapy treated group. Conclusion: MHC I expression and TIL levels might be affected by ER pathway modulation and IFN treatment. Further studies elucidating the mechanism of MHC I regulation could suggest a way to boost TIL influx in cancer in a clinical setting.

The Association of Estrogen Receptor Activity, Interferon Signaling, and MHC Class I Expression in Breast Cancer

Background: The expression of major histocompatibility complex class I (MHC I) has previously been reported to be negatively associated with estrogen receptor (ER) expression. Furthermore, the expression of MHC I, the level of tumor-infiltrating lymphocytes (TILs), and the expression of interferon (IFN) mediator MxA are positively associated with one another in human breast cancers. This study aimed to investigate the mechanisms of association of MHC I with estrogen and IFN signaling. Methods: The human leukocyte antigen (HLA)-ABC protein expression was analyzed in breast cancer cell lines. The expressions of HLA-A and MxA mRNAs were analyzed in MCF-7 cells in Gene Expression Omnibus (GEO) data. ER and HLA-ABC expressions and TIL levels in tumor tissue were also analyzed in ER+/ human epidermal growth factor receptor 2 (HER2)- breast cancer patients who randomly received either neoadjuvant chemotherapy or estrogen modulator treatment following surgical resection. Results: HLA-ABC protein expression was decreased after β-estradiol treatment or hESR-GFP transfection and increased after fulvestrant or IFN-γ treatment in breast cancer cell lines. In GEO data, HLA-A and MxA expression was increased after ESR1 shRNA transfection. When comparing the two patient groups, ER Allred score was significantly lower and the HLA-ABC expression and TIL levels were significantly higher in the estrogen modulator treated group than the chemotherapy treated group. Conclusion: MHC I expression and TIL levels might be affected by ER pathway modulation and IFN treatment. Further studies elucidating the mechanism of MHC I regulation could suggest a way to boost TIL influx in cancer in a clinical setting.

Structure of Ycf1p reveals the transmembrane domain TMD0 and the regulatory region of ABCC transporters

ATP binding cassette (ABC) proteins typically function in active transport of solutes across membranes. The ABC core structure is composed of two transmembrane domains (TMD1 and TMD2) and two cytosolic nucleotide binding domains (NBD1 and NBD2). Some members of the C-subfamily of ABC (ABCC) proteins, including human multidrug resistance proteins (MRPs), also possess an N-terminal transmembrane domain (TMD0) that contains five transmembrane α-helices and is connected to the ABC core by the L0 linker. While TMD0 was resolved in SUR1, the atypical ABCC protein that is part of the hetero-octameric ATP-sensitive K+ channel, little is known about the structure of TMD0 in monomeric ABC transporters. Here, we present the structure of yeast cadmium factor 1 protein (Ycf1p), a homolog of human MRP1, determined by electron cryo-microscopy (cryo-EM). A comparison of Ycf1p, SUR1, and a structure of MRP1 that showed TMD0 at low resolution demonstrates that TMD0 can adopt different orientations relative to the ABC core, including a ∼145° rotation between Ycf1p and SUR1. The cryo-EM map also reveals that segments of the regulatory (R) region, which links NBD1 to TMD2 and was poorly resolved in earlier ABCC structures, interacts with the L0 linker, NBD1, and TMD2. These interactions, combined with fluorescence quenching experiments of isolated NBD1 with and without the R region, suggest how posttranslational modifications of the R region modulate ABC protein activity. Mapping known mutations from MRP2 and MRP6 onto the Ycf1p structure explains how mutations involving TMD0 and the R region of these proteins lead to disease.

Structure of Ycf1p reveals the transmembrane domain TMD0 and the regulatory R region of ABCC transporters

ATP binding cassette (ABC) proteins typically function in active transport of solutes across membranes. The ABC core structure is comprised of two transmembrane domains (TMD1 and TMD2) and two cytosolic nucleotide binding domains (NBD1 and NBD2). Some members of the C-subfamily of ABC (ABCC) proteins, including human multidrug resistance proteins (MRPs), also possess an N-terminal transmembrane domain (TMD0) that contains five transmembrane α-helices and is connected to the ABC core by the L0 linker. While TMD0 was resolved in SUR1, the atypical ABCC protein that is part of the hetero-octameric ATP-sensitive K+ channel, little is known about the structure of TMD0 in monomeric ABC transporters. Here, we present the structure of yeast cadmium factor 1 protein (Ycf1p), a homologue of human MRP1, determined by electron cryomicroscopy (cryo-EM). Comparison of Ycf1p, SUR1, and a structure of MRP1 that showed TMD0 at low resolution demonstrates that TMD0 can adopt different orientations relative to the ABC core, including a 145° rotation between Ycf1p and SUR1. The cryo-EM map also reveals that segments of the regulatory (R) region, which links NBD1 to TMD2 and was poorly resolved in earlier ABCC structures, interacts with the L0 linker, NBD1, and TMD2. These interactions, combined with fluorescence quenching experiments of isolated NBD1 with and without the R region, suggests how post-translational modifications of the R region modulate ABC protein activity. Mapping known mutations from MRP2 and MRP6 onto the Ycf1p structure explains how mutations involving TMD0 and the R region of these proteins lead to disease.Statement of SignificanceThe Ycf1p structure provides an atomic model for the TMD0 domain of ABCC transporters and for two segments of the regulatory (R) region that links NBD1 to TMD2. The orientation of TMD0 in Ycf1p differs from that seen in SUR1, the regulatory ABCC protein in KATP channels, demonstrating flexibility in TMD0/ABC core contacts. The structure suggests how post-translational modifications of the R region modulate ABC protein activity and provides a mechanistic understanding of several diseases that occur due to mutation of human homologues of Ycf1p.

Acyl-CoA thioesterase activity of peroxisomal ABC protein ABCD1 is required for the transport of very long-chain acyl-CoA into peroxisomes

The ABCD1 protein, one of the four ATP-binding cassette (ABC) proteins in subfamily D, is located on the peroxisomal membrane and is involved in the transport of very long chain fatty acid (VLCFA)-CoA into peroxisomes. Its mutation causes X-linked adrenoleukodystophy (X-ALD): an inborn error of peroxisomal β-oxidation of VLCFA. Whether ABCD1 transports VLCFA-CoA as a CoA ester or free fatty acid is controversial. Recently, Comatose (CTS), a plant homologue of human ABCD1, has been shown to possess acyl-CoA thioesterase (ACOT) activity, and it is suggested that this activity is required for transport of acyl-CoA into peroxisomes. However, the precise transport mechanism is unknown. Here, we expressed human His-tagged ABCD1 in methylotrophic yeast, and characterized its ACOT activity and transport mechanism. The expressed ABCD1 possessed both ATPase and ACOT activities. The ACOT activity of ABCD1 was inhibited by p-chloromercuribenzoic acid (pCMB), a cysteine-reactive compound. Furthermore, we performed a transport assay with ABCD1-containing liposomes using 7-nitro-2–1,3-benzoxadiazol-4-yl (NBD)-labeled acyl-CoA as the substrate. The results showed that the fatty acid produced from VLCFA-CoA by ABCD1 is transported into liposomes and that ACOT activity is essential during this transport process. We propose a detailed mechanism of VLCFA-CoA transport by ABCD1.