An investigation into a cardiolipin acyl chain insertion site in cytochrome c

Mitochondrial cytochrome c associates with the phosphoplipid cardiolipin (CL) through a combination of electrostatic and hydrophobic interactions. The latter occurs by insertion into cytochrome c of an acyl chain, resulting in the dissociation of the axial Met-80 heme-iron ligand. The resulting five coordinate cytochrome c/CL complex has peroxidatic properties leading to peroxidation of CL and dissociation of the complex. These events are considered to be pre-apoptotic and culminate with release of cytochrome c from the mitochondria into the cytoplasm. Two distinct surface regions on cytochrome c have been suggested to mediate CL acyl chain insertion and this study has probed one of these regions. We have constructed a series of alanine mutants aimed at disrupting a surface cleft formed between residues 67-71 and 82-85. The physicochemical properties, peroxidase activity, CL binding, and kinetics of carbon monoxide (CO) binding to the ferrous cytochrome c/CL complex have been assessed for the individual mutants. Our findings reveal that the majority of mutants are capable of binding CL in the same apparent stoichiometry as the wild-type protein, with the extent to which the Met-80 ligand is bound in the ferrous cytochrome c/CL complex being mutant specific at neutral pH. Mutation of the species conserved Arg-91 residue, that anchors the cleft, results in the greatest changes to physicochemical properties of the protein leading to a change in the CL binding ratio required to effect structural changes and to the ligand-exchange properties of the ferrous cytochrome c/CL complex.     

Recognition and fluorescence sensing of specific amino acid residue on protein surface using designed molecules

Many biological processes are mediated by surface recognition between proteins. Small molecules that recognize and bind a specific region of a protein surface may be promising agents for disrupting certain protein-protein surface interactions, which consequently leads to regulation of cellar functions. This article describes our recent efforts toward the development of the designed small molecules, which can recognize histidine or phosphorylated amino acid residues on peptide surfaces in a sequence-selective manner. These results demonstrate that cooperative metal-ligand interaction is powerful for tight and selective binding to the specific amino acid residues of proteins in aqueous medium.     

Surface plasmon resonance studies of complex formation between cytochrome c and bovine cytochrome c oxidase incorporated into a supported planar lipid bilayer. II. Binding of cytochrome c to oxidase-containing cardiolipin/phosphatidylcholine membranes.

Complex formation between horse heart cytochrome c (cyt c) and bovine cytochrome c oxidase (cco) incorporated into a supported planar egg phosphatidylcholine membrane containing varying amounts of cardiolipin (CL) (0-20 mol%) has been studied under low (10 mM) and medium (160 mM) ionic strength conditions by surface plasmon resonance (SPR) spectroscopy. Both specific and nonspecific modes of cyt c binding are observed. The dissociation constant of the specific interaction between cyt c and cco increases from approximately 6.5 microM at low ionic strength to 18 microM at medium ionic strength, whereas the final saturation level of bound protein is independent of salt concentration and corresponds to approximately 53% of the total cco molecules present in the membrane. This suggests a 1:1 binding stoichiometry between the two proteins. The nonspecific binding component is governed by electrostatic interactions between cyt c and the membrane lipids and results in a partially ionic strength-reversible protein-membrane association. Thus, hydrophobic interactions between cyt c and the membrane, which are the predominant mode of binding in the absence of cco, are greatly suppressed. Both the amount of nonspecifically bound protein and the binding affinity can be varied over a broad range by changing the ionic strength and the extent of CL incorporation into the membrane. Under conditions approximating the physiological state in the mitochondrion (i.e., 20 mol% CL and medium ionic strength), 1-1.5 cyt c molecules are bound to the lipid phase per molecule of cco, with a dissociation constant of 0.1 microM. The possible physiological significance of these observations is discussed.     

Bilayer structure in phospholipid-cytochrome c model membranes.

