Theoretical Studies on Coupled Electron and Proton Transfer in Cytochrome c Oxidase

The complexity of life is based on an effective energy transduction machinery, which has evolved during the last 3.5 billion years. In aerobic life, the utilization of the high oxidizing potential of molecular oxygen powers this machinery. Oxygen is safely reduced by a membrane bound enzyme, cytochrome c oxidase (CcO), to produce an electrochemical proton gradient over the mitochondrial or bacterial membrane. This gradient is used for energy-requiring reactions such as synthesis of ATP by F0F1-ATPase and active transport. In this thesis, the molecular mechanism by which CcO couples the oxygen reduction chemistry to proton-pumping has been studied by theoretical computer simulations. By building both classical and quantum mechanical model systems based on the X-ray structure of CcO from Bos taurus, the dynamics and energetics of the system were studied in different intermediate states of the enzyme. As a result of this work, a mechanism was suggested by which CcO can prevent protons from leaking backwards in proton-pumping. The use and activation of two proton conducting channels were also enlightened together with a mechanism by which CcO sorts the chemical protons from pumped protons. The latter problem is referred to as the gating mechanism of CcO, and has remained a challenge in the bioenergetics field for more than three decades. Furthermore, a new method for deriving charge parameters for classical simulations of complex metalloenzymes was developed.

Proton translocation coupled to electron transfer reactions in terminal oxidases

Terminal oxidases are the final proteins of the respiratory chain in eukaryotes and some bacteria. They catalyze most of the biological oxygen consumption on Earth done by aerobic organisms. During the catalytic reaction terminal oxidases reduce dioxygen to water and use the energy released in this process to maintain the electrochemical proton gradient by functioning as a redox-driven proton pump. This membrane gradient of protons is extremely important for cells as it is used for many cellular processes, such as transportation of substrates and ATP synthesis. Even though the structures of several terminal oxidases are known, they are not sufficient in themselves to explain the molecular mechanism of proton pumping. In this work we have applied a complex approach using a variety of different techniques to address the properties and the mechanism of proton translocation by the terminal oxidases. The combination of direct measurements of pH changes during catalytic turnover, time-resolved potentiometric electrometry and optical spectroscopy, made it possible to obtain valuable information about various aspects of oxidase functioning. We compared oxygen binding properties of terminal oxidases from the distinct heme-copper (CcO) and cytochrome bd families and found that cytochrome bd has a high affinity for oxygen, which is 3 orders of magnitude higher than that of CcO. Interestingly, the difference between CcO and cytochrome bd is not only in higher affinity of the latter to oxygen, but also in the way that each of these enzymes traps oxygen during catalysis. CcO traps oxygen kinetically - the molecule of bound dioxygen is rapidly reduced before it can dissociate. Alternatively, cytochrome bd employs an alternative mechanism of oxygen trapping - part of the redox energy is invested into tight oxygen binding, and the price paid for this is the lack of proton pumping. A single cycle of oxygen reduction to water is characterized by translocation of four protons across the membrane. Our results make it possible to assign the pumping steps to discrete transitions of the catalytic cycle and indicate that during in vivo turnover of the oxidase these four protons are transferred, one at a time, during the P→F, F→OH, Oh→Eh, and Eh→R transitions. At the same time, each individual proton translocation step in the catalytic cycle is not just a single reaction catalyzed by CcO, but rather a complicated sequence of interdependent electron and proton transfers. We assume that each single proton translocation cycle of CcO is assured by internal proton transfer from the conserved Glu-278 to an as yet unidentified pump site above the hemes. Delivery of a proton to the pump site serves as a driving reaction that forces the proton translocation cycle to continue.

Structure, function and intracellular dynamics of alphavirus replication complexes

