A kinetically trapped intermediate of FK506 binding protein forms in vitro: chaperone machinery dominates protein folding in vivo

Single-Molecule-Based, Label-Free Monitoring of Molecular Glue Efficacies for Promoting Protein-Protein Interactions Using YaxAB Nanopores

Modulating protein-protein interactions (PPIs) is an attractive strategy in drug discovery. Molecular glues, bifunctional small-molecule drugs that promote PPIs, offer an approach to targeting traditionally undruggable targets. However, the efficient discovery of molecular glues has been hampered by the current limitations of conventional ensemble-averaging-based methods. In this study, we present a YaxAB nanopore for probing the efficacy of molecular glues in inducing PPIs. Using YaxAB nanopores, we demonstrate single-molecule-based, label-free monitoring of protein complex formation between mammalian target of rapamycin (mTOR) and FK506-binding proteins (FKBPs) triggered by the molecular glue, rapamycin. Owing to its wide entrance and adjustable pore size, in combination with potent electro-osmotic flow (EOF), a single funnel-shaped YaxAB nanopore enables the simultaneous detection and single-molecule-level quantification of multiprotein states, including single proteins, binary complexes, and ternary complexes induced by rapamycin. Notably, YaxAB nanopores could sensitively discriminate between the binary complexes or ternary complexes induced by rapamycin and its analogues, despite the subtle size differences of ∼122 or ∼116 Da, respectively. Taken together, our results provide proof-of-concept for single-molecule-based, label-free, and ultrasensitive screening and structure-activity relationship (SAR) analysis of molecular glues, which will contribute to low-cost, highly efficient discovery, and rational design of bifunctional modality of drugs, such as molecular glues.

19F chemical library and 19F-NMR for a weakly bound complex structure

Fragment-based drug discovery (FBDD), which involves small compounds <300 Da, has been recognized as one of the most powerful tools for drug discovery. In FBDD, the affinity of hit compounds tends to be low, and the analysis of protein-compound interactions becomes difficult. In an effort to overcome such difficulty, we developed a ¹⁹F-NMR screening method optimizing a ¹⁹F chemical library focusing on highly soluble monomeric molecules. Our method was successfully applied to four proteins, including protein kinases and a membrane protein. For FKBP12, hit compounds were carefully validated by protein thermal shift analysis, ¹H-¹⁵N HSQC NMR spectroscopy, and isothermal titration calorimetry to determine dissociation constants and model complex structures. It should be noted that the ¹H and ¹⁹F saturation transfer difference experiments were crucial to obtaining highly precise model structures. The combination of ¹⁹F-NMR analysis and the optimized ¹⁹F chemical library enables the modeling of the complex structure made up of a weak binder and its target protein.

A Curvilinear-Path Umbrella Sampling Approach to Characterizing the Interactions Between Rapamycin and Three FKBP12 Variants

Rapamycin is an immunosuppressant macrolide that exhibits anti-proliferative properties through inhibiting the mTOR kinase. In fact, the drug first associates with the FKBP12 enzyme before interacting with the FRB domain of its target. Despite the availability of structural and thermodynamic information on the interaction of FKBP12 with rapamycin, the energetic and mechanistic understanding of this process is still incomplete. We recently reported a multiple-walker umbrella sampling simulation approach to characterizing the protein–protein interaction energetics along curvilinear paths. In the present paper, we extend our investigations to a protein-small molecule duo, the FKBP12•rapamycin complex. We estimate the binding free energies of rapamycin with wild-type FKBP12 and two mutants in which a hydrogen bond has been removed, D37V and Y82F. Furthermore, the underlying mechanistic details are analyzed. The calculated standard free energies of binding agree well with the experimental data, and the roles of the hydrogen bonds are shown to be quite different for each of these two mutated residues. On one hand, removing the carboxylate group of D37 strongly destabilizes the association; on the other hand, the hydroxyl group of Y82 is nearly unnecessary for the stability of the complex because some nonconventional, cryptic, indirect interaction mechanisms seem to be at work.

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.

