Computational simulation of the docking of Prochlorothrix hollandica plastocyanin to potosystem I: modeling the electron transfer complex

Fundamental limits on wavelength, efficiency and yield of the charge separation triad

In an attempt to optimize a high yield, high efficiency artificial photosynthetic protein we have discovered unique energy and spatial architecture limits which apply to all light-activated photosynthetic systems. We have generated an analytical solution for the time behavior of the core three cofactor charge separation element in photosynthesis, the photosynthetic cofactor triad, and explored the functional consequences of its makeup including its architecture, the reduction potentials of its components, and the absorption energy of the light absorbing primary-donor cofactor. Our primary findings are two: First, that a high efficiency, high yield triad will have an absorption frequency more than twice the reorganization energy of the first electron transfer, and second, that the relative distance of the acceptor and the donor from the primary-donor plays an important role in determining the yields, with the highest efficiency, highest yield architecture having the light absorbing cofactor closest to the acceptor. Surprisingly, despite the increased complexity found in natural solar energy conversion proteins, we find that the construction of this central triad in natural systems matches these predictions. Our analysis thus not only suggests explanations for some aspects of the makeup of natural photosynthetic systems, it also provides specific design criteria necessary to create high efficiency, high yield artificial protein-based triads.

Modeling the interaction between the N-terminal domain of the tumor suppressor p53 and azurin

It is known that the half life of the tumor suppressor p53 can be increased by the interaction with the bacterial protein azurin, resulting in an enhanced anti-tumoral activity. The understanding of the molecular mechanisms on the basis of this phenomenon can open the way to new anti-cancer strategies. Some experimental works have given evidence of an interaction between p53 and azurin (AZ); however the binding regions of the proteins are still unknown. Recently, fluorescence studies have shown that p53 partakes in the binding with the bacterial protein by its N-terminal (NT) domain. Here we have used a computational method to get insight into this interacting mode. The model that we propose for the best complex between AZ and p53 has been obtained from a rigid-body docking, coupled with a molecular dynamics (MD) simulation, a free energy calculation, and validated by mutagenesis analysis. We have found a high degree of geometric fit between the two proteins that are kept together by several hydrophobic interactions and numerous hydrogen bonds. Interestingly, it has emerged that AZ binds essentially to the helices H(I) and H(III) of the p53 NT domain, which are also interacting regions for the foremost inhibitor of p53, MDM2.

High quality binding modes in docking ligands to proteins

Multiple near-optimal conformations of protein-ligand complexes provide a better chance for accurate representation of biomolecular interactions, compared with a single structure. We present ISE-dock--a docking program which is based on the iterative stochastic elimination (ISE) algorithm. ISE eliminates values that consistently lead to the worst results, thus optimizing the search for docking poses. It constructs large sets of such poses with no additional computational cost compared with single poses. ISE-dock is validated using 81 protein-ligand complexes from the PDB and its performance was compared with those of Glide, GOLD, and AutoDock. ISE-dock has a better chance than the other three to find more than 60% top single poses under RMSD = 2.0 A and more than 80% under RMSD = 3.0 A from experimental. ISE alone produced at least one 3.0 A or better solutions among the top 20 poses in the entire test set. In 98% of the examined molecules, ISE produced solutions that are closer than 2.0 A from experimental. Paired t-tests (PTT) were used throughout to assess the significance of comparisons between the performances of the different programs. ISE-dock provides more than 100-fold docking solutions in a similar time frame as LGA in AutoDock. We demonstrate the usefulness of the large near optimal populations of ligand poses by showing a correlation between the docking results and experiments that support multiple binding modes in p38 MAP kinase (Pargellis et al., Nat Struct Biol 2002;9:268-272] and in Human Transthyretin (Hamilton, Benson, Cell Mol Life Sci 2001;58:1491-1521).

Correlation of electron transfer rate in photosynthetic reaction centers with intraprotein dielectric properties

