Evolution of the 12 kDa FK506-binding protein gene

Background information. The FKBPs (FK506‐binding proteins) belong to a ubiquitous family of proteins that are found in a wide range of taxonomic groups. These proteins participate in a variety of pathways, including protein folding, down‐regulation of T‐cell activation and inhibition of cell‐cycle progression. Results. A cDNA encoding the 12 kDa FKBP gene orthologue (FKBP12) in Bombyx mori was been isolated from both Bm‐5 cultured cells and silk‐gland tissue. Using the FKBP12 cDNA in combination with the B. mori 6× whole‐genome shotgun database, we were able to identify the FKBP12 gene, as well as the positions of its intron—exon junctions. Conclusions. FKBP12 exon sizes and intronic positions are highly conserved among FKBP12 orthologues in 24 diverse genomes. Comparison of 41 FKBP12 genes revealed several intronic insertion and deletion events throughout evolution. In addition, paralogous FKBP12 isoforms were identified in all 12 vertebrate genomes. Both structural and phylogenetics analyses suggest that the isoforms may be evolving independently, possibly due to the distinct functional roles played by each paralogue.     

The FK506-binding protein, Fpr4, is an acidic histone chaperone

Fpr4, a FK506-binding protein (FKBP), is a recently identified novel histone chaperone. How it interacts with histones and facilitates their deposition onto DNA, however, are not understood. Here, we report a functional analysis that shows Fpr4 forms complexes with histones and facilitates nucleosome assembly like previously characterized acidic histone chaperones. We also show that the chaperone activity of Fpr4 resides solely in an acidic domain, while the peptidylprolyl isomerase domain conserved among all FKBPs inhibits the chaperone activity. These observations argue that Fpr4, while unique structurally, deposits histones onto DNA for nucleosome assembly through the well-established mechanism shared by other chaperones.     

Solution structure of FK506-binding protein 12 from Aedes aegypti

Dengue remains one of the major public concerns as the virus eludes the immune response. Currently, no vaccines or antiviral therapeutics are available for dengue prevention or treatment. Immunosuppressive drug FK506 shows an antimalarial activity, and its molecular target, FK506‐binding protein (FKBP), was identified in human Plasmodium parasites. Likewise, a conserved FKBP family protein has also been identified in Aedes aegypti (AaFKBP12), which is expected to play a similar role in the life cycle of Aedes aegypti, the primary vector of dengue virus infection. As FKBPs belong to a highly conserved class of immunophilin family and are involved in key biological regulations, they are considered as attractive pharmacological targets. In this study, we have determined the nuclear magnetic resonance solution structure of AaFKBP12, a novel FKBP member from Aedes aegypti, and presented its structural features, which may facilitate the design of potential inhibitory ligands against the dengue‐transmitting mosquitoes. Proteins 2012;. © 2012 Wiley Periodicals, Inc.     

Crystal structure of the FK506 binding domain of human FKBP25 in complex with FK506

Human FKBP25 (hFKBP25) is a nuclear immunophilin and interacts with several nuclear proteins, hence involving in many nuclear events. Similar to other FKBPs, FK506 binding domain (FKBD) of hFKBP25 also binds to immunosuppressive drugs such as rapamycin and FK506, albeit with a lower affinity for the latter. The molecular basis underlying this difference in affinity could not be addressed due to the lack of the crystal structure of hFKBD25 in complex with FK506. Here, we report the crystal structure of hFKBD25 in complex with FK506 determined at 1.8 Å resolution and its comparison with the hFKBD25–rapamycin complex, bringing out the microheterogeneity in the mode of interaction of these drugs, which could possibly explain the lower affinity for FK506.     

Structure-based prediction of DNA-binding sites on proteins using the empirical preference of electrostatic potential and the shape of molecular surfaces

Protein-DNA interactions play an essential role in the genetic activities of life. Many structures of protein-DNA complexes are already known, but the common rules on how and where proteins bind to DNA have not emerged. Many attempts have been made to predict protein-DNA interactions using structural information, but the success rate is still about 80%. We analyzed 63 protein-DNA complexes by focusing our attention on the shape of the molecular surface of the protein and DNA, along with the electrostatic potential on the surface, and constructed a new statistical evaluation function to make predictions of DNA interaction sites on protein molecular surfaces. The shape of the molecular surface was described by a combination of local and global average curvature, which are intended to describe the small convex and concave and the large-scale concave curvatures of the protein surface preferentially appearing at DNA-binding sites. Using these structural features, along with the electrostatic potential obtained by solving the Poisson-Boltzmann equation numerically, we have developed prediction schemes with 86% and 96% accuracy for DNA-binding and non-DNA-binding proteins, respectively.     

FKBP133: a novel mouse FK506-binding protein homolog alters growth cone morphology

