Modeling conformational changes in cyclosporin A.

NMR and X-ray structures for the immunosuppressant cyclosporin A (CsA) reveal a remarkable difference between the unbound (free) conformation in organic solvents and the conformation bound to cyclophilin. We have performed computer simulations of the molecular dynamics of CsA under a variety of conditions and confirmed the stability of these two conformations at room temperature in water and in vacuum. However, when the free conformation was modeled in vacuum at 600 K, a transition pathway leading to the bound conformation was observed. This involved a change in the cis MeLeu-9 peptide bond to a trans conformation and the movement of the side chains forming the dominant hydrophobic cluster (residues MeBmt-1, MeLeu-4, MeLeu-6, and MeLeu-10) to the opposite side of the plane formed by the backbone atoms in the molecular ring. The final conformation had a backbone RMS deviation from the bound conformation of 0.53 A and was as stable in dynamics simulations as the bound conformation. Our calculations allowed us to make a detailed analysis of a transition pathway between the free and the bound conformations of CsA and to identify two distinct regions of coordinated movement in CsA, both of which underwent transitions independently.     

Solution structure of cyclosporin A and a nonimmunosuppressive analog bound to fully deuterated cyclophilin.

A simple strategy involving 1H nuclear magnetic resonance (NMR) spectroscopy and complete protein deuteration was used to determine the structures of two receptor-bound drugs. A potent immunosuppressive, cyclosporin A (CsA) binds tightly to the ubiquitous and highly conserved 17.7-kDa immunophilin, cyclophilin (CyP). Fully deuterated CyP was produced by overexpressing the human CyP gene in Escherichia coli grown on deuterated algal hydrolysate in 98% D2O. As only the CsA molecule is protonated in the CsA-CyP complex, we were able to make a complete sequential assignment of the bound drug using standard two-dimensional proton NMR experiments. The structure determination was accomplished using dynamical simulated annealing calculations with a total of 124 NMR-derived distance and torsion angle restraints. Aside from binding CsA, CyP also acts as a peptidyl-prolyl cis-trans isomerase. Thus, much importance had been ascribed to the cis peptide bond present in the structures reported for free CsA in organic solvents and in crystal studies. Interestingly, CyP-bound CsA exists in an all-trans conformation with no detectable elements of regular secondary structure and no intramolecular hydrogen bonds. A nonactive CsA analog, MeAla6-CsA, was studied using the same CyP deuteration strategy. In addition to structural elucidation of the two bound drugs, we were able to differentiate between the bound and surface-exposed residues of the drugs and also validate our previous hypothesis that the single CyP tryptophan is located in the CsA-binding site.(ABSTRACT TRUNCATED AT 250 WORDS)     

Improved modeling of side-chains in proteins with rotamer-based methods: a flexible rotamer model

Side-chain modeling has a widespread application in many current methods for protein tertiary structure determination, prediction, and design. Of the existing side-chain modeling methods, rotamer-based methods are the fastest and most efficient. Classically, a rotamer is conceived as a single, rigid conformation of an amino acid sidechain. Here, we present a flexible rotamer model in which a rotamer is a continuous ensemble of conformations that cluster around the classic rigid rotamer. We have developed a thermodynamically based method for calculating effective energies for the flexible rotamer. These energies have a one-to-one correspondence with the potential energies of the rigid rotamer. Therefore, the flexible rotamer model is completely general and may be used with any rotamer-based method in substitution of the rigid rotamer model. We have compared the performance of the flexible and rigid rotamer models with one side-chain modeling method in particular (the self-consistent mean field theory method) on a set of 20 high quality crystallographic protein structures. For the flexible rotamer model, we obtained average predictions of 85.8% for chi1, 76.5% for chi1+2 and 1.34 A for root-mean-square deviation (RMSD); the corresponding values for core residues were 93.0%, 87.7% and 0.70 A, respectively. These values represent improvements of 7.3% for chi1, 8.1% for chi1+2 and 0.23 A for RMSD over the predictions obtained with the rigid rotamer model under otherwise identical conditions; the corresponding improvements for core residues were 6.9%, 10.5% and 0.43 A, respectively. We found that the predictions obtained with the flexible rotamer model were also significantly better than those obtained for the same set of proteins with another state-of-the-art side-chain placement method in the literature, especially for core residues. The flexible rotamer model represents a considerable improvement over the classic rigid rotamer model. It can, therefore, be used with considerable advantage in all rotamer-based methods commonly applied to protein tertiary structure determination, prediction, and design and also in predictions of free energies in mutational studies.     