Lipid protein interactions in biological membranes differ markedly depending on whether the protein is intrinsic or extrinsic. These interactions are studied using lipid spin labels diffused into model systems consisting of phospholipid bilayers and a specific protein. Recently, an intrinsic protein complex, cytochrome oxidase, was examined and the data suggest there is a boundary layer of immobilized lipid between the hydrophobic protein surfaces and adjacent fluid bilayer regions. In the present study, a typical extrinsic protein, cytochrome c, was complexed with a cardiolipin/lecithin (1:4 by weight) mixture. The phospholipids in the presence and absence of cytochrome c exhibit typical bilayer behavior as jedged by four spin-labeling criteria: fluidity gradient, spectral anisotropy of oriented bilayers, response to hydration and the polarity profile. Any effects of cytochrome c on the ESR spectra of lipid spin labels are small, in contrast to the effects of intrinsic proteins. These data are consistent with electrostatic binding of cytochrome c to the charged groups of the phospholipids, and indicate that the presence of extrinsic proteins will not interfere with measurements of boundary lipid in intact biological membranes.     

Definition of the interaction domain for cytochrome c on cytochrome c oxidase. III. Prediction of the docked complex by a complete, systematic search

The electron transfer complex between bovine cytochrome c oxidase and horse cytochrome c has been predicted with the docking program DOT, which performs a complete, systematic search over all six rotational and translational degrees of freedom. Energies for over 36 billion configurations were calculated, providing a free-energy landscape showing guidance of positively charged cytochrome c to the negative region on the cytochrome c oxidase surface formed by subunit II. In a representative configuration, the solvent-exposed cytochrome c heme edge is within 4 A of the indole ring of subunit II residue Trp(104), indicating a likely electron transfer path. These two groups are surrounded by a small, hydrophobic contact region, which is surrounded by electrostatically complementary hydrophilic interactions. Cytochrome c/cytochrome c oxidase interactions of Lys(13) with Asp(119) and Lys(72) with Gln(103) and Asp(158) are the most critical polar interactions due to their proximity to the hydrophobic region and exclusion from bulk solvent. The predicted complex matches previous mutagenesis, binding, and time-resolved kinetics studies that implicate Trp(104) in electron transfer and show the importance of specific charged residues to protein affinity. Electrostatic forces not only enhance long range protein/protein association; they also predominate in short range alignment, creating the transient interaction needed for rapid turnover.     

Desulfovibrio desulfuricans G20 tetraheme cytochrome structure at 1.5 Angstrom and cytochrome interaction with metal complexes

The structure of the type I tetraheme cytochrome c(3) from Desulfovibrio desulfuricans G20 was determined to 1.5 Angstrom by X-ray crystallography. In addition to the oxidized form, the structure of the molybdate-bound form of the protein was determined from oxidized crystals soaked in sodium molybdate. Only small structural shifts were obtained with metal binding, consistent with the remarkable structural stability of this protein. In vitro experiments with pure cytochrome showed that molybdate could oxidize the reduced cytochrome, although not as rapidly as U(VI) present as uranyl acetate. Alterations in the overall conformation and thermostability of the metal-oxidized protein were investigated by circular dichroism studies. Again, only small changes in protein structure were documented. The location of the molybdate ion near heme IV in the crystal structure suggested heme IV as the site of electron exit from the reduced cytochrome and implicated Lys14 and Lys56 in binding. Analysis of structurally conserved water molecules in type I cytochrome c(3) crystal structures identified interactions predicted to be important for protein stability and possibly for intramolecular electron transfer among heme molecules.     

How does a small molecule bind at a cryptic binding site?

