Chromophore-protein interactions and the function of the photosynthetic reaction center: a molecular dynamics study.

The coupling between electron transfer and protein structure and dynamics in the photosynthetic reaction center of Rhodopseudomonas viridis is investigated. For this purpose molecular dynamics simulations of the essential portions (a segment of 5797 atoms) of this protein complex have been carried out. Electron transfer in the primary event is modeled by altering the charge distributions of the chromophores according to quantum chemical calculations. The simulations show (i) that fluctuations of the protein matrix, which are coupled electrostatically to electron transfer, play an important role in controlling the electron transfer rates and (ii) that the protein matrix stabilizes the separated electron pair state through rapid (200 fs) and temperature-independent dielectric relaxation. The photosynthetic reaction center resembles a polar liquid in that the internal motions of the whole protein complex, rather than only those of specific side groups, contribute to i and ii. The solvent reorganization energy is about 4.5 kcal/mol. The simulations indicate that rather small structural rearrangements and changes in motional amplitudes accompany the primary electron transfer.     

Photoexcited breathers in conjugated polyenes: an excited-state molecular dynamics study

pi-conjugated polymers have become an important class of materials for electronic devices. Design of these devices requires understanding such processes as photochemical reactions, spatial dynamics of photoexcitations, and energy and charge transport, which in turn involve complex coupled electron-vibrational dynamics. Here we study nonlinear photoexcitation dynamics in the polyene oligomers by using a quantum-chemical method suitable for the simulation of excited-state molecular dynamics in extended molecular systems with sizes up to hundreds of atoms. The method is based on the adiabatic propagation of the ground-state and transition single-electron density matrices along the trajectory. The simulations reveal formation of a self-localized vibronic excitation ("breather" or multiquanta bound state) with a typical period of 34 fs and allows us to identify specific slow and fast nuclear motions strongly coupled to the electronic degrees of freedom. The effect of chain imperfections and chemical defects on the dynamics is also investigated. A complementary two-dimensional analysis of corresponding transition density matrices provides an efficient way to monitor time-dependent real-space localization of the photoexcitation by identifying the underlying changes in charge densities and bond orders. Possible correlated electronic and vibrational spectroscopic signatures of photoexcited breathers are predicted, and generalizations to energy localization in complex macromolecules are discussed.     

Femtosecond direct observation of charge transfer between bases in DNA

Charge transfer in supramolecular assemblies of DNA is unique because of the notion that the pi-stacked bases within the duplex may mediate the transport, possibly leading to damage and/or repair. The phenomenon of transport through pi-stacked arrays over a long distance has an analogy to conduction in molecular electronics, but the mechanism still needs to be determined. To decipher the elementary steps and the mechanism, one has to directly measure the dynamics in real time and in suitably designed, structurally well characterized DNA assemblies. Here, we report our first observation of the femtosecond dynamics of charge transport processes occurring between bases within duplex DNA. By monitoring the population of an initially excited 2-aminopurine, an isomer of adenine, we can follow the charge transfer process and measure its rate. We then study the effect of different bases next to the donor (acceptor), the base sequence, and the distance dependence between the donor and acceptor. We find that the charge injection to a nearest neighbor base is crucial and the time scale is vastly different: 10 ps for guanine and up to 512 ps for inosine. Depending on the base sequence the transfer can be slowed down or inhibited, and the distance dependence is dramatic over the range of 14 A. These observations provide the time scale, and the range and efficiency of the transfer. The results suggest the invalidity of an efficient wire-type behavior and indicate that long-range transport is a slow process of a different mechanism.     

Femtosecond coherent transient infrared spectroscopy of reaction centers from Rhodobacter sphaeroides.

Protein and cofactor vibrational dynamics associated with photoexcitation and charge separation in the photosynthetic reaction center were investigated with femto-second (300-400 fs) time-resolved infrared (1560-1960 cm-1) spectroscopy. The experiments are in the coherent transient limit where the quantum uncertainty principle governs the evolution of the protein vibrational changes. No significant protein relaxation accompanies charge separation, although the electric field resulting from charge separation modifies the polypeptide carbonyl spectra. The potential energy surfaces of the "special pair" P and the photoexcited singlet state P* and environmental perturbations on them are similar as judged from coherence transfer measurements. The vibrational dephasing time of P* modes in this region is 600 fs. A subpicosecond transient at 1665 cm-1 was found to have the kinetics expected for a sequential electron transfer process. Kinetic signatures of all other transient intermediates, P, P*, and P+, participating in the primary steps of photosynthesis were identified in the difference infrared spectra.     

