Interplay of self-association and conformational flexibility in regulating protein function

Oligomerization of drug transporters: Forms, functions, and mechanisms

Drug transporters are essential players in the transmembrane transport of a wide variety of clinical drugs. The broad substrate spectra and versatile distribution pattern of these membrane proteins infer their pharmacological and clinical significance. With our accumulating knowledge on the three-dimensional structure of drug transporters, their oligomerization status has become a topic of intense study due to the possible functional roles carried out by such kind of post-translational modification (PTM). In-depth studies of oligomeric complexes formed among drug transporters as well as their interactions with other regulatory proteins can help us better understand the regulatory mechanisms of these membrane proteins, provide clues for the development of novel drugs, and improve the therapeutic efficacy. In this review, we describe different oligomerization forms as well as their structural basis of major drug transporters in the ATP-binding cassette and solute carrier superfamilies, summarize our current knowledge on the influence of oligomerization for protein expression level and transport function of these membrane proteins, and discuss the regulatory mechanisms of oligomerization. Finally, we highlight the challenges associated with the current oligomerization studies and propose some thoughts on the pharmaceutical application of this important drug transporter PTM.

Unveiling clinically significant PPARγ mutations for thiazolidinedione treatment responsiveness through atomistic simulations

In Type 2 diabetes, increased insulin sensitivity is induced by thiazolidinedione activation of the peroxisome proliferator-activated receptor gamma (PPARγ). Recent data indicate a relationship between SNPs in PPARγ and poor drug response. Therefore, understanding the pathogenic consequences of mutations in PPARγ-mediated protein-drug interactions will be prima-facie for establishing personalized medicine. The PPARG gene has 197 missense SNPs, 22 of which were determined to be both deleterious and destabilizing, employing in silico approaches. Molecular docking analysis suggested that the mutation influenced the binding energy of at least seven of the variants. The mutant R316H was identified as the most damaging and deleterious from the observed results. For a better understanding of the dynamic variation upon mutation at the atomic level, molecular dynamics simulations of the wild-type and R316H mutant PPARγ structure were performed. The analysis indicates that the mutation increased protein structural compactness while decreasing flexibility. The reduced dynamics in the mutant structure was further validated by principal component analysis. This mechanistic evaluation of the PPARγ protein variants provides insight into the relationship between genetic variation and interindividual variability of drug responsiveness and will facilitate the future studies for the development of tailored treatment regime for precision medicine.

Unraveling the folding and dimerization properties of the human FoxP subfamily of transcription factors

Human FoxP proteins share a highly conserved DNA-binding domain that dimerizes via three-dimensional domain swapping, although showing varying oligomerization propensities among its members. Here, we present an experimental and computational characterization of all human FoxP proteins to unravel how their amino acid substitutions impact their folding and dimerization mechanism. We solved the crystal structure of the forkhead domain of FoxP4 to then perform a comparison across all members, finding that their sequence changes impact not only the structural heterogeneity of their forkhead domains but also the protein-protein association energy barrier. Lastly, we demonstrate that the accumulation of a monomeric intermediate is an oligomerization-dependent feature rather than a common aspect of monomers and dimers in this protein subfamily.

Soft disorder modulates the assembly path of protein complexes

The relationship between interactions, flexibility and disorder in proteins has been explored from many angles over the years: folding upon binding, flexibility of the core relative to the periphery, entropy changes, etc. In this work, we provide statistical evidence for the involvement of highly mobile and disordered regions in complex assembly. We ordered the entire set of X-ray crystallographic structures in the Protein Data Bank into hierarchies of progressive interactions involving identical or very similar protein chains, yielding 40205 hierarchies of protein complexes with increasing numbers of partners. We then examine them as proxies for the assembly pathways. Using this database, we show that upon oligomerisation, the new interfaces tend to be observed at residues that were characterised as softly disordered (flexible, amorphous or missing residues) in the complexes preceding them in the hierarchy. We also rule out the possibility that this correlation is just a surface effect by restricting the analysis to residues on the surface of the complexes. Interestingly, we find that the location of soft disordered residues in the sequence changes as the number of partners increases. Our results show that there is a general mechanism for protein assembly that involves soft disorder and modulates the way protein complexes are assembled. This work highlights the difficulty of predicting the structure of large protein complexes from sequence and emphasises the importance of linking predictors of soft disorder to the next generation of predictors of complex structure. Finally, we investigate the relationship between the Alphafold2’s confidence metric pLDDT for structure prediction in unbound versus bound structures, and soft disorder. We show a strong correlation between Alphafold2 low confidence residues and the union of all regions of soft disorder observed in the hierarchy. This paves the way for using the pLDDT metric as a proxy for predicting interfaces and assembly paths.

