Structural basis of conformational transitions in the active site and 80's loop in the FK506-binding protein FKBP12

Propagation of conformational instability in FK506-binding protein FKBP12

FKBP12 is the archetype of the FK506 binding domains that define the family of FKBP proteins which participate in the regulation of various distinct physiological signaling processes. As the drugs FK506 and rapamycin inhibit many of these FKBP proteins, there is need to develop therapeutics which exhibit selectivity within this family. The long β4-β5 loop of the FKBP domain is known to regulate transcriptional activity for the steroid hormone receptors and appears to participate in regulating calcium channel activity for the cardiac and skeletal muscle ryanodine receptors. The β4-β5 loop of FKBP12 has been shown to undergo extensive conformational dynamics, and here we report hydrogen exchange measurements for a series of mutational variants in that loop which indicate deviations from a two-state kinetics for those dynamics. In addition to a previously characterized local transition near the tip of this loop, evidence is presented for a second site of conformational dynamics in the stem of this loop. These mutation-dependent hydrogen exchange effects extend beyond the β4-β5 loop, primarily by disrupting the hydrogen bond between the Gly 58 amide and the Tyr 80 carbonyl oxygen which links the two halves of the structural rim that surrounds the active site cleft. Mutationally-induced opening of the cleft between Gly 58 and Tyr 80 not only modulates the global stability of the protein, it promotes a conformational transition in the distant β2-β3a hairpin that modulates the binding affinity for a FKBP51-selective inhibitor previously designed to exploit a localized conformational transition at the homologous site.

Transient conformations in the unliganded FK506 binding domain of FKBP51 correspond to two distinct inhibitor-bound states

Members of the FK506-binding protein (FKBP) family regulate a range of important physiological processes. Unfortunately, current therapeutics such as FK506 and rapamycin exhibit only modest selectivity among these functionally distinct proteins. Recent progress in developing selective inhibitors has been reported for FKBP51 and FKBP52, which act as mutual antagonists in the regulation of steroid hormone signaling. Two structurally similar inhibitors yield distinct protein conformations at the binding site. Localized conformational transition in the binding site of the unliganded FK1 domain of FKBP51 is suppressed by a K58T mutation that also suppresses the binding of these inhibitors. Here, it is shown that the changes in amide hydrogen exchange kinetics arising from this K58T substitution are largely localized to this structural region. Accurate determination of the hydroxide-catalyzed exchange rate constants in both the wildtype and K58T variant proteins impose strong constraints upon the pattern of amide exchange reactivities within either a single or a pair of transient conformations that could give rise to the differences between these two sets of measured rate constants. Poisson-Boltzmann continuum dielectric calculations provide moderately accurate predictions of the structure-dependent hydrogen exchange reactivity for solvent-exposed protein backbone amides. Applying such calculations to the local protein conformations observed in the two inhibitor-bound FKBP51 domains demonstrated that the experimentally determined exchange rate constants for the wildtype domain are robustly predicted by a population-weighted sum of the experimental hydrogen exchange reactivity of the K58T variant and the predicted exchange reactivities in model conformations derived from the two inhibitor-bound protein structures.

Impact of distant peptide substrate residues on enzymatic activity of SlyD

Peptidyl-prolyl isomerases (PPIases) catalyze intrinsically slow and often rate-limiting isomerization of prolyl-peptide bonds in unfolded or partially folded proteins, thereby speeding up the folding process and preventing misfolding. They often possess binding and chaperone domains in addition to the domain carrying the isomerization activity. Although generally, their substrates display no identity in their amino acid sequence upstream and downstream of the proline with 20 possibilities for each residue, PPIases are efficient enzymes. SlyD is a highly efficient PPIase consisting of an isomerase domain and an additional chaperone domain. The binding of peptide substrates to SlyD and its enzymatic activity depend to some extend on the proline-proximal residues, however, the impact of proline-distant residues has not been investigated so far. Here, we introduce a label-free NMR-based method to measure SlyD activity on different peptide substrates and analysed the data in the context of obtained binding affinities and several co-crystal structures. We show that especially charged and aromatic residues up to eight positions downstream and three positions upstream of the proline and outside the canonical region of similar conformations affect the activity and binding, although they rarely display distinct conformations in our crystal structures. We hypothesize that these positions primarily influence the association reaction. In the absence of the chaperone domain the isomerase activity strongly correlates with substrate affinity, whereas additional factors play a role in its presence. The mutual orientation of isomerase and chaperone domains depends on the presence of substrates in both binding sites, implying allosteric regulation of enzymatic activity.