Femtosecond dynamics of dative bonding: concepts of reversible and dissociative electron transfer reactions

With fs time, speed, and angular resolution of the elementary steps in electron transfer reactions, we report direct observation of reversible and dissociative processes for dative bonding involving covalent and ionic characters. For bimolecular reactions of various donors and acceptors we find strong correlation between the structure and the dynamics. The dynamics from the transition state to final products involve two elementary processes, with different reaction times, speed, and angular distributions. For example, for the R2S.I2 (R = C2H5) system, it is shown that after charge separation, the reversible electron transfer occurs in less than 150 fs (fastest trajectory) and is followed by the rupture of the I-I bond with the release of the first I-atom in 510 fs. However, the second process of the remaining and trapped I-atom takes 1.15 ps with its speed (500 m/s) being much smaller than the first one (1,030 m/s). The S-I-I average angle is 130 degrees. These findings, on this and the other systems reported here, elucidate the mechanism and address some concepts of nonconcertedness, caging, and restricted energy redistribution.     

Effect of protein dynamics on biological electron transfer

Computer simulations of the effect of protein dynamics on the long distance tunneling mediated by the protein matrix have been carried out for a Ru-modified (His 126) azurin molecule. We find that the tunneling matrix element is a sensitive function of the atomic configuration of the part of the protein matrix in which tunneling currents (pathways) are localized. Molecular dynamics simulations show that fluctuations of the matrix element can occur on a time scale as short as 10 fs. These short time fluctuations are an indication of a strong dynamic coupling of a tunneling electron to vibrational motions of the protein nuclear coordinates. The latter results in a modification of the conventional Marcus picture of electron transfer in proteins. The new element in the modified theory is that the tunneling electron is capable of emitting or absorbing vibrational energy (phonons) from the medium. As a result, some biological reactions may occur in an activationless fashion. An analytical theoretical model is proposed to account for thermal fluctuations of the medium in long distance electron transfer reactions. The model shows that, at long distances, the phonon-modified inelastic tunneling always dominates over the conventional elastic tunneling.     

Paired moving charges in mitochondrial energy coupling. II. Universality of the principles for energy coupling in biological systems.

The thesis is developed that an acceptable model of biological energy coupling must have universal application. The paired moving charge model of mitochondrial energy coupling is examined from the standpoint of this thesis. Fundamental to this model is the notion that energy coupling involves interaction between paired uncompensated charged species in two vectorially aligned and spatially separated reaction centers. The two charge-separating devices are assumed to be the electron transfer chain (in chloroplast and mitochondria) and intrinsic ionophores (in all transducing organelles and kinases). The universality of the ionophore principle becomes then the crucial test of the validity of the paired moving charge model. The multiple facets of ionophore-mediated couples processes are explored, e.g., coupled hydrolysis of ATP, hormonal control of ion movements, and active transport.     

Femtosecond dynamics of flavoproteins: charge separation and recombination in riboflavine (vitamin B2)-binding protein and in glucose oxidase enzyme

Flavoproteins can function as hydrophobic sites for vitamin B(2) (riboflavin) or, in other structures, with cofactors for catalytic reactions such as glucose oxidation. In this contribution, we report direct observation of charge separation and recombination in two flavoproteins: riboflavin-binding protein and glucose oxidase. With femtosecond resolution, we observed the ultrafast electron transfer from tryptophan(s) to riboflavin in the riboflavin-binding protein, with two reaction times: approximately 100 fs (86% component) and 700 fs (14%). The charge recombination was observed to take place in 8 ps, as probed by the decay of the charge-separated state and the recovery of the ground state. The time scale for charge separation and recombination indicates the local structural tightness for the dynamics to occur that fast and with efficiency of more than 99%. In contrast, in glucose oxidase, electron transfer between flavin-adenine-dinucleotide and tryptophan(s)/tyrosine(s) takes much longer times, 1.8 ps (75%) and 10 ps (25%); the corresponding charge recombination occurs on two time scales, 30 ps and nanoseconds, and the efficiency is still more than 97%. The contrast in time scales for the two structurally different proteins (of the same family) correlates with the distinction in function: hydrophobic recognition of the vitamin in the former requires a tightly bound structure (ultrafast dynamics), and oxidation-reduction reactions in the latter prefer the formation of a charge-separated state that lives long enough for chemistry to occur efficiently. Finally, we also studied the influence on the dynamics of protein conformations at different ionic strengths and denaturant concentrations and observed the sharp collapse of the hydrophobic cleft and, in contrast, the gradual change of glucose oxidase.     

