The Pauli-Lubanski Vector in a Group-Theoretical Approach to Relativistic Wave Equations

154866-Thumbnail Image.png
Description
Chapter 1 introduces some key elements of important topics such as; quantum mechanics,

representation theory of the Lorentz and Poincare groups, and a review of some basic rela- ´

tivistic wave equations that will play an important role in the work to

Chapter 1 introduces some key elements of important topics such as; quantum mechanics,

representation theory of the Lorentz and Poincare groups, and a review of some basic rela- ´

tivistic wave equations that will play an important role in the work to follow. In Chapter 2,

a complex covariant form of the classical Maxwell’s equations in a moving medium or at

rest is introduced. In addition, a compact, Lorentz invariant, form of the energy-momentum

tensor is derived. In chapter 3, the concept of photon helicity is critically analyzed and its

connection with the Pauli-Lubanski vector from the viewpoint of the complex electromag- ´

netic field, E+ iH. To this end, a complex covariant form of Maxwell’s equations is used.

Chapter 4 analyzes basic relativistic wave equations for the classical fields, such as Dirac’s

equation, Weyl’s two-component equation for massless neutrinos and the Proca, Maxwell

and Fierz-Pauli equations, from the viewpoint of the Pauli-Lubanski vector and the Casimir ´

operators of the Poincare group. A connection between the spin of a particle/field and ´

consistency of the corresponding overdetermined system is emphasized in the massless

case. Chapter 5 focuses on the so-called generalized quantum harmonic oscillator, which

is a Schrodinger equation with a time-varying quadratic Hamiltonian operator. The time ¨

evolution of exact wave functions of the generalized harmonic oscillators is determined

in terms of the solutions of certain Ermakov and Riccati-type systems. In addition, it is

shown that the classical Arnold transform is naturally connected with Ehrenfest’s theorem

for generalized harmonic oscillators. In Chapter 6, as an example of the usefulness of the

methods introduced in Chapter 5 a model for the quantization of an electromagnetic field

in a variable media is analyzed. The concept of quantization of an electromagnetic field

in factorizable media is discussed via the Caldirola-Kanai Hamiltonian. A single mode

of radiation for this model is used to find time-dependent photon amplitudes in relation

to Fock states. A multi-parameter family of the squeezed states, photon statistics, and the

uncertainty relation, are explicitly given in terms of the Ermakov-type system.
Date Created
2016
Agent

Non-Gaussian Lineshapes and Dynamics of Time-Resolved Linear and Nonlinear (Correlation) Spectra

129578-Thumbnail Image.png
Description