Conformational Insights into the Control of CNF1 Toxin Activity by Peptidyl-Prolyl Isomerization: A Molecular Dynamics Perspective

The cytotoxic necrotizing factor 1 (CNF1) toxin from uropathogenic Escherichia coli constitutively activates Rho GTPases by catalyzing the deamidation of a critical glutamine residue located in the switch II (SWII). In crystallographic structures of the CNF1 catalytic domain (CNF1^CD), surface-exposed P768 and P968 peptidyl-prolyl imide bonds (X-Pro) adopt an unusual cis conformation. Here, we show that mutation of each proline residue into glycine abrogates CNF1^CD in vitro deamidase activity, while mutant forms of CNF1 remain functional on RhoA in cells. Using molecular dynamics simulations coupled to protein-peptide docking, we highlight the long-distance impact of peptidyl-prolyl cis-trans isomerization on the network of interactions between the loops bordering the entrance of the catalytic cleft. The energetically favorable isomerization of P768 compared with P968, induces an enlargement of loop L1 that fosters the invasion of CNF1^CD catalytic cleft by a peptide encompassing SWII of RhoA. The connection of the P968 cis isomer to the catalytic cysteine C866 via a ladder of stacking interactions is alleviated along the cis-trans isomerization. Finally, the cis-trans conversion of P768 favors a switch of the thiol side chain of C866 from a resting to an active orientation. The long-distance impact of peptidyl-prolyl cis-trans isomerizations is expected to have implications for target modification.

Genomic discovery of an evolutionarily programmed modality for small-molecule targeting of an intractable protein surface

The vast majority of intracellular protein targets are refractory toward small-molecule therapeutic engagement, and additional therapeutic modalities are needed to overcome this deficiency. Here, the identification and characterization of a natural product, WDB002, reveals a therapeutic modality that dramatically expands the currently accepted limits of druggability. WDB002, in complex with the FK506-binding protein (FKBP12), potently and selectively binds the human centrosomal protein 250 (CEP250), resulting in disruption of CEP250 function in cells. The recognition mode is unprecedented in that the targeted domain of CEP250 is a coiled coil and is topologically featureless, embodying both a structural motif and surface topology previously considered on the extreme limits of "undruggability" for an intracellular target. Structural studies reveal extensive protein-WDB002 and protein-protein contacts, with the latter being distinct from those seen in FKBP12 ternary complexes formed by FK506 and rapamycin. Outward-facing structural changes in a bound small molecule can thus reprogram FKBP12 to engage diverse, otherwise "undruggable" targets. The flat-targeting modality demonstrated here has the potential to expand the druggable target range of small-molecule therapeutics. As CEP250 was recently found to be an interaction partner with the Nsp13 protein of the SARS-CoV-2 virus that causes COVID-19 disease, it is possible that WDB002 or an analog may exert useful antiviral activity through its ability to form high-affinity ternary complexes containing CEP250 and FKBP12.

Crystal structure and transient dimerization for the FKBP12 protein from the pathogenic fungus Candida auris

