Structural coupling between FKBP12 and buried water

Map conformational landscapes of intrinsically disordered proteins with polymer physics quantities

Disordered proteins are conformationally flexible proteins that are biologically important and have been implicated in devastating diseases such as Alzheimer's disease and cancer. Unlike stably folded structured proteins, disordered proteins sample a range of different conformations that needs to be accounted for. Here, we treat disordered proteins as polymer chains, and compute a dimensionless quantity called instantaneous shape ratio (Rs), as Rs = Ree2/Rg2, where Ree is end-to-end distance and Rg is radius of gyration. Extended protein conformations tend to have high Ree compared with Rg, and thus have high Rs values, whereas compact conformations have smaller Rs values. We use a scatter plot of Rs (representing shape) against Rg (representing size) as a simple map of conformational landscapes. We first examine the conformational landscape of simple polymer models such as Random Walk, Self-Avoiding Walk, and Gaussian Walk (GW), and we notice that all protein/polymer maps lie within the boundaries of the GW map. We thus use the GW map as a reference and, to assess conformational diversity, we compute the fraction of the GW conformations (fC) covered by each protein/polymer. Disordered proteins all have high fC scores, consistent with their disordered nature. Each disordered protein accesses a different region of the reference map, revealing differences in their conformational ensembles. We additionally examine the conformational maps of the nonviral gene delivery vector polyethyleneimine at various protonation states, and find that they resemble disordered proteins, with coverage of the reference map decreasing with increasing protonation state, indicating decreasing conformational diversity. We propose that our method of combining Rs and Rg in a scatter plot generates a simple, meaningful map of the conformational landscape of a disordered protein, which in turn can be used to assess conformational diversity of disordered proteins.

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.

Chemical remodeling of a cellular chaperone to target the active state of mutant KRAS

The discovery of small-molecule inhibitors requires suitable binding pockets on protein surfaces. Proteins that lack this feature are considered undruggable and require innovative strategies for therapeutic targeting. KRAS is the most frequently activated oncogene in cancer, and the active state of mutant KRAS is such a recalcitrant target. We designed a natural product-inspired small molecule that remodels the surface of cyclophilin A (CYPA) to create a neomorphic interface with high affinity and selectivity for the active state of KRASG12C (in which glycine-12 is mutated to cysteine). The resulting CYPA:drug:KRASG12C tricomplex inactivated oncogenic signaling and led to tumor regressions in multiple human cancer models. This inhibitory strategy can be used to target additional KRAS mutants and other undruggable cancer drivers. Tricomplex inhibitors that selectively target active KRASG12C or multiple RAS mutants are in clinical trials now (NCT05462717 and NCT05379985).

Single-molecule fingerprinting of protein-drug interaction using a funneled biological nanopore

In drug discovery, efficient screening of protein-drug interactions (PDIs) is hampered by the limitations of current biophysical approaches. Here, we develop a biological nanopore sensor for single-molecule detection of proteins and PDIs using the pore-forming toxin YaxAB. Using this YaxAB nanopore, we demonstrate label-free, single-molecule detection of interactions between the anticancer Bcl-xL protein and small-molecule drugs as well as the Bak-BH3 peptide. The long funnel-shaped structure and nanofluidic characteristics of the YaxAB nanopore enable the electro-osmotic trapping of diverse folded proteins and high-resolution monitoring of PDIs. Distinctive nanopore event distributions observed in the two-dimensional (ΔI/I_o-versus-I_N) plot illustrate the ability of the YaxAB nanopore to discriminate individual small-molecule drugs bound to Bcl-xL from non-binders. Taken together, our results present the YaxAB nanopore as a robust platform for label-free, ultrasensitive, single-molecule detection of PDIs, opening up a possibility for low-cost, highly efficient drug discovery against diverse drug targets.

RosettaDDGPrediction for high-throughput mutational scans: From stability to binding

Reliable prediction of free energy changes upon amino acid substitutions (ΔΔGs) is crucial to investigate their impact on protein stability and protein-protein interaction. Advances in experimental mutational scans allow high-throughput studies thanks to multiplex techniques. On the other hand, genomics initiatives provide a large amount of data on disease-related variants that can benefit from analyses with structure-based methods. Therefore, the computational field should keep the same pace and provide new tools for fast and accurate high-throughput ΔΔG calculations. In this context, the Rosetta modeling suite implements effective approaches to predict folding/unfolding ΔΔGs in a protein monomer upon amino acid substitutions and calculate the changes in binding free energy in protein complexes. However, their application can be challenging to users without extensive experience with Rosetta. Furthermore, Rosetta protocols for ΔΔG prediction are designed considering one variant at a time, making the setup of high-throughput screenings cumbersome. For these reasons, we devised RosettaDDGPrediction, a customizable Python wrapper designed to run free energy calculations on a set of amino acid substitutions using Rosetta protocols with little intervention from the user. Moreover, RosettaDDGPrediction assists with checking completed runs and aggregates raw data for multiple variants, as well as generates publication-ready graphics. We showed the potential of the tool in four case studies, including variants of uncertain significance in childhood cancer, proteins with known experimental unfolding ΔΔGs values, interactions between target proteins and disordered motifs, and phosphomimetics. RosettaDDGPrediction is available, free of charge and under GNU General Public License v3.0, at https://github.com/ELELAB/RosettaDDGPrediction.