Signatures of nonlinear and non-Gaussian dynamics in time-resolved linear and nonlinear (correlation) 2D spectra are analyzed in a model considering a linear plus quadratic dependence of the spectroscopic transition frequency on a Gaussian nuclear coordinate of the thermal bath (quadratic

Signatures of nonlinear and non-Gaussian dynamics in time-resolved linear and nonlinear (correlation) 2D spectra are analyzed in a model considering a linear plus quadratic dependence of the spectroscopic transition frequency on a Gaussian nuclear coordinate of the thermal bath (quadratic coupling). This new model is contrasted to the commonly assumed linear dependence of the transition frequency on the medium nuclear coordinates (linear coupling). The linear coupling model predicts equality between the Stokes shift and equilibrium correlation functions of the transition frequency and time-independent spectral width. Both predictions are often violated, and we are asking here the question of whether a nonlinear solvent response and/or non-Gaussian dynamics are required to explain these observations. We find that correlation functions of spectroscopic observables calculated in the quadratic coupling model depend on the chromophore’s electronic state and the spectral width gains time dependence, all in violation of the predictions of the linear coupling models. Lineshape functions of 2D spectra are derived assuming Ornstein–Uhlenbeck dynamics of the bath nuclear modes. The model predicts asymmetry of 2D correlation plots and bending of the center line. The latter is often used to extract two-point correlation functions from 2D spectra. The dynamics of the transition frequency are non-Gaussian. However, the effect of non-Gaussian dynamics is limited to the third-order (skewness) time correlation function, without affecting the time correlation functions of higher order. The theory is tested against molecular dynamics simulations of a model polar–polarizable chromophore dissolved in a force field water.

Date Created
2014-07-17
Agent

Protein Electron Transfer: Dynamics and Statistics

129615-Thumbnail Image.png
Description

Electron transfer between redox proteins participating in energy chains of biology is required to proceed with high energetic efficiency, minimizing losses of redox energy to heat. Within the standard models of electron transfer, this requirement, combined with the need for

Electron transfer between redox proteins participating in energy chains of biology is required to proceed with high energetic efficiency, minimizing losses of redox energy to heat. Within the standard models of electron transfer, this requirement, combined with the need for unidirectional (preferably activationless) transitions, is translated into the need to minimize the reorganization energy of electron transfer. This design program is, however, unrealistic for proteins whose active sites are typically positioned close to the polar and flexible protein-water interface to allow inter-protein electron tunneling. The high flexibility of the interfacial region makes both the hydration water and the surface protein layer act as highly polar solvents. The reorganization energy, as measured by fluctuations, is not minimized, but rather maximized in this region. Natural systems in fact utilize the broad breadth of interfacial electrostatic fluctuations, but in the ways not anticipated by the standard models based on equilibrium thermodynamics.

The combination of the broad spectrum of static fluctuations with their dispersive dynamics offers the mechanism of dynamical freezing (ergodicity breaking) of subsets of nuclear modes on the time of reaction/residence of the electron at a redox cofactor. The separation of time-scales of nuclear modes coupled to electron transfer allows dynamical freezing. In particular, the separation between the relaxation time of electro-elastic fluctuations of the interface and the time of conformational transitions of the protein caused by changing redox state results in dynamical freezing of the latter for sufficiently fast electron transfer. The observable consequence of this dynamical freezing is significantly different reorganization energies describing the curvature at the bottom of electron-transfer free energy surfaces (large) and the distance between their minima (Stokes shift, small). The ratio of the two reorganization energies establishes the parameter by which the energetic efficiency of protein electron transfer is increased relative to the standard expectations, thus minimizing losses of energy to heat. Energetically efficient electron transfer occurs in a chain of conformationally quenched cofactors and is characterized by flattened free energy surfaces, reminiscent of the flat and rugged landscape at the stability basin of a folded protein.

Date Created
2013
Agent

Network models for materials and biological systems

150164-Thumbnail Image.png
Description
The properties of materials depend heavily on the spatial distribution and connectivity of their constituent parts. This applies equally to materials such as diamond and glasses as it does to biomolecules that are the product of billions of years of

The properties of materials depend heavily on the spatial distribution and connectivity of their constituent parts. This applies equally to materials such as diamond and glasses as it does to biomolecules that are the product of billions of years of evolution. In science, insight is often gained through simple models with characteristics that are the result of the few features that have purposely been retained. Common to all research within in this thesis is the use of network-based models to describe the properties of materials. This work begins with the description of a technique for decoupling boundary effects from intrinsic properties of nanomaterials that maps the atomic distribution of nanomaterials of diverse shape and size but common atomic geometry onto a universal curve. This is followed by an investigation of correlated density fluctuations in the large length scale limit in amorphous materials through the analysis of large continuous random network models. The difficulty of estimating this limit from finite models is overcome by the development of a technique that uses the variance in the number of atoms in finite subregions to perform the extrapolation to large length scales. The technique is applied to models of amorphous silicon and vitreous silica and compared with results from recent experiments. The latter part this work applies network-based models to biological systems. The first application models force-induced protein unfolding as crack propagation on a constraint network consisting of interactions such as hydrogen bonds that cross-link and stabilize a folded polypeptide chain. Unfolding pathways generated by the model are compared with molecular dynamics simulation and experiment for a diverse set of proteins, demonstrating that the model is able to capture not only native state behavior but also partially unfolded intermediates far from the native state. This study concludes with the extension of the latter model in the development of an efficient algorithm for predicting protein structure through the flexible fitting of atomic models to low-resolution cryo-electron microscopy data. By optimizing the fit to synthetic data through directed sampling and context-dependent constraint removal, predictions are made with accuracies within the expected variability of the native state.
Date Created
2011
Agent

Applications of adaptive umbrella sampling in biomolecular simulation

150118-Thumbnail Image.png
Description
Conformational changes in biomolecules often take place on longer timescales than are easily accessible with unbiased molecular dynamics simulations, necessitating the use of enhanced sampling techniques, such as adaptive umbrella sampling. In this technique, the conformational free energy is calculated

Conformational changes in biomolecules often take place on longer timescales than are easily accessible with unbiased molecular dynamics simulations, necessitating the use of enhanced sampling techniques, such as adaptive umbrella sampling. In this technique, the conformational free energy is calculated in terms of a designated set of reaction coordinates. At the same time, estimates of this free energy are subtracted from the potential energy in order to remove free energy barriers and cause conformational changes to take place more rapidly. This dissertation presents applications of adaptive umbrella sampling to a variety of biomolecular systems. The first study investigated the effects of glycosylation in GalNAc2-MM1, an analog of glycosylated macrophage activating factor. It was found that glycosylation destabilizes the protein by increasing the solvent exposure of hydrophobic residues. The second study examined the role of bound calcium ions in promoting the isomerization of a cis peptide bond in the collagen-binding domain of Clostridium histolyticum collagenase. This study determined that the bound calcium ions reduced the barrier to the isomerization of this peptide bond as well as stabilizing the cis conformation thermodynamically, and identified some of the reasons for this. The third study represents the application of GAMUS (Gaussian mixture adaptive umbrella sampling) to on the conformational dynamics of the fluorescent dye Cy3 attached to the 5' end of DNA, and made predictions concerning the affinity of Cy3 for different base pairs, which were subsequently verified experimentally. Finally, the adaptive umbrella sampling method is extended to make use of the roll angle between adjacent base pairs as a reaction coordinate in order to examine the bending both of free DNA and of DNA bound to the archaeal protein Sac7d. It is found that when DNA bends significantly, cations from the surrounding solution congregate on the concave side, which increases the flexibility of the DNA by screening the repulsion between phosphate backbones. The flexibility of DNA on short length scales is compared to the worm-like chain model, and the contribution of cooperativity in DNA bending to protein-DNA binding is assessed.
Date Created
2011
Agent

Rationalization of protein conformational dynamics by molecular simulations: studies of the ERK2 kinase and the LAC repressor headpiece-O1 operator complex

150032-Thumbnail Image.png
Description
Molecular dynamics (MD) simulations provide a particularly useful approach to understanding conformational change in biomolecular systems. MD simulations provide an atomistic, physics-based description of the motions accessible to biomolecular systems on the pico- to micro-second timescale, yielding important insight into

Molecular dynamics (MD) simulations provide a particularly useful approach to understanding conformational change in biomolecular systems. MD simulations provide an atomistic, physics-based description of the motions accessible to biomolecular systems on the pico- to micro-second timescale, yielding important insight into the free energy of the system, the dynamical stability of contacts and the role of correlated motions in directing the motions of the system. In this thesis, I use molecular dynamics simulations to provide molecular mechanisms that rationalize structural, thermodynamic, and mutation data on the interactions between the lac repressor headpiece and its O1 operator DNA as well as the ERK2 protein kinase. I performed molecular dynamics simulations of the lac repressor headpiece - O1 operator complex at the natural angle as well as at under- and overbent angles to assess the factors that determine the natural DNA bending angle. I find both energetic and entropic factors contribute to recognition of the natural angle. At the natural angle the energy of the system is minimized by optimization of protein-DNA contacts and the entropy of the system is maximized by release of water from the protein-DNA interface and decorrelation of protein motions. To identify the mechanism by which mutations lead to auto-activation of ERK2, I performed a series of molecular dynamics simulations of ERK1/2 in various stages of activation as well as the constitutively active Q103A, I84A, L73P and R65S ERK2 mutants. My simulations indicate the importance of domain closure for auto-activation and activity regulation. My results enable me to predict two loss-of-function mutants of ERK2, G83A and Q64C, that have been confirmed in experiments by collaborators. One of the powerful capabilities of MD simulations in biochemistry is the ability to find low free energy pathways that connect and explain disparate structural data on biomolecular systems. An extention of the targeted molecular dynamics technique using constraints on internal coordinates will be presented and evaluated. The method gives good results for the alanine dipeptide, but breaks down when applied to study conformational changes in GroEL and adenylate kinase.
Date Created
2011
Agent