Molecular dynamics simulations of site point mutations in the TPR domain of cyclophilin 40 identify conformational states with distinct dynamic and enzymatic properties

Cyclophilin 40 (Cyp40) is a member of the immunophilin family that acts as a peptidyl-prolyl-isomerase enzyme and binds to the heat shock protein 90 (Hsp90). Its structure comprises an N-terminal cyclophilin domain and a C-terminal tetratricopeptide (TPR) domain. Cyp40 is overexpressed in prostate cancer and certain T-cell lymphomas. The groove for Hsp90 binding on the TPR domain includes residues Lys227 and Lys308, referred to as the carboxylate clamp, and is essential for Cyp40-Hsp90 binding. In this study, the effect of two mutations, K227A and K308A, and their combinative mutant was investigated by performing a total of 5.76 μs of all-atom molecular dynamics (MD) simulations in explicit solvent. All simulations, except the K308A mutant, were found to adopt two distinct (extended or compact) conformers defined by different cyclophilin-TPR interdomain distances. The K308A mutant was only observed in the extended form which is observed in the Cyp40 X-ray structure. The wild-type, K227A, and combined mutant also showed bimodal distributions. The experimental melting temperature, Tm, values of the mutants correlate with the degree of compactness with the K308A extended mutant having a marginally lower melting temperature. Another novel measure of compactness determined from the MD data, the "coordination shell volume," also shows a direct correlation with Tm. In addition, the MD simulations show an allosteric effect with the mutations in the remote TPR domain having a pronounced effect on the molecular motions of the enzymatic cyclophilin domain which helps rationalise the experimentally observed increase in enzyme activity measured for all three mutations.

A Streamlined, Automated Protocol for the Production of Milligram Quantities of Untagged Recombinant Rat Lactate Dehydrogenase A Using ÄKTAxpressTM

We developed an efficient, automated 2-step purification protocol for the production of milligram quantities of untagged recombinant rat lactate dehydrogenase A (rLDHA) from E. coli, using the ÄKTAxpress™ chromatography system. Cation exchange followed by size exclusion results in average final purity in excess of 93% and yields ~ 14 milligrams per 50 ml of original cell culture in EnPresso B media, in under 8 hrs, including all primary sample processing and column equilibration steps. The protein is highly active and coherent biophysically and a viable alternative to the more problematic human homolog for structural and ligand-binding studies; an apo structure of untagged rLDHA was solved to a resolution 2.29 Å (PDB ID 5ES3). Our automated methodology uses generic commercially available pre-packed columns and simple buffers, and represents a robust standard method for the production of milligram amounts of untagged rLDHA, facilitating a novel fragment screening approach for new inhibitors.

Cyclophilin40 isomerase activity is regulated by a temperature-dependent allosteric interaction with Hsp90

Cyclophilin 40 (Cyp40) comprises an N-terminal cyclophilin domain with peptidyl-prolyl isomerase (PPIase) activity and a C-terminal tetratricopeptide repeat (TPR) domain that binds to the C-terminal-EEVD sequence common to both heat shock protein 70 (Hsp70) and Hsp90. We show in the present study that binding of peptides containing the MEEVD motif reduces the PPIase activity by ∼30%. CD and fluorescence assays show that the TPR domain is less stable than the cyclophilin domain and is stabilized by peptide binding. Isothermal titration calorimetry (ITC) shows that the affinity for the-MEEVD peptide is temperature sensitive in the physiological temperature range. Results from these biophysical studies fit with the MD simulations of the apo and holo (peptide-bound) structures which show a significant reduction in root mean square (RMS) fluctuation in both TPR and cyclophilin domains when-MEEVD is bound. The MD simulations of the apo-protein also highlight strong anti-correlated motions between residues around the PPIase-active site and a band of residues running across four of the seven helices in the TPR domain. Peptide binding leads to a distortion in the shape of the active site and a significant reduction in these strongly anti-correlated motions, providing an explanation for the allosteric effect of ligand binding and loss of PPIase activity. Together the experimental and MD results suggest that on heat shock, dissociation of Cyp40 from complexes mediated by the TPR domain leads to an increased pool of free Cyp40 capable of acting as an isomerase/chaperone in conditions of cellular stress.