Intracellular membrane alterations are hallmarks of positive-sense RNA (+RNA) virus replication. Strong evidence indicates that within these exotic compartments, viral replicase proteins engage in RNA genome replication and transcription. To date, fundamental questions such as the origin of altered membranes, mechanisms of membrane deformation and topological distribution and function of viral components, are still waiting for comprehensive answers. This study addressed some of the above mentioned questions for the membrane alterations induced during Semliki Forest virus (SFV) infection of mammalian cells. With the aid of electron and fluorescence microscopy coupled with radioactive labelling and immuno-cytochemistry techniques, our group and others showed that few hours after infection the four non structural proteins (nsP1-4) and newly synthesized RNAs of SFV colocalized in close proximity of small membrane invaginations, designated as spherules . These 50-70 nm structures were mainly detected in the perinuclear area, at the limiting membrane of modified endosomes and lysosomes, named CPV-I (cytopathic vacuoles type I). More rarely, spherules were also found at the plasma membrane (PM). In the first part of this study I present the first three-dimensional reconstruction of the CPV-I and the spherules, obtained by electron tomography after chemical or cryo-fixation. Different approaches for imaging these macromolecular assemblies to obtain better structure preservation and higher resolution are presented as unpublished data. This study provides insights into spherule organization and distribution of viral components. The results of this and other experiments presented in this thesis will challenge currently accepted models for virus replication complex formation and function. In a revisitation of our previous models, the second part of this work provides the first complete description of the biogenesis of the CPV-I. The results demonstrate that these virus-induced vacuoles, where hundreds of spherules accumulate at late stages during infection, represent the final phase of a journey initiated at the PM, which apparently serves as a platform for spherule formation. From the PM spherules were internalized by an endocytic event that required the activity of the class I PI3K, caveolin-1, cellular cholesterol and functional actin-myosin network. The resulting neutral endocytic carrier vesicle delivered the spherules to the membrane of pre-existing acidic endosomes via multiple fusion events. Microtubule based transport supported the vectorial transfer of these intermediates to the pericentriolar area where further fusions generated the CPV-I. A signal for spherule internalization was identified in one of the replicase proteins, nsP3. Infections of cells with viruses harbouring a deletion in a highly phosphorylated region of nsP3 did not result in the formation of CPV-Is. Instead, thousands of spherules remained at the PM throughout the infection cycle. Finally, the role of the replicase protein nsP2 during viral RNA replication and transcription was investigated. Three enzymatic activities, protease, NTPase and RNA-triphosphatase were studied with the aid of temperature sensitive mutants in vitro and, when possible, in vivo. The results highlighted the interplay of the different nsP2 functions during different steps of RNA replication and sub-genomic promoter regulation, and suggest that the protein could have different activities when participating in the replication complex or as a free enzyme.

Hyperexpression of chicken riboflavin carrier protein: Antibodies to the recombinant protein curtail pregnancy in rodents

Chicken riboflavin carrier protein (RCP) is a phosphoglycoprotein present in the egg white and yolk of egg-laying animals and in the sera of laying hens and of estrogenized chicks. The RCP cDNA, encoding a protein of predictedMr27,000, has been cloned into a T7 polymerase-driven vector, and high-level expression was observed on induction with IPTG inEscherichia coli.The protein was largely localized in inclusion bodies when expressed at 37°C but was present in the cytosolic fraction when induced at 22°C. At 37°C, two major bands were detected in whole-cell lysates of the strain expressing the protein. N-terminal sequence analysis indicated that the two proteins represented translated products with and without the pelB leader sequence encoded in the pET20b vector, but both included an additional 10 amino acids generated during cloning procedures. The inclusion body obtained at 37°C, on extraction with detergent, led to preferential solubilization of the protein without the pelB signal sequence. The solubilized recombinant RCP was recognized by polyclonal antisera to native RCP but radioimmunoassay revealed quantitative differences in the epitopes exhibited by the recombinant protein. Thus, sequence-specific monoclonal antibodies to chicken RCP also cross-reacted with the recombinant protein with almost equal efficiency, but antibodies which recognize conformation-dependent epitopes showed relatively reduced cross-reactivity with the recombinant protein. Polyclonal antibodies to recombinant RCP were able to recognize both the native and the denatured RCP. Administration of recombinant RCP antisera to pregnant mice led to embryonic resorption leading to early pregnancy termination. These findings reveal that the recombinant protein will be useful for investigations related to the mechanism of pregnancy termination on immunoneutralization of RCP in mammals, as well as in unraveling folding properties of RCP in terms of its ligand binding and antigenetic determinants exposed at its surface.

