Biomedical subjects

R I Cukier

Publications and source records attributed to R I Cukier.

8 recordsLinked to original sources

PubMed2005

A targeted reweighting method for accelerating the exploration of high-dimensional configuration space.

Time scales available to biomolecular simulations are limited by barriers among states in a high-dimensional configuration space. If equilibrium averages are to be computed, methods that accelerate barrier passage can be carried out by non-Boltzmann sampling. Barriers can be reduced by modifying the potential-energy function and running dynamics on the modified surface. The Boltzmann average can be restored by reweighting each point along the trajectory. We introduce a targeted reweighting scheme where some barriers are reduced, while others are not modified. If only equilibrium properties are desired, trajectories in configuration space can be generated by Langevin dynamics. Once past a transient time, these trajectories guarantee equilibrium sampling when reweighted. A relatively high-order stochastic integration method can be used to generate trajectories. The targeted reweighting scheme is illustrated by a series of double-well models with varying degrees of freedom and shown to be a very efficient method to provide the correct equilibrium distributions, in comparison with analytic results. The scheme is applied to a protein model consisting of a chain of connected beads characterized by dihedral angles and the van der Waals interactions among the beads. We investigate the sampling of configuration space for a model of a helix-turn-helix motif. The targeted reweighting is found to be essential to permit the original all-helical conformation to bend and generate turn structures while still maintaining the alpha-helical segments.

Chemistry, Physical↗

PubMed2005

Structural character and energetics of tyrosyl radical formation by electron/proton transfers of a covalently linked histidine-tyrosine: a model for cytochrome C oxidase.

The structural, energetic, and electronic and IR spectroscopic properties for a model of the cross-linked histidine-tyrosine (His-Tyr) residues as found in cytochrome c oxidase (CcO) are investigated by ab initio methods. The formation of a His-Tyr radical is studied by two paths: proton release followed by electron release and vice versa. The energetics for the proton/electron releases of the Tyr depend modestly on the cross-linked His substituent and, more sensitively, on the charge of the cation attached to the imino N site of the His residue. Protonation of the imino N site significantly increases the electron ionization potential and decreases the proton dissociation energy, making them competitive processes. A positive charge placed at the imino N site, whose value is scanned from zero to one, shows a continuous increase in ionization potential and a decrease in proton dissociation energy, with the +1 limit agreeing well with the protonated imino N site result, indicating a dominant electrostatic effect. The charge populations and the spin density distributions of the His-Tyr model, the radical cation formed by electron ionization, the anion formed by proton dissociation, and the final His-Tyr radical depend sensitively on the substituents, implying a modulation role on the charge transfer between the phenol and imidazole rings, especially for the charged species. His-Tyr and protonated His-Tyr exhibit differences among their respective structural isomers with consequences on their IR absorptions. Small barriers between their pseudo-cis and pseudo-trans rotamers demonstrate the relative flexibility between the two rings, and these may facilitate proton release and charge transfer. The cation effect demonstrates that the cationized cross-linked His-Tyr should be the best candidate to mimic the covalently ring-linked histidine-tyrosine structure in CcO.

Electron Transport Complex IV↗

PubMed2005

A molecular dynamics study of water chain formation in the proton-conducting K channel of cytochrome c oxidase.

The formation of water chains in cytochrome c oxidase (CcO) is studied by molecular dynamics (MD). Focus is on water chains in the K channel that can supply a proton to the binuclear center (the heme a3 Fe/CuB region), the site of O2 reduction. By assessing the presence of chains of any length on a short time scale (0.1 ps), a view of the kinds of chains and their persistence is obtained. Chains from the entry of the channel on the inner membrane to Thr359 (Rhodobacter sphaeroides numbering) are often present but are blocked at that point until a rotation of the Thr359 side chain occurs, permitting formation of chains from Thr359 towards the binuclear center. No continuous hydrogen-bonded water chains are found connecting Thr359 and the binuclear center. Instead, waters hydrogen bond from Thr359 to the hydroxyl of the heme a3 farnesyl and then continue to the binuclear center via Tyr288, which has been identified as a source of a proton for O2 reduction. Three hydrogen-bonded waters are found to be present in the binuclear center after a sufficiently long simulation time. One is ligated to the CuB and could be associated with a water (or hydroxyl) identified in the crystal structure as the fourth ligand of CuB. The water hydrogen-bonded to the hydroxyl of Tyr288 is extremely persistent and well positioned to participate in O2 reduction. The third water is located where O2 is often suggested to reside in mechanistic studies of O2 reduction.

Algorithms↗

PubMed2004

Quantum molecular dynamics simulation of proton transfer in cytochrome c oxidase.

