Non-Traditional Biochemistry: Disordered Proteins and Educational Pathways

193587-Thumbnail Image.png
Description
Since understanding the nature of proteins, it has been a long held belief that protein sequence dictated structure which then determined function. As such, all proteins contained structure and those that did not must not serve a purpose. For the

Since understanding the nature of proteins, it has been a long held belief that protein sequence dictated structure which then determined function. As such, all proteins contained structure and those that did not must not serve a purpose. For the last 25 years, scientists have begun to understand that disordered proteins, lacking structure, did not lack function. Their unique ability to undergo liquid-liquid phase separation served a cellular purpose, most involving nucleic acids. As more is uncovered, these unique proteins are being used to build new systems. Phase separated disordered proteins were used to design a functional organelle using the enzyme horseradish peroxidase and its chromatic substrate ABTS. Upon doing so, it was discovered that disordered proteins are highly susceptible to chemical modification through radical reactions with tyrosine. The increased frequency of tyrosine in disordered proteins provides multiple sites of conjugation by the ABTS radical and other substrates. These modifications then alter the physical properties of the proteins. The phase separated system was also incorporated with shell proteins from bacterial microcompartments in an attempt to limit access to the droplets. Through expression with truncations of the disordered sequence, shell proteins were able to interact with the droplets. Despite the appearance of complete coatings, they were found to be permeable to their surroundings, though much more stable than uncoated droplets. Just as disordered proteins are considered outside the traditional structures, so too are many students entering higher education. Non-traditional students are becoming more prevalent in the undergraduate population, though they are woefully underrepresented in the natural sciences. The benefits these students bring to their programs is highlighted and the circumstances that drive them away from STEM is explored. Non-traditional students contribute to the diversity of the scientific population, though many pursue education in non-STEM fields. To support these students, focus is put on andragogy (the teaching of adults), rather than pedagogy (the teaching of children). Non-traditional students face isolation and discrimination that is not being addressed by higher education institutions, hindering their ability to succeed. Through infrastructure designed for adult learners, STEM fields can be diversified in non-traditional ways.
Date Created
2024
Agent

Comparative Analysis of Molecular Simulations

Description
The purpose of this project was to compare the different physical models behind four algorithms in computational chemistry: Molecular dynamics with a thermostat (specifically simple velocity rescaling, Berendsen, and Nosé-Hoover), Langevin dynamics, Brownian dynamics, and Monte Carlo. These algorithms were

The purpose of this project was to compare the different physical models behind four algorithms in computational chemistry: Molecular dynamics with a thermostat (specifically simple velocity rescaling, Berendsen, and Nosé-Hoover), Langevin dynamics, Brownian dynamics, and Monte Carlo. These algorithms were programmed in C and the impact of specific parameters, such as the coupling parameter and time step, were studied. Their results were compared based on their radial distribution functions and, when the thermostats were in use, fluctuations in temperature.
Date Created
2022-12
Agent

Single-Focus Confocal Data Analysis with Bayesian Nonparametrics

Description
The cell is a dense environment composes of proteins, nucleic acids, as well as other small molecules, which are constantly bombarding each other and interacting. These interactions and the diffusive motions are driven by internal thermal fluctuations. Upon collision, molecules

The cell is a dense environment composes of proteins, nucleic acids, as well as other small molecules, which are constantly bombarding each other and interacting. These interactions and the diffusive motions are driven by internal thermal fluctuations. Upon collision, molecules can interact and form complexes. It is of interest to learn kinetic parameters such as reaction rates of one molecule converting to different species or two molecules colliding and form a new species as well as to learn diffusion coefficients.