International concern over the recent emergence of Candida auris infections reflects not only its comparative ease of transmission and substantial mortality but the increasing level of resistance observed to all three major classes of antifungal drugs. Diminution in virulence has been reported for a wide range of fungal pathogens when the FK506-binding protein FKBP12 binds to that immunosuppressant drug and the binary complex then inhibits the fungal calcineurin signaling pathway. Structure-based drug design efforts have described modifications of FK506 which modestly reduce virulence for a number of fungal pathogens while also lessening the side effect of suppressing the tissue immunity response in the patient. To aid in such studies, we report the crystal structure of Candida auris FKBP12. As physiological relevance has been proposed for transient homodimerization interactions of distantly related fungal FKBP12 proteins, we report the solution NMR characterization of the homodimerization interactions of the FKBP12 proteins from both Candida auris and Candida glabrata.

Dynamics expression of DmFKBP12/Calstabin during embryonic early development of Drosophila melanogaster

Background

Calcium signaling are conserved from invertebrates to vertebrates and plays critical roles in many molecular mechanisms of embryogenesis and postnatal development. As a critical component of the signaling pathway, the RyR medicated calcium-induced calcium release signaling system, has been well studied along with their regulator FK506-binding protein 12 (FKBP12/Calstabin). Lack of FKBP12 is known to result in lethal cardiac dysfunction in mouse. However, precisely how FKBP12 is regulated and effects calcium signaling in Drosophila melanogaster remains largely unknown.

Results

In this study, we identified both temporal and localization changes in expression of DmFKBP12, a translational and transcriptional regulator of Drosophila RyR (DmRyR) and FKBP12, through embryonic development. DmFKBP12 is first expressed at the syncytial blastoderm stage and undergoes increased expression during the cellular blastoderm and early gastrulation stages. At late gastrulation, DmFKBP12 expression begins to decline until it reaches homeostasis, which it then maintains throughout the rest of development. Throughout these described changes in expression, DmFKBP12 mRNA remain stable, which indicates that protein dynamics are attributed to regulation at the mRNA to protein translation level. In addition to temporal changes in expression, dynamic expression profiles during Drosophila development also revealed DmFKBP12 localization. Although DmFKBP12 is distributed evenly between the anterior to posterior poles of the blastoderm egg, the protein is expressed more strongly in the cortex of the early Drosophila gastrula with the highest concentration found in the basement membrane of the cellular blastoderm. Fertilized egg, through the profile as under-membrane cortex distribution concentering onto basement at cellular blastoderm, to the profile as three-gem layer localization in primitive neuronal and digestion architecture of early Drosophila gastrula. By late gastrulation, DmFKBP12 is no longer identified in the yolk or lumen of duct structures and has relocated to the future brain (suboesophageal and supraesophageal ganglions), ventral nervous system, and muscular system. Throughout these changes in distribution, in situ DmFKBP12 mRNA monitoring detected equal distribution of DmFKBP12 mRNA, once again indicating that regulation of DmFKBP12 occurs at the translational level in Drosophila development.

Conclusion

As a critical regulator of the DmRyR-FKBP complex, DmFKBP12 expression in Drosophila fluctuates temporally and geographically with the formation of organ systems. These finding indicate that DmFKBP12 and RyR associated calcium signaling plays an essential role in the successful development of Drosophila melanogaster. Further study on the differences between mammalian RyR-FKBP12 and Drosophila DmRyR-FKBP12 can be exploited to develop safe pesticides.

Electronic supplementary material

The online version of this article (10.1186/s13578-019-0270-6) contains supplementary material, which is available to authorized users.

Quantifying protein dynamics in the ps-ns time regime by NMR relaxation