Structural features of chloroplast trigger factor determined at 2.6 Å resolution

The folding of newly synthesized polypeptides requires the coordinated action of molecular chaperones. Prokaryotic cells and the chloroplasts of plant cells possess the ribosome-associated chaperone trigger factor, which binds nascent polypeptides at their exit stage from the ribosomal tunnel. The structure of bacterial trigger factor has been well characterized and it has a dragon-shaped conformation, with flexible domains responsible for ribosome binding, peptidyl-prolyl cis-trans isomerization (PPIase) activity and substrate protein binding. Chloroplast trigger-factor sequences have diversified from those of their bacterial orthologs and their molecular mechanism in plant organelles has been little investigated to date. Here, the crystal structure of the plastidic trigger factor from the green alga Chlamydomonas reinhardtii is presented at 2.6 Å resolution. Due to the high intramolecular flexibility of the protein, diffraction to this resolution was only achieved using a protein that lacked the N-terminal ribosome-binding domain. The eukaryotic trigger factor from C. reinhardtii exhibits a comparable dragon-shaped conformation to its bacterial counterpart. However, the C-terminal chaperone domain displays distinct charge distributions, with altered positioning of the helical arms and a specifically altered charge distribution along the surface responsible for substrate binding. While the PPIase domain shows a highly conserved structure compared with other PPIases, its rather weak activity and an unusual orientation towards the C-terminal domain points to specific adaptations of eukaryotic trigger factor for function in chloroplasts.

Leveraging Fungal and Human Calcineurin-Inhibitor Structures, Biophysical Data, and Dynamics To Design Selective and Nonimmunosuppressive FK506 Analogs

Calcineurin is a critical enzyme in fungal pathogenesis and antifungal drug tolerance and, therefore, an attractive antifungal target. Current clinically accessible calcineurin inhibitors, such as FK506, are immunosuppressive to humans, so exploiting calcineurin inhibition as an antifungal strategy necessitates fungal specificity in order to avoid inhibiting the human pathway. Harnessing fungal calcineurin-inhibitor crystal structures, we recently developed a less immunosuppressive FK506 analog, APX879, with broad-spectrum antifungal activity and demonstrable efficacy in a murine model of invasive fungal infection. Our overarching goal is to better understand, at a molecular level, the interaction determinants of the human and fungal FK506-binding proteins (FKBP12) required for calcineurin inhibition in order to guide the design of fungus-selective, nonimmunosuppressive FK506 analogs. To this end, we characterized high-resolution structures of the Mucor circinelloides FKBP12 bound to FK506 and of the Aspergillus fumigatus, M. circinelloides, and human FKBP12 proteins bound to the FK506 analog APX879, which exhibits enhanced selectivity for fungal pathogens. Combining structural, genetic, and biophysical methodologies with molecular dynamics simulations, we identify critical variations in these structurally similar FKBP12-ligand complexes. The work presented here, aimed at the rational design of more effective calcineurin inhibitors, indeed suggests that modifications to the APX879 scaffold centered around the C₁₅, C₁₆, C₁₈, C₃₆, and C₃₇ positions provide the potential to significantly enhance fungal selectivity. IMPORTANCE Invasive fungal infections are a leading cause of death in the immunocompromised patient population. The rise in drug resistance to current antifungals highlights the urgent need to develop more efficacious and highly selective agents. Numerous investigations of major fungal pathogens have confirmed the critical role of the calcineurin pathway for fungal virulence, making it an attractive target for antifungal development. Although FK506 inhibits calcineurin, it is immunosuppressive in humans and cannot be used as an antifungal. By combining structural, genetic, biophysical, and in silico methodologies, we pinpoint regions of the FK506 scaffold and a less immunosuppressive analog, APX879, centered around the C₁₅ to C₁₈ and C₃₆ to C₃₇ positions that could be altered with selective extensions and/or deletions to enhance fungal selectivity. This work represents a significant advancement toward realizing calcineurin as a viable target for antifungal drug discovery.

Consequences of Stability-Induced Epistasis for Substitution Rates

Do interactions between residues in a protein (i.e., epistasis) significantly alter evolutionary dynamics? If so, what consequences might they have on inference from traditional codon substitution models which assume site-independence for the sake of computational tractability? To investigate the effects of epistasis on substitution rates, we employed a mechanistic mutation-selection model in conjunction with a fitness framework derived from protein stability. We refer to this as the stability-informed site-dependent (S-SD) model and developed a new stability-informed site-independent (S-SI) model that captures the average effect of stability constraints on individual sites of a protein. Comparison of S-SI and S-SD offers a novel and direct method for investigating the consequences of stability-induced epistasis on protein evolution. We developed S-SI and S-SD models for three natural proteins and showed that they generate sequences consistent with real alignments. Our analyses revealed that epistasis tends to increase substitution rates compared with the rates under site-independent evolution. We then assessed the epistatic sensitivity of individual site and discovered a counterintuitive effect: Highly connected sites were less influenced by epistasis relative to exposed sites. Lastly, we show that, despite the unrealistic assumptions, traditional models perform comparably well in the presence and absence of epistasis and provide reasonable summaries of average selection intensities. We conclude that epistatic models are critical to understanding protein evolutionary dynamics, but epistasis might not be required for reasonable inference of selection pressure when averaging over time and sites.