Several experimental measurements can probe diffusion coefficients at the single-molecule and bulk level. The target of this thesis is on single-molecule methods, which can assess diffusion coefficients at the individual molecular level. For instance, super resolution methods like stochastic optical reconstruction microscopy (STORM) and photo activated localization microscopy (PALM), have a high spatial resolution with the cost of lower temporal resolution. Also, there is a different group of methods, such as MINFLUX, multi-detector tracking, which can track a single molecule with high spatio-temporal resolution. The problem with these methods is that they are only applicable to very diluted samples since they need to ensure existence of a single molecule in the region of interest (ROI).

In this thesis, the goal is to have the best of both worlds by achieving high spatio-temporal resolutions without being limited to a few molecules. To do so, one needs to refocus on fluorescence correlation spectroscopy (FCS) as a method that applies to both in vivo and in vitro systems with a high temporal resolution and relies on multiple molecules traversing a confocal volume for an extended period of time. The difficulty here is that the interpretation of the signal leads to different estimates for the kinetic parameters such as diffusion coefficients based on a different number of molecules we consider in the model. It is for this reason that the focus of this thesis is now on using Bayesian nonparametrics (BNPs) as a way to solve this model selection problem and extract kinetic parameters such as diffusion coefficients at the single-molecule level from a few photons, and thus with the highest temporal resolution as possible.
Date Created
2020
Agent

Temperature and polarizability effects on electron transfer in biology and artificial photosynthesis

Description
This study aims to address the deficiencies of the Marcus model of electron transfer

(ET) and then provide modifications to the model. A confirmation of the inverted energy

gap law, which is the cleanest verification so far, is presented for donor-acceptor complexes.

In

This study aims to address the deficiencies of the Marcus model of electron transfer

(ET) and then provide modifications to the model. A confirmation of the inverted energy

gap law, which is the cleanest verification so far, is presented for donor-acceptor complexes.

In addition to the macroscopic properties of the solvent, the physical properties of the solvent

are incorporated in the model via the microscopic solvation model. For the molecules

studied in this dissertation, the rate constant first increases with cooling, in contrast to the

prediction of the Arrhenius law, and then decreases at lower temperatures. Additionally,

the polarizability of solute, which was not considered in the original Marcus theory, is included

by the Q-model of ET. Through accounting for the polarizability of the reactants, the

Q-model offers an important design principle for achieving high performance solar energy

conversion materials. By means of the analytical Q-model of ET, it is shown that including

molecular polarizability of C60 affects the reorganization energy and the activation barrier

of ET reaction.

The theory and Electrochemistry of Ferredoxin and Cytochrome c are also investigated.

By providing a new formulation for reaction reorganization energy, a long-standing disconnect

between the results of atomistic simulations and cyclic voltametery experiments is

resolved. The significant role of polarizability of enzymes in reducing the activation energy

of ET is discussed. The binding/unbinding of waters to the active site of Ferredoxin leads

to non-Gaussian statistics of energy gap and result in a smaller activation energy of ET.

Furthermore, the dielectric constant of water at the interface of neutral and charged

C60 is studied. The dielectric constant is found to be in the range of 10 to 22 which is

remarkably smaller compared to bulk water( 80). Moreover, the interfacial structural

crossover and hydration thermodynamic of charged C60 in water is studied. Increasing the

charge of the C60 molecule result in a dramatic structural transition in the hydration shell,

which lead to increase in the population of dangling O-H bonds at the interface.
Date Created
2019
Agent

Computational approaches to simulation and analysis of large conformational transitions in proteins

156046-Thumbnail Image.png
Description
In a typical living cell, millions to billions of proteins—nanomachines that fluctuate and cycle among many conformational states—convert available free energy into mechanochemical work. A fundamental goal of biophysics is to ascertain how 3D protein structures encode specific functions, such