Both ¹⁵N chemical shift anisotropy (CSA) and sufficiently rapid exchange linebroadening transitions exhibit relaxation contributions that are proportional to the square of the magnetic field. Deconvoluting these contributions is further complicated by residue-dependent variations in protein amide ¹⁵N CSA values which have proven difficult to accurately measure. Exploiting recently reported improvements for the implementation of T₁ and T_1ρ experiments, field strength-dependent studies have been carried out on the B3 domain of protein G (GB3) as well as on the immunophilin FKBP12 and a H87V variant of that protein in which the major conformational exchange linebroadening transition is suppressed. By applying a zero frequency spectral density rescaling analysis to the relaxation data collected at magnetic fields from 500 to 900 MHz ¹H, differential residue-specific ¹⁵N CSA values have been obtained for GB3 which correlate with those derived from solid state and liquid crystalline NMR measurements to a level similar to the correlation among those previously reported studies. Application of this analysis protocol to FKBP12 demonstrated an efficient quantitation of both weak exchange linebroadening contributions and differential residue-specific ¹⁵N CSA values. Experimental access to such differential residue-specific ¹⁵N CSA values should significantly facilitate more accurate comparisons with molecular dynamics simulations of protein motion that occurs within the timeframe of global molecular tumbling.

Conformational Dynamics in FKBP Domains: Relevance to Molecular Signaling and Drug Design

Among the 22 FKBP domains in the human genome, FKBP12.6 and the first FKBP domains (FK1) of FKBP51 and FKBP52 are evolutionarily and structurally most similar to the archetypical FKBP12. As such, the development of inhibitors with selectivity among these four FKBP domains poses a significant challenge for structure-based design. The pleiotropic effects of these FKBP domains in a range of signaling processes such as the regulation of ryanodine receptor calcium channels by FKBP12 and FKBP12.6 and steroid receptor regulation by the FK1 domains of FKBP51 and FKBP52 amply justify the efforts to develop selective therapies. In contrast to their close structural similarities, these four FKBP domains exhibit a substantial diversity in their conformational flexibility. A number of distinct conformational transitions have been characterized for FKBP12 spanning timeframes from 20 s to 10 ns and in each case these dynamics have been shown to markedly differ from the conformational behavior for one or more of the other three FKBP domains. Protein flexibilitybased inhibitor design could draw upon the transitions that are significantly populated in only one of the targeted proteins. Both the similarities and differences among these four proteins valuably inform the understanding of how dynamical effects propagate across the FKBP domains as well as potentially how such intramolecular transitions might couple to the larger scale transitions that are central to the signaling complexes in which these FKBP domains function.

Coupling of Conformational Transitions in the N-terminal Domain of the 51-kDa FK506-binding Protein (FKBP51) Near Its Site of Interaction with the Steroid Receptor Proteins

Interchanging Leu-119 for Pro-119 at the tip of the β4-β5 loop in the first FK506 binding domain (FK1) of the FKBP51 and FKBP52 proteins, respectively, has been reported to largely reverse the inhibitory (FKBP51) or stimulatory (FKBP52) effects of these co-chaperones on the transcriptional activity of glucocorticoid and androgen receptor-protein complexes. Previous NMR relaxation studies have identified exchange line broadening, indicative of submillisecond conformational motion, throughout the β4-β5 loop in the FK1 domain of FKBP51, which are suppressed by the FKBP52-like L119P substitution. This substitution also attenuates exchange line broadening in the underlying β2 and β3a strands that is centered near a bifurcated main chain hydrogen bond interaction between these two strands. The present study demonstrates that these exchange line broadening effects arise from two distinct coupled conformational transitions, and the transition within the β2 and β3a strands samples a transient conformation that resembles the crystal structures of the selectively inhibited FK1 domain of FKBP51 recently reported. Although the crystal structures for their series of inhibitors were interpreted as evidence for an induced fit mechanism of association, the presence of a similar conformation being significantly populated in the unliganded FKBP51 domain is more consistent with a conformational selection binding process. The contrastingly reduced conformational plasticity of the corresponding FK1 domain of FKBP52 is consistent with the current model in which FKBP51 binds to both the apo- and hormone-bound forms of the steroid receptor to modulate its affinity for ligand, whereas FKBP52 binds selectively to the latter state.

