Docking and electron transfer studies between rubredoxin and rubredoxin:oxygen oxidoreductase

The interaction and electron transfer (ET) between rubredoxin (Rd) and rubredoxin:oxygen oxidoreductase (ROO) from Desulfovibrio gigas is studied by molecular modelling techniques. Experimental kinetic assays using recombinant proteins show that the Rd reoxidation by ROO displays a bell-shaped dependence on ionic strength, suggesting a non-trivial electrostatic dependence of the interaction between these two proteins. Rigid docking studies reveal a prevalence for Rd to interact, in a very specific way, with the surface of the ROO dimer near its FMN cofactors. The optimization of the lowest energy complexes, using molecular dynamics simulation, shows a very tight interaction between the surface of the two proteins, with a high probability for Rd residues (but not the iron centre directly) to be in direct contact with the FMN cofactors of ROO. Both electrostatics and van der Waals interactions contribute to the final energy of the complex. In these complexes, the major contributions for complex formation are polar interactions between acidic residues of Rd and basic residues of ROO, plus substantial non-polar interactions between different groups. Important residues for this process are identified. ET estimates (using the Pathways model), in the optimized lowest energy complexes, suggest that these configurations are efficient for transferring electrons. The experimental bell-shaped dependence of kinetics on ionic strength is analysed in view of the molecular modelling results, and hypotheses for the molecular basis of this phenomenon are discussed.     

Classical molecular dynamics simulation of the photoinduced electron transfer dynamics of plastocyanin.

Classical molecular dynamics simulations are used to investigate the nuclear motions associated with photoinduced electron transfer in plastocyanin. The blue copper protein is modeled using a molecular mechanics potential; potential parameters for the copper-protein interactions are determined using an x-ray crystallographic structure and absorption and resonance Raman spectra. Molecular dynamics simulations yield a variety of information about the ground (oxidized) and optically excited (charge-transfer) states: 1) The probability distribution of the potential difference between the states, which is used to determine the coordinate and energy displacements, places the states well within the Marcus inverted region. 2) The two-time autocorrelation function of the difference potential in the ground state and the average of the difference potential after instantaneous excitation to the excited state are very similar (confirming linear response in this system); their decay indicates that vibrational relaxation occurs in about 1 ps in both states. 3) The spectral densities of various internal coordinates begin to identify the vibrations that affect the optical transition; the spectral density of the difference potential correlation function should also prove useful in quantum simulations of the back electron transfer. 4) Correlation functions of the protein atomic motions with the difference potential show that the nuclear motions are correlated over a distance of more than 20 A, especially along proposed electron transport paths.     

Protein dynamics enhance electronic coupling in electron transfer complexes.

Electron-transferring flavoproteins (ETFs) from human and Paracoccus denitrificans have been analyzed by small angle x-ray scattering, showing that neither molecule exists in a rigid conformation in solution. Both ETFs sample a range of conformations corresponding to a large rotation of domain II with respect to domains I and III. A model of the human ETF.medium chain acyl-CoA dehydrogenase complex, consistent with x-ray scattering data, indicates that optimal electron transfer requires domain II of ETF to rotate by approximately 30 to 50 degrees toward domain I relative to its position in the x-ray structure. Domain motion establishes a new "robust engineering principle" for electron transfer complexes, tolerating multiple configurations of the complex while retaining efficient electron transfer.     

Protein structure and dynamics in nonaqueous solvents: insights from molecular dynamics simulation studies