FK506-binding proteins are the peptidyl prolyl cis-trans isomerases that are involved in various intracellular events. We characterized a novel mouse FK506-binding protein homolog, FKBP133/KIAA0674, in the developing nervous system. FKBP133 contains a domain similar to Wiskott-Aldrich syndrome protein homology region 1 (WH1) and a domain homologous to FK506-binding protein motif. FKBP133 was predominantly expressed in cerebral cortex, hippocampus, and peripheral ganglia at embryonic day 18.5. FKBP133 protein was distributed in the axonal shafts and was partially co-localized with F-actin in the growth cones of dorsal root ganglion neurons (DRG). The number of filopodia was increased in the DRG neurons overexpressing FKBP133. In contrast, the overexpression of a mutant deleted the WH1 domain reduced the growth cone size and the number of filopodia. Furthermore, the neurons overexpressing FKBP133 became significantly resistant to Semaphorin-3A induced collapse response. These results suggest that FKBP133 modulates growth cone behavior with the WH1 domain.     

N-terminal extension changes the folding mechanism of the FK506-binding protein

Many of the protein fusion systems used to enhance the yield of recombinant proteins result in the addition of a small number of amino acid residues onto the desired protein. Here, we investigate the effect of short (three amino acid) N-terminal extensions on the equilibrium denaturation and kinetic folding and unfolding reactions of the FK506-binding protein (FKBP) and compare the results obtained with data collected on an FKBP variant lacking this extension. Isothermal equilibrium denaturation experiments demonstrated that the N-terminal extension had a slight destabilizing effect. NMR investigations showed that the N-terminal extension slightly perturbed the protein structure near the site of the extension, with lesser effects being propagated into the single alpha-helix of FKBP. These structural perturbations probably account for the differential stability. In contrast to the relatively minor equilibrium effects, the N-terminal extension generated a kinetic-folding intermediate that is not observed in the shorter construct. Kinetic experiments performed on a construct with a different amino acid sequence in the extension showed that the length and the sequence of the extension both contribute to the observed equilibrium and kinetic effects. These results point to an important role for the N terminus in the folding of FKBP and suggest that a biological consequence of N-terminal methionine removal observed in many eukaryotic and prokaryotic proteins is to increase the folding efficiency of the polypeptide chain.     

Effects of pH, salt, and macromolecular crowding on the stability of FK506-binding protein: an integrated experimental and theoretical study

Environmental variables can exert significant influences on the folding stability of a protein, and elucidating these influences provides insight on the determinants of protein stability. Here, experimental data on the stability of FKBP12 are reported for the effects of three environmental variables: pH, salt, and macromolecular crowding. In the pH range of 5-9, contribution to the pH dependence of the unfolding free energy from residual charge-charge interactions in the unfolded state was found to be negligible. The negligible contribution was attributed to the lack of sequentially nearest neighboring charged residues around groups that titrate in the pH range. KCl lowered the stability of FKBP12 and the E31Q/D32N double mutant at small salt concentrations but raised stability after approximately 0.5 M salt. Such a turnover behavior was accounted for by the balance of two opposing types of protein-salt interactions: the Debye-Hückel type, modeling the response of the ions to protein charges, favors the unfolded state while the Kirkwood type, accounting for the disadvantage of the ions moving toward the low-dielectric protein cavity from the bulk solvent, disfavors the unfolded state. Ficoll 70 as a crowding agent was found to have a modest effect on protein stability, in qualitative agreement with a simple model suggesting that the folded and unfolded states are nearly equally adversely affected by macromolecular crowding. For any environmental variable, it is the balance of its effects on the folded and unfolded states that determines the outcome on the folding stability.     

Escherichia coli and other species of the enterobacteriaceae encode a protein similar to the family of Mip-like FK506-binding proteins

A newly identified gene in Escherichia coli, fkpA, encodes a protein with extensive similarity to the macrophage infectivity potentiator (Mip) proteins of Legionella pneumophila and Chlamydia trachomatis. The FkpA protein may be a new member of the family of FK506-binding proteins (FKBPs) because its carboxyl domain includes a sequence that matches the consensus FK506-binding motif in 40 of 48 positions. including those amino acids at the active site that form hydrogen bonds with the drug FK506. The amino acid sequence of the 29kDa FkpA protein is 30–35% identical to the Mip proteins of L. pneumophila, L. micdadei, and C. trachomatis. Of the 270 amino acids of FkpA, 113 (42%) are identical to the sequence of one or another of these Mip proteins. Overexpression of FkpA or deletion of fkpA from the E. coli chromosome had no detrimental effect on bacterial growth, indicating that fkpA is not an essential gene. Hybridization of fkpA-specific DNA probes to genomic blots révealed that similar genes exist in several representatives of the Enterobacteriaceae. Thus, mip-like genes are not found exelusively in bacteria having a predominately intracellular life style, but instead appear to be a new FKBP subfamily that is a common constituent of many bacteria.     