Environment-Specific Force Field for Intrinsically Disordered and Ordered Proteins

The need for accurate and efficient force fields for modeling 3D structures of macrobiomolecules and in particular intrinsically disordered proteins (IDPs) has increased with recent findings to associate IDPs and human diseases. However, most conventional protein force fields and recent IDP-specific force fields are limited in reproducing accurate structural features of IDPs. Here, we present an environmental specific precise force field (ESFF1) based on CMAP corrections of 71 different sequence environments to improve the accuracy and efficiency of MD simulation for both IDPs and folded proteins. MD simulations of 84 different short peptides, IDPs, and structured proteins show that ESFF1 can accurately reproduce spectroscopic properties for different peptides and proteins whether they are disordered or ordered. The successful ab initio folding of five fast-folding proteins further supports the reliability of ESFF1. The extensive analysis documented here shows that ESFF1 is able to achieve a reasonable balance between ordered and disordered states in protein simulations.

Identification of Protein Functional Regions

Protein sequence stores the information relative to both functionality and stability, thus making it difficult to disentangle the two contributions. However, the identification of critical residues for function and stability has important implications for the mapping of the proteome interactions, as well as for many pharmaceutical applications, e. g. the identification of ligand binding regions for targeted pharmaceutical protein design. In this work, we propose a computational method to identify critical residues for protein functionality and stability and to further categorise them in strictly functional, structural and intermediate. We evaluate single site conservation and use Direct Coupling Analysis (DCA) to identify co-evolved residues both in natural and artificial evolution processes. We reproduce artificial evolution using protein design and base our approach on the hypothesis that artificial evolution in the absence of any functional constraint would exclusively lead to site conservation and co-evolution events of the structural type. Conversely, natural evolution intrinsically embeds both functional and structural information. By comparing the lists of conserved and co-evolved residues, outcomes of the analysis on natural and artificial evolution, we identify the functional residues without the need of any a priori knowledge of the biological role of the analysed protein.

Enhanced multiple stress tolerance in Arabidopsis by overexpression of the polar moss peptidyl prolyl isomerase FKBP12 gene

PaFKBP12 overexpression in Arabidopsis resulted in stress tolerance to heat, ABA, drought, and salt stress, in addition to growth promotion under normal conditions. Polytrichastrum alpinum (alpine haircap moss) is one of polar organisms that can withstand the severe conditions of the Antarctic. In this study, we report the isolation of a peptidyl prolyl isomerase FKBP12 gene (PaFKBP12) from P. alpinum collected in the Antarctic and its functional implications in development and stress responses in plants. In P. alpinum, PaFKBP12 expression was induced by heat and ABA. Overexpression of PaFKBP12 in Arabidopsis increased the plant size, which appeared to result from increased rates of cell cycle. Under heat stress conditions, PaFKBP12-overexpressing lines (PaFKBP12-OE) showed better growth and survival than the wild type. PaFKBP12-OE also showed higher root elongation rates, better shoot growth and enhanced survival at higher concentrations of ABA in comparison to the wild type. In addition, PaFKBP12-OE were more tolerant to drought and salt stress than the wild type. All these phenotypes were accompanied with higher induction of the stress responsive genes in PaFKBP12-OE than in the wild type. Taken together, our findings revealed important functions of PaFKBP12 in plant development and abiotic stress responses.

The Roles of Water in the Protein Matrix: A Largely Untapped Resource for Drug Discovery

The value of thoroughly understanding the thermodynamics specific to a drug discovery/design study is well known. Over the past decade, the crucial roles of water molecules in protein structure, function, and dynamics have also become increasingly appreciated. This Perspective explores water in the biological environment by adopting its point of view in such phenomena. The prevailing thermodynamic models of the past, where water was seen largely in terms of an entropic gain after its displacement by a ligand, are now known to be much too simplistic. We adopt a set of terminology that describes water molecules as being "hot" and "cold", which we have defined as being easy and difficult to displace, respectively. The basis of these designations, which involve both enthalpic and entropic water contributions, are explored in several classes of biomolecules and structural motifs. The hallmarks for characterizing water molecules are examined, and computational tools for evaluating water-centric thermodynamics are reviewed. This Perspective's summary features guidelines for exploiting water molecules in drug discovery.

Ubiquitylation Directly Induces Fold Destabilization of Proteins

Ubiquitin is a common post-translational modifier and its conjugation is a key signal for proteolysis by the proteasome. Because the molecular mass of ubiquitin is larger than that of other modifiers such as phosphate, acetyl, or methyl groups, ubiquitylation not only influences biochemical signaling, but also may exert physical effects on its substrate proteins by increasing molecular volume and altering shape anisotropy. Here we show that ubiquitylation destabilizes the fold of two proteins, FKBP12 and FABP4, and that elongation of the conjugated ubiquitin chains further enhances this destabilization effect. Moreover, NMR relaxation analysis shows that ubiquitylation induces characteristic structural fluctuations in the backbone of both proteins. These results suggest that the ubiquitylation-driven structural fluctuations lead to fold destabilization of its substrate proteins. Thus, physical destabilization by ubiquitylation may facilitate protein degradation by the proteasome.

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.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Prediction of Water Binding to Protein Hydration Sites with a Discrete, Semiexplicit Solvent Model