A number of the electrogenic reactions in photosystem I, photosystem II, and bacterial reaction centers (RC) were comparatively analyzed, and the variation of the dielectric permittivity (epsilon) in the vicinity of electron carriers along the membrane normal was calculated. The value of epsilon was minimal at the core of the complexes and gradually increased towards the periphery. We found that the rate of electron transfer (ET) correlated with the value of the dielectric permittivity: the fastest primary ET reactions occur in the low-polarity core of the complexes within the picosecond time range, whereas slower secondary reactions take place at the high-polarity periphery of the complexes within micro- to millisecond time range. The observed correlation was quantitatively interpreted in the framework of the Marcus theory. We calculated the reorganization energy of ET carriers using their van der Waals volumes and experimentally determined epsilon values. The electronic coupling was calculated by the empirical Moser-Dutton rule for the distance-dependent electron tunneling rate in nonadiabatic ET reactions. We concluded that the local dielectric permittivity inferred from the electrometric measurements could be quantitatively used to estimate the rate constant of ET reactions in membrane proteins with resolved atomic structure with the accuracy of less than one order of magnitude.

Docking and molecular dynamics simulation of the Azurin–Cytochrome c551 electron transfer complex

We coupled protein-protein docking procedure with molecular dynamics (MD) simulation to investigate the electron transfer (ET) complex Azurin-Cytochrome c551 whose transient character makes difficult a direct experimental investigation. The ensemble of complexes generated by the docking algorithm are filtered according to both the distance between the metal ions in the redox centres of the two proteins and to the involvement of suitable residues at the interface. The resulting best complex (BC) is characterized by a distance of 1.59 nm and involves Val23 and Ile59 of Cytochrome c551. The ET properties have been evaluated in the framework of the Pathways model and compared with experimental data. A 60 ns long MD simulation, carried on at full hydration, evidenced that the two protein molecules retain their mutual spatial positions upon forming the complex. An analysis of the ET properties of the complex, monitored at regular time intervals, has revealed that several different ET paths are possible, with the occasional intervening of water molecules. Furthermore, the temporal evolution of the geometric distance between the two redox centres is characterized by very fast fluctuations around an average value of 1.6 nm, with periodic jumps at 2 nm with a frequency of about 70 MHz. Such a behaviour is discussed in connection with a nonlinear dynamics of protein systems and its possible implications in the ET process are explored.

Oligomerization states of Bowman-Birk inhibitor by atomic force microscopy and computational approaches

Several methods have been applied to study protein-protein interaction from structural and thermodynamic point of view. The present study reveals that atomic force microscopy (AFM), molecular modeling, and docking approaches represent alternative methods offering new strategy to investigate structural aspects in oligomerization process of proteinase inhibitors. The topography of the black-eyed pea trypsin/chymotrypsin inhibitor (BTCI) was recorded by AFM and compared with computational rigid-bodies docking approaches. Multimeric states of BTCI identified from AFM analysis showed globular-ellipsoidal shapes. Monomers, dimers, trimers, and hexamers were the most prominent molecular arrays observed in AFM images as evaluated by molecular volume calculations and corroborated by in silico docking and theoretical approaches. We therefore propose that BTCI adopts stable and well-packed self-assembled states in monomer-dimer-trimer-hexamer equilibrium. Although there are no correlation between specificity and packing efficiency among proteinases and proteinase inhibitors, the AFM and docked BTCI analyses suggest that these assemblies may exist in situ to play their potential function in oligomerization process.

Structure of the Intermolecular Complex between Plastocyanin and Cytochrome f from Spinach

In oxygenic photosynthesis, plastocyanin shuttles electrons between the membrane-bound complexes cytochrome b6f and photosystem I. The homologous complex between cytochrome f and plastocyanin, both from spinach, is the object of this study. The solution structure of the reduced spinach plastocyanin was determined using high field NMR spectroscopy, whereas the model structure of oxidized cytochrome f was obtained by homology modeling calculations and molecular dynamics. The model structure of the intermolecular complex was calculated using the program AUTODOCK, taking into account biological information obtained from mutagenesis experiments. The best electron transfer pathway from the heme group of cytochrome f to the copper ion of plastocyanin was calculated using the program HARLEM, obtaining a coupling decay value of 1.8 x 10(-4). Possible mechanisms of interaction and electron transfer between plastocyanin and cytochrome f were discussed considering the possible formation of a supercomplex that associates one cytochrome b6f, one photosystem I, and one plastocyanin.

