Electron and proton transfer in NADH:ubiquinone oxidoreductase (Complex I) from Escherichia coli

Cells of every living organism on our planet − bacterium, plant or animal − are organized in such a way that despite differences in structure and function they utilize the same metabolic energy represented by electrochemical proton gradient across a membrane. This gradient of protons is generated by the series of membrane bound multisubunit proteins, Complex I, II, III and IV, organized in so-called respiratory or electron transport chain. In the eukaryotic cell it locates in the inner mitochondrial membrane while in the bacterial cell it locates in the cytoplasmic membrane. The function of the respiratory chain is to accept electrons from NADH and ubiquinol and transfer them to oxygen resulting in the formation of water. The free energy released upon these redox reactions is converted by respiratory enzymes into an electrochemical proton gradient, which is used for synthesis of ATP as well as for many other energy dependent processes. This thesis is focused on studies of the first member of the respiratory chain − NADH:ubiquinone oxidoreductase or Complex I. This enzyme has a boot-shape structure with hydrophilic and hydrophobic domains, the former of which has all redox groups of the protein, the flavin and eight to nine iron-sulfur clusters. Complex I serves as a proton pump coupling transfer of two electrons from NADH to ubiquinone to the translocation of four protons across the membrane. So far the mechanism of energy transduction by Complex I is unknown. In the present study we applied a set of different methods to study the electron and proton transfer reactions in Complex I from Escherichia coli. The main achievement was the experiment that showed that the electron transfer through the hydrophilic domain of Complex I is unlikely to be coupled to proton transfer directly or to conformational changes in the protein. In this work for the first time properties of all redox centers of Complex I were characterized in the intact purified bacterial enzyme. We also probed the role of several conserved amino acid residues in the electron transfer of Complex I. Finally, we found that highly conserved amino acid residues in several membrane subunits form a common pattern with a very prominent feature – the presence of a few lysines within the membrane. Based on the experimental data, we suggested a tentative principle which may govern the redox-coupled proton pumping in Complex I.

Monte Carlo Simulations of Molecular Clusters in Nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.

Monte Carlo simulations of molecular clusters in nucleation

A better understanding of the limiting step in a first order phase transition, the nucleation process, is of major importance to a variety of scientific fields ranging from atmospheric sciences to nanotechnology and even to cosmology. This is due to the fact that in most phase transitions the new phase is separated from the mother phase by a free energy barrier. This barrier is crossed in a process called nucleation. Nowadays it is considered that a significant fraction of all atmospheric particles is produced by vapor-to liquid nucleation. In atmospheric sciences, as well as in other scientific fields, the theoretical treatment of nucleation is mostly based on a theory known as the Classical Nucleation Theory. However, the Classical Nucleation Theory is known to have only a limited success in predicting the rate at which vapor-to-liquid nucleation takes place at given conditions. This thesis studies the unary homogeneous vapor-to-liquid nucleation from a statistical mechanics viewpoint. We apply Monte Carlo simulations of molecular clusters to calculate the free energy barrier separating the vapor and liquid phases and compare our results against the laboratory measurements and Classical Nucleation Theory predictions. According to our results, the work of adding a monomer to a cluster in equilibrium vapour is accurately described by the liquid drop model applied by the Classical Nucleation Theory, once the clusters are larger than some threshold size. The threshold cluster sizes contain only a few or some tens of molecules depending on the interaction potential and temperature. However, the error made in modeling the smallest of clusters as liquid drops results in an erroneous absolute value for the cluster work of formation throughout the size range, as predicted by the McGraw-Laaksonen scaling law. By calculating correction factors to Classical Nucleation Theory predictions for the nucleation barriers of argon and water, we show that the corrected predictions produce nucleation rates that are in good comparison with experiments. For the smallest clusters, the deviation between the simulation results and the liquid drop values are accurately modelled by the low order virial coefficients at modest temperatures and vapour densities, or in other words, in the validity range of the non-interacting cluster theory by Frenkel, Band and Bilj. Our results do not indicate a need for a size dependent replacement free energy correction. The results also indicate that Classical Nucleation Theory predicts the size of the critical cluster correctly. We also presents a new method for the calculation of the equilibrium vapour density, surface tension size dependence and planar surface tension directly from cluster simulations. We also show how the size dependence of the cluster surface tension in equimolar surface is a function of virial coefficients, a result confirmed by our cluster simulations.