In a typical living cell, millions to billions of proteins—nanomachines that fluctuate and cycle among many conformational states—convert available free energy into mechanochemical work. A fundamental goal of biophysics is to ascertain how 3D protein structures encode specific functions, such as catalyzing chemical reactions or transporting nutrients into a cell. Protein dynamics span femtosecond timescales (i.e., covalent bond oscillations) to large conformational transition timescales in, and beyond, the millisecond regime (e.g., glucose transport across a phospholipid bilayer). Actual transition events are fast but rare, occurring orders of magnitude faster than typical metastable equilibrium waiting times. Equilibrium molecular dynamics (EqMD) can capture atomistic detail and solute-solvent interactions, but even microseconds of sampling attainable nowadays still falls orders of magnitude short of transition timescales, especially for large systems, rendering observations of such "rare events" difficult or effectively impossible.

Advanced path-sampling methods exploit reduced physical models or biasing to produce plausible transitions while balancing accuracy and efficiency, but quantifying their accuracy relative to other numerical and experimental data has been challenging. Indeed, new horizons in elucidating protein function necessitate that present methodologies be revised to more seamlessly and quantitatively integrate a spectrum of methods, both numerical and experimental. In this dissertation, experimental and computational methods are put into perspective using the enzyme adenylate kinase (AdK) as an illustrative example. We introduce Path Similarity Analysis (PSA)—an integrative computational framework developed to quantify transition path similarity. PSA not only reliably distinguished AdK transitions by the originating method, but also traced pathway differences between two methods back to charge-charge interactions (neglected by the stereochemical model, but not the all-atom force field) in several conserved salt bridges. Cryo-electron microscopy maps of the transporter Bor1p are directly incorporated into EqMD simulations using MD flexible fitting to produce viable structural models and infer a plausible transport mechanism. Conforming to the theme of integration, a short compendium of an exploratory project—developing a hybrid atomistic-continuum method—is presented, including initial results and a novel fluctuating hydrodynamics model and corresponding numerical code.
Date Created
2017
Agent

Quantifying Solvent Kinetics in Molecular Dynamics Simulations of Biomolecules

135992-Thumbnail Image.png
Description
The relation between water and protein physics is a topic of much interest. Molecular dynamics (MD) simulations of biomolecules are a common computational technique to obtain atomistic insight into the physical behavior of biomolecules, including the nature of the interaction

The relation between water and protein physics is a topic of much interest. Molecular dynamics (MD) simulations of biomolecules are a common computational technique to obtain atomistic insight into the physical behavior of biomolecules, including the nature of the interaction between water and the protein. In order to model biomolecules at the highest level of accuracy, an explicit, atomistic representation of the water is typically necessary. The number of water molecules that need to be simulated is normally on the order of thousands. The high dimensional MD dataset is then expanded with considerably more dimensions. We describe here a set of tools which can be used to extract general features of the water behavior, which can then be utilized to build simplified models of the water kinetics which make quantitative predictions, such as the flux rate through a pore.
Date Created
2015-12
Agent

Dielectric Constant of Water in the Interface

128051-Thumbnail Image.png
Description

We define the dielectric constant (susceptibility) that should enter the Maxwell boundary value problem when applied to microscopic dielectric interfaces polarized by external fields. The dielectric constant (susceptibility) of the interface is defined by exact linear-response equations involving correlations of

We define the dielectric constant (susceptibility) that should enter the Maxwell boundary value problem when applied to microscopic dielectric interfaces polarized by external fields. The dielectric constant (susceptibility) of the interface is defined by exact linear-response equations involving correlations of statistically fluctuating interface polarization and the Coulomb interaction energy of external charges with the dielectric. The theory is applied to the interface between water and spherical solutes of altering size studied by molecular dynamics (MD) simulations. The effective dielectric constant of interfacial water is found to be significantly lower than its bulk value, and it also depends on the solute size. For TIP3P water used in MD simulations, the interface dielectric constant changes from 9 to 4 when the solute radius is increased from ∼5 to 18 Å.

Date Created
2016-07-06
Agent

Configurational Entropy of Polar Glass Formers and the Effect of Electric Field on Glass Transition

128062-Thumbnail Image.png
Description