Buried water molecules are ubiquitous in protein structures and are found at the interface of most protein-ligand complexes. Determining their distribution and thermodynamic effect is a challenging yet important task, of great of practical value for the modeling of biomolecular structures and their interactions. In this study, we present a novel method aimed at the prediction of buried water molecules in protein structures and estimation of their binding free energies. It is based on a semiexplicit, discrete solvation model, which we previously introduced in the context of small molecule hydration. The method is applicable to all macromolecular structures described by a standard all-atom force field, and predicts complete solvent distribution within a single run with modest computational cost. We demonstrate that it indicates positions of buried hydration sites, including those filled by more than one water molecule, and accurately differentiates them from sterically accessible to water but void regions. The obtained estimates of water binding free energies are in fair agreement with reference results determined with the double decoupling method.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

Statistical survey of the buried waters in the Protein Data Bank

The structures of buried water molecules were studied in an ensemble of high-quality and non-redundant protein crystal structures. Buried water molecules were clustered and classified in lake-like clusters, which are completely isolated from the bulk solvent, and bay-like clusters, which are in contact with the bulk solvent through a surface water molecule. Buried water molecules are extremely common: lake-like clusters are found in 89 % of the protein crystal structures and bay-like clusters in 93 %. Clusters with only one water molecule are much more common than larger clusters. Both cluster types incline to be surrounded by loop residues, and to a minor extent by residues in extended secondary structure. Helical residues on the contrary do not tend to surround clusters of buried water molecules. One buried water molecule is found every 30-50 amino acid residues, depending on the secondary structures that are more abundant in the protein. Both main- and side-chain atoms are in contact with buried waters; they form four hydrogen bonds with the first water and 1-1.5 additional hydrogen bond for each additional water in the cluster. Consequently, buried water molecules appear to be firmly packed and rigid like the protein atoms. In this regard, it is remarkable to observe that prolines often surround water molecules buried in the protein interior. Interestingly, clusters of buried water molecules tend to be just beneath the protein surface. Moreover, water molecules tend to form a one-dimensional wire rather than more compact arrangements. This agrees with recent evidence of the mechanisms of solvent exchange between internal cavities and bulk solvent.

A chemical switch for controlling viral infectivity

Chemically triggered molecular switches for controlling the fate and function of biological systems are fundamental to the emergence of synthetic biology and the development of biomedical applications. We here present the first chemically triggered switch for controlling the infectivity of adeno-associated viral (AAV) vectors.

Kinetics of exchange between zero-, one-, and two-hydrogen-bonded states of methyl and ethyl acetate in methanol

It has recently been shown that the ester carbonyl stretching vibration can be used as a sensitive probe of local electrostatic field in molecular systems. To further characterize this vibrational probe and extend its potential applications, we studied the kinetics of chemical exchange between differently hydrogen-bonded (H-bonded) ester carbonyl groups of methyl acetate (MA) and ethyl acetate (EA) in methanol. We found that, while both MA and EA can form zero, one, or two H-bonds with the solvent, the population of the 2hb state in MA is significantly smaller than that in EA. Using a combination of linear and nonlinear infrared measurements and numerical simulations, we further determined the rate constants for the exchange between these differently H-bonded states. We found that for MA the chemical exchange reaction between the two dominant states (i.e., 0hb and 1hb states) has a relaxation rate constant of 0.14 ps(-1), whereas for EA the three-state chemical exchange reaction occurs in a predominantly sequential manner with the following relaxation rate constants: 0.11 ps(-1) for exchange between 0hb and 1hb states and 0.12 ps(-1) for exchange between 1hb and 2hb states.

Computational perspective and evaluation of plausible catalytic mechanisms of peptidyl-prolyl cis-trans isomerases

BACKGROUND

Peptidyl prolyl cis-trans isomerization of the protein backbone is involved in the regulation of many biological processes. Cis-trans isomerization is notoriously slow and is catalyzed by a family of cis-trans peptidyl prolyl isomerases (PPIases) that have been implicated in many diseases. A general consensus on how these enzymes speed up prolyl isomerization has not been reached after decades of both experimental and computational studies.

SCOPE OF REVIEW

Computational studies carried out to understand the catalytic mechanism of the prototypical FK506 binding protein 12, Cyclophilin A and peptidyl-prolyl cis-trans isomerase NIMA-interacting 1 (Pin1) are reviewed. A summary and an evaluation of the implications of the proposed mechanisms from computational studies are presented.

MAJOR CONCLUSIONS

The analysis of computational studies and evaluation of the proposed mechanisms provide a general consensus and a better understanding of PPIase catalysis. The speedup of the rate of peptidyl-prolyl isomerization by PPIases can be best described by a catalytic mechanism in which the substrate in transition state configuration is stabilized. The enzymes preferentially bind the transition state configuration of the substrate relative to the cis conformation, which in most cases is bound better than the trans conformation of the substrate. Stabilization of the transition state configuration of the substrate leads to a lower free energy barrier and a faster rate of isomerization when compared to the uncatalyzed isomerization reaction.

GENERAL SIGNIFICANCE