Coupled Solution for Forced Recharge in Confined Aquifers

Analytical solutions for forced well recharge currently in use were initially developed for pumping scenarios and applied for recharge cases assuming that radial flow in the recharge well replicates a mirror image of that in to a pumping well. Moreover these solutions were not extended to multiaquifer systems. Well bore numerical solutions were generally not considering the effect of well bore interaction, which has a significant effect in the case of a recharge well. In the present paper, improved analytical solutions are developed for a well fully penetrating either single or multiaquifers in respect.to of well storage, well loss, and interactions between the individual aquifers through well bore. The solution developed for constant and variable rates of injection and well loss is applied to the experimental data of the Hansol well injection project near the city of Ahmedabad in the Gujarat state in India. The paper also discusses the difference in well hydraulics of recharge and recovery wells.

A strategy for scheduling tightly coupled parallel applications on clusters

Although various strategies have been developed for scheduling parallel applications with independent tasks, very little work exists for scheduling tightly coupled parallel applications on cluster environments. In this paper, we compare four different strategies based on performance models of tightly coupled parallel applications for scheduling the applications on clusters. In addition to algorithms based on existing popular optimization techniques, we also propose a new algorithm called Box Elimination that searches the space of performance model parameters to determine the best schedule of machines. By means of real and simulation experiments, we evaluated the algorithms on single cluster and multi-cluster setups. We show that our Box Elimination algorithm generates up to 80% more efficient schedule than other algorithms. We also show that the execution times of the schedules produced by our algorithm are more robust against the performance modeling errors.

Femtosecond carrier dynamics and saturable absorption in graphene suspensions

Nonlinear optical properties and carrier relaxation dynamics in graphene, suspended in three different solvents, are investigated using femtosecond (80 fs pulses) Z-scan and degenerate pump-probe spectroscopy at 790 nm. The results demonstrate saturable absorption property of graphene with a nonlinear absorption coefficient, beta of (similar to 2-9) x 10(-8) cm/W. Two distinct time scales associated with the relaxation of photoexcited carriers, a fast one in the range of 130-330 fs (related to carrier-carrier scattering) followed by it slower one in 3.5-4.9 ps range (associated with carrier-phonon scattering) are observed. (C) 2009 American Institute of Physics.

Thermal anomalies in ternary Ge42-xPbxSe58 glasses near the charge carrier reversal threshold

The carrier type reversal (CTR) from p- to n-type in semiconducting chalcogenide glasses is an important and a long standing problem in glass science. Ge-Se glasses exhibit CTR when the metallic elements Bi and Pb are added. For example, bulk Ge42-xSe58Pbx glasses exhibit CTR around 8-9 at. % of Pb. These glasses have been prepared by melt quenching method. Glass transition temperature (T-g), Specific heat change between the liquid and the glassy states (Delta C-p) at T-g and the nonreversing heat flow (Delta H-nr) measured by modulated differential scanning calorimetry exhibit anomalies at 9 at. % of Pb. These observed anomalies are interpreted on the basis of the nano scale phase separation occurring in these glasses.

Asymptotic analysis for the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell: The beam mode