How does nitrous oxide reductase interact with its electron donors?-A docking study

Electron transfer reactions are crucial for respiration and denitrification. In this article, we analyze the interaction of nitrous oxide reductase with its electron donors cytochrome c550 and pseudoazurin. Our docking protocol comprises generation of candidate complexes followed by a selection step based on the distance of the donor and acceptor groups in each partner protein. Finally, the structures of the candidate complexes were optimized using a force field calculation, together with a second distance filtering step. The prediction power of this protocol was studied using the crystal structure of the cytochrome c2/photosynthetic reaction center of Rhodobacter sphaeroides as a reference. The results suggest that both cytochrome c550 and pseudoazurin bind at the same hydrophobic surface patch residing near the CuA center of nitrous oxide reductase. The central, well-conserved interaction surface of the donors is hydrophobic, but it is surrounded by numerous lysine side-chains, which interact electrostatically with analogously positioned side-chain carboxylates of the acceptor. The prediction output is an ensemble of energetically similar structures that are rotationally related to each other. While such an ensemble may reflect incomplete prediction power of the docking protocol, it may also manifest a biological situation where there are multiple ways of forming a productive electron transfer complex. Analyses of the predicted structures and the conservation pattern of the amino acid residues suggest the existence of specific electron transfer pathways to and from the CuA center of nitrous oxide reductase.

Photoelectric studies of the transmembrane charge transfer reactions in photosystem I pigment-protein complexes

The results of studies of charge transfer in cyanobacterial photosystem I (PS I) using the photoelectric method are reviewed. The electrogenicity in the PS I complex and its interaction with natural donors (plastocyanin, cytochrome c(6)), natural acceptors (ferredoxin, flavodoxin), or artificial acceptors and donors (methyl viologen and other redox dyes) were studied. The operating dielectric constant values in the vicinity of the charge transfer carriers in situ were calculated. The profile of distribution of the dielectric constant along the PS I pigment-protein complex (from plastocyanin or cytochrome c(6) through the chlorophyll dimer P700 to the acceptor complex) was estimated, and possible mechanisms of correlation between the local dielectric constant and electron transfer rate constant were discussed.

Mutagenesis of prochlorothrix plastocyanin reveals additional features in photosystem I interactions

Three surface residues of plastocyanin from Prochlorothrix hollandica have been modified by site-directed mutagenesis. Changes have been made in methionine 33, located in the hydrophobic patch of the copper protein, and in arginine 86 and proline 53, both located in the eastern hydrophilic area. The reactivity toward photosystem I of single mutants M33N, P53A, P53E, R86Q, R86E, and the double mutant M33N/P14L has been studied by laser flash absorption spectroscopy. All the mutations yield increased reactivity of plastocyanin toward photosystem I as compared with wild type plastocyanin, thus indicating that in Prochlorothrix electron donation to photosystem I is not optimized. The most drastic increases in the intracomplex electron transfer rate are obtained with mutants in methionine 33, whereas replacing arginine 86 only modestly affects the plastocyanin-photosystem I equilibrium constant for complex formation. Mutations at position 53 also promote major changes in the association of plastocyanin with photosystem I, yielding a change from a mechanism involving complex formation to a simpler collisional interaction. Molecular dynamics calculations indicate that mutations at position 33 promote changes in the H-bond network around the copper center. The comparative kinetic analysis of the reactivity of Prochlorothrix plastocyanin mutants toward photosystem I from other cyanobacteria reveals that mutations M33N, P53A, and P53E result in enhanced general reactivity.