Weak interactions and molecular recognition in systems involving electron transfer proteins

Electrostatic interactions and other weak interactions between amino acid side chains on protein surfaces play important roles in molecular recognition, and the mechanism of their intermolecular interactions has gained much interest. We established that charged peptides are useful for investigating the molecular recognition character of proteins and their molecular interaction induced structural changes. Positively charged lysine peptides competitively inhibited electron transfer from reduced cytochrome f (cyt f or cytochrome c (cyt c) to oxidized plastocyanin (PC), due to neutralization of the negatively charged site of PC by formation of PC-lysine peptide complexes. Lysine peptides also inhibited electron transfer from cyt c to cytochrome c peroxidase. Likewise, negatively charged aspartic acid peptides interacted with the positively charged sites of cytfand cyt c, and competitively inhibited electron transfer from reduced cytfor cyt c to oxidized PC and from [Fe(CN)6]4- to oxidized cyt c. Changes in the geometry and a shift to a higher redox potential of the active site Cu of PC on oligolysine binding were detected by spectroscopic and electrochemical measurements, owing to the absence of absorption in the visible region for lysine peptides. Structural and redox potential changes were also observed for cyt f and cyt c by interaction with aspartic acid peptides.     

Simulation of electron transfer between cytochrome C2 and the bacterial photosynthetic reaction center: Brownian dynamics analysis of the native proteins and double mutants

Electron transfer is essential for bacterial photosynthesis which converts light energy into chemical energy. This paper theoretically studies the interprotein electron transfer from cytochrome c(2) of Rhodobacter capsulatus to the photosynthetic reaction center of Rhodobacter sphaeroides in native and mutated systems. Brownian dynamics is used with an exponential distance-dependent electron-transfer rate model to compute bimolecular rate constants, which are consistent with experimental data when reasonable prefactors and decay constants are used. Interestingly, switching of the reaction mechanism from the diffusion-controlled limit in the native proteins to the activation-controlled limit in one of the mutants (DK(L261)/KE(C99)) was found. We also predict that the second-order rate for the native reaction center/cytochrome c(2) system will decrease with increasing ionic strength, a characteristic of electrostatically controlled docking.     

表面带负电荷氨基酸残基的突变对细胞色素b5与细胞色素c结合和识别的影响

利用核磁共振方法研究表面带不同负电荷氨基酸残基突变后的细胞色素b₅与细胞色素c的结合与识别.结果表明,静电作用在细胞色素b₅与细胞色素c的结合过程中有着重要的贡献,而且这些静电贡献在一定程度上具有累加性,E48的贡献略大于E44.同时还证明Browniandynamicssimulations优化出的Glu48-Lys13,Glu56-Lys₈₇,Asp60-Lys86和heme6-propionate-Tml₇₂(细胞色素b₅的残基排在前面)的结合方式在溶液中的确存在.细胞色素b₅突变体(E48,E56/A,D60/A)及[Cr(oxalate)₃]^3-对细胞色素c的表面结合竞争实验表明,细胞色素c表面结合区Site仍然同细胞色素b₅突变体(E48,E56/A,D60/A)有结合作用,只是结合强度上相对于野生细胞色素b₅同细胞色素c的结合有所降低.这表明除上述的Brownian dynamics simulations模型外,尚有其它如Salemme模型等的结合方式,这也揭示出细胞色素b₅和细胞色素c之间的结合是比较动态的.     

Context-based identification of protein-protein interfaces and "hot-spot" residues

Reliable determination of protein-protein interaction sites is of critical importance for structure-based design of small molecules modulating protein function through macromolecular interfaces. We present an alignment-free computational method for prediction of protein-protein interface residues. The method ("iPred") is based on a knowledge-based scoring function adapted from the field of protein folding and small molecule docking. Based on a training set of 394 hetero-dimeric proteins iPred achieves sustained accuracy on an external unbound test set. Prediction robustness was assessed from more than 1500 diverse complexes containing homo- and hetero-dimers. The technique does not rely on sequence conservation, so that rapid interface identification is possible even for proteins for which homologs are unknown or lack conserved residue patterns in interface region. Functional "hot-spot" residues are enriched among the predicted interface residues, rendering the method predestined for macromolecular binding site identification and drug design studies aiming at modulating protein-protein interaction that might influence protein function. For a comparative structural model of peptidase HtrA from Helicobacter pylori, we performed mutation studies for predicted hot-spot residues, which were confirmed as functionally relevant for HtrA activity or oligomerization.     

Cytochrome c(2) Exit Strategy: Dissociation Studies and Evolutionary Implications

Small, water-soluble, type c cytochromes form a transient network connecting major bioenergetic membrane protein complexes in both photosynthesis and respiration. In the photosynthesis cycle of Rhodobacter sphaeroides, cytochrome c2 (cyt c2) docks to the reaction center (RC), undergoes electron transfer, and exits for the cytochrome bc1 complex. Translations of cyt c2 about the RC-cyt c2 docking interface and surrounding membrane reveal possible exit pathways. A pathway at a minimal elevation allowed by the architecture of the RC is analyzed using both an all-atom steered molecular dynamics simulation of the RC-cyt c2 complex and a bioinformatic analysis of the structures and sequences of cyt c. The structure-based phylogenetic analysis allows for the identification of structural elements that have evolved to satisfy the requirements of having multiple functional partners. The patterns of evolutionary variation obtained from the phylogenetic analysis of both docking partners of cyt c2 reveal conservation of key residues involved in the interaction interfaces that would be candidates for further experimental studies. Additionally, using the molecular mechanics Poisson-Boltzmann surface area method we calculate that the binding free energy of reduced cyt c2 to the RC is nearly 6 kcal/mol more favorable than with oxidized cyt c2. The redox-dependent variations lead to changes in structural flexibility, behavior of the interfacial water molecules, and eventually changes in the binding free energy of the complex.     

Metalloprotein association,self-association,and dynamics governed by hydrophobic interactions: simultaneous occurrence of gated and true electron-transfer reactions between cytochrome f and cytochrome c(6) from Chlamydomonas reinhardtii

Noninvasive reconstitution of the heme in cytochrome c(6) with zinc(II) ions allowed us to study the photoinduced electron-transfer reaction (3)Zncyt c(6) + cyt f(III) --> Zncyt c(6)(+) + cyt f(II) between physiological partners cytochrome c(6) and cytochrome f, both from Chlamydomonas reinhardtii. The reaction kinetics was analyzed in terms of protein docking and electron transfer. In contrast to various protein pairs studied before, both the unimolecular and the bimolecular reactions of this oxidative quenching take place at all ionic strengths from 2.5 through 700 mM. The respective intracomplex rate constants are k(uni) (1.2 +/- 0.1) x 10(4) s(-1) for persistent and k(bi) (9 +/- 4) x 10(2) s(-1) for the transient protein complex. The former reaction seems to be true electron transfer, and the latter seems to be electron transfer gated by a structural rearrangement. Remarkably, these reactions occur simultaneously, and both rate constants are invariant with ionic strength. The association constant K(a) for zinc cytochrome c(6) and cytochrome f(III) remains (5 +/- 3) x 10(5) M(-1) in the ionic strength range from 700 to 10 mM and then rises slightly to (7 +/- 2) x 10(6) M(-1), as ionic strength is lowered to 2.5 mM. Evidently, docking of these proteins from C. reinhardtii is due to hydrophobic interaction, slightly augmented by weak electrostatic attraction. Kinetics, chromatography, and cross-linking consistently show that cytochrome f self-dimerizes at ionic strengths of 200 mM and higher. Cytochrome f(III) quenches triplet state (3)Zncyt c(6), but its dimer does not. Formation of this unreactive dimer is an important step in the mechanism of electron transfer. Not only association between the reacting proteins, but also their self-association, should be considered when analyzing reaction mechanisms.     

Electrostatic docking of a supramolecular host-guest assembly to cytochrome c probed by bidirectional photoinduced electron transfer

A water-soluble octacarboxyhemicarcerand was used as a shuttle to transport redox-active substrates across the aqueous medium and deliver them to the target protein. The results show that weak multivalent interactions and conformational flexibility can be exploited to reversibly bind complex supramolecular assemblies to biological molecules. Hydrophobic electron donors and acceptors were encapsulated within the hemicarcerand, and photoinduced electron transfer (ET) between the Zn-substituted cytochrome c (MW = 12.3 kD) and the host-guest complexes (MW = 2.2 kD) was used to probe the association between the negatively charged hemicarceplex and the positively charged protein. The behavior of the resulting ternary protein-hemicarcerand-guest assembly was investigated in two binding limits: (1) when K(encaps) ? K(assoc), the hemicarcerand transports the ligand to the protein while protecting it from the aqueous medium; and (2) when K(assoc) > K(encaps), the hemicarcerand-protein complex is formed first, and the hemicarcerand acts as an artificial receptor site that intercepts ligands from solution and positions them close to the active site of the metalloenzyme. In both cases, ET mediated by the protein-bound hemicarcerand is much faster than that due to diffusional encounters with the respective free donor or acceptor in solution. The measured ET rates suggest that the dominant binding region of the host-guest complex on the surface of the protein is consistent with the docking area of the native redox partner of cytochrome c. The strong association with the protein is attributed to the flexible conformation and adaptable charge distribution of the hemicarcerand, which allow for surface-matching with the cytochrome.     

Proteo-dendrimers designed for complementary recognition of cytochrome c: dendrimer architecture toward nanoscale protein complexation

"Proteo-dendrimers" in which polyanionic hepta(glutamic acids), fluorescent zinc porphyrinate cores, hydrophilic polyether surfaces, and nonpeptide hydrophobic dendrons are combined, were developed as a new series of synthetic receptors for protein recognition. They have polyanionic "patch" structures on their surfaces and undergo complementary electrostatic interactions with a positively charged cytochrome c patch, as observed in biological protein-protein recognition systems. Stability constants of the resulting supramolecular complexes were determined in phosphate buffer (pH 7) by monitoring the fluorescence quenching of the zinc porphyrinates. These proteo-dendrimer receptors exhibited higher affinities with cytochrome c proteins in aqueous solutions than with biological cytochrome b5. Furthermore, they effectively blocked complexation of biological cytochrome b5 with cytochrome c, indicating that the proteo-dendrimers and cytochrome b5 similarly occupy the polycationic patch of cytochrome c.     

鳞状上皮细胞癌抗原2与组蛋白G的柔性对接研究

蛋白质与蛋白质间的相互作用在生物体内具有极其重要的意义 .确定蛋白质间的相互作用的机制对生物调节和药物设计都是关键的步骤 .近年来 ,随着计算机和软件技术的发展 ,使建立蛋白质 -蛋白质复合物模型成为可能 ,为研究蛋白质间的相互作用提供了直观和比较准确的参考 .目前的蛋白质 -蛋白质对接软件将对接的两个蛋白质看成刚性 ,并不考虑蛋白质对接时的柔性 ,即没有考虑由于界面氨基酸间的相互作用而导致的界面结构的变化 .本文在直接的蛋白质 -蛋白质刚性对接的基础上 ,应用动力学模拟的方法 ,对刚性对接的结构进行蛋白质界面的构象搜索 ,…     

Redox-coupled dynamics and folding in cytochrome c

Cytochrome c functions as an electron carrier in the mitochondrial electron-transport chain using the Fe(II)-Fe(III) redox couple of a covalently attached heme prosthetic group, and it has served as a paradigm for both biological redox activity and protein folding. On the basis of a wide variety of biophysical techniques, it has been suggested that the protein is more flexible in the oxidized state than in the reduced state, which has led to speculation that it is the dynamics of the protein that has been evolved to control the cofactor's redox properties. To test this hypothesis, we incorporated carbon-deuterium bonds throughout cytochrome c and characterized their absorption frequencies and line widths using IR spectroscopy. The absorption frequencies of several residues on the proximal side of the heme show redox-dependent changes, but none show changes in line width, implying that the flexibility of the oxidized and reduced proteins is not different. However, the spectra demonstrate that folded protein is in equilibrium with a surprisingly large amount of locally unfolded protein, which increases with oxidation for residues localized to the proximal side of the heme. The data suggest that while the oxidized protein is not more flexible than the reduced protein, it is more locally unfolded. Local unfolding of cytochrome c might be one mechanism whereby the protein evolved to control electron transfer.     

SDA 7: A modular and parallel implementation of the simulation of diffusional association software

The simulation of diffusional association (SDA) Brownian dynamics software package has been widely used in the study of biomacromolecular association. Initially developed to calculate bimolecular protein–protein association rate constants, it has since been extended to study electron transfer rates, to predict the structures of biomacromolecular complexes, to investigate the adsorption of proteins to inorganic surfaces, and to simulate the dynamics of large systems containing many biomacromolecular solutes, allowing the study of concentration‐dependent effects. These extensions have led to a number of divergent versions of the software. In this article, we report the development of the latest version of the software (SDA 7). This release was developed to consolidate the existing codes into a single framework, while improving the parallelization of the code to better exploit modern multicore shared memory computer architectures. It is built using a modular object‐oriented programming scheme, to allow for easy maintenance and extension of the software, and includes new features, such as adding flexible solute representations. We discuss a number of application examples, which describe some of the methods available in the release, and provide benchmarking data to demonstrate the parallel performance. © 2015 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.     

H-bond refinement for electron transfer membrane-bound protein–protein complexes: Cytochrome c oxidase and cytochrome c552

In this study we propose a protocol to evaluate membrane-bound cytochrome c oxidase–cytochrome c552 docking candidates. An initial rigid docking algorithm generates docking poses of the cytochrome c oxidase–cytochrome c552, candidates are then aggregated into a 512-DPPC membrane model and solvated in explicit solvent. Molecular dynamic simulations are performed to induce conformational changes to membrane-bound protein complexes. Lastly each protein–protein complex is optimized in terms of its hydrogen bond network, evaluated energetically and ranked. The protocol is directly applicable to other membrane-protein complexes, such as protein–ligand systems.     

Electrostatics in protein binding and function

Protein electrostatic properties stem from the proportion and distribution of polar and charged residues. Polar and charged residues regulate the electrostatic properties by forming short-range interactions, like salt-bridges and hydrogen-bonds, and by defining the over-all electrostatic environment in the protein. Electrostatics play a major role in defining the mechanisms of protein-protein complex formation, molecular recognitions, thermal stabilities, conformational adaptabilities and protein movements. For example:- Functional hinges, or flexible regions of the protein, lack short-range electrostatic interactions; Thermophilic proteins have higher electrostatic interactions than their mesophilic counter parts; Increase in binding specificity and affinity involve optimization of electrostatics; High affinity antibodies have higher, and stronger, electrostatic interactions with their antigens; Rigid parts of proteins have higher and stronger electrostatic interactions. In this review we address the significance of electrostatics in protein folding, binding and function. We discuss that the electrostatic properties are evolutionally selected by a protein to perform an specific function. We also provide bona fide examples to illustrate this. Additionally, using continuum electrostatic and molecular dynamics approaches we show that the "hot-spot" inter-molecular interactions in a very specific antibody-antigen binding are mainly established through charged residues. These "hot-spot" molecular interactions stay intact even during high temperature molecular dynamics simulations, while the other inter-molecular interactions, of lesser functional significance, disappear. This further corroborates the significance of charge-charge interactions in defining binding mechanisms. High affinity binding frequently involves "electrostatic steering". The forces emerge from over-all electrostatic complementarities and by the formation of charged and polar interactions. We demonstrate that although the high affinity binding of barnase-barstar and anti-hen egg white lysozyme (HEL) antibody-HEL complexes involve different molecular mechanisms, it is electrostatically regulated in both the cases. These observations, and several other studies, suggest that a fine tuning of local and global electrostatic properties are essential for protein binding and function.     

Structures of Cytochrome b_5 Mutated at the Charged Surface-Residues and Their Interactions with Cytochrome c

ＩｎｔｒｏｄｕｃｔｉｏｎＭｉｃｒｏｓｏｍａｌｃｙｔｏｃｈｒｏｍｅｂ5(Ｃｙｔｂ5)ｉｓａｍｅｍｂｅｒｏｆｃｙｔｏｃｈｒｏｍｅｂ5ｆａｍｉｌｙ ,ａｎｄｉｔｓｅｒｖｅｓａｓａｎｅｌｅｃｔｒｏｎｃａｒｒｉｅｒｉｎａｓｅｒｉｅｓｏｆｅｌｅｃｔｒｏｎ ｔｒａｎｓｆｅｒｐｒｏｃｅｓｓｅｓｉｎｂｉｏｌｏｇｉｃａｌｓｙｓ ｔｅｍｓ .1 3 Ｃｙｔｂ5ｉｓａｍｅｍｂｒａｎｅｐｒｏｔｅｉｎｗｉｔｈＭｒ～ 16ｋＤａ ,ｃｏｎｓｉｓｔｉｎｇｏｆｔｗｏｄｏｍａｉｎｓ ,ｏｎｅｈｙｄｒｏｐｈｏｂｉｃｄｏｍａｉｎｗｈｉｃｈａｎｃｈｏｒｓｔｈ…     

HADDOCK: a protein-protein docking approach based on biochemical or biophysical information

The structure determination of protein-protein complexes is a rather tedious and lengthy process, by both NMR and X-ray crystallography. Several methods based on docking to study protein complexes have also been well developed over the past few years. Most of these approaches are not driven by experimental data but are based on a combination of energetics and shape complementarity. Here, we present an approach called HADDOCK (High Ambiguity Driven protein-protein Docking) that makes use of biochemical and/or biophysical interaction data such as chemical shift perturbation data resulting from NMR titration experiments or mutagenesis data. This information is introduced as Ambiguous Interaction Restraints (AIRs) to drive the docking process. An AIR is defined as an ambiguous distance between all residues shown to be involved in the interaction. The accuracy of our approach is demonstrated with three molecular complexes. For two of these complexes, for which both the complex and the free protein structures have been solved, NMR titration data were available. Mutagenesis data were used in the last example. In all cases, the best structures generated by HADDOCK, that is, the structures with the lowest intermolecular energies, were the closest to the published structure of the respective complexes (within 2.0 A backbone RMSD).     

细菌化学趋向性受体聚集体的相互作用

细菌化学趋向性受体的最小结构单元为二聚体,在细胞膜上这些二聚体会聚集成大团簇. X射线晶体结构和低分辨电镜结构测定表明,这些团簇有两类不同的形式,一种是在晶体结构中观察到的倒金字塔式二聚体的三聚体重复形成的聚集,另一种为由二聚体尾部相互盘绕形成的拉链状聚集. 有关拉链状聚集的详细分子模型目前尚不清楚. 本文使用蛋白质-蛋白质对接的方法研究了大肠杆菌丝氨酸化学趋向性受体Tsr 二聚体之间的相互作用. 分子对接计算表明,倒金字塔式聚集和拉链状聚集的基本复合物都是可以出现的,相应复合物的分子动力学模拟表明这些结构都具有一定的稳定性. 对于所获得的拉链状聚集体的基本复合物结构模型进行了详细的二聚体作用界面分析,发现二聚体间主要通过静电和疏水作用形成复合物,其中Arg388、Phe373 和Ile377是形成拉链状聚集的关键作用残基. 所建立的Tsr 拉链状聚集的结构模型有助于揭示细菌化学趋向性受体在细胞膜上聚集的分子机制,为进一步的聚集理论及模拟研究提供了基础.     

细菌化学趋向性受体聚集体的相互作用

细菌化学趋向性受体的最小结构单元为二聚体,在细胞膜上这些二聚体会聚集成大团簇.X射线晶体结构和低分辨电镜结构测定表明,这些团簇有两类不同的形式,一种是在晶体结构中观察到的倒金字塔式二聚体的三聚体重复形成的聚集,另一种为由二聚体尾部相互盘绕形成的拉链状聚集.有关拉链状聚集的详细分子模型目前尚不清楚.本文使用蛋白质-蛋白质对接的方法研究了大肠杆菌丝氨酸化学趋向性受体Tsr二聚体之间的相互作用.分子对接计算表明,倒金字塔式聚集和拉链状聚集的基本复合物都是可以出现的,相应复合物的分子动力学模拟表明这些结构都具有一定的稳定性.对于所获得的拉链状聚集体的基本复合物结构模型进行了详细的二聚体作用界面分析,发现二聚体间主要通过静电和疏水作用形成复合物,其中Arg388、Phe373和Ile377是形成拉链状聚集的关键作用残基.所建立的Tsr拉链状聚集的结构模型有助于揭示细菌化学趋向性受体在细胞膜上聚集的分子机制,为进一步的聚集理论及模拟研究提供了基础.     

Recent progress and future directions in protein-protein docking

This article gives an overview of recent progress in protein-protein docking and it identifies several directions for future research. Recent results from the CAPRI blind docking experiments show that docking algorithms are steadily improving in both reliability and accuracy. Current docking algorithms employ a range of efficient search and scoring strategies, including e.g. fast Fourier transform correlations, geometric hashing, and Monte Carlo techniques. These approaches can often produce a relatively small list of up to a few thousand orientations, amongst which a near-native binding mode is often observed. However, despite the use of improved scoring functions which typically include models of desolvation, hydrophobicity, and electrostatics, current algorithms still have difficulty in identifying the correct solution from the list of false positives, or decoys. Nonetheless, significant progress is being made through better use of bioinformatics, biochemical, and biophysical information such as e.g. sequence conservation analysis, protein interaction databases, alanine scanning, and NMR residual dipolar coupling restraints to help identify key binding residues. Promising new approaches to incorporate models of protein flexibility during docking are being developed, including the use of molecular dynamics snapshots, rotameric and off-rotamer searches, internal coordinate mechanics, and principal component analysis based techniques. Some investigators now use explicit solvent models in their docking protocols. Many of these approaches can be computationally intensive, although new silicon chip technologies such as programmable graphics processor units are beginning to offer competitive alternatives to conventional high performance computer systems. As cryo-EM techniques improve apace, docking NMR and X-ray protein structures into low resolution EM density maps is helping to bridge the resolution gap between these complementary techniques. The use of symmetry and fragment assembly constraints are also helping to make possible docking-based predictions of large multimeric protein complexes. In the near future, the closer integration of docking algorithms with protein interface prediction software, structural databases, and sequence analysis techniques should help produce better predictions of protein interaction networks and more accurate structural models of the fundamental molecular interactions within the cell.