Using asymptotics, the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell vibrating in the beam mode (viz. circumferential wave order n = 1) are studied. Initially, the uncoupled wavenumbers of the acoustic fluid and the cylindrical shell structure are discussed. Simple closed form expressions for the structural wavenumbers (longitudinal, torsional and bending) are derived using asymptotic methods for low- and high-frequencies. It is found that at low frequencies the cylinder in the beam mode behaves like a Timoshenko beam. Next, the coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation of the structure and the acoustic fluid, with an added fluid-loading term involving a parameter mu due to the coupling. An asymptotic expansion involving mu is substituted in this equation. Analytical expressions are derived for the coupled wavenumbers (as modifications to the uncoupled wavenumbers) separately for low- and high-frequency ranges and further, within each frequency range, for large and small values of mu. Only the flexural wavenumber, the first rigid duct acoustic cut-on wavenumber and the first pressure-release acoustic cut-on wavenumber are considered. The general trend found is that for small mu, the coupled wavenumbers are close to the in vacuo structural wavenumber and the wavenumbers of the rigid-acoustic duct. With increasing mu, the perturbations increase, until the coupled wavenumbers are better identified as perturbations to the pressure-release wavenumbers. The systematic derivation for the separate cases of small and large mu gives more insight into the physics and helps to continuously track the wavenumber solutions as the fluid-loading parameter is varied from small to large values. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an intersection in the coupled case, but a gap is created at that frequency. This method of asymptotics is simple to implement using a symbolic computation package (like Maple). (C) 2008 Elsevier Ltd. All rights reserved.

An asymptotic analysis for the coupled dispersion characteristics of a structural acoustic waveguide

Analytical expressions are derived, using asymptotics, for the fluid-structure coupled wavenumbers in a one-dimensional (1-D) structural acoustic waveguide. The coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation with an added term due to the fluid-structure coupling. As a result of this coupling, the prior uncoupled structural and acoustic wavenumbers, now become coupled structural and acoustic wavenumbers. A fluid-loading parameter e, defined as the ratio of mass of fluid to mass of the structure per unit area, is introduced which when set to zero yields the uncoupled dispersion equation. The coupled wavenumber is then expressed in terms of an asymptotic series in e. Analytical expressions are found as e is varied from small to large values. Different asymptotic expansions are used for different frequency ranges with continuous transitions occurring between them. This systematic derivation helps to continuously track the wavenumber solutions as the fluid-loading parameter is varied from small to large values. Though the asymptotic expansion used is limited to the first-order correction factor, the results are close to the numerical results. A general trend is that a given wavenumber branch transits from a rigid-walled solution to a pressure-release solution with increasing E. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an-intersection in the coupled case, but a gap is created at that frequency. (c) 2007 Elsevier Ltd. All rights reserved.

Asymptotic analysis for the coupled wavenumbers in an infinite fluid-filled flexible cylindrical shell: The axisymmetric mode

The coupled wavenumbers of a fluid-filled flexible cylindrical shell vibrating in the axisymmetric mode are studied. The coupled dispersion equation of the system is rewritten in the form of the uncoupled dispersion equation of the structure and the acoustic fluid, with an added fluid-loading term involving a parameter e due to the coupling. Using the smallness of Poisson's ratio (v), a double-asymptotic expansion involving e and v 2 is substituted in this equation. Analytical expressions are derived for the coupled wavenumbers (for large and small values of E). Different asymptotic expansions are used for different frequency ranges with continuous transitions occurring between them. The wavenumber solutions are continuously tracked as e varies from small to large values. A general trend observed is that a given wavenumber branch transits from a rigidwalled solution to a pressure-release solution with increasing E. Also, it is found that at any frequency where two wavenumbers intersect in the uncoupled analysis, there is no more an intersection in the coupled case, but a gap is created at that frequency. Only the axisymmetric mode is considered. However, the method can be extended to the higher order modes.

Evolutionary Significance of Self-acylation Property in Acyl Carrier Proteins