Protein-protein interactions (PPIs) are ubiquitous biomolecular processes that are central to virtually all aspects of cellular function. Identifying small molecules that modulate specific disease-related PPIs is a strategy with enormous promise for drug discovery. The design of drugs to disrupt PPIs is challenging, however, because many potential drug-binding sites at PPI interfaces are “cryptic”: When unoccupied by a ligand, cryptic sites are often flat and featureless, and thus not readily recognizable in crystal structures, with the geometric and chemical characteristics of typical small-molecule binding sites only emerging upon ligand binding. The rational design of small molecules to inhibit specific PPIs would benefit from a better understanding of how such molecules bind at PPI interfaces. To this end, we have conducted unbiased, all-atom MD simulations of the binding of four small-molecule inhibitors (SP4206 and three SP4206 analogs) to interleukin 2 (IL2)—which performs its function by forming a PPI with its receptor—without incorporating any prior structural information about the ligands’ binding. In multiple binding events, a small molecule settled into a stable binding pose at the PPI interface of IL2, resulting in a protein–small-molecule binding site and pose virtually identical to that observed in an existing crystal structure of the IL2-SP4206 complex. Binding of the small molecule stabilized the IL2 binding groove, which when the small molecule was not bound emerged only transiently and incompletely. Moreover, free energy perturbation (FEP) calculations successfully distinguished between the native and non-native IL2–small-molecule binding poses found in the simulations, suggesting that binding simulations in combination with FEP may provide an effective tool for identifying cryptic binding sites and determining the binding poses of small molecules designed to disrupt PPI interfaces by binding to such sites.     

The interaction of NADH-cytochrome b5 reductase and cytochrome b5 bound to egg lecithin liposomes.

Incubation of liposomes prepared by sonication of egg lecithin with the amphipathic form of cytochrome b5 results in the binding of a maximum of 244 molecules of cytochrome b5 per liposomal vesicle. Interactions of the phospholipid with the hydrophobic segment of cytochrome b5 are involved in this binding which does not disrupt the liposome. When a small amount of NADH-cytochrome b5 reductase is bound liposomes simultaneously with cytochrome b5, the two proteins catalyze the reduction of cytochrome c by NADH. A qualitative kinetic analysis reveals that all of the cytochrome b5 interacts with reductase, a result consistent with these protein undergoing translational diffusion in the plane of the membrane. This system and the purified stearyl coenzyme A desaturase provide a model to study the dynamics of protein andlipid interactions in this membrane-bound oxidative sequence.     

The Overlap of Small Molecule and Protein Binding Sites within Families of Protein Structures

Protein–protein interactions are challenging targets for modulation by small molecules. Here, we propose an approach that harnesses the increasing structural coverage of protein complexes to identify small molecules that may target protein interactions. Specifically, we identify ligand and protein binding sites that overlap upon alignment of homologous proteins. Of the 2,619 protein structure families observed to bind proteins, 1,028 also bind small molecules (250–1000 Da), and 197 exhibit a statistically significant (p<0.01) overlap between ligand and protein binding positions. These “bi-functional positions”, which bind both ligands and proteins, are particularly enriched in tyrosine and tryptophan residues, similar to “energetic hotspots” described previously, and are significantly less conserved than mono-functional and solvent exposed positions. Homology transfer identifies ligands whose binding sites overlap at least 20% of the protein interface for 35% of domain–domain and 45% of domain–peptide mediated interactions. The analysis recovered known small-molecule modulators of protein interactions as well as predicted new interaction targets based on the sequence similarity of ligand binding sites. We illustrate the predictive utility of the method by suggesting structural mechanisms for the effects of sanglifehrin A on HIV virion production, bepridil on the cellular entry of anthrax edema factor, and fusicoccin on vertebrate developmental pathways. The results, available at http://pibase.janelia.org, represent a comprehensive collection of structurally characterized modulators of protein interactions, and suggest that homologous structures are a useful resource for the rational design of interaction modulators.     

Discovery and characterization of cooperative ligand binding in the adaptive region of interleukin-2

The cytokine hormone interleukin-2 (IL-2) contains a highly adaptive region that binds small, druglike molecules. The binding properties of this adaptive region have been explored using a "tethering" method that relies on the formation of a disulfide bond between the protein and small-molecule ligands. Using tethering, surface plasmon resonance (SPR), and X-ray crystallography, we have discovered that the IL-2 adaptive region contains at least two cooperative binding sites where the binding of a first ligand to one site promotes or antagonizes the binding of a second ligand to the second site. Cooperative energies of interaction of -2 kcal/mol are observed. The observation that the adaptive region contains two adjacent sites may lead to the development of tight-binding antagonists of a protein-protein interaction. Cooperative ligand binding in the adaptive region of IL-2 underscores the importance of protein dynamics in molecular recognition. The tethering approach provides a novel and general strategy for discovering such cooperative binding interactions in specific, flexible regions of protein structure.     