Cloud computing approaches for prediction of ligand binding poses and pathways

We describe an innovative protocol for ab initio prediction of ligand crystallographic binding poses and highly effective analysis of large datasets generated for protein-ligand dynamics. We include a procedure for setup and performance of distributed molecular dynamics simulations on cloud computing architectures, a model for efficient analysis of simulation data, and a metric for evaluation of model convergence. We give accurate binding pose predictions for five ligands ranging in affinity from 7 nM to > 200 μM for the immunophilin protein FKBP12, for expedited results in cases where experimental structures are difficult to produce. Our approach goes beyond single, low energy ligand poses to give quantitative kinetic information that can inform protein engineering and ligand design.

Differential conformational dynamics in the closely homologous FK506-binding domains of FKBP51 and FKBP52

As co-chaperones of Hsp90 (heat-shock protein 90), FKBP51 (FK506-binding protein of 51 kDa) and FKBP52 (FK506-binding protein of 52 kDa) act as antagonists in regulating the hormone affinity and nuclear transport of steroid receptor complexes. Exchange of Leu¹¹⁹ in FKBP51 for Pro¹¹⁹ in FKBP52 has been shown to largely reverse the steroid receptor activities of FKBP51 and FKBP52. To examine whether differences in conformational dynamics/plasticity might correlate with changes in the reported receptor activities, ¹⁵N-NMR relaxation measurements were carried out on the N-terminal FKBP domains of FKBP51 and FKBP52 as well as their residue-swapped variants. Both proteins exhibit a similar pattern of motion in the picosecond–nanosecond timeframe as well as a small degree of ¹⁵N line-broadening, indicative of motion in the microsecond–millisecond timeframe, in the β_3a strand of the central sheet. Only the FKBP51 domain exhibits much larger line-broadening in the adjacent β₃ bulge (40′s loop of FKBP12) and throughout the long β₄–β₅ loop (80′s loop of FKBP12). The L119P mutation at the tip of the β₄–β₅ loop completely suppressed the line-broadening in this loop while partially suppressing the line-broadening in the neighbouring β₂ and β_3a strands. The complementary P119L and P119L/P124S variants of FKBP52 yielded similar patterns of line-broadening for the β₄–β₅ loop as that for FKBP51, although only 20% and 60% as intense respectively. However, despite the close structural similarity in the packing interactions between the β₄–β₅ loop and the β_3a strand for FKBP51 and FKBP52, the line-broadening in the β_3a strand is unaffected by the P119L or P119L/P124S mutations in FKBP52.

Unlike FKBP52, the FK1 domain of FKBP51 exhibits microsecond–millisecond conformational dynamics in the β₃ bulge and the β₄–β₅ loop, known sites of protein signalling interactions. Swapping residue 119 yields altered conformational dynamics in a pattern reminiscent of reported modulations in steroid receptor activity.

Statistical allosteric coupling to the active site indole ring flip equilibria in the FK506-binding domain

In solution, the Trp 59 indole ring at the base of the active site cleft in the FKBP domain protein FKBP12 is rotated by ~90° at a population level of 20%, relative to its canonical crystallographic orientation. NMR measurements on the homologous FK1 domains of human FKBP51 and FKBP52 indicate no observable indole ring flip conformation, while the V101I variant of FKBP12 decreases the population having a perpendicular indole orientation by 10-fold. A set of three parallel 400 ns CHARMM27 molecular simulations for both wild type FKBP12 and the V101I variant examined how this ring flip might be energetically coupled to a transition of the Glu 60 sidechain which interacts with the backbone of the 50's loop located ~12 Å from the indole nitrogen. Analysis of the transition matrix for the local dynamics of the Glu 60 sidechain, the Trp 59 sidechain, and of the structurally interposed α-helix hydrogen bonding pattern yielded a statistical allosteric coupling of 10 kJ/mol with negligible concerted dynamical coupling for the transitions of the two sidechains.