Accounting for global protein deformability during protein-protein and protein-ligand docking

Computational docking methods are valuable tools aimed to simplify the costly process of drug development and improvement. Most current approaches assume a rigid receptor structure to allow virtual screening of large numbers of possible ligands and putative binding sites on a receptor molecule. However, inclusion of receptor flexibility can be of critical importance since binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a ligand. Recent approaches to efficiently account for receptor flexibility during docking simulations are reviewed. In particular, accounting efficiently for global conformational changes of the protein backbone during docking is a still challenging unsolved problem. An approximate method has recently been suggested that is based on relaxing the receptor conformation during docking in pre-calculated soft collective degrees of freedom (M. Zacharias, Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP, Proteins: Struct., Funct., Genet. 54 (2004) 759-767). Test applications on protein-protein docking and on docking the inhibitor staurosporine to the apo-form of cAMP-dependent protein kinase A catalytic domain indicate significant improvement of docking results compared to rigid docking at a very modest computational demand. Accounting for receptor conformational changes in pre-calculated global degrees of freedom might offer a promising route to improve systematic docking screening simulations.     

Docking Studies on the Complexed and Uncomplexed FKBP12 Structures with Bound and Unbound Ligands: An Implication of a Conformational Selection Mechanism for Binding

Docking of FK506, rapamycin, and L-685,818 into their receptor, FKBP12, suggests that unlike the respective structures determined by X-ray crystallography, the uncomplexed FKBP12 structures determined by NMR may not be directly usable to identify high affinity ligands by docking studies for computational drug screening. In view of the resolution of the experimentally determined structures of FKBP12 and relatively small difference of the receptor binding sites between the complexed and uncomplexed states, it is unclear if the conformational induction mechanism is relevant to the binding of FKBP12 with its ligands. Alternatively, we advocate a conformation selection mechanism fundamentally akin to a mechanism proposed by Burgen. This mechanism better explains the experimental and calculated results for the binding of FKBP12 with FK506. It emphasizes that both guest and host select their most compatible preformed conformers to effect binding, and that the observed free energy of binding is a sum of the free energy change in complexation of the two most compatible conformers and the free energy changes in conversion of the Boltzmann-weighted principal conformers to the most compatible conformers. Conceptually, this mechanism represents one physical or nonphysical path of a thermodynamic cycle that is closed by the other path represented by the conformational induction mechanism, which can also be physical or nonphysical; it provides a theoretical means to estimate the affinity of the guest to the host with the experimentally available 3D structures of the two partners.     

Selective ligand purification using high-performance affinity beads

Since the development of affinity chromatography, affinity purification technology has been applied to many aspects of biological research, becoming an indispensable tool. Efficient strategies for the identification of biologically active compounds based on biochemical specificity have not yet been established, despite widespread interest in identifying chemicals that directly alter biomolecular functions. Here, we report a novel method for purifying chemicals that specifically interact with a target biomolecule using reverse affinity beads, a receptor-immobilized high-performance solid-phase matrix. When FK506-binding protein 12 (FKBP12) immobilized beads were used in this process, FK506 was efficiently purified in one step either from a mixture of chemical compounds or from fermented broth extract. The reverse affinity beads facilitated identification of drug/receptor complex binding proteins by reconstitution of immobilized ligand/receptor complexes on the beads. When FKBP12/FK506 and FKBP12/rapamycin complexes were immobilized, calcineurin and FKBP/rapamycin-associated protein were purified from a crude cell extract, respectively. These data indicate that reverse affinity beads are powerful tools for identification of both specific ligands and proteins that interact with receptor/ligand complexes.     

Solution structure of FK506 bound to FKBP-12.

The complex of the immunosuppressant FK506 bound to FKBP-12 has been studied in solution using 1H and inverse-detected 13C NMR methods. The resonances of bound, 13C-labelled FK506 were assigned and a set of 66 intraligand NOE distance restraints were used to calculate the structure of the bound ligand by distance geometry and restrained molecular dynamics methods. The structure of bound FK506 in solution closely resembles that seen in the X-ray structure [17], except for the allyl region. The differences reflect the influence of intermolecular crystal contacts and have implications for interpretation of the interaction of the FK506/FKBP complex with its putative biological receptor.     

Differential control of glucocorticoid receptor hormone-binding function by tetratricopeptide repeat (TPR) proteins and the immunosuppressive ligand FK506