Computational design of protein-small molecule interfaces

The computational design of proteins that bind small molecule ligands is one of the unsolved challenges in protein engineering. It is complicated by the relatively small size of the ligand which limits the number of intermolecular interactions. Furthermore, near-perfect geometries between interacting partners are required to achieve high binding affinities. For apolar, rigid small molecules the interactions are dominated by short-range van der Waals forces. As the number of polar groups in the ligand increases, hydrogen bonds, salt bridges, cation–π, and π–π interactions gain importance. These partial covalent interactions are longer ranged, and additionally, their strength depends on the environment (e.g. solvent exposure). To assess the current state of protein-small molecule interface design, we benchmark the popular computer algorithm Rosetta on a diverse set of 43 protein–ligand complexes. On average, we achieve sequence recoveries in the binding site of 59% when the ligand is allowed limited reorientation, and 48% when the ligand is allowed full reorientation. When simulating the redesign of a protein binding site, sequence recovery among residues that contribute most to binding was 52% when slight ligand reorientation was allowed, and 27% when full ligand reorientation was allowed. As expected, sequence recovery correlates with ligand displacement.     

Vectorially oriented monolayers of the cytochrome c/cytochrome oxidase bimolecular complex.

Vectorially oriented monolayers of yeast cytochrome c and its bimolecular complex with bovine heart cytochrome c oxidase have been formed by self-assembly from solution. Both quartz and Ge/Si multilayer substrates were chemical vapor deposited with an amine-terminated alkylsiloxane monolayer that was then reacted with a hetero-bifunctional cross-linking reagent, and the resulting maleimide endgroup surface then provided for covalent interactions with the naturally occurring single surface cysteine 102 of the yeast cytochrome c. The bimolecular complex was formed by further incubating these cytochrome c monolayers in detergent-solubilized cytochrome oxidase. The sequential formation of such monolayers and the vectorially oriented nature of the cytochrome oxidase was studied via meridional x-ray diffraction, which directly provided electron density profiles of the protein(s) along the axis normal to the substrate plane. The nature of these profiles is consistent with previous work performed on vectorially oriented monolayers of either cytochrome c or cytochrome oxidase alone. Furthermore, optical spectroscopy has indicated that the rate of binding of cytochrome oxidase to the cytochrome c monolayer is an order of magnitude faster than the binding of cytochrome oxidase to an amine-terminated surface that was meant to mimic the ring of lysine residues around the heme edge of cytochrome c, which are known to be involved in the binding of this protein to cytochrome oxidase.     

Structure of ferricytochrome c' from Rhodospirillum molischianum at 1.67 A resolution

The structure of ferricytochrome c' from Rhodospirillum molischianum has been crystallographically refined to 1.67 A resolution using a combination of reciprocal space and restrained least-squares refinement methods. The final crystallographic R-factor for 30,533 reflections measured with I greater than sigma (I) between infinity and 1.67 A is 0.188. The final model incorporates 1944 unique protein atoms (of a total of 1972) together with 194 bound solvent molecules. The structure has been analysed with respect to its detailed conformational properties, secondary structural features, temperature factor behavior, bound solvent sites, and heme geometry. The asymmetric unit of the cytochrome c' crystal contains a dimer composed of chemically identical 128-residue polypeptide chains. Although the refined structure shows the monomers to be very similar, examination of the differences that do occur allows an evaluation of how different lattice contacts affect protein conformation and solvent binding. In particular, comparison of solvent binding sites in the two subunits allows identification of a common set that are not altered by lattice interactions. The preservation of these solvent interactions in different lattice environments suggests that they play a structural role in protein stabilization in solution. The refined structure additionally reveals some new features that relate to the ligand binding properties and unusual mixed-spin state character of cytochrome c'. Finally, comparison of the heme binding geometry in cytochrome c' and other structurally unrelated c-type cytochromes shows that two alternative, but sterically favorable, conformational variants occur among the seven examples examined.     

FGO: a novel ontology for identification of ligand functional group