Proton transfer/translocation is studied in cytochrome c oxidase (CcO) by a combination of quantum mechanics (QM) for the transferring protons and classical molecular dynamics (MD) for the protein and solvent. The possibility of a glutamate, Glu286 in the Rhodobacter sphaeroides numbering scheme, acting as a rely point for proton translocation is investigated. The MD finds a hydrogen-bonded cycle of two waters and the carboxylate oxygens of Glu286. The possibility of protonating Glu286 to form neutral GluH is studied and we find that, as experimentally inferred, this glutamate can spend most of its time as GluH. Since translocation relies on the presence of water chains within CcO channels, MD is used to assess their formation. Glu286 and Mg(2+) can be connected by continuous hydrogen-bonded chains that are robust, though transient, and the protein appears spongy above (toward the outer membrane) the Mg(2+). In contrast, the D-channel spanning Asp132, close to the inner membrane surface, to Glu286, forms water chains that are much sparser and do not continuously connect these residues. Rather, there are chains spanning Glu286 to the vicinity of Asn140, and other more robust and ramified water structures that connect Asp132 with waters close to Asn140.

Asparagine↗

PubMed2004

Theory and simulation of proton-coupled electron transfer, hydrogen-atom transfer, and proton translocation in proteins.

A theory of proton coupled electron transfer (PCET) is reviewed with application to charge transfer steps in the photosystem II oxygen-evolving complex (PSII/OEC). The relation between PCET when it is a concerted electron proton transfer (ETPT) process and hydrogen-atom transfer (HAT) reactions is discussed. Signatures expected for HAT reactions in terms of the size of the kinetic isotope effect and overall magnitude of the rate constant are discussed in the context of PSII/OEC. The formal similarity of ETPT to proton transfer and translocation is used to introduce a combined quantum mechanical (for the transferring protons) and molecular dynamics for the heavy-atom degrees of freedom approach. The method is used to examine double proton transfer in cytochrome c oxidase where two waters and a glutamate (Glu286) that is implicated in the proton translocation mechanism form a cyclic hydrogen bonded structure. Protonation of the glutamate is found to occur in agreement with experimental results.

Electron Transport↗

PubMed2000

Effects of Sr2+-substitution on the reduction rates of Yz* in PSII membranes--evidence for concerted hydrogen-atom transfer in oxygen evolution.

Several groups have recently investigated the kinetic effects of biochemical treatments, site-directed mutagenesis, or substitution of essential cofactors on the stepwise, water-oxidizing chemistry catalyzed by Photosystem II. Consistently, these studies show evidence for a slowing of the final, oxygen-releasing step, S(3) --> S(0), of the catalytic cycle. To a degree, some of this work also shows a slowing of the earlier S-state transitions. To study these processes in more detail, we have investigated the effect of replacing Ca(2+) with Sr(2+)on the rates of the S-state transitions by using time-resolved electron paramagnetic resonance. The results show a slowdown of the last transition in the cycle, consistent with a report from Boussac et al. [Boussac, A., Sétif, P., and Rutherford, A. W. (1992) Biochemistry 31, 1224-1234], and of the earlier S-state transitions as well, which suggests that a common molecular mechanism is at work and that Sr(2+) is less effective than Ca(2+) in supporting it. While the oxidation of Y(z) by P(680)(+) has been extensively studied and can be understood within the context of nonadiabatic electron tunneling combined with rapid, non-rate-limiting proton transfer in the holo-system [Tommos, C., and Babcock, G. T. (2000) Biochim. Biophys. Acta 1458, 199], the reduction of Y(z*) by the Mn cluster cannot be described effectively by a nonadiabatic electron-transfer formalism. This indicates that this reaction is rate limited by processes other than electron tunneling. We discuss our results for Y(z*) reduction and those of others for the activation parameters (E(a), A, KIE, and rates) associated with this process, in terms of both sequential and concerted proton-coupled, electron transfer. Our analysis indicates that concerted hydrogen-atom transfer processes best explain the observed characteristics of the S-state advances.

Arabidopsis Proteins↗

PubMed2000

Concerted hydrogen-atom abstraction in photosynthetic water oxidation.

Photosystem II evolves oxygen by using water in the unlikely role of a reductant. The absorption of sunlight by chlorophyll produces highly oxidizing equivalents that are filled with electrons stripped from water. This proton-coupled redox chemistry occurs at the oxygen-evolving complex, which contains a tetramanganese cluster, a redox-active tyrosine amino acid hydrogen-bonded to a histidine amino acid, a calcium ion and chloride. Hydrogen-atom abstraction by the tyrosyl radical from water bound to the manganese cluster is now widely held to occur in this process, at least for some of the steps in the catalytic cycle. We discuss kinetic and energetic constraints on the hydrogen-atom abstraction process.

Hydrogen↗

PubMed1998

Proton-coupled electron transfer.

Proton-coupled electron transfer (PCET) is an important mechanism for charge transfer in a wide variety of systems including biology- and materials-oriented venues. We review several areas where the transfer of an electron and proton is tightly coupled and discuss model systems that can provide an experimental basis for a test of PCET theory. In a PCET reaction, the electron and proton may transfer consecutively (ET/PT) or concertedly (ETPT). The distinction between these processes is formulated, and rate-constant expressions for the two reaction channels are presented. Methods for the evaluation of these rate constants are discussed that are based on dielectric continuum theory. Electron donor hydrogen-bonded-interface electron acceptor systems displaying PCET reactivity are presented, and the rate-constant expressions corresponding to the ETPT and ET/PT channels for several model reaction complexes are evaluated.

Electron Transport↗