Many laboratories have documented the existence of tetratricopeptide repeat (TPR) proteins (also known as immunophilins) in hormone-free steroid receptor complexes. Yet, the distinct roles of these proteins in steroid receptor action are poorly understood. In this work, we have investigated the effects of four TPR proteins (FKBP52, FKBP51, Cyp40, and PP5) on hormone-binding function of glucocorticoid receptor (GR) endogenously expressed in mammalian L929 cells. As a first step, we treated L929 cells with select immunophilin ligands [FK506, rapamycin, cyclosporin A (CsA), and cyclosporin H (CsH)], which are commonly thought to increase the GR response to hormone by inhibiting membrane-based steroid exporters. As expected, all four immunophilin ligands increased both the intracellular concentration of dexamethasone and GR activity at the MMTV-CAT reporter. To determine whether these ligands could target GR function independent of steroid export mechanisms, we performed GR reporter gene assays under conditions of immunophilin ligand and dexamethasone treatment that yielded equal intracellular hormone concentrations. FK506 was found to stimulate GR transactivity beyond the effect of this ligand on hormone retention. In contrast, CsA only affected the GR through upregulation of hormone retention. By Scatchard analysis, FK506 was found to increase GR hormone-binding affinity while decreasing total binding sites for hormone. This result correlated with loss of GR-associated FKBP51 and replacement with PP5. Interestingly, no GR-associated Cyp40 was found in these cells, consistent with the ability of CsA ligand to only affect GR through the hormone export mechanism. To test the role of FKBP52 independent of FK506, FKBP52 was placed under the control of a tetracycline-inducible promoter. Upregulation of FKBP52 caused an increase in both GR hormone-binding affinity and transactivity, even in the absence of FK506. These results show that immunosuppressive ligands can alter GR hormone-binding function by changing the TPR protein composition of receptor complexes and that TPR proteins exert a hierarchical effect on this GR function in the following order: FKBP52 > PP5 > FKBP51.     

Molecular cloning of a 25-kDa high affinity rapamycin binding protein, FKBP25.

Two FK506 binding proteins of molecular mass 12 kDa (FKBP12) and 13 kDa (FKBP13) have been identified as common cellular receptors of the immunosuppressants FK506 and rapamycin. Here we report the molecular cloning and overexpression of a 25-kDa rapamycin and FK506 binding protein (termed FKBP25) with peptidylprolyl cis-trans-isomerase (PPIase) activity. The amino acid sequence, predicted from the FKBP25 cDNA, shares identity with FKBP12 (44%) and FKBP13 (47%) in the C-terminal 97 amino acids. Unlike either FKBP12 or FKBP13, the nucleotide sequence of FKBP25 contains a number of putative nuclear localization sequences. The PPIase activity of recombinant FKBP25 was comparable with that of FKBP12. The PPIase activity of FKBP25 was far more sensitive to inhibition by rapamycin (IC50 = 50 nM) than FK506 (IC50 = 400 nM). PPIase activity of 100 nM FKBP25 was almost completely inhibited by 150 nM rapamycin while only 90% inhibition was achieved by 4 microM FK506. These data demonstrate that FKBP25 has a higher affinity for rapamycin than for FK506 and suggest that this cellular receptor may be an important target molecule for immunosuppression by rapamycin.     

Template-based docking of a prolactin receptor proline-rich motif octapeptide to FKBP12: implications for cytokine receptor signaling.

A conserved proline-rich motif (PRM) in the cytoplasmic domain of cytokine receptors has been suggested to be a signaling switch regulated by the action of the FK506 binding protein (FKBP) family of peptidylprolyl isomerases (O'Neal KD, Yu-Lee LY, Shearer WT, 1995, Ann NY Acad Sci 766:282-284). We have docked the prolactin receptor PRM (Ile1-Phe2-Pro3-Pro4-Val5-Pro6-Gly7-Pro8) to the ligand binding site of FKBP12. The procedure involved conformational search restricted by NMR restraints (O'Neal KD et al., 1996, Biochem J 315:833-844), energy minimization of the octapeptide conformers so obtained, template-based docking of a selected conformer to FKBP12, and energy refinement of the resulting complex. The template used was the crystal structure of a cyclic FK506-peptide hybrid bound to FKBP12. Val5-Pro6 of the PRM was taken to be the biologically relevant Xaa-Pro bond. The docked conformer is stabilized by two intramolecular hydrogen bonds, N7H7-->O4 and N2H2-->O8, and two intermolecular ones, Ile56; N-H-->O = C:Pro6 and Tyr82:O-H-->O = C:Gly7. This conformer features a Type I beta-turn and has extensive hydrophobic contacts with the FKBP12 binding surface. The observed interactions support the hypothesis that FKBP12 catalyzes cis-trans isomerization in the PRM when it is part of the longer cytoplasmic domain of a cytokine receptor, and suggest a significant role for the PRM in signal transduction.     