Small molecules play crucial role in the modulation of biological functions by interacting with specific macromolecules. Hence small molecule interactions are captured by a variety of experimental methods to estimate and propose correlations between molecular structures to their biological activities. The tremendous expanse in publicly available small molecules is also driving new efforts to better understand interactions involving small molecules particularly in area of drug docking and pharmacogenomics. We have studied and designed a functional group identification system with the associated ontology for it. The functional group identification system can detect the functional group components from given ligand structure with specific coordinate information. Functional group ontology (FGO) proposed by us is a structured classification of chemical functional group which acts as an important source of prior knowledge that may be automatically integrated to support identification, categorization and predictive data analysis tasks. We have used a new annotation method which can be used to construct the original structure from given ontological expression using exact coordinate information. Here, we also discuss about ontology-driven similarity measure of functional groups and uses of such novel ontology for pharmacophore searching and de-novo ligand designing.     

Properties of small molecules affecting insulin receptor function.

Small molecules with insulin mimetic effects and oral availability are of interest for potential substitution of insulin injections in the treatment of diabetes. We have searched databases for compounds capable of mimicking one epitope of the insulin molecule known to be involved in binding to the insulin receptor (IR). This approach identifies thymolphthalein, which is an apparent weak agonist that displaces insulin from its receptor, stimulates auto- and substrate phosphorylation of IR, and potentiates lipogenesis in adipocytes in the presence of submaximal concentrations of insulin. The various effects are observed in the 10(-5)-10(-3) M range of ligand concentration and result in partial insulin activity. Furthermore, analogues of the related phenol red and fluorescein molecules fully displace insulin from the IR ectodomain, however, without insulin agonistic effects. The interactions are further characterized by NMR, UV-vis, and fluorescence spectroscopies. It is shown that both fluorescence and UV-vis changes in the ligand spectra induced by IR fragments occur with Kd values similar to those obtained in the displacement assay. Nevertheless, insulin itself cannot completely abolish binding of the small molecules. Determination of the binding stoichiometry reveals multiple binding sites for ligands of which one overlaps with the insulin binding site on the receptor.     

Structural parameterization of the binding enthalpy of small ligands

A major goal in ligand and drug design is the optimization of the binding affinity of selected lead molecules. However, the binding affinity is defined by the free energy of binding, which, in turn, is determined by the enthalpy and entropy changes. Because the binding enthalpy is the term that predominantly reflects the strength of the interactions of the ligand with its target relative to those with the solvent, it is desirable to develop ways of predicting enthalpy changes from structural considerations. The application of structure/enthalpy correlations derived from protein stability data has yielded inconsistent results when applied to small ligands of pharmaceutical interest (MW < 800). Here we present a first attempt at an empirical parameterization of the binding enthalpy for small ligands in terms of structural information. We find that at least three terms need to be considered: (1) the intrinsic enthalpy change that reflects the nature of the interactions between ligand, target, and solvent; (2) the enthalpy associated with any possible conformational change in the protein or ligand upon binding; and, (3) the enthalpy associated with protonation/deprotonation events, if present. As in the case of protein stability, the intrinsic binding enthalpy scales with changes in solvent accessible surface areas. However, an accurate estimation of the intrinsic binding enthalpy requires explicit consideration of long-lived water molecules at the binding interface. The best statistical structure/enthalpy correlation is obtained when buried water molecules within 5-7 A of the ligand are included in the calculations. For all seven protein systems considered (HIV-1 protease, dihydrodipicolinate reductase, Rnase T1, streptavidin, pp60c-Src SH2 domain, Hsp90 molecular chaperone, and bovine beta-trypsin) the binding enthalpy of 25 small molecular weight peptide and nonpeptide ligands can be accounted for with a standard error of +/-0.3 kcal x mol(-1).     

Druggable Protein Interaction Sites Are More Predisposed to Surface Pocket Formation than the Rest of the Protein Surface