NMR studies of [U-13C]cyclosporin A bound to cyclophilin: bound conformation and portions of cyclosporin involved in binding

Cyclosporin A (CsA), a potent immunosuppressant, is known to bind with high specificity to cyclophilin (CyP), a 17.7 kDa protein with peptidyl-prolyl isomerase activity. In order to investigate the three-dimensional structure of the CsA/CyP complex, we have applied a variety of multidimensional NMR methods in the study of uniformly 13C-labeled CsA bound to cyclophilin. The 1H and 13C NMR signals of cyclosporin A in the bound state have been assigned, and from a quantitative interpretation of the 3D NOE data, the bound conformation of CsA has been determined. Three-dimensional structures of CsA calculated from the NOE data by using a distance geometry/simulated appealing protocol were found to be very different from previously determined crystalline and solution conformations of uncomplexed CsA. In addition, from CsA/CyP NOEs, the portions of CsA that interact with cyclophilin were identified. For the most part, those CsA residues with NOEs to cyclophilin were the same residues important for cyclophilin binding and immunosuppressive activity as determined from structure/activity relationships. The structural information derived in this study together with the known structure/activity relationships for CsA analogues may prove useful in the design of improved immunosuppressants. Moreover, the approach that is described for obtaining the structural information is widely applicable to the study of small molecule/large molecule interactions.     

The cyclophilin multigene family of peptidyl-prolyl isomerases. Characterization of three separate human isoforms.

Cyclophilin (CyP), a major cytosolic protein possessing peptidyl-prolyl cis-trans isomerase activity, has been implicated as the specific receptor of the immunosuppressive drug cyclosporin A (CsA). To identify other potential CsA receptors related to CyP, two human cDNA libraries were screened under low stringency conditions using human CyP cDNA (encoding hCyP1) as a probe. Two cDNAs were identified which encode distinct proteins related to human hCyP1. These two novel proteins, designated hCyP2 and hCyP3, share 65 and 76% amino acid sequence homology with hCyP1, respectively. Both hCyP2 and hCyP3 contain NH2-terminal hydrophobic extensions of 32 and 42 amino acids, respectively. Protein-specific antibodies revealed the predominant association of hCyP2 and hCyP3 with membranes and subcellular organelles, which suggests that the amino-terminal leader sequences of the two CyP isoforms may act as signal peptides. In contrast to the results with hCyP1, Southern blot analysis indicated that both hCyP2 and hCyP3 gene sequences are represented infrequently in the human genome. Northern and Western blot analysis showed that the distribution of mRNA and proteins of the three hCyPs in differing tissues and cell types was similar. Each hCyP protein was expressed in Escherichia coli, purified, and shown to be an active peptidyl-prolyl isomerase. Substrate specificity was examined with 11 synthetic peptides (Suc-Xaa-Yaa-Pro-Phe-4-nitroanilide), and inhibition of the peptidyl-prolyl isomerase activities associated with hCyP1, hCyP2, and hCyP3 was studied with CsA, MeAla6-CsA and MeBm2t1-CsA. From both equilibrium considerations and the results of kinetic characterizations it is proposed that of these three CyP proteins, hCyP1 is the most likely intracellular target for CsA.     