A fast protein-ligand docking algorithm based on hydrogen bond matching and surface shape complementarity

With the rapid development of structural determination of target proteins for human diseases, high throughout virtual screening based drug discovery is gaining popularity gradually. In this paper, a fast docking algorithm (H-DOCK) based on hydrogen bond matching and surface shape complementarity was developed. In H-DOCK, firstly a divide-and-conquer strategy based enumeration approach is applied to rank the intermolecular modes between protein and ligand by maximizing their hydrogen bonds matching, then each docked conformation of the ligand is calculated according to the matched hydrogen bonding geometry, finally a simple but effective scoring function reflecting mainly the van der Waals interaction is used to evaluate the docked conformations of the ligand. H-DOCK is tested for rigid ligand docking and flexible one, the latter is implemented by repeating rigid docking for multiple conformations of a small molecule and ranking all together. For rigid ligands, H-DOCK was tested on a set of 271 complexes where there is at least one intermolecular hydrogen bond, and H-DOCK achieved success rate (RMSD<2.0?Å) of 91.1%. For flexible ligands, H-DOCK was tested on another set of 93 complexes, where each case was a conformation ensemble containing native ligand conformation as well as 100 decoy ones generated by AutoDock [1], and the success rate reached 81.7%. The high success rate of H-DOCK indicates that the hydrogen bonding and steric hindrance can grasp the key interaction between protein and ligand. H-DOCK is quite efficient compared with the conventional docking algorithms, and it takes only about 0.14 seconds for a rigid ligand docking and about 8.25 seconds for a flexible one on average. According to the preliminary docking results, it implies that H-DOCK can be potentially used for large scale virtual screening as a pre-filter for a more accurate but less efficient docking algorithm.     

FLIPDock: docking flexible ligands into flexible receptors

Conformational changes of biological macromolecules when binding with ligands have long been observed and remain a challenge for automated docking methods. Here we present a novel protein-ligand docking software called FLIPDock (Flexible LIgand-Protein Docking) allowing the automated docking of flexible ligand molecules into active sites of flexible receptor molecules. In FLIPDock, conformational spaces of molecules are encoded using a data structure that we have developed recently called the Flexibility Tree (FT). While the FT can represent fully flexible ligands, it was initially designed as a hierarchical and multiresolution data structure for the selective encoding of conformational subspaces of large biological macromolecules. These conformational subspaces can be built to span a range of conformations important for the biological activity of a protein. A variety of motions can be combined, ranging from domains moving as rigid bodies or backbone atoms undergoing normal mode-based deformations, to side chains assuming rotameric conformations. In addition, these conformational subspaces are parameterized by a small number of variables which can be searched during the docking process, thus effectively modeling the conformational changes in a flexible receptor. FLIPDock searches the variables using genetic algorithm-based search techniques and evaluates putative docking complexes with a scoring function based on the AutoDock3.05 force-field. In this paper, we describe the concepts behind FLIPDock and the overall architecture of the program. We demonstrate FLIPDock's ability to solve docking problems in which the assumption of a rigid receptor previously prevented the successful docking of known ligands. In particular, we repeat an earlier cross docking experiment and demonstrate an increased success rate of 93.5%, compared to original 72% success rate achieved by AutoDock over the 400 cross-docking calculations. We also demonstrate FLIPDock's ability to handle conformational changes involving backbone motion by docking balanol to an adenosine-binding pocket of protein kinase A.     

The 12 kD FK 506 binding protein FKBP12 is released in the male reproductive tract and stimulates sperm motility.