A model of low-temperature polar liquids is constructed that accounts for the configurational heat capacity, entropy, and the effect of a strong electric field on the glass transition. The model is based on the Padé-truncated perturbation expansions of the liquid

A model of low-temperature polar liquids is constructed that accounts for the configurational heat capacity, entropy, and the effect of a strong electric field on the glass transition. The model is based on the Padé-truncated perturbation expansions of the liquid state theory. Depending on parameters, it accommodates an ideal glass transition of vanishing configurational entropy and its avoidance, with a square-root divergent enumeration function at the point of its termination. A composite density-temperature parameter ργ/T, often used to represent combined pressure and temperature data, follows from the model. The theory is in good agreement with the experimental data for excess (over the crystal state) thermodynamics of molecular glass formers. We suggest that the Kauzmann entropy crisis might be a signature of vanishing configurational entropy of a subset of degrees of freedom, multipolar rotations in our model. This scenario has observable consequences: (i) a dynamical crossover of the relaxation time and (ii) the fragility index defined by the ratio of the excess heat capacity and excess entropy at the glass transition. The Kauzmann temperature of vanishing configurational entropy and the corresponding glass transition temperature shift upward when the electric field is applied. The temperature shift scales quadratically with the field strength.

Date Created
2016-07-20
Agent

Electron-Transfer Chain in Respiratory Complex I

128139-Thumbnail Image.png
Description

Complex I is a part of the respiration energy chain converting the redox energy into the cross-membrane proton gradient. The electron-transfer chain of iron-sulfur cofactors within the water-soluble peripheral part of the complex is responsible for the delivery of electrons

Complex I is a part of the respiration energy chain converting the redox energy into the cross-membrane proton gradient. The electron-transfer chain of iron-sulfur cofactors within the water-soluble peripheral part of the complex is responsible for the delivery of electrons to the proton pumping subunit. The protein is porous to water penetration and the hydration level of the cofactors changes when the electron is transferred along the chain. High reaction barriers and trapping of the electrons at the iron-sulfur cofactors are prevented by the combination of intense electrostatic noise produced by the protein-water interface with the high density of quantum states in the iron-sulfur clusters caused by spin interactions between paramagnetic iron atoms. The combination of these factors substantially lowers the activation barrier for electron transfer compared to the prediction of the Marcus theory, bringing the rate to the experimentally established range. The unique role of iron-sulfur clusters as electron-transfer cofactors is in merging protein-water fluctuations with quantum-state multiplicity to allow low activation barriers and robust operation. Water plays a vital role in electron transport energetics by electrowetting the cofactors in the chain upon arrival of the electron. A general property of a protein is to violate the fluctuation-dissipation relation through nonergodic sampling of its landscape. High functional efficiency of redox enzymes is a direct consequence of nonergodicity.

Date Created
2017-07-14
Agent

Polarizability of the Active Site of Cytochrome C Reduces the Activation Barrier for Electron Transfer

128515-Thumbnail Image.png
Description

Enzymes in biology’s energy chains operate with low energy input distributed through multiple electron transfer steps between protein active sites. The general challenge of biological design is how to lower the activation barrier without sacrificing a large negative reaction free

Enzymes in biology’s energy chains operate with low energy input distributed through multiple electron transfer steps between protein active sites. The general challenge of biological design is how to lower the activation barrier without sacrificing a large negative reaction free energy. We show that this goal is achieved through a large polarizability of the active site. It is polarized by allowing a large number of excited states, which are populated quantum mechanically by electrostatic fluctuations of the protein and hydration water shells. This perspective is achieved by extensive mixed quantum mechanical/molecular dynamics simulations of the half reaction of reduction of cytochrome c. The barrier for electron transfer is consistently lowered by increasing the number of excited states included in the Hamiltonian of the active site diagonalized along the classical trajectory. We suggest that molecular polarizability, in addition to much studied electrostatics of permanent charges, is a key parameter to consider in order to understand how enzymes work.

Date Created
2016-06-16
Agent