Skip to main content
eScholarship
Open Access Publications from the University of California

Scholarly Works (5 results)

  • Thesis
  • Peer Reviewed
UC San Diego Electronic Theses and Dissertations (2015)

As human understandings of physics, chemistry and biology converge and the development of computers proceeds, computational chemistry or computational biophysics has become a substantial field of research. It serves to explore the fundamentals of life and also has extended applications in the field of medicine. Among the many aspects of computational chemistry, this Ph. D. work focuses on the numerical methods for studying diffusion and electrostatics of biomolecules at the nanoscale. Diffusion and electrostatics are two independent subjects in terms of their physics, but closely related in applications. In living cells, the mechanism of diffusion powers a ligand to move towards its binding target. And electrostatic forces between the ligand and the target or the ligand and the environment guide the direction of the diffusion, the correct binding orientation and, together with other molecular forces, ensure the stability of the bound complex. More abstractly, diffusion describes the stochastic manner biomolecules move on their energy landscape and electrostatic forces are a major contributor to the shape of the energy landscape. This Ph. D. work aims to acquire a good understanding of both biomolecular diffusion and electrostatics and how the two are used together in numerical calculations. Three projects are presented. The first project is a proof of concept of the bead-model approach to calculate the diffusion tensor. The second project is the benchmark for a new electrostatics method, the size-modified Poisson-Boltzmann equation. The third project is an application that combines diffusion and electrostatics to calculate the substrate channeling efficiency between the human thymidylate synthase and dihydrofolate reductase.

    Cover page: Computational Studies on Biomolecular Diffusion and Electrostatics
    • Article
    • Peer Reviewed
    UC San Diego Previously Published Works (2016)
    In vivo, as an advanced catalytic strategy, transient non-covalently bound multi-enzyme complexes can be formed to facilitate the relay of substrates, i. e. substrate channeling, between sequential enzymatic reactions and to enhance the throughput of multi-step enzymatic pathways. The human thymidylate synthase and dihydrofolate reductase catalyze two consecutive reactions in the folate metabolism pathway, and experiments have shown that they are very likely to bind in the same multi-enzyme complex in vivo. While reports on the protozoa thymidylate synthase-dihydrofolate reductase bifunctional enzyme give substantial evidences of substrate channeling along a surface "electrostatic highway," attention has not been paid to whether the human thymidylate synthase and dihydrofolate reductase, if they are in contact with each other in the multi-enzyme complex, are capable of substrate channeling employing surface electrostatics. This work utilizes protein-protein docking, electrostatics calculations, and Brownian dynamics to explore the existence and mechanism of the substrate channeling between the human thymidylate synthase and dihydrofolate reductase. The results show that the bound human thymidylate synthase and dihydrofolate reductase are capable of substrate channeling and the formation of the surface "electrostatic highway." The substrate channeling efficiency between the two can be reasonably high and comparable to that of the protozoa.
      • Article
      • Peer Reviewed
      UC San Diego Previously Published Works (2018)
      Malate dehydrogenase (MDH) and citrate synthase (CS) are two pacemaking enzymes involved in the tricarboxylic acid (TCA) cycle. Oxaloacetate (OAA) molecules are the intermediate substrates that are transferred from the MDH to CS to carry out sequential catalysis. It is known that, to achieve a high flux of intermediate transport and reduce the probability of substrate leaking, a MDH-CS metabolon forms to enhance the OAA substrate channeling. In this study, we aim to understand the OAA channeling within possible MDH-CS metabolons that have different structural orientations in their complexes. Three MDH-CS metabolons from native bovine, wild-type porcine, and recombinant sources, published in recent work, were selected to calculate OAA transfer efficiency by Brownian dynamics (BD) simulations and to study, through electrostatic potential calculations, a possible role of charges that drive the substrate channeling. Our results show that an electrostatic channel is formed in the metabolons of native bovine and recombinant porcine enzymes, which guides the oppositely charged OAA molecules passing through the channel and enhances the transfer efficiency. However, the channeling probability in a suggested wild-type porcine metabolon conformation is reduced due to an extended diffusion length between the MDH and CS active sites, implying that the corresponding arrangements of MDH and CS result in the decrease of electrostatic steering between substrates and protein surface and then reduce the substrate transfer efficiency from one active site to another.
        • Article
        • Peer Reviewed
        UC San Diego Previously Published Works (2014)
        Mean-field methods, such as the Poisson-Boltzmann equation (PBE), are often used to calculate the electrostatic properties of molecular systems. In the past two decades, an enhancement of the PBE, the size-modified Poisson-Boltzmann equation (SMPBE), has been reported. Here, the PBE and the SMPBE are reevaluated for realistic molecular systems, namely, lipid bilayers, under eight different sets of input parameters. The SMPBE appears to reproduce the molecular dynamics simulation results better than the PBE only under specific parameter sets, but in general, it performs no better than the Stern layer correction of the PBE. These results emphasize the need for careful discussions of the accuracy of mean-field calculations on realistic systems with respect to the choice of parameters and call for reconsideration of the cost-efficiency and the significance of the current SMPBE formulation.
          Cover page: Poisson–Boltzmann versus Size-Modified Poisson–Boltzmann Electrostatics Applied to Lipid Bilayers
          • Article
          • Peer Reviewed
          UC San Diego Previously Published Works (2018)
          Increases in CO2 concentration in plant leaves due to respiration in the dark and the continuing atmospheric [CO2] rise cause closing of stomatal pores, thus affecting plant-water relations globally. However, the underlying CO2/bicarbonate (CO2/HCO3-) sensing mechanisms remain unknown. [CO2] elevation in leaves triggers stomatal closure by anion efflux mediated via the SLAC1 anion channel localized in the plasma membrane of guard cells. Previous reconstitution analysis has suggested that intracellular bicarbonate ions might directly up-regulate SLAC1 channel activity. However, whether such a CO2/HCO3- regulation of SLAC1 is relevant for CO2 control of stomatal movements in planta remains unknown. Here, we computationally probe for candidate bicarbonate-interacting sites within the SLAC1 anion channel via long-timescale Gaussian accelerated molecular dynamics (GaMD) simulations. Mutations of two putative bicarbonate-interacting residues, R256 and R321, impaired the enhancement of the SLAC1 anion channel activity by CO2/HCO3- in Xenopus oocytes. Mutations of the neighboring charged amino acid K255 and residue R432 and the predicted gate residue F450 did not affect HCO3- regulation of SLAC1. Notably, gas-exchange experiments with slac1-transformed plants expressing mutated SLAC1 proteins revealed that the SLAC1 residue R256 is required for CO2 regulation of stomatal movements in planta, but not for abscisic acid (ABA)-induced stomatal closing. Patch clamp analyses of guard cells show that activation of S-type anion channels by CO2/HCO3-, but not by ABA, was impaired, indicating the relevance of R256 for CO2 signal transduction. Together, these analyses suggest that the SLAC1 anion channel is one of the physiologically relevant CO2/HCO3- sensors in guard cells.
            Top