Fully understanding the catalytic mechanism of PPIases has broad implications for drug design, elucidation of the molecular basis of many diseases, protein engineering, and enzyme catalysis in general. This article is part of a Special Issue entitled Proline-directed Foldases: Cell Signaling Catalysts and Drug Targets.

Structure of the rabbit ryanodine receptor RyR1 at near-atomic resolution

The ryanodine receptors (RyRs) are high-conductance intracellular Ca²⁺ channels that play a pivotal role in the excitation-contraction coupling of skeletal and cardiac muscles. RyRs are the largest known ion channels, with a homotetrameric organization and approximately 5000 residues in each protomer. Here we report the structure of the rabbit RyR1 in complex with its modulator FKBP12 at an overall resolution of 3.8 Å, determined by single-particle electron cryo-microscopy. Three previously uncharacterized domains, named Central, Handle, and Helical domains, display the armadillo repeat fold. These domains, together with the amino-terminal domain, constitute a network of superhelical scaffold for binding and propagation of conformational changes. The channel domain exhibits the voltage-gated ion channel superfamily fold with distinct features. A negative charge-enriched hairpin loop connecting S5 and the pore helix is positioned above the entrance to the selectivity filter vestibule. The four elongated S6 segments form a right-handed helical bundle that closes the pore at the cytoplasmic border of the membrane. Allosteric regulation of the pore by the cytoplasmic domains is mediated through extensive interactions between the Central domains and the channel domain. These structural features explain high ion conductance by RyRs and the long-range allosteric regulation of channel activities.

Structure-affinity properties of a high-affinity ligand of FKBP12 studied by molecular simulations of a binding intermediate

With a view to explaining the structure-affinity properties of the ligands of the protein FKBP12, we characterized a binding intermediate state between this protein and a high-affinity ligand. Indeed, the nature and extent of the intermolecular contacts developed in such a species may play a role on its stability and, hence, on the overall association rate. To find the binding intermediate, a molecular simulation protocol was used to unbind the ligand by gradually decreasing the biasing forces introduced. The intermediate was subsequently refined with 17 independent stochastic boundary molecular dynamics simulations that provide a consistent picture of the intermediate state. In this state, the core region of the ligand remains stable, notably because of the two anchoring oxygen atoms that correspond to recurrent motifs found in all FKBP12 ligand core structures. Besides, the non-core regions participate in numerous transient intermolecular and intramolecular contacts. The dynamic aspect of most of the contacts seems important both for the ligand to retain at least a part of its configurational entropy and for avoiding a trapped state along the binding pathway. Since the transient and anchoring contacts contribute to increasing the stability of the intermediate, as a corollary, the dissociation rate constant of this intermediate should be decreased, resulting in an increase of the affinity constant . The present results support our previous conclusions and provide a coherent rationale for explaining the prevalence in high-affinity ligands of (i) the two oxygen atoms found in carbonyl or sulfonyl groups of dissimilar core structures and of (ii) symmetric or pseudo-symmetric mobile groups of atoms found as non-core moieties. Another interesting aspect of the intermediate is the distortion of the flexible 80 s loop of the protein, mainly in its tip region, that promotes the accessibility to the bound state.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

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.

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

The extensive set of NMR doublings exhibited by the immunophilin FKBP12 (FK506-binding protein 12) arose from a slow transition to the cis-peptide configuration at Gly⁸⁹ near the tip of the 80′s loop, the site for numerous protein-recognition interactions for both FKBP12 and other FKBP domain proteins. The 80′s loop also exhibited linebroadening, indicative of microsecond to millisecond conformational dynamics, but only in the trans-peptide state. The G89A variant shifted the trans–cis peptide equilibrium from 88:12 to 33:67, whereas a proline residue substitution induced fully the cis-peptide configuration. The 80′s loop conformation in the G89P crystal structure at 1.50 Å resolution differed from wild-type FKBP12 primarily at residues 88, 89 and 90, and it closely resembled that reported for FKBP52. Structure-based chemical-shift predictions indicated that the microsecond to millisecond dynamics in the 80′s loop probably arose from a concerted main chain (ψ₈₈ and ϕ₈₉) torsion angle transition. The indole side chain of Trp⁵⁹ at the base of the active-site cleft was reoriented ~90^o and the adjacent backbone was shifted in the G89P crystal structure. NOE analysis of wild-type FKBP12 demonstrated that this indole populates the perpendicular orientation at 20%. The ¹⁵N relaxation analysis was consistent with the indole reorientation occurring in the nanosecond timeframe. Recollection of the G89P crystal data at 1.20 Å resolution revealed a weaker wild-type-like orientation for the indole ring. Differences in the residues that underlie the Trp⁵⁹ indole ring and altered interactions linking the 50′s loop to the active site suggested that reorientation of this ring may be disfavoured in the other six members of the FKBP domain family that bear this active-site tryptophan residue.

Extensive resonance doubling arises from a cis–trans peptide transition at Gly⁸⁹, whereas linebroadening appears due to a concerted shift in the neighbouring torsion angles. The active site Trp⁵⁹ ring adopts a perpendicular orientation at a population of 20%.

Crystal structure and conformational flexibility of the unligated FK506-binding protein FKBP12.6