UFSRAT: Ultra-fast Shape Recognition with Atom Types--the discovery of novel bioactive small molecular scaffolds for FKBP12 and 11βHSD1

MOTIVATION

Using molecular similarity to discover bioactive small molecules with novel chemical scaffolds can be computationally demanding. We describe Ultra-fast Shape Recognition with Atom Types (UFSRAT), an efficient algorithm that considers both the 3D distribution (shape) and electrostatics of atoms to score and retrieve molecules capable of making similar interactions to those of the supplied query.

RESULTS

Computational optimization and pre-calculation of molecular descriptors enables a query molecule to be run against a database containing 3.8 million molecules and results returned in under 10 seconds on modest hardware. UFSRAT has been used in pipelines to identify bioactive molecules for two clinically relevant drug targets; FK506-Binding Protein 12 and 11β-hydroxysteroid dehydrogenase type 1. In the case of FK506-Binding Protein 12, UFSRAT was used as the first step in a structure-based virtual screening pipeline, yielding many actives, of which the most active shows a KD, app of 281 µM and contains a substructure present in the query compound. Success was also achieved running solely the UFSRAT technique to identify new actives for 11β-hydroxysteroid dehydrogenase type 1, for which the most active displays an IC50 of 67 nM in a cell based assay and contains a substructure radically different to the query. This demonstrates the valuable ability of the UFSRAT algorithm to perform scaffold hops.

AVAILABILITY AND IMPLEMENTATION

A web-based implementation of the algorithm is freely available at http://opus.bch.ed.ac.uk/ufsrat/.

Engineering the acceptor specificity of trehalose phosphorylase for the production of trehalose analogs

Trehalose (α-D-glucopyranosyl-(1,1)-α-D-glucopyranoside) is widely used in the food industry, thanks to its protective effect against freezing and dehydration. Analogs of trehalose have the additional benefit that they are not digested and thus do not contribute to our caloric intake. Such trehalose analogs can be produced with the enzyme trehalose phosphorylase, when it is applied in the reverse, synthetic mode. Despite the enzyme's broad acceptor specificity, its catalytic efficiency for alternative monosaccharides is much lower than for glucose. For galactose, this difference is shown here to be caused by a lower K(m) whereas the k(cat) for both substrates is equal. Consequently, increasing the affinity was attempted by enzyme engineering of the trehalose phosphorylase from Thermoanaerobacter brockii, using both semirational and random mutagenesis. While a semirational approach proved unsuccessful, high-throughput screening of an error-prone PCR library resulted in the discovery of three beneficial mutations that lowered K(m) two- to three-fold. In addition, it was found that mutation of these positions also leads to an improved catalytic efficiency for mannose and fructose, suggesting their involvement in acceptor promiscuity. Combining the beneficial mutations did not further improve the affinity, and even resulted in a decreased catalytic activity and thermostability. Therefore, enzyme variant R448S is proposed as new biocatalyst for the industrial production of lactotrehalose (α-D-glucopyranosyl-(1,1)-α-D-galactopyranoside).

Streamlined, automated protocols for the production of milligram quantities of untagged recombinant human cyclophilin-A (hCypA) and untagged human proliferating cell nuclear antigen (hPCNA) using AKTAxpress

We developed streamlined, automated purification protocols for the production of milligram quantities of untagged recombinant human cyclophilin-A (hCypA) and untagged human proliferating cell nuclear antigen (hPCNA) from Escherichia coli, using the AKTAxpress chromatography system. The automated 2-step (cation exchange and size exclusion) purification protocol for untagged hCypA results in final purity and yields of 93% and approximately 5 mg L(-1) of original cell culture, respectively, in under 12h, including all primary sample processing and column equilibration steps. The novel automated 4-step (anion exchange, desalt, heparin-affinity and size exclusion, in linear sequence) purification protocol for untagged hPCNA results in final purity and yields of 87% and approximately 4 mg L(-1) of original cell culture, respectively, in under 24h, including all primary sample processing and column equilibration steps. This saves in excess of four full working days when compared to the traditional protocol, producing protein with similar final yield, purity and activity. Furthermore, it limits a time-dependent protein aggregation, a problem with the traditional protocol that results in a loss of final yield. Both automated protocols were developed to use generic commercially available pre-packed columns and automatically prepared minimal buffers, designed to eliminate user and system variations, maximize run reproducibility, standardize yield and purity between batches, increase throughput and reduce user input to a minimum. Both protocols represent robust generic methods for the automated production of untagged hCypA and hPCNA.