BACKGROUND: The 12 kD FK506 binding protein FKBP12 is a cytosolic receptor for the immunosuppressant drugs FK506 and rapamycin. In addition to its critical role in drug-induced T-cell immunosuppression, FKBP12 associates physiologically with ryanodine and inositol 1,4,5-trisphosphate (IP3) receptors, regulating their ability to flux calcium. We investigated a role for FKBP12 in male reproductive physiology on the basis of our identification of extremely high levels of [3H]FK506 binding in male reproductive tissues. MATERIALS AND METHODS: [3H]FK506 binding studies were performed to identify tissues enriched with FK506 binding sites. The abundant [3H]FK506 binding sites identified in the male reproductive tract were localized by [3H]FK506 autoradiography. FK506 affinity chromatography was employed to purify FKBP from epididymal fluid. Anti-FKBP12 Western analysis was used to confirm the identity of the purified FKBP. The binding of exogenous FKBP12 to sperm was evaluated by [32P]FKBP12 binding studies and [33P]FKBP12 autoradiography. The effect of recombinant FKBP12 on sperm motility was investigated using a Hamilton Thorne motility analyzer. RESULTS: Male reproductive tissues contained high levels of [3H]FK506 binding. The localization of [3H]FK506 binding sites to the tubular epithelium of the caput epididymis and the lumen of the cauda and vas deferens suggested that FKBP is released in the male reproductive tract. FKBP12 was purified from epididymal plasma by FK506 affinity chromatography. Radiolabeled FKBP12 specifically bound to immature but not mature sperm. In sperm motility studies, FKBP12-treated caput sperm exhibited double the curvilinear velocity of untreated controls. CONCLUSIONS: High levels of FKBP12 are released in the male reproductive tract and specifically associate with maturing sperm. Recombinant FKBP12 enhances the curvilinear velocity of immature sperm, suggesting a role for FKBP12 in motility initiation. The highest concentrations of soluble FKBP12 in the male reproductive tract occur in the lumen of the vas deferens, a site of sperm storage and the conduit for ejaculated sperm. Preservation of mammalian sperm for reproductive technologies may be optimized by supplementing incubation or storage media with FKBP12.     

Structural characterization of the RyR1-FKBP12 interaction

The 12 kDa FK506-binding protein (FKBP12) constitutively binds to the calcium release channel RyR1. Removal of FKBP12 using FK506 or rapamycin causes an increased open probability and an increase in the frequency of sub-conductance states in RyR1. Using cryo-electron microscopy and single-particle image processing, we have determined the 3D difference map of FKBP12 associated with RyR1 at 16 A resolution that can be fitted with the atomic model of FKBP12 in a unique orientation. This has allowed us to better define the surfaces of close apposition between FKBP12 and RyR1. Our results shed light on the role of several FKBP12 residues that had been found critical for the specificity of the RyR1-FKBP12 interaction. As predicted from previous immunoprecipitation studies, our results suggest that Gln3 participates directly in this interaction. The orientation of RyR1-bound FKBP12, with part of its FK506 binding site facing towards RyR1, allows us to propose how FK506 is involved in the dissociation of FKBP12 from RyR1.     

Crystal Structures of the Free and Ligand-Bound FK1–FK2 Domain Segment of FKBP52 Reveal a Flexible Inter-Domain Hinge

The human Hsp90 co-chaperone FKBP52 belongs to the family of FK506-binding proteins, which act as peptidyl-prolyl isomerases. FKBP52 specifically enhances the signaling of steroid hormone receptors, modulates ion channels and regulates neuronal outgrowth dynamics. In turn, small-molecule ligands of FKBP52 have been suggested as potential neurotrophic or anti-prostate cancer agents. The usefulness of available ligands is however limited by a lack of selectivity. The immunophilin FKBP52 is composed of three domains, an FK506-binding domain with peptidyl-prolyl isomerase activity, an FKBP-like domain of unknown function and a TPR-clamp domain, which recognizes the C-terminal peptide of Hsp90 with high affinity. The herein reported crystal structures of FKBP52 reveal that the short linker connecting the FK506-binding domain and the FKBP-like domain acts as a flexible hinge. This enhanced flexibility and its modulation by phosphorylation might explain some of the functional antagonism between the closely related homologs FKBP51 and FKBP52. We further present two co-crystal structures of FKBP52 in complex with the prototypic ligand FK506 and a synthetic analog thereof. These structures revealed the molecular interactions in great detail, which enabled in-depth comparison with the corresponding complexes of the other cytosolic FKBPs, FKBP51 and FKBP12. The observed subtle differences provide crucial insights for the rational design of ligands with improved selectivity for FKBP52.     