Despite intense interest and considerable effort via high-throughput screening, there are few examples of small molecules that directly inhibit protein-protein interactions. This suggests that many protein interaction surfaces may not be intrinsically “druggable” by small molecules, and elevates in importance the few successful examples as model systems for improving our fundamental understanding of druggability. Here we describe an approach for exploring protein fluctuations enriched in conformations containing surface pockets suitable for small molecule binding. Starting from a set of seven unbound protein structures, we find that the presence of low-energy pocket-containing conformations is indeed a signature of druggable protein interaction sites and that analogous surface pockets are not formed elsewhere on the protein. We further find that ensembles of conformations generated with this biased approach structurally resemble known inhibitor-bound structures more closely than equivalent ensembles of unbiased conformations. Collectively these results suggest that “druggability” is a property encoded on a protein surface through its propensity to form pockets, and inspire a model in which the crude features of the predisposed pocket(s) restrict the range of complementary ligands; additional smaller conformational changes then respond to details of a particular ligand. We anticipate that the insights described here will prove useful in selecting protein targets for therapeutic intervention.     

The B10.A mouse B cell response to pigeon cytochrome c is directed against the same area of the protein that is recognized by B10.A T cells in association with the Ek beta:Ek alpha Ia molecule

An analysis of the fine specificities of the primary and hyperimmune antibody responses of B10.A mice to pigeon cytochrome c showed that both were qualitatively very similar. Small amounts of antibody appeared to be directed against the regions of serine 15 and/or glutamic acid 44. The remaining antibodies (greater than 70%) bound to the same complex topographic determinant (including residues 3, 103, and 104) on the back surface of pigeon cytochrome c which had been found to dominate the rabbit antibody response to this protein, and to be involved in Ia-restricted T cell stimulation. The mouse antibodies reacted very poorly with fragmented forms of the immunogen or with tobacco hornworm moth cytochrome c, even though both of these antigens had been shown previously to strongly stimulate pigeon cytochrome c-primed T cells. The specificities of the primary IgG responses of seven other mouse strains were found to be very similar, but not identical, to that of B10.A mice. The cytochrome c-specific antibodies in the hyperimmune serum were shown to bind to determinants involving residues that vary between pigeon and mouse cytochromes c. Comparison of the binding of the antibodies to the immunogen and to the corresponding host protein enabled the calculation of the proportion of the overall binding energy contributed by the variant residues. This was as low as 19 to 35% for the primary response, rose to 25 to 46% for the hyperimmune mouse antibodies, and reached 40 to 63% for hyperimmune rabbit antibodies. The remaining energy of interaction (37 to 81%) was necessarily contributed by the surface of the protein surrounding the variant residues, which is the same for the immunogen and the host protein. These results illustrate the relatively subtle differences in binding affinities which can distinguish self from non-self recognition by antibody molecules.     

Excision of a proposed electron transfer pathway in cytochrome c peroxidase and its replacement by a ligand-binding channel

A previously proposed electron transfer (ET) pathway in the heme enzyme cytochrome c peroxidase has been excised from the structure, leaving an open ligand-binding channel in its place. Earlier studies on cavity mutants of this enzyme have revealed structural plasticity in this region of the molecule. Analysis of these structures has allowed the design of a variant in which the specific section of protein backbone representing a previously proposed ET pathway is accurately extracted from the protein. A crystal structure verified the creation of an open channel that overlays the removed segment, extending from the surface of the protein to the heme at the core of the protein. A number of heterocyclic cations were found to bind to the proximal-channel mutant with affinities that can be rationalized based on the structures. It is proposed that small ligands bind more weakly to the proximal-channel mutant than to the W191G cavity due to an increased off rate of the open channel, whereas larger ligands are able to bind to the channel mutant without inducing large conformational changes. The structure of benzimidazole bound to the proximal-channel mutant shows that the ligand accurately overlays the position of the tryptophan radical center that was removed from the wild-type enzyme and displaces four of the eight ordered solvent molecules seen in the empty cavity. Ligand binding also caused a small rearrangement of the redesigned protein loop, perhaps as a result of improved electrostatic interactions with the ligand. The engineered channel offers the potential for introducing synthetic replacements for the removed structure, such as sensitizer-linked substrates. These installed "molecular wires" could be used to rapidly initiate reactions, trap reactive intermediates, or answer unresolved questions about ET pathways.