Fluorescent probes to characterise FK506-binding proteins

Talented all-rounders: Fluorescence polarisation assays were developed for members of the FK506-binding protein family by using fluorescent rapamycin analogues (demonstrated in the figure). These tracers retain medium to high affinity to all tested proteins (FKBP12, -12.6, -13, -25, -51, -52). They can be used for active-site titrations, competition assays with unlabelled ligands and enable a robust, miniaturized assay adequate for high-throughput screening.FK506-binding proteins (FKBPs) convey the immunosuppressive action of FK506 and rapamycin and mediate the neuroprotective properties of these compounds, and participate in the regulation of calcium channels. In addition, the larger homologues FKBP51 and FKBP52 act as cochaperones for Hsp90 and regulate the transactivational activity of steroid hormone receptors. To further characterize these FKBPs, we have synthesized fluorescein-coupled rapamycin analogues. In fluorescence polarization assays one of these compounds retained high affinity to all tested proteins (K(d): 0.1-20 nM) and could be used for active-site titrations. To adapt the fluorescence polarization assay for high-throughput purposes, a simplified rapamycin derivative was synthesized and labelled with fluorescein. This probe showed moderate affinity for the FK1 domains of FKBP51 (177 nM) and FKBP52 (469 nM) and allowed a highly robust, optimized, miniaturized assay (Z'>0.7) sufficient for high-throughput screening of large compound libraries.

A survey of the year 2007 literature on applications of isothermal titration calorimetry

Elucidation of the energetic principles of binding affinity and specificity is a central task in many branches of current sciences: biology, medicine, pharmacology, chemistry, material sciences, etc. In biomedical research, integral approaches combining structural information with in-solution biophysical data have proved to be a powerful way toward understanding the physical basis of vital cellular phenomena. Isothermal titration calorimetry (ITC) is a valuable experimental tool facilitating quantification of the thermodynamic parameters that characterize recognition processes involving biomacromolecules. The method provides access to all relevant thermodynamic information by performing a few experiments. In particular, ITC experiments allow to by-pass tedious and (rarely precise) procedures aimed at determining the changes in enthalpy and entropy upon binding by van't Hoff analysis. Notwithstanding limitations, ITC has now the reputation of being the "gold standard" and ITC data are widely used to validate theoretical predictions of thermodynamic parameters, as well as to benchmark the results of novel binding assays. In this paper, we discuss several publications from 2007 reporting ITC results. The focus is on applications in biologically oriented fields. We do not intend a comprehensive coverage of all newly accumulated information. Rather, we emphasize work which has captured our attention with originality and far-reaching analysis, or else has provided ideas for expanding the potential of the method.

Unfolding pathways of goat alpha-lactalbumin as revealed in multiple alignment of molecular dynamics trajectories

Molecular dynamics simulations of protein unfolding were performed at an elevated temperature for the authentic and recombinant forms of goat alpha-lactalbumin. Despite very similar three-dimensional structures, the two forms have significantly different unfolding rates due to an extra N-terminal methionine in the recombinant protein. To identify subtle differences between the two forms in the highly stochastic kinetics of unfolding, we classified the unfolding trajectories using the multiple alignment method based on the analogy between the biological sequences and the molecular dynamics trajectories. A dendrogram derived from the multiple trajectory alignment revealed a clear difference in the unfolding pathways of the authentic and recombinant proteins, i.e. the former reached the transition state in an all-or-none manner while the latter unfolded less cooperatively. It was also found in the classification that the two forms of the protein shared a common transition state structure, which was in excellent agreement with the transition state structure observed experimentally in the Phi-value analysis.