Electrically driven motor in the outer hair cell: Effect of a mechanical constraint

The outer hair cell has a unique voltage-dependent motility associated with charge transfer across the plasma membrane. To examine mechanical changes in the membrane that are coupled with such charge movements, we digested the undercoating of the membrane with trypsin. We inflated the cell into a sphere and constrained the surface area by not allowing volume changes. We found that this constraint on the membrane area sharply reduced motor-associated charge movement across the membrane, demonstrating that charge transfer is directly coupled with membrane area change. This electromechanical coupling in the plasma membrane must be the key element for the motile mechanism of the outer hair cell.     

Electron tunneling in respiratory complex I

NADH:ubiquinone oxidoreductase (complex I) plays a central role in the respiratory electron transport chain by coupling the transfer of electrons from NADH to ubiquinone to the creation of the proton gradient across the membrane necessary for ATP synthesis. Here the atomistic details of electronic wiring of all Fe/S clusters in complex I are revealed by using the tunneling current theory and computer simulations; both density functional theory and semiempirical electronic structure methods were used to examine antiferromagnetically coupled spin states and corresponding tunneling wave functions. Distinct electron tunneling pathways between neighboring Fe/S clusters are identified; the pathways primarily consist of two cysteine ligands and one additional key residue. Internal water between protein subunits is identified as an essential mediator enhancing the overall electron transfer rate by almost three orders of magnitude to achieve a physiologically significant value. The identified key residues are further characterized by sensitivity of electron transfer rates to their mutations, examined in simulations, and their conservation among complex I homologues. The unusual electronic structure properties of Fe(4)S(4) clusters in complex I explain their remarkable efficiency of electron transfer.     

Long-time quantum simulation of the primary charge separation in bacterial photosynthesis.

Accurate quantum mechanical simulations of the primary charge transfer in photosynthetic reaction centers are reported. The process is modeled by three coupled electronic states corresponding to the photoexcited chlorophyll special pair (donor), the reduced bacteriopheophytin (acceptor), and the reduced accessory chlorophyll (bridge) that interact with a dissipative medium of protein and solvent degrees of freedom. The time evolution of the excited special pair is followed over 17 ps by using a fully quantum mechanical path integral scheme. We find that a free energy of the reduced accessory chlorophyll state approximately equal to 400 cm(-1) lower than that of the excited special pair state yields state populations in agreement with experimental results on wild-type and modified reaction centers. For this energetic configuration electron transfer is a two-step process.     

Molecular electronics: some views on transport junctions and beyond

The field of molecular electronics comprises a fundamental set of issues concerning the electronic response of molecules as parts of a mesoscopic structure and a technology-facing area of science. We will overview some important aspects of these subfields. The most advanced ideas in the field involve the use of molecules as individual logic or memory units and are broadly based on using the quantum state space of the molecule. Current work in molecular electronics usually addresses molecular junction transport, where the molecule acts as a barrier for incoming electrons: This is the fundamental Landauer idea of "conduction as scattering" generalized to molecular junction structures. Another point of view in terms of superexchange as a guiding mechanism for coherent electron transfer through the molecular bridge is discussed. Molecules generally exhibit relatively strong vibronic coupling. The last section of this overview focuses on vibronic effects, including inelastic electron tunneling spectroscopy, hysteresis in junction charge transport, and negative differential resistance in molecular transport junctions.     

The generation of oscillations in networks of electrically coupled cells

In several biological systems, the electrical coupling of nonoscillating cells generates synchronized membrane potential oscillations. Because the isolated cell is nonoscillating and electrical coupling tends to equalize the membrane potentials of the coupled cells, the mechanism underlying these oscillations is unclear. Here we present a dynamic mechanism by which the electrical coupling of identical nonoscillating cells can generate synchronous membrane potential oscillations. We demonstrate this mechanism by constructing a biologically feasible model of electrically coupled cells, characterized by an excitable membrane and calcium dynamics. We show that strong electrical coupling in this network generates multiple oscillatory states with different spatio-temporal patterns and discuss their possible role in the cooperative computations performed by the system.     