Effects of cyclosporine A and its immunosuppressive or non-immunosuppressive derivatives [D-Ser]8-CsA and Cs9 on mitochondria from different brain regions

We studied the functional properties of isolated brain mitochondria (BM) prepared from total rat brain (BM(total)) or from cerebral subregions under basal and Ca(2+) overload conditions in order to evaluate the effects of cyclosporine A (CsA) in a regiospecific manner. CsA-induced effects were compared with those of two derivatives-the none-immunosuppressive [O-(NH(2)(CH2)(5)NHC(O)CH(2))-D-Ser](8)-CsA (Cs9) and its congener, the immunosuppressive [D-Ser](8)-CsA. The glutamate/malate-dependent state 3 respiration of mitochondria (state 3(glu/mal)) differed in region-specific manner (cortex > striatum = cerebellum > substantia nigra > hippocampus), but was significantly increased by 1μM CsA (+21±5%) in all regions. Ca(2+) overload induced by addition of 20μM Ca(2+) caused a significant decrease of state 3(glu/mal) (-45 to -55%) which was almost completely prevented in the presence of 1μM CsA, 1μM Cs9 or 1μM [D-Ser](8)-CsA. Mitochondrial Ca(2+) accumulation thresholds linked to permeability transition (PT) as well as the rate and completeness of mitochondrial Ca(2+) accumulation differed between different brain regions. For the first time, we provide a detailed, regiospecific analysis of Ca(2+)-dependent properties of brain mitochondria. Regardless of their immunosuppressive impact, CsA and its analogues improved mitochondrial functional properties under control conditions. They also preserved brain mitochondria against Ca(2+) overload-mediated PT and functional impairments. Since Cs9 does not mediate immunosuppression, it might be used as a more specific PT inhibitor than CsA.     

The backrub motion: how protein backbone shrugs when a sidechain dances

Surprisingly, the frozen structures from ultra-high-resolution protein crystallography reveal a prevalent, but subtle, mode of local backbone motion coupled to much larger, two-state changes of sidechain conformation. This "backrub" motion provides an influential and common type of local plasticity in protein backbone. Concerted reorientation of two adjacent peptides swings the central sidechain perpendicular to the chain direction, changing accessible sidechain conformations while leaving flanking structure undisturbed. Alternate conformations in sub-1 angstroms crystal structures show backrub motions for two-thirds of the significant Cbeta shifts and 3% of the total residues in these proteins (126/3882), accompanied by two-state changes in sidechain rotamer. The Backrub modeling tool is effective in crystallographic rebuilding. For homology modeling or protein redesign, backrubs can provide realistic, small perturbations to rigid backbones. For large sidechain changes in protein dynamics or for single mutations, backrubs allow backbone accommodation while maintaining H bonds and ideal geometry.     

Identification of actin and HSP 70 as cyclosporin A binding proteins by photoaffinity labeling and fluorescence displacement assays.

A novel family of cyclosporin A (CsA) binding proteins was identified by using the biologically active, radioiodinated photoaffinity probe [D-Lys-N epsilon-(4-azido-3-[125I]iodophenyl)propionyl)]8-CsA. In addition to cyclophilin, proteins with molecular masses of 43 kDa and approximately 50-55 kDa were labeled in Jurkat extracts and bovine calf thymus. Sequence analysis of the 43-kDa protein purified from calf thymus and subsequent Western analysis of CsA affinity-purified material from Jurkat extracts identified the 43-kDa component as actin. [D-Lys-N epsilon-(5-dimethylamino-1-naphthalenesulfonyl)]8-CsA, a fluorescent analogue of CsA, was prepared and used to measure the binding constants of cyclosporin derivatives to actin by means of a new fluorescence displacement assay. [D-Lys-N epsilon-(5-dimethylamino-1-naphthalenesulfonyl)]8-CsA and [N delta-t-butoxycarbonyl diaminobutyryl)]8-CsA bind to bovine actin at physiologically relevant concentrations, with dissociation constants of 60 +/- 33 and 570 +/- 380 nM, respectively. Because the ATPase fragment of heat shock cognate 70 (HSC 70) is structurally related to actin, the yeast homologue SSA1 was tested and found to be radiolabeled by the cyclosporin A photoaffinity reagent. The binding constant for [D-Lys-N epsilon-(5-dimethylamino-1-naphthalenesulfonyl)]8-CsA to SSA1 was determined and is 53 +/- 48 nM. These results indicate that actin and the 70-kDa heat shock protein family contain a structurally related domain for binding of cyclosporin A-related peptides.     