Molecular dynamics simulation,binding free energy calculation and unbinding pathway analysis on selectivity difference between FKBP51 and FKBP52: Insight into the molecular mechanism of isoform selectivity

As co‐chaperones of the 90‐kDa heat shock protein(HSP90), FK506 binding protein 51 (FKBP51) and FK506 binding protein 52 (FKBP52) modulate the maturation of steroid hormone receptor through their specific FK1 domains (FKBP12‐like domain 1). The inhibitors targeting FK1 domains are potential therapies for endocrine‐related physiological disorders. However, the structural conservation of the FK1 domains between FKBP51 and FKBP52 make it difficult to obtain satisfactory selectivity in FK506‐based drug design. Fortunately, a series of iFit ligands synthesized by Hausch et al exhibited excellent selectivity for FKBP51, providing new opportunity for design selective inhibitors. We performed molecular dynamics simulation, binding free energy calculation and unbinding pathway analysis to reveal selective mechanism for the inhibitor iFit4 binding with FKBP51 and FKBP52. The conformational stability evaluation of the “Phe67‐in” and “Phe67‐out” states implies that FKBP51 and FKBP52 have different preferences for “Phe67‐in” and “Phe67‐out” states, which we suggest as the determinant factor for the selectivity for FKBP51. The binding free energy calculations demonstrate that nonpolar interaction is favorable for the inhibitors binding, while the polar interaction and entropy contribution are adverse for the inhibitors binding. According to the results from binding free energy decomposition, the electrostatic difference of residue 85 causes the most significant thermodynamics effects on the binding of iFit4 to FKBP51 and FKBP52. Furthermore, the importance of substructure units on iFit4 were further evaluated by unbinding pathway analysis and residue‐residue contact analysis between iFit4 and the proteins. The results will provide new clues for the design of selective inhibitors for FKBP51.     

Mutation of FKBP associated protein 48 (FAP48) at proline 219 disrupts the interaction with FKBP12 and FKBP52

Immunophilins are known as intracellular receptors for the immunosuppressant drugs, cyclosporin A, FK506, and rapamycin. They can be divided into two groups, cyclophilins that bind cyclosporin A and FK506 binding proteins (FKBPs) that bind FK506 and rapamycin. Many efforts were made to elucidate the physiological role of the immunophilins. Many of them are involved in intracellular signalling as they bind to calcium channels or to steroid receptor complexes. A yeast two-hybrid screen was used to identify further target proteins that interact with known proteins. Recently, a 48-kDa FKBP associated protein (FAP48) was isolated that binds to FKBP12 and FKBP52. Binding of FAP48 to FKBPs is inhibited by the macrolide FK506 indicating that the binding sites on the immunophilins coincide with the binding site for FK506. A peptidyl-prolyl motif on FAP48 should be responsible for the binding of the protein to FKBPs. We sequentially point mutated proline sites on FAP48 and checked the mutant proteins for interaction with FKBP12 and FKBP52. Mutation of proline 219 to alanine leads to a loss of interaction indicating that a cysteinyl prolyl site might be responsible for the binding of FAP48 to FKBPs. Thus we identified proline 219 being essential for the interaction.     

Solution structure of a neurotrophic ligand bound to FKBP12 and its effects on protein dynamics.

The structure of a recently reported neurotrophic ligand, 3-(3-pyridyl)-1-propyl(2S)-1-(3,3-dimethyl-1, 2-dioxopentyl)-2-pyrrolidinecarboxylate, in complex with FKBP12 was determined using heteronuclear NMR spectroscopy. The inhibitor exhibits a binding mode analogous to that observed for the macrocycle FK506, used widely as an immunosuppressant, with the prolyl ring replacing the pipecolyl moiety and the amide bond in a trans conformation. However, fewer favourable protein-ligand interactions are detected in the structure of the complex, suggesting weaker binding compared with the immunosuppressant drug. Indeed, a micromolar dissociation constant was estimated from the NMR ligand titration profile, in contrast to the previously published nanomolar inhibition activity. Although the inhibitor possesses a remarkable structural simplicity with respect to FK506, 15N relaxation studies show that it induces similar effects on the protein dynamics, stabilizing the conformation of solvent-exposed residues which are important for mediating the interaction of immunophilin/ligand complexes with molecular targets and potentially for the transmission of the neurotrophic action of FKBP12 inhibitors.     