Evidence for a third sodium-binding site in glutamate transporters suggests an ion/substrate coupling model

Excitatory amino acid transporters (EAATs) remove glutamate from synapses. They maintain an efficient synaptic transmission and prevent glutamate from reaching neurotoxic levels. Glutamate transporters couple the uptake of one glutamate to the cotransport of three sodium ions and one proton and the countertransport of one potassium ion. The molecular mechanism for this coupled uptake of glutamate and its co- and counter-transported ions is not known. In a crystal structure of the bacterial glutamate transporter homolog, Glt_Ph, only two cations are bound to the transporter, and there is no indication of the location of the third sodium site. In experiments using voltage clamp fluorometry and simulations based on molecular dynamics combined with grand canonical Monte Carlo and free energy simulations performed on different isoforms of Glt_Ph as well on a homology model of EAAT3, we sought to locate the third sodium-binding site in EAAT3. Both experiments and computer simulations suggest that T370 and N451 (T314 and N401 in Glt_Ph) form part of the third sodium-binding site. Interestingly, the sodium bound at T370 forms part of the binding site for the amino acid substrate, perhaps explaining both the strict coupling of sodium transport to uptake of glutamate and the ion selectivity of the affinity for the transported amino acid in EAATs.     

Contactless experiments on individual DNA molecules show no evidence for molecular wire behavior

A fundamental requirement for a molecule to be considered a molecular wire (MW) is the ability to transport electrical charge with a reasonably low resistance. We have carried out two experiments that measure first, the charge transfer from an electrode to the molecule, and second, the dielectric response of the MW. The latter experiment requires no contacts to either end of the molecule. From our experiments we conclude that adsorbed individual DNA molecules have a resistivity similar to mica, glass, and silicon oxide substrates. Therefore adsorbed DNA is not a conductor, and it should not be considered as a viable candidate for MW applications. Parallel studies on other nanowires, including single-walled carbon nanotubes, showed conductivity as expected.     

Uncouplers and the molecular mechanism of uncoupling in mitochondria.

Uncouplers are molecules with protonophoric and ionophoric capabilities that mediate coupled cyclical transport of cations--a transport that takes precedence over all other coupled processes. Uncouplers form cation-containing complexes with electrogenic ionophores that potentiate cyclical transport of cations. The molecular mechanism of uncoupling sheds strong light on the mechanism of coupling.     

Electrolytic transport through a synthetic nanometer-diameter pore

We have produced single, synthetic nanometer-diameter pores by using a tightly focused, high-energy electron beam to sputter atoms in 10-nm-thick silicon nitride membranes. Subsequently, we measured the ionic conductance as a function of time, bath concentration, and pore diameter to infer the conductivity and ionic mobility through the pores. The pore conductivity is found to be much larger than the bulk conductivity for dilute bath concentrations, where the Debye length is larger than the pore radius, whereas it is comparable with or less than the bulk for high bath concentrations. We interpret these observations by using multiscale simulations of the ion transport through the pores. Molecular dynamics is used to estimate the ion mobility, and ion transport in the pore is described by the coupled Poisson-Nernst-Planck and the Stokes equations that are solved self-consistently for the ion concentration and velocity and electrical potential. We find that the measurements are consistent with the presence of fixed negative charge in the pore wall and a reduction of the ion mobility because of the fixed charge and the ion proximity to the pore wall.     

Impact of distal mutations on the network of coupled motions correlated to hydride transfer in dihydrofolate reductase

A comprehensive analysis of the network of coupled motions correlated to hydride transfer in dihydrofolate reductase is presented. Hybrid quantum/classical molecular dynamics simulations are combined with a rank correlation analysis method to extract thermally averaged properties that vary along the collective reaction coordinate according to a prescribed target model. Coupled motions correlated to hydride transfer are identified throughout the enzyme. Calculations for wild-type dihydrofolate reductase and a triple mutant, along with the associated single and double mutants, indicate that each enzyme system samples a unique distribution of coupled motions correlated to hydride transfer. These coupled motions provide an explanation for the experimentally measured nonadditivity effects in the hydride transfer rates for these mutants. This analysis illustrates that mutations distal to the active site can introduce nonlocal structural perturbations and significantly impact the catalytic rate by altering the conformational motions of the entire enzyme and the probability of sampling conformations conducive to the catalyzed reaction.