FK-506, a potent novel immunosuppressive agent, binds to a cytosolic protein which is distinct from the cyclosporin A-binding protein, cyclophilin

A novel macrolide antibiotic, FK-506, isolated from Streptomyces tsukubaensis, has been shown to be a potent immunosuppressive agent in vivo and in vitro. FK-506 shares a number of immunosuppressive properties with the cyclic peptide, cyclosporin A (CsA), although 10 to 100 times more potent in this regard. These similarities suggest that both agents may share a similar mechanism(s) of action at the biochemical level. We have identified a cytoplasmic binding protein for FK-506 in the human T cell line, JURKAT, using [3H]FK-506. The FK-506 binding protein has a mr of 10 to 12 kDa (as determined by gel filtration), is heat stable and does not bind CsA. This contrasts with the CsA binding protein, cyclophilin, in that cyclophilin is heat labile and has a mr of 15 to 17 kDa. Our data suggest that FK-506 binds to a low m.w. protein(s) in JURKAT cells, which is distinct from cyclophilin. This protein may mediate the immunosuppressive effects of FK-506 in T cells. In addition, our results suggest that the immunosuppressive activity of FK-506, as with CsA, is mediated by an intracellular mechanism.     

A computational study of the closed and open states of the influenza a M2 proton channel

In this study, four possible conformations of the His-37 and Trp-41 residues for the closed state of the influenza M2 ion channel were identified by a conformation scan based on a solid-state NMR restraint. In the four conformations, the His-37 residue can be of either the t-160 or t60 rotamer, whereas Trp-41 can be of either the t-105 or t90 rotamer. These conformations were further analyzed by density functional theory calculations and molecular dynamics simulations, and the data indicate that the His-37 residue most likely adopts the t60 rotamer and should be monoprotonated at the delta-nitrogen site, whereas Trp-41 adopts the t90 rotamer. This result is consistent with published experimental data and points to a simple gating mechanism: in the closed state, the His-37 and Trp-41 residues adopt the (t60, t90) conformation, which nearly occludes the pore, preventing nonproton ions from passing through due to the steric and desolvation effects. Moreover, the His-37 tetrad interrupts the otherwise continuous hydrogen-bonding network of the pore water by forcing the water molecules above and below it to adopt opposite orientations, thus adding to the blockage of proton shuttling. The channel can be easily opened by rotating the His-37 chi2 angle from 60 to 0 degrees . This open structure allows pore water to penetrate the constrictive region and to form a continuous water wire for protons to shuttle through, while being still narrow enough to exclude other ions.     

Conformational Heterogeneity of Cyclosporin A in Cyclophilin 18 Binding

The immunosuppressive drug cyclosporin A (CsA) binds to its receptor protein cyclophilin 18 (Cyp18) in two distinct kinetic phases, while the mechanism remains elusive. Stopped-flow measurements coupled with titration and competition experiments were used to investigate the puzzling two-phase process of CsA and Cyp18 interaction. This study leads to the dissection of different conformational fractions of either direct fast binding or slow binding with rate-limiting conformational inter-conversion and the real-time measurement of k_on value (8.34 ± 0.22 x10⁶ M^-1s^-1) in solution. Furthermore, our study indicates that the structure of CsA during dissociation from the protein possesses a distribution of conformations different from those in solution under equilibrium condition.     