Protein structure and dynamics in nonaqueous solvents are here investigated using molecular dynamics simulation studies, by considering two model proteins (ubiquitin and cutinase) in hexane, under varying hydration conditions. Ionization of the protein groups is treated assuming "pH memory," i.e., using the ionization states characteristic of aqueous solution. Neutralization of charged groups by counterions is done by considering a counterion for each charged group that cannot be made neutral by establishing a salt bridge with another charged group; this treatment is more physically reasonable for the nonaqueous situation, contrasting with the usual procedures. Our studies show that hydration has a profound effect on protein stability and flexibility in nonaqueous solvents. The structure becomes more nativelike with increasing values of hydration, up to a certain point, when further increases render it unstable and unfolding starts to occur. There is an optimal amount of water, approximately 10% (w/w), where the protein structure and flexibility are closer to the ones found in aqueous solution. This behavior can explain the experimentally known bell-shaped dependence of enzyme catalysis on hydration, and the molecular reasons for it are examined here. Water and counterions play a fundamental and dynamic role on protein stabilization, but they also seem to be important for protein unfolding at high percentages of bound water.     

Investigations of the thermostability of rubredoxin models using molecular dynamics simulations.

The affects of differences in amino acid sequence on the temperature stability of the three-dimensional structure of the small beta-sheet protein, rubredoxin (Rd), was revealed when a set of homology models was subjected to molecular dynamics simulations at relatively high temperatures. Models of Rd from the hyperthermophile, Pyrococcus furiosus (Pf), an organism that grows optimally at 100 degrees C, were compared to three mesophilic Rds of known X-ray crystal structure. Simulations covering the limits of known Rd thermostabilities were carried out at temperatures of 300 K, 343 K, 373 K, and 413 K. They suggest that Rd stability is correlated with structural dynamics. Because the dynamic behavior of three Pf Rd models was consistently different from the dynamic behavior of the three mesophilic Rd structures, detailed analysis of the temperature-dependent dynamic behavior was carried out. The major differences between the models of the protein from the hyperthermophile and the others were: (1) an obvious temperature-dependent transition in the mesophilic structures not seen with the Pf Rd models, (2) consistent AMBER energy for the Pf Rd due to differences in nonbonded interaction terms, (3) less variation in the average conformations for the Pf Rd models with temperature, and (4) the presence of more extensive secondary structure for the Pf Rd models. These unsolvated dynamics simulations support a simple, general hypothesis to explain the hyperthermostability of Pf Rd. Its structure simplifies the conformational space to give a single minimum accessible over an extreme range of temperatures, whereas the mesophilic proteins sample a more complex conformational space with two or more minima over the same temperature range.     

OmpA: Gating and dynamics via molecular dynamics simulations

Outer membrane proteins (OMPs) of Gram-negative bacteria have a variety of functions including passive transport, active transport, catalysis, pathogenesis and signal transduction. Whilst the structures of ∼ 25 OMPs are currently known, there is relatively little known about their dynamics in different environments. The outer membrane protein, OmpA from Escherichia coli has been studied extensively in different environments both experimentally and computationally, and thus provides an ideal test case for the study of the dynamics and environmental interactions of outer membrane proteins. We review molecular dynamics simulations of OmpA and its homologues in a variety of different environments and discuss possible mechanisms of pore gating. The transmembrane domain of E. coli OmpA shows subtle differences in dynamics and interactions between a detergent micelle and a lipid bilayer environment. Simulations of the crystallographic unit cell reveal a micelle-like network of detergent molecules interacting with the protein monomers. Simulation and modelling studies emphasise the role of an electrostatic-switch mechanism in the pore-gating mechanism. Simulation studies have been extended to comparative models of OmpA homologues from Pseudomonas aeruginosa (OprF) and Pasteurella multocida (PmOmpA), the latter model including the periplasmic C-terminal domain.     

Seven-helix bundles: molecular modeling via restrained molecular dynamics.