Acyl carrier protein is an integral component of many cellular metabolic processes. A number of studies have reported self-acylation behavior in acyl carrier proteins. Although AM exhibit high levels of similarity in their primary and tertiary structures, self-acylation behavior is restricted to only some ACPs that can be classified into two major families based on their function. The first family of ACPs is involved in polyketide biosynthesis, whereas the second family participates in fatty acid synthesis. Facilitated by the growing number of genome sequences available for analyses, large-scale phylogenetic studies were used in these studies to uncover as to how self-acylation behavior of acyl carrier proteins is linked with the evolution of metabolic pathways in organisms. These studies show that self-acylation behavior in acyl carrier proteins was lost during the course of evolution, with certain organisms and organelles viz. plastids, retaining it for specified functions. (C) 2009 IUBMB IUBMB Life, 61(8): 853-859, 2009

Developing an Unmanned Aerial Vehicle (UAV) multi-sensor payload carrier for air quality monitoring

Emissions of gases and particles from sea-faring ships have been shown to impact on the atmospheric chemistry and climate. To efficiently monitor and report these emissions found from a ship’s plume, the concept of using a multi-rotor or UAV to hover inside or near the exhaust of the ship to actively record the data in real time is being developed. However, for the required sensors obtain the data; their sensors must face into the airflow of the ships plume. This report presents an approach to have sensors able to read in the chemicals and particles emitted from the ship without affecting the flight dynamics of the multi-rotor UAV by building a sealed chamber in which a pump can take in the surrounding air (outside the downwash effect of the multi-rotor) where the sensors are placed and can analyse the gases safely. Results show that the system is small, lightweight and air-sealed and ready for flight test.

Structural Insights into the acyl-intermediates of plasmodium falciparum fatty acid synthesis pathway; the mechanism of expansion of acyl carrier protein core

Acyl carrier protein (ACP) plays a central role in fatty acid biosynthesis. However, the molecular machinery that mediates its function is not yet fully understood. Therefore, structural studies were carried out on the acyl-ACP intermediates of Plasmodium falciparum using NMR as a spectroscopic probe. Chemical shift perturbation studies put forth a new picture of the interaction of ACP molecule with the acyl chain, namely, the hydrophobic core can protect up to 12 carbon units, and additional carbons protrude out from the top of the hydrophobic cavity. The latter hypothesis stems from chemical shift changes observed in C-alpha and C-beta of Ser-37 in tetradecanoyl-ACP. C-13, N-15-Double-filtered nuclear Overhauser effect (NOE) spectroscopy experiments further substantiate the concept; in octanoyl (C-8)- and dodecanoyl (C-12)-ACP, a long range NOE is observed within the phosphopantetheine arm, suggesting an arch-like conformation. This NOE is nearly invisible in tetradecanoyl (C-14)-ACP, indicating a change in conformation of the prosthetic group. Furthermore, the present study provides insights into the molecular mechanism of ACP expansion, as revealed from a unique side chain-to-backbone hydrogen bond between two fairly conserved residues, Ile-55 HN and Glu-48 O. The backbone amide of Ile-55 HN reports a pK(a) value for the carboxylate, similar to 1.9 pH units higher than model compound value, suggesting strong electrostatic repulsion between helix II and helix III. Charge-charge repulsion between the helices in combination with thrust from inside due to acyl chain would energetically favor the separation of the two helices. Helix III has fewer structural restraints and, hence, undergoes major conformational change without altering the overall-fold of P. falciparum ACP.

Critical lattice size limit for synchronized chaotic state in one- and two-dimensional diffusively coupled map lattices

We consider diffusively coupled map lattices with P neighbors (where P is arbitrary) and study the stability of the synchronized state. We show that there exists a critical lattice size beyond which the synchronized state is unstable. This generalizes earlier results for nearest neighbor coupling. We confirm the analytical results by performing numerical simulations on coupled map lattices with logistic map at each node. The above analysis is also extended to two-dimensional P-neighbor diffusively coupled map lattices.