Extending the accuracy limits of prediction for side-chain conformations

Current techniques for the prediction of side-chain conformations on a fixed backbone have an accuracy limit of about 1.0-1.5 A rmsd for core residues. We have carried out a detailed and systematic analysis of the factors that influence the prediction of side-chain conformation and, on this basis, have succeeded in extending the limits of side-chain prediction for core residues to about 0.7 A rmsd from native, and 94 % and 89 % of chi(1) and chi(1+2 ) dihedral angles correctly predicted to within 20 degrees of native, respectively. These results are obtained using a force-field that accounts for only van der Waals interactions and torsional potentials. Prediction accuracy is strongly dependent on the rotamer library used. That is, a complete and detailed rotamer library is essential. The greatest accuracy was obtained with an extensive rotamer library, containing over 7560 members, in which bond lengths and bond angles were taken from the database rather than simply assuming idealized values. Perhaps the most surprising finding is that the combinatorial problem normally associated with the prediction of the side-chain conformation does not appear to be important. This conclusion is based on the fact that the prediction of the conformation of a single side-chain with all others fixed in their native conformations is only slightly more accurate than the simultaneous prediction of all side-chain dihedral angles.     

Side-chain modeling with an optimized scoring function

Modeling side-chain conformations on a fixed protein backbone has a wide application in structure prediction and molecular design. Each effort in this field requires decisions about a rotamer set, scoring function, and search strategy. We have developed a new and simple scoring function, which operates on side-chain rotamers and consists of the following energy terms: contact surface, volume overlap, backbone dependency, electrostatic interactions, and desolvation energy. The weights of these energy terms were optimized to achieve the minimal average root mean square (rms) deviation between the lowest energy rotamer and real side-chain conformation on a training set of high-resolution protein structures. In the course of optimization, for every residue, its side chain was replaced by varying rotamers, whereas conformations for all other residues were kept as they appeared in the crystal structure. We obtained prediction accuracy of 90.4% for chi(1), 78.3% for chi(1 + 2), and 1.18 A overall rms deviation. Furthermore, the derived scoring function combined with a Monte Carlo search algorithm was used to place all side chains onto a protein backbone simultaneously. The average prediction accuracy was 87.9% for chi(1), 73.2% for chi(1 + 2), and 1.34 A rms deviation for 30 protein structures. Our approach was compared with available side-chain construction methods and showed improvement over the best among them: 4.4% for chi(1), 4.7% for chi(1 + 2), and 0.21 A for rms deviation. We hypothesize that the scoring function instead of the search strategy is the main obstacle in side-chain modeling. Additionally, we show that a more detailed rotamer library is expected to increase chi(1 + 2) prediction accuracy but may have little effect on chi(1) prediction accuracy.     

Cyclosporin A—cyclophilin complex formation A model based on X-ray and NMR data

The previously determined 3D NMR solution structure of cyclophilin-bound cyclosporin A (CsA) was docked onto the X-ray crystal structure of cyclophilin. Intermolecular nuclear Overhauser effects (NOE) between CsA and cyclophilin were used as constraints in a restrained energy minimization to generate a model of the complex which satisfied all the NOE distance constraints. The model shows that the residues 9 to 11 and 1 to 5 of the cyclic CsA molecule are in contact with cyclophilin. Comparing the model of the CsA—cyclophilin complex to the X-ray crystal structure of a complex of cyclophilin with a substrate for peptidyl-proline cis-trans isomerase activity, i.e. the linear tetrapeptide substrate ae-Ala-Ala-Pro-Ala-amc (ac. acetyl; amc. amidomethylcoumarin), one notices that the contacting peptide segments in the two ligands are oriented in opposite directions, and that the side chain or MeVal-11 of CsA superposes rather precisely with the position of the prolyl residue in ae-Ala-Ala-Pro-Ala-amc.     