Structure-based optimization of designed Armadillo-repeat proteins

The armadillo domain is a right‐handed super‐helix of repeating units composed of three α‐helices each. Armadillo repeat proteins (ArmRPs) are frequently involved in protein–protein interactions, and because of their modular recognition of extended peptide regions they can serve as templates for the design of artificial peptide binding scaffolds. On the basis of sequential and structural analyses, different consensus‐designed ArmRPs were synthesized and show high thermodynamic stabilities, compared to naturally occurring ArmRPs. We determined the crystal structures of four full‐consensus ArmRPs with three or four identical internal repeats and two different designs for the N‐ and C‐caps. The crystal structures were refined at resolutions ranging from 1.80 to 2.50 Å for the above mentioned designs. A redesign of our initial caps was required to obtain well diffracting crystals. However, the structures with the redesigned caps caused domain swapping events between the N‐caps. To prevent this domain swap, 9 and 6 point mutations were introduced in the N‐ and C‐caps, respectively. Structural and biophysical analysis showed that this subsequent redesign of the N‐cap prevented domain swapping and improved the thermodynamic stability of the proteins. We systematically investigated the best cap combinations. We conclude that designed ArmRPs with optimized caps are intrinsically stable and well‐expressed monomeric proteins and that the high‐resolution structures provide excellent structural templates for the continuation of the design of sequence‐specific modular peptide recognition units based on armadillo repeats.     

C‐H…O hydrogen bonds in FK506‐binding protein–ligand interactions

Hydrogen bonds are important interaction forces observed in protein structures. They can be classified as stronger or weaker depending on their energy, thereby reflecting on the type of donor. The contribution of weak hydrogen bonds is deemed as an important factor toward structure stability along with the stronger bonds. One such bond, the C‐H…O type hydrogen bond, is shown to make a contribution in maintaining three dimensional structures of proteins. Apart from their presence within protein structures, the role of these bonds in protein–ligand interactions is also noteworthy. In this study, we present a statistical analysis on the presence of C‐H…O hydrogen bonds observed between FKBPs and their cognate ligands. The FK506‐binding proteins (FKBPs) carry peptidyl cis–trans isomerase activity apart from the immunosuppressive property by binding to the immunosuppressive drugs FK506 or rapamycin. Because the active site of FKBPs is lined up by many hydrophobic residues, we speculated that the prevalence of C‐H…O hydrogen bonds will be considerable. In a total of 25 structures analyzed, a higher frequency of C‐H…O hydrogen bonds is observed in comparison with the stronger hydrogen bonds. These C‐H…O hydrogen bonds are dominated by a highly conserved donor, the C^α/β of Val55 and an acceptor, the backbone oxygen of Glu54. Both these residues are positioned in the β4‐α1 loop, whereas the other residues Tyr26, Phe36 and Phe99 with higher frequencies are lined up at the opposite face of the active site. These preferences could be implicated in FKBP pharmacophore models toward enhancing the ligand affinity. This study could be a prelude to studying other proteins with hydrophobic pockets to gain better insights into ligand recognition. Copyright © 2013 John Wiley & Sons, Ltd.     

FK-506结合蛋白对钙释放通道的调控

细胞内自由钙作为一种重要的细胞信使广泛地参与细胞生理功能调控.胞内钙库(内质网系和肌浆网系)对调节细胞内自由钙水平起着重要的作用.钙库膜上的钙释放通道(ryanodine受体和三磷酸肌醇受体)受许多因素调控,其中之一就是新近研究得相当多的FK506结合蛋白.免疫抑制剂FK506能特异地结合钙库上一种分子质量为12 ku左右的蛋白,这种FK506结合蛋白与钙释放通道形成一种紧密连接的复合体,在正常生理情况下对钙释放通道起着十分重要的调控作用.     

FK506-binding proteins 12 and 12.6 (FKBPs) as regulators of cardiac Ryanodine Receptors: Insights from new functional and structural knowledge

Ryanodine Receptors (RyRs) are intracellular Ca²⁺ channels that mediate Ca²⁺ flux from the sarco(endo)plasmic reticulum in many cell types. The interaction of RyRs with FK506-binding proteins (FKBPs) has been proposed as an important regulatory mechanism, where the loss of this interaction leads to channel dysfunction. In the heart, phosphorylation of RyR has been suggested to disrupt the RyR-FKBP interaction promoting altered Ca²⁺ signaling, heart failure and arrhythmias. However, the functional result of FKBP interaction with RyR and how this interaction is regulated remains highly controversial. Recently, high resolution structures of RyR have provided novel aspects to the ongoing debate. This review will discuss the most recent functional data in light of these new structures.     