Dimerization of PHGDH via the catalytic unit is essential for its enzymatic function

Human D-3-phosphoglycerate dehydrogenase (PHGDH), a key enzyme in de novo serine biosynthesis, is amplified in various cancers and serves as a potential target for anticancer drug development. To facilitate this process, more information is needed on the basic biochemistry of this enzyme. For example, PHGDH was found to form tetramers in solution and the structure of its catalytic unit (sPHGDH) was solved as a dimer. However, how the oligomeric states affect PHGDH enzyme activity remains elusive. We studied the dependence of PHGDH enzymatic activity on its oligomeric states. We found that sPHGDH forms a mixture of monomers and dimers in solution with a dimer dissociation constant of ∼0.58 μM, with the enzyme activity depending on the dimer content. We computationally identified hotspot residues at the sPHGDH dimer interface. Single-point mutants at these sites disrupt dimer formation and abolish enzyme activity. Molecular dynamics simulations showed that dimer formation facilitates substrate binding and maintains the correct conformation required for enzyme catalysis. We further showed that the full-length PHGDH exists as a dynamic mixture of monomers, dimers, and tetramers in solution with enzyme concentration-dependent activity. Mutations that can completely disrupt the sPHGDH dimer show different abilities to interrupt the full-length PHGDH tetramer. Among them, E108A and I121A can also disrupt the oligomeric structures of the full-length PHGDH and abolish its enzyme activity. Our study indicates that disrupting the oligomeric structure of PHGDH serves as a novel strategy for PHGDH drug design and the hotspot residues identified can guide the design process.

Homodimerization of a glycoside hydrolase family GH1 β-glucosidase suggests distinct activity of enzyme different states

In this work, we investigated how activity and oligomeric state are related in a purified GH1 β-glucosidase from Spodoptera frugiperda (Sfβgly). Gel filtration chromatography coupled to a multiple angle light scattering detector allowed separation of the homodimer and monomer states and determination of the dimer dissociation constant (K_D ), which was in the micromolar range. Enzyme kinetic parameters showed that the dimer is on average 2.5-fold more active. Later, we evaluated the kinetics of homodimerization, scanning the changes in the Sfβgly intrinsic fluorescence over time when the dimer dissociates into the monomer after a large dilution. We described how the rate constant of monomerization (k_off ) is affected by temperature, revealing the enthalpic and entropic contributions to the process. We also evaluated how the rate constant (k_obs ) by which equilibrium is reached after dimer dilution behaves when varying the initial Sfβgly concentration. These data indicated that Sfβgly dimerizes through the conformational selection mechanism, in which the monomer undergoes a conformational exchange and then binds to a similar monomer, forming a more active homodimer. Finally, we noted that conformational selection reports and experiments usually rely on a ligand whose concentration is in excess, but for homodimerization, this approach does not hold. Hence, since our approach overcomes this limitation, this study not only is a new contribution to the comprehension of GH1 β-glucosidases, but it can also help to elucidate protein interaction pathways.

Mechanism of conditional partner selectivity in MITF/TFE family transcription factors with a conserved coiled coil stammer motif

Interrupted dimeric coiled coil segments are found in a broad range of proteins and generally confer selective functional properties such as binding to specific ligands. However, there is only one documented case of a basic-helix–loop–helix leucine zipper transcription factor—microphthalmia-associated transcription factor (MITF)—in which an insertion of a three-residue stammer serves as a determinant of conditional partner selectivity. To unravel the molecular principles of this selectivity, we have analyzed the high-resolution structures of stammer-containing MITF and an engineered stammer-less MITF variant, which comprises an uninterrupted symmetric coiled coil. Despite this fundamental difference, both MITF structures reveal identical flanking in-phase coiled coil arrangements, gained by helical over-winding and local asymmetry in wild-type MITF across the stammer region. These conserved structural properties allow the maintenance of a proper functional readout in terms of nuclear localization and binding to specific DNA-response motifs regardless of the presence of the stammer. By contrast, MITF heterodimer formation with other bHLH-Zip transcription factors is only permissive when both factors contain either the same type of inserted stammer or no insert. Our data illustrate a unique principle of conditional partner selectivity within the wide arsenal of transcription factors with specific partner-dependent functional readouts.