Functional analysis of Leishmania major cyclophilin

A potent immunosuppressive drug cyclosporin A (CsA) is known to inhibit human cell infection by the pathogenic protozoan parasite Leishmania major both in vitro and in vivo. The proposed mechanism of action involves CsA binding to Leishmania major-expressed cyclophilin and subsequent down-regulation of signaling events necessary for establishing productive infection. Recently, we identified a ubiquitously expressed membrane protein, CD147, as a signaling receptor for extracellular cyclophilins in mammalian cells. Here we demonstrate that, while being enzymatically active, the Leishmania cyclophilin, unlike its human homologue, does not interact with CD147 on the cell surface of target cells. CD147 facilitates neither Leishmania binding nor infection. Primary structure and biochemical analyses revealed that the parasite's cyclophilin is defective in heparan binding, an event required for signaling interaction between CD147 and human cyclophilin. When the heparan-binding motif was reconstituted in Leishmania cyclophilin, it regained the CD147-dependent signaling activity. These results underscore a critical role of cyclophilin-heparan interactions in CD147-mediated signaling events and argue against the role of Leishmania cyclophilin in parasite binding to target cells.     

Cloning and characterization of Giardia intestinalis cyclophilin

The cyclophilins (Cyps) are family members of proteins that exhibit peptidylprolyl cis-trans isomerase (PPIase, EC 5.2.1.8) activity and bind the immunosuppressive agent cyclosprin A (CsA) in varying degrees. During the process of random sequencing of a cDNA library made from Giardia intestinalis WB strain, the cyclophilin gene (gicyp 1) was isolated. An open reading frame of gicyp 1 gene was 576 nucleotides, which corresponded to a translation product of 176 amino acids (Gicyp 1). The identity with other Cyps was about 58-71%. The 13 residues that constituted the CsA binding site of human cyclophilin were also detected in the amino acid sequence of Gicyp 1, including tryptophan residue essential for the drug binding. The single copy of the gicyp 1 gene was detected in the G. intestinalis chromosome by southern hybridization analysis. Recombinant Gicyp 1 protein clearly accelerated the rate of cis-->trans isomerization of the peptide substrate and the catalysis was completely inhibited by the addition of 0.5 microM CsA.     

Anti-Trypanosoma cruzi effects of cyclosporin A derivatives: possible role of a P-glycoprotein and parasite cyclophilins

Cyclophilins are target molecules for cyclosporin A (CsA), an immunosuppressive antimicrobial drug. We have previously reported the in vitro anti-Trypanosoma cruzi activity of H-7-94 and F-7-62 non-immunosuppressive CsA analogues. In this work, we continue the study of the parasiticidal effect of H-7-94 and F-7-62 CsA analogues in vitro and in vivo and we analyse 3 new CsA derivatives: MeIle-4-CsA (NIM 811), MeVal-4-CsA (MeVal-4) and D-MeAla-3-EtVal-4-CsA, (EtVal-4). The most efficient anti-T. cruzi effect was observed with H-7-94, F-7-62 and MeVal-4 CsA analogues evidenced as inhibition of epimastigote proliferation, trypomastigote penetration, intracellular amastigote development and in vivo T. cruzi infection. This trypanocidal activity could be due to inhibition of the peptidyl prolyl cis-trans isomerase activity on the T. cruzi recombinant cyclophilins tested. Furthermore, CsA and F-7-62 derivative inhibited the efflux of rhodamine 123 from T. cruzi epimastigotes, suggesting an interference with a P-glycoprotein activity. Moreover, H-7-94 and F-7-62 CsA analogues were not toxic as shown by cell viability and by aminopyrine-N-demethylase activity on mammalian cells. Our results show that H-7-94, F-7-62 and MeVal-4 CsA analogues expressed the highest inhibiting effects on T. cruzi, being promissory parasiticidal drugs worthy of further studies.