NMR solution structure and dynamics of the peptidyl-prolyl cis-trans isomerase domain of the trigger factor from Mycoplasma genitalium compared to FK506-binding protein

We have solved the solution structure of the peptidyl-prolyl cis-trans isomerase (PPIase) domain of the trigger factor from Mycoplasma genitalium by homo- and heteronuclear NMR spectroscopy. Our results lead to a well-defined structure with a backbone rmsd of 0.23 A. As predicted, the PPIase domain of the trigger factor adopts the FK506 binding protein (FKBP) fold. Furthermore, our NMR relaxation data indicate that the dynamic behavior of the trigger factor PPIase domain and of FKBP are similar. Structural variations when compared to FKBP exist in the flap region and within the bulges of strand 5 of the beta sheet. Although the active-site crevice is similar to that of FKBP, subtle steric variations in this region can explain why FK506 does not bind to the trigger factor. Sequence variability (27% identity) between trigger factor and FKBP results in significant differences in surface charge distribution and the absence of the first strand of the central beta sheet. Our data indicate, however, that this strand may be partially structured as "nascent" beta strand. This makes the trigger factor PPIase domain the most minimal representative of the FKBP like protein family of PPIases.     

New classification of supersecondary structures of sandwich-like proteins uncovers strict patterns of strand assemblage

To describe the supersecondary structure (SSS) of beta sandwich-like proteins (SPs), we introduce a structural unit called the "strandon." A strandon is defined as a set of sequentially consecutive strands connected by hydrogen bonds in 3D structures. Representing beta-proteins as the assembly of strandons exposes the underlying similarities in their SSS and enables us to construct a novel classification scheme of SPs. Classification of all known SPs is based on shared supersecondary structural features and is presented in the SSS database (http://binfs.umdnj.edu/sssdb/). Analysis of the SSS reveals two common specific patterns. The first pattern defines the arrangement of strandons and was found in 95% of all examined SPs. The second pattern establishes the ordering of strands in the protein domain and was observed in 82% of the analyzed SPs. Knowledge of these two patterns that uncover the spatial arrangement of strands will likely prove useful in protein structure prediction.     

The full electron structure of the FKBP12/FK506 complex

We present a study of FKBP12/FK506 using an electron structure calculation. These calculations employ a novel technique called eCADD on the protein’s full electron structure along with its hydrophobic pocket and the frontier-orbital-perturbation theory. We first obtain the energy bands and orbital coefficients of protein FKBP12. On this basis, we found that the activity atoms and activity residues of FKBP12 were in good agreement with X-ray crystallography experiments. The results indicate that the interactions occur only between the LUMOs of FKBP12 and the HOMO of FK506, not between the HOMOs of FKBP12 and the LUMO of FK506. In other words, the activity sites of protein FKBP12 are located on its LUMOs, not HOMOs. The electron structures of FKBP12/FK506 give us a clearer understanding of their interaction mechanism and will help us design new ligands of FKBP12.     

Sequence-based study of two related proteins with different folding behaviors

Z(SPA-1) is an engineered protein that binds to its parent, the three-helix-bundle Z domain of staphylococcal protein A. Uncomplexed Z(SPA-1) shows a reduced helix content and a melting behavior that is less cooperative, compared with the wild-type Z domain. Here we show that the difference in folding behavior between these two sequences can be partly understood in terms of an off-lattice model with 5-6 atoms per amino acid and a minimalistic potential, in which folding is driven by backbone hydrogen bonding and effective hydrophobic attraction.     

Solution structure of an alternate conformation of helix27 from Escherichia coli16S rRNA

Helix (H)27 of 16S ribosomal (r)RNA from Escherichia coli was dubbed the "switch helix" when mutagenesis suggested that two alternative base pair registers may have distinct functional roles in the bacterial ribosome. Although more recent genetic analyses suggest that H27 conformational switching is not required for translation, previous solution studies demonstrated that the isolated E. coli H27 can dynamically convert between the 885 and 888 conformations. Here, we have solved the nuclear magnetic resonance solution structure of a locked 888 conformation. NOE and residual dipolar coupling restraints reveal an architecture that markedly differs from that of the 885 conformation found in crystal structures of the bacterial ribosome. In place of the loop E motif that characterizes the 885 conformer and that the 888 conformer cannot adopt, we find evidence for an asymmetrical A-rich internal loop stabilized by stacking interactions among the unpaired A's. Comparison of the isolated H27 888 solution structure with the 885 crystal structure within the context of the ribosome suggests a difference in overall length of H27 that presents one plausible reason for the absence of H27 conformational switching within the sterically confining ribosome.