The primary known physiological function of FKBP12.6 involves its role in regulating the RyR2 isoform of ryanodine receptor Ca(2+) channels in cardiac muscle, pancreatic β islets and the central nervous system. With only a single previously reported X-ray structure of FKBP12.6, bound to the immunosuppressant rapamycin, structural inferences for this protein have been drawn from the more extensive studies of the homologous FKBP12. X-ray structures at 1.70 and 1.90 Å resolution from P2₁ and P3₁21 crystal forms are reported for an unligated cysteine-free variant of FKBP12.6 which exhibit a notable diversity of conformations. In one monomer from the P3₁21 crystal form, the aromatic ring of Phe59 at the base of the active site is rotated perpendicular to its typical orientation, generating a steric conflict for the immunosuppressant-binding mode. The peptide unit linking Gly89 and Val90 at the tip of the protein-recognition `80s loop' is flipped in the P2₁ crystal form. Unlike the >30 reported FKBP12 structures, the backbone conformation of this loop closely follows that of the first FKBP domain of FKBP51. The NMR resonances for 21 backbone amides of FKBP12.6 are doubled, corresponding to a slow conformational transition centered near the tip of the 80s loop, as recently reported for 31 amides of FKBP12. The comparative absence of doubling for residues along the opposite face of the active-site pocket in FKBP12.6 may in part reflect attenuated structural coupling owing to increased conformational plasticity around the Phe59 ring.

Analysing the visible conformational substates of the FK506-binding protein FKBP12

The ¹H-¹⁵N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80′s loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50′s loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~10⁵-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80′s loop. Line-broadening was also decreased for the residues in the α-helix and 50′s loop, whereas line-broadening in the 40′s loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40′s loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.

The Polarizable Atomic Multipole-based AMOEBA Force Field for Proteins

Development of the AMOEBA (Atomic Multipole Optimized Energetics for Biomolecular Simulation) force field for proteins is presented. The current version (AMOEBA-2013) utilizes permanent electrostatic multipole moments through the quadrupole at each atom, and explicitly treats polarization effects in various chemical and physical environments. The atomic multipole electrostatic parameters for each amino acid residue type are derived from high-level gas phase quantum mechanical calculations via a consistent and extensible protocol. Molecular polarizability is modeled via a Thole-style damped interactive induction model based upon distributed atomic polarizabilities. Inter- and intramolecular polarization is treated in a consistent fashion via the Thole model. The intramolecular polarization model ensures transferability of electrostatic parameters among different conformations, as demonstrated by the agreement between QM and AMOEBA electrostatic potentials, and dipole moments of dipeptides. The backbone and side chain torsional parameters were determined by comparing to gas-phase QM (RI-TRIM MP2/CBS) conformational energies of dipeptides and to statistical distributions from the Protein Data Bank. Molecular dynamics simulations are reported for short peptides in explicit water to examine their conformational properties in solution. Overall the calculated conformational free energies and J-coupling constants are consistent with PDB statistics and experimental NMR results, respectively. In addition, the experimental crystal structures of a number of proteins are well maintained during molecular dynamics (MD) simulation. While further calculations are necessary to fully validate the force field, initial results suggest the AMOEBA polarizable multipole force field is able to describe the structure and energetics of peptides and proteins, in both gas-phase and solution environments.

Analysing the visible conformational substates of the FK506-binding protein FKBP12

The 1H-15N 2D NMR correlation spectrum of the widely studied FK506-binding protein FKBP12 (FK506-binding protein of 12 kDa) contains previously unreported peak doublings for at least 31 residues that arise from a minor conformational state (12% of total) which exchanges with the major conformation with a time constant of 3.0 s at 43°C. The largest differences in chemical shift occur for the 80's loop that forms critical recognition interactions with many of the protein partners for the FKBP family. The residues exhibiting doubling extend into the adjacent strands of the β-sheet, across the active site to the α-helix and into the 50's loop. Each of the seven proline residues adopts a trans-peptide linkage in both the major and minor conformations, indicating that this slow transition is not the result of prolyl isomerization. Many of the residues exhibiting resonance doubling also participate in conformational line-broadening transition(s) that occur ~105-fold more rapidly, proposed previously to arise from a single global process. The 1.70 Å (1 Å=0.1 nm) resolution X-ray structure of the H87V variant is strikingly similar to that of FKBP12, yet this substitution quenches the slow conformational transition throughout the protein while quenching the line-broadening transition for residues near the 80's loop. Line-broadening was also decreased for the residues in the α-helix and 50's loop, whereas line-broadening in the 40's loop was unaffected. The K44V mutation selectively reduces the line-broadening in the 40's loop, verifying that at least three distinct conformational transitions underlie the line-broadening processes of FKBP12.

Structure-based engineering of streptavidin monomer with a reduced biotin dissociation rate

We recently reported the engineering of monomeric streptavidin, mSA, corresponding to one subunit of wild type (wt) streptavidin tetramer. The monomer was designed by homology modeling, in which the streptavidin and rhizavidin sequences were combined to engineer a high affinity binding pocket containing residues from a single subunit only. Although mSA is stable and binds biotin with nanomolar affinity, its fast off rate (koff ) creates practical challenges during applications. We obtained a 1.9 Å crystal structure of mSA bound to biotin to understand their interaction in detail, and used the structure to introduce targeted mutations to improve its binding kinetics. To this end, we compared mSA to shwanavidin, which contains a hydrophobic lid containing F43 in the binding pocket and binds biotin tightly. However, the T48F mutation in mSA, which introduces a comparable hydrophobic lid, only resulted in a modest 20-40% improvement in the measured koff . On the other hand, introducing the S25H mutation near the bicyclic ring of bound biotin increased the dissociation half life (t½ ) from 11 to 83 min at 20°C. Molecular dynamics (MD) simulations suggest that H25 stabilizes the binding loop L3,4 by interacting with A47, and protects key intermolecular hydrogen bonds by limiting solvent entry into the binding pocket. Concurrent T48F or T48W mutation clashes with H25 and partially abrogates the beneficial effects of H25. Taken together, this study suggests that stabilization of the binding loop and solvation of the binding pocket are important determinants of the dissociation kinetics in mSA.

Small-molecule binding sites on proteins established by paramagnetic NMR spectroscopy

Determining the three-dimensional structure of a small molecule-protein complex with weak affinity can be a significant challenge. We present a paramagnetic NMR method to determine intermolecular structure restraints based on pseudocontact shifts (PCSs). Since the ligand must be in fast exchange between free and bound states and the fraction bound can be as low as a few percent, the method is ideal for ligands with high micromolar to millimolar dissociation constants. Paramagnetic tags are attached, one at a time, in a well-defined way via two arms at several sites on the protein surface. The ligand PCSs were measured from simple 1D (1)H spectra and used as docking restraints. An independent confirmation of the complex structure was carried out using intermolecular NOEs. The results show that structures derived from these two approaches are similar. The best results are obtained if the magnetic susceptibility tensors of the tags are known, but it is demonstrated that with two-armed probes, the magnetic susceptibility tensor can be predicted with sufficient accuracy to provide a low-resolution model of the ligand orientation and the location of the binding site in the absence of isotope-labeled protein. This approach can facilitate fragment-based drug discovery in obtaining structural information on the initial fragment hits.

Assessing the chemical accuracy of protein structures via peptide acidity

Although the protein native state is a Boltzmann conformational ensemble, practical applications often require a representative model from the most populated region of that distribution. The acidity of the backbone amides, as reflected in hydrogen exchange rates, is exquisitely sensitive to the surrounding charge and dielectric volume distribution. For each of four proteins, three independently determined X-ray structures of differing crystallographic resolution were used to predict exchange for the static solvent-exposed amide hydrogens. The average correlation coefficients range from 0.74 for ubiquitin to 0.93 for Pyrococcus furiosus rubredoxin, reflecting the larger range of experimental exchange rates exhibited by the latter protein. The exchange prediction errors modestly correlate with the crystallographic resolution. MODELLER 9v6-derived homology models at ~60% sequence identity (36% identity for chymotrypsin inhibitor CI2) yielded correlation coefficients that are ~0.1 smaller than for the cognate X-ray structures. The most recently deposited NOE-based ubiquitin structure and the original NMR structure of CI2 fail to provide statistically significant predictions of hydrogen exchange. However, the more recent RECOORD refinement study of CI2 yielded predictions comparable to the X-ray and homology model-based analyses.

The FKBP families of higher plants: Exploring the structures and functions of protein interaction specialists

The FK506-binding proteins (FKBPs) are known both as the receptors for immunosuppressant drugs and as prolyl isomerase (PPIase) enzymes that catalyse rotation of prolyl bonds. FKBPs are characterised by the inclusion of at least one FK506-binding domain (FKBd), the receptor site for proline and the active site for PPIase catalysis. The FKBPs form large and diverse families in most organisms, with the largest FKBP families occurring in higher plants. Plant FKBPs are molecular chaperones that interact with specific protein partners to regulate a diversity of cellular processes. Recent studies have found that plant FKBPs operate in intricate and coordinated mechanisms for regulating stress response and development processes, and discoveries of new interaction partners expand their cellular influences to gene expression and photosynthetic adaptations. This review presents an examination of the molecular and structural features and functional roles of the higher plant FKBP family within the context of these recent findings, and discusses the significance of domain conservation and variation for the development of a diverse, versatile and complex chaperone family.

A role for Fkbp6 and the chaperone machinery in piRNA amplification and transposon silencing

Epigenetic silencing of transposons by Piwi-interacting RNAs (piRNAs) constitutes an RNA-based genome defense mechanism. Piwi endonuclease action amplifies the piRNA pool by generating new piRNAs from target transcripts by a poorly understood mechanism. Here, we identified mouse Fkbp6 as a factor in this biogenesis pathway delivering piRNAs to the Piwi protein Miwi2. Mice lacking Fkbp6 derepress LINE1 (L1) retrotransposon and display reduced DNA methylation due to deficient nuclear accumulation of Miwi2. Like other cochaperones, Fkbp6 associates with the molecular chaperone Hsp90 via its tetratricopeptide repeat (TPR) domain. Inhibition of the ATP-dependent Hsp90 activity in an insect cell culture model results in the accumulation of short antisense RNAs in Piwi complexes. We identify these to be byproducts of piRNA amplification that accumulate only in nuage-localized Piwi proteins. We propose that the chaperone machinery normally ejects these inhibitory RNAs, allowing turnover of Piwi complexes for their continued participation in piRNA amplification.

Molecular interactions of c-ABL mutants in complex with imatinib/nilotinib: a computational study using linear interaction energy (LIE) calculations

In spite of the effectiveness of Imatinib for chronic myeloid leukemia (CML) treatment, resistance has repeatedly been reported and is associated with point mutations in the BCR-ABL chimeric gene. To overcome this resistance, several inhibitors of BCR-ABL tyrosine kinase activity were developed. In this context, computational simulations have become a powerful tool for understanding drug-protein interactions. Herein, we report a comparative molecular dynamics analysis of the interaction between two tyrosine kinase inhibitors (imatinib or nilotinib) against wild type c-ABL protein and 12 mutants, using the semi-empirical linear interaction energy (LIE) method, to assess the feasibility of this approach for studying resistance against the inhibitory activity of these drugs. In addition, to understand the structural changes that are associated with resistance, we describe the behavior of water molecules that interact simultaneously with specific residues (Glu286, Lys271 and Asp381) of c-ABL (wild type or mutant) and their relationship with drug resistance. Experimental IC50 values for the interaction between imatinib, wild type c-ABL, and 12 mutants were used to obtain the proper LIE coefficients (α, β and γ) to estimate the free energy of the binding of imatinib with wild-type and mutant proteins, and values were extrapolated for the analysis of the nilotinib/c-ABL interaction. Our results indicate that LIE was suitable to predict the superior inhibitory activity of nilotinib and the resistance to inhibition that was observed in c-ABL mutants. Additionally, for c-ABL mutants, the observed number of water molecules being turned over while interacting with amino acids Glu286, Lys271 and Asp381 was associated with resistance to imatinib, resulting in a less effective inhibition of the kinase activity.

Exploration of pipecolate sulfonamides as binders of the FK506-binding proteins 51 and 52

FK506-binding proteins (FKBP) 51 and 52 are cochaperones that modulate the signal transduction of steroid hormone receptors. Single nucleotide polymorphisms in the gene encoding FKBP51 have been associated with a variety of psychiatric disorders. Rapamycin and FK506 are two macrocyclic natural products, which tightly bind to most FKBP family members, including FKBP51 and FKBP52. A bioisosteric replacement of the α-ketoamide moiety of rapamycin and FK506 with a sulfonamide was envisaged with the retention of the conserved hydrogen bonds. A focused solid support-based synthesis protocol was developed, which led to ligands with submicromolar affinity for FKBP51 and FKBP52. The molecular binding mode for one sulfonamide analogue was confirmed by X-ray crystallography.

Towards automated crystallographic structure refinement with phenix.refine

phenix.refine is a program within the PHENIX package that supports crystallographic structure refinement against experimental data with a wide range of upper resolution limits using a large repertoire of model parameterizations. It has several automation features and is also highly flexible. Several hundred parameters enable extensive customizations for complex use cases. Multiple user-defined refinement strategies can be applied to specific parts of the model in a single refinement run. An intuitive graphical user interface is available to guide novice users and to assist advanced users in managing refinement projects. X-ray or neutron diffraction data can be used separately or jointly in refinement. phenix.refine is tightly integrated into the PHENIX suite, where it serves as a critical component in automated model building, final structure refinement, structure validation and deposition to the wwPDB. This paper presents an overview of the major phenix.refine features, with extensive literature references for readers interested in more detailed discussions of the methods.

The structure of a Burkholderia pseudomallei immunophilin-inhibitor complex reveals new approaches to antimicrobial development

Mips (macrophage infectivity potentiators) are a subset of immunophilins associated with virulence in a range of micro-organisms. These proteins possess peptidylprolyl isomerase activity and are inhibited by drugs including rapamycin and tacrolimus. We determined the structure of the Mip homologue [BpML1 (Burkholderia pseudomallei Mip-like protein 1)] from the human pathogen and biowarfare threat B. pseudomallei by NMR and X-ray crystallography. The crystal structure suggests that key catalytic residues in the BpML1 active site have unexpected conformational flexibility consistent with a role in catalysis. The structure further revealed BpML1 binding to a helical peptide, in a manner resembling the physiological interaction of human TGFβRI (transforming growth factor β receptor I) with the human immunophilin FKBP12 (FK506-binding protein 12). Furthermore, the structure of BpML1 bound to the class inhibitor cycloheximide N-ethylethanoate showed that this inhibitor mimics such a helical peptide, in contrast with the extended prolyl-peptide mimicking shown by inhibitors such as tacrolimus. We suggest that Mips, and potentially other bacterial immunophilins, participate in protein-protein interactions in addition to their peptidylprolyl isomerase activity, and that some roles of Mip proteins in virulence are independent of their peptidylprolyl isomerase activity.

Targeting FKBP isoforms with small-molecule ligands

The FK506 binding protein (FKBP) family of proteins provide an interesting series of drug targets since different isoforms modulate diverse cellular pathways. There are therapeutic opportunities in the fields of cancer therapy, neurodegenerative conditions and psychiatric disorders. X-ray crystallographic or NMR data are available for eight of fourteen human FKBPs covering ten of the twenty-two different FKBP domains. We have made a detailed sequence and structural comparison of human FKBP domains. These data show that the chemical scaffolds common to the immunosuppressive inhibitors FK506 and rapamycin bind to the most conserved region of the binding site. This observation opens the way to the design of isoform specific inhibitors.