Mean field analysis of FKBP12 complexes with FK506 and rapamycin: Implications for a role of crystallographic water molecules in molecular recognition and specificity

Mean field analysis of FKBP12 complexes with FK506 and rapamycin has been performed by using structures obtained from molecular docking simulations on a simple, yet robust molecular recognition energy landscape. When crystallographic water molecules are included in the simulations as an extension of the FKBP12 protein surface, there is an appreciable stability gap between the energy of the native FKBP12–FK506 complex and energies of conformations with the “native-like” binding mode. By contrast, the energy spectrum of the FKBP12–rapamycin complex is dense regardless of the presence of the water molecules. The stability gap in the FKBP12–FK506 system is determined by two critical water molecules from the effector region that participate in a network of specific hydrogen bond interactions. This interaction pattern protects the integrity and precision of the composite ligand-protein effector surface in the binary FKBP12–FK506 complex and is preserved in the crystal structure of the FKBP12–FK506–calcineurin ternary complex. These features of the binding energy landscapes provide useful insights into specific and nonspecific aspects of FK506 and rapamycin recognition. Proteins 28:313–324, 1997. © 1997 Wiley-Liss, Inc.     

NMR and crystallographic structures of the FK506 binding domain of human malarial parasite Plasmodium vivax FKBP35

The emergence of drug‐resistant malaria parasites is the major threat to effective malaria control, prompting a search for novel compounds with mechanisms of action that are different from the traditionally used drugs. The immunosuppressive drug FK506 shows an antimalarial activity. The mechanism of the drug action involves the molecular interaction with the parasite target proteins PfFKBP35 and PvFKBP35, which are novel FK506 binding protein family (FKBP) members from Plasmodium falciparum and Plasmodium vivax, respectively. Currently, molecular mechanisms of the FKBP family proteins in the parasites still remain elusive. To understand their functions, here we have determined the structures of the FK506 binding domain of Plasmodium vivax (PvFKBD) in unliganded form by NMR spectroscopy and in complex with FK506 by X‐ray crystallography. We found out that PvFKBP35 exhibits a canonical FKBD fold and shares kinetic profiles similar to those of PfFKBP35, the homologous protein in P. falciparum, indicating that the parasite FKBP family members play similar biological roles in their life cycles. Despite the similarity, differences were observed in the ligand binding modes between PvFKBD and HsFKBP12, a human FKBP homolog, which could provide insightful information into designing selective antimalarial drug against the parasites.     

An FKBP12 binding assay based upon biotinylated FKBP12.

A binding assay was developed for measuring the affinity of FKBP12 ligands. A biotinylation signal sequence was fused to the 5' end of the human FKBP12 gene, and the fusion protein was expressed in Escherichia coli with biotin ligase. The fusion protein was immobilized in avidin-coated multiwell plates, and varying concentrations of test ligands were allowed to compete with [3H]FK506 for FKBP12 sites on the plate. The assay provided Kd values for FK520, 32-hydroxyethyl indolyl FK520, and 18-ene, 20-oxa FK520 that are in agreement with previously reported values. The assay provides a convenient and rapid method for the assessment of FKBP12 binding by small molecules.     

Mechanism of osteogenic induction by FK506 via BMP/Smad pathways

FK506 is an immunosuppressant that exerts effects by binding to FK506-binding protein 12 (FKBP12). Recently, FK506 has also been reported to promote osteogenic differentiation when administered locally or in vitro in combination with bone morphogenetic proteins (BMPs), although the underlying mechanism remains unclarified. The present study initially showed that FK506 alone at a higher concentration (1muM) induced osteogenic differentiation of mesenchymal cell lines, which was suppressed by adenoviral introduction of Smad6. FK506 rapidly activates the BMP-dependent Smads in the absence of BMPs, and the activation was blocked by Smad6. Overexpression of FKBP12, which was reported to block the ligand-independent activation of BMP type I receptor A (BMPRIA), suppressed Smad signaling induced by FK506, but not that induced by BMP2. BMPRIA and FKBP12 bound to each other, and this binding was suppressed by FK506. These data suggest that FK506 promotes osteogenic differentiation by activating BMP receptors through interacting with FKBP12.