Simulated annealing via restrained molecular dynamics (SA/MD) has been used to model compact bundles of seven approximately (anti)parallel alpha-helices. Seven such helix bundles occur, e.g., in bacteriorhodopsin, in rhodopsin, and in the channel-forming N-terminal domain of Bacillus thuringiensis delta-endotoxin. Two classes of model are considered: (a) those consisting of seven Ala20 peptide chains; and (b) those containing a single polypeptide chain, made up of seven Ala20 helices linked by GlyN interhelix loops (where N = 5 or 10). Three different starting C alpha templates for SA/MD are used, in which the seven helices are arranged (a) on a left-handed circular template, (b) on a bacteriorhodopsin-like template, or (c) on a zig-zag template. The ensembles of models generated by SA/MD are analyzed in terms of their geometry and energetics, and the most stable structures from each ensemble are examined in greater detail. Structures resembling bacteriorhodopsin and structures resembling delta-endotoxin are both represented among the most stable structures. delta-Endotoxin-like structures arise from both circular and bacteriorhodopsin-like C alpha templates. A third helix-packing mode occurs several times among the stable structures, regardless of the C alpha template and of the presence or absence of interhelix loops. It is characterized by a "4 + 1" core, in which four helices form a distorted left-handed supercoil around a central, buried helix. The remaining two helices pack onto the outside of the core. This packing mode is comparable with that proposed for rhodopsin on the basis of two-dimensional electron crystallographic and sequence analysis studies.     

Protein film voltammetry and co-factor electron transfer dynamics in spinach photosystem II core complex

Direct protein film voltammetry (PFV) was used to investigate the redox properties of the photosystem II (PSII) core complex from spinach. The complex was isolated using an improved protocol not used previously for PFV. The PSII core complex had high oxygen-evolving capacity and was incorporated into thin lipid and polyion films. Three well-defined reversible pairs of reduction and oxidation voltammetry peaks were observed at 4 °C in the dark. Results were similar in both types of films, indicating that the environment of the PSII-bound cofactors was not influenced by film type. Based on comparison with various control samples including Mn-depleted PSII, peaks were assigned to chlorophyll a (Chl a) (E _m = ?0.47 V, all vs. NHE, at pH 6), quinones (?0.12 V), and the manganese (Mn) cluster (E _m = 0.18 V). PFV of purified iron heme protein cytochrome b-559 (Cyt b-559), a component of PSII, gave a partly reversible peak pair at 0.004 V that did not have a potential similar to any peaks observed from the intact PSII core complex. The closest peak in PSII to 0.004 V is the 0.18 V peak that was found to be associated with a two-electron process, and thus is inconsistent with iron heme protein voltammetry. The ?0.47 V peak had a peak potential and peak potential-pH dependence similar to that found for purified Chl a incorporated into DMPC films. The midpoint potentials reported here may differ to various extents from previously reported redox titration data due to the influence of electrode double-layer effects. Heterogeneous electron transfer (hET) rate constants were estimated by theoretical fitting and digital simulations for the ?0.47 and 0.18 V peaks. Data for the Chl a peaks were best fit to a one-electron model, while the peak assigned to the Mn cluster was best fit by a two-electron/one-proton model.     

Redox reactions via vanadium-induced electron transfer

One-electron reduction and oxidation induced by vanadium complexes are demonstrated to be useful in oxidative and reductive transformations of carbonyl compounds. The redox interaction between vanadium complexes and redox-active ligands is achieved with coenzyme PQQ and polyanilines that afford the corresponding redox systems.     

Protein structure to function via dynamics

Protein folding, binding, catalytic activity and molecular recognition all involve molecular movements, with varying extents. The molecular movements are brought upon via flexible regions. Stemming from sequence, a fine tuning of electrostatic and hydrophobic properties of the protein fold determine flexible and rigid regions. Studies show flexible regions usually lack electrostatic interactions, such as salt-bridges and hydrogen-bonds, while the rigid regions often have larger number of such electrostatic interactions. Protein flexible regions are not simply an outcome of looser packing or instability, rather they are evolutionally selected. In this review article we highlight the significance of protein flexibilities in folding, binding and function, and their structural and thermodynamic determinants. Our electrostatic calculations and molecular dynamic simulations on an antibody-antigen complex further illustrate the importance of protein flexibilities in binding and function.     

Protein dynamics under light control

A stochastic view of allostery is providing quantitative estimates of the energy made available through protein photoswitches.     

Chemically gated electron transfer. A means of accelerating and regulating rates of biological electron transfer

Long-range protein electron transfer [ET] reactions may be relatively slow because of long ET distance and low driving force. It is possible to dramatically increase the rate of such nonadiabatic reactions by using an adiabatic chemical reaction to activate the system for rapid ET. Three such examples are discussed; nitrogenase, pyruvate:ferredoxin oxidoreductase, and the methylamine dehydrogenase-amicyanin complex. In each example, the faster activated ET reaction is gated (i.e., rate-limited) by the chemical reaction. However, the reaction rate is still orders of magnitude greater than that of the ungated true ET reaction in the absence of chemical activation. Models are presented to describe the mechanisms of activation in the context of ET theory, and the relevance of such chemically gated ET to the regulation of metabolism is discussed.     

Theoretical studies of proton-coupled electron transfer reactions

A theoretical formulation for proton-coupled electron transfer (PCET) is described. This theory allows the calculation of rates and kinetic isotope effects and provides insight into the underlying fundamental principles of PCET reactions. Applications of this theory to PCET reactions in iron bi-imidazoline complexes, oxoruthenium polypyridyl complexes, osmium-benzoquinone systems, amidinium-carboxylate salt bridges, DNA-acrylamide complexes, and ruthenium polypyridyl-tyrosine systems are summarized. The mechanistic insight gained from theoretical calculations on these model systems is relevant to PCET in more complex biological processes such as photosynthesis and respiration.     

The molecular dynamics of pain control

Pain is necessary for survival, but persistent pain can result in anxiety, depression and a reduction in the quality of life. The discriminative and affective qualities of pain are both thought to be regulated in an activity-dependent fashion. Recent studies have identified cells and molecules that regulate pain sensitivity and the parallel pathways that distribute nociceptive information to limbic or sensory areas of the forebrain. Here, we emphasize the cellular and neurobiological consequences of pain, especially those that are involved in the generation and maintenance of chronic pain. These new insights into pain processing will significantly alter our approach to pain control and the development of new analgesics.     

Protein dynamics in solution and in a crystalline environment: a molecular dynamics study

The effect of a solvent and a crystalline environment on the dynamics of proteins is investigated by the method of computer simulation. Three 25-ps molecular dynamics simulations at 300 K of the bovine pancreatic trypsin inhibitor (BPTI), consisting of 454 heavy atoms, are compared: one of BPTI in vacuo, one of BPTI in a box with 2647 spherical nonpolar solvent atoms, and one of BPTI surrounded by fixed crystal image atoms. Both average and time-dependent molecular properties are examined to determine the effect of the environment on the behavior of the protein. The dynamics of BPTI in solution or in the crystal environment are found to be very similar to that found in the vacuum calculation. The primary difference in the average properties is that the equilibrium structure in the presence of solvent or the crystal field is significantly closer to the X-ray structure than is the vacuum result; concomitantly, the more realistic environment leads to a number density closer to experiment. The presence of solvent has a negligible effect on the overall magnitude of the positional or dihedral angle fluctuations in the interior of the protein; however, there are changes in the decay times of the fluctuations of interior atoms. For surface residues, both the magnitude and the time course of the motions are significantly altered by the solvent. There tends to be an increase in the displacements of long side chains and the flexible parts of the main chain that protrude into the solvent. Further, these motions tend to have a more diffusive character with longer relaxation times than in vacuo. The crystal environment has a specific effect on a number of side chains which are held in relatively fixed positions through hydrogen-bond and electric interactions with the neighboring protein atoms. Most of the effects of the solution environment seem to be sufficiently nonspecific that it may be possible to model them by applying a mean field and stochastic dynamic methods.