skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Phase space trajectories and Lyapunov exponents in the dynamics of an alpha-helical protein lattice with intra- and inter-spine interactions

Abstract

The nonlinear dynamics of intra- and inter-spine interaction models of alpha-helical proteins is investigated by proposing a Hamiltonian using the first quantized operators. Hamilton's equations of motion are derived, and the dynamics is studied by constructing the trajectories and phase space plots in both cases. The phase space plots display a chaotic behaviour in the dynamics, which opens questions about the relationship between the chaos and exciton-exciton and exciton-phonon interactions. This is verified by plotting the Lyapunov characteristic exponent curves.

Authors:
;  [1];  [2]
  1. Department of Physics, Women's Christian College, Nagercoil 629 001 (India)
  2. Nuclear Physics Division, Bhabha Atomic Research Centre, Mumbai 400085 (India)
Publication Date:
OSTI Identifier:
22482296
Resource Type:
Journal Article
Resource Relation:
Journal Name: Chaos (Woodbury, N. Y.); Journal Volume: 25; Journal Issue: 11; Other Information: (c) 2015 AIP Publishing LLC; Country of input: International Atomic Energy Agency (IAEA)
Country of Publication:
United States
Language:
English
Subject:
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSICS; CHAOS THEORY; DIAGRAMS; EQUATIONS OF MOTION; HAMILTONIANS; INTERACTIONS; LYAPUNOV METHOD; NONLINEAR PROBLEMS; PHASE SPACE; PHONONS; PROTEINS; TRAJECTORIES

Citation Formats

Angelin Jeba, K., Latha, M. M., E-mail: lathaisaac@yahoo.com, and Jain, Sudhir R.. Phase space trajectories and Lyapunov exponents in the dynamics of an alpha-helical protein lattice with intra- and inter-spine interactions. United States: N. p., 2015. Web. doi:10.1063/1.4935957.
Angelin Jeba, K., Latha, M. M., E-mail: lathaisaac@yahoo.com, & Jain, Sudhir R.. Phase space trajectories and Lyapunov exponents in the dynamics of an alpha-helical protein lattice with intra- and inter-spine interactions. United States. doi:10.1063/1.4935957.
Angelin Jeba, K., Latha, M. M., E-mail: lathaisaac@yahoo.com, and Jain, Sudhir R.. 2015. "Phase space trajectories and Lyapunov exponents in the dynamics of an alpha-helical protein lattice with intra- and inter-spine interactions". United States. doi:10.1063/1.4935957.
@article{osti_22482296,
title = {Phase space trajectories and Lyapunov exponents in the dynamics of an alpha-helical protein lattice with intra- and inter-spine interactions},
author = {Angelin Jeba, K. and Latha, M. M., E-mail: lathaisaac@yahoo.com and Jain, Sudhir R.},
abstractNote = {The nonlinear dynamics of intra- and inter-spine interaction models of alpha-helical proteins is investigated by proposing a Hamiltonian using the first quantized operators. Hamilton's equations of motion are derived, and the dynamics is studied by constructing the trajectories and phase space plots in both cases. The phase space plots display a chaotic behaviour in the dynamics, which opens questions about the relationship between the chaos and exciton-exciton and exciton-phonon interactions. This is verified by plotting the Lyapunov characteristic exponent curves.},
doi = {10.1063/1.4935957},
journal = {Chaos (Woodbury, N. Y.)},
number = 11,
volume = 25,
place = {United States},
year = 2015,
month =
}
  • We study the behavior of the generalized Lyapunov exponents for chaotic symplectic dynamical systems and products of random matrices in the limit of large dimensions D. For products of random matrices without any particular structure the generalized Lyapunov exponents become equal in this limit and the value of one of the generalized Lyapunov exponents is obtained by simple arguments. On the contrary, for random symplectic matrices with peculiar structures and for chaotic symplectic maps the generalized Lyapunov exponents remains different for D ..-->.. infinity, indicating that high dimensionality cannot always destroy intermittency.
  • Purpose: The intra- and inter-observer variability in delineation of the parotids on the kilo-voltage computed tomography (kVCT) and mega-voltage computed tomography (MVCT) were examined to establish their impact on the dose calculation during adaptive head and neck helical tomotherapy (HT). Methods: Three observers delineated left and right parotids for ten randomly selected patients with oropharynx cancer treated on HT. The pre-treatment kVCT and the MVCT from the first fraction of irradiation were selected to delineation. The delineation procedure was repeated three times by each observer. The parotids were delineated according to the institutional protocol. The analyses included intra-observer reproducibility andmore » inter-structure, -observer and -modality variability of the volume and dose. Results: The differences between the left and right parotid outlines were not statistically significant (p>0.3). The reproducibility of the delineation was confirmed for each observer on the kVCT (p>0.2) and on the MVCT (p>0.1). The inter-observer variability of the outlines was significant (p<0.001) as well as the inter-modality variability (p<0.006). The parotids delineated on the MVCT were 10% smaller than on the kVCT. The inter-observer variability of the parotids delineation did not affect the average dose (p=0.096 on the kVCT and p=0.176 on the MVCT). The dose calculated on the MVCT was higher by 3.3% than dose from the kVCT (p=0.009). Conclusion: Usage of the institutional protocols for the parotids delineation reduces intra-observer variability and increases reproducibility of the outlines. These protocols do not eliminate delineation differences between the observers, but these differences are not clinically significant and do not affect average doses in the parotids. The volumes of the parotids delineated on the MVCT are smaller than on the kVCT, which affects the differences in the calculated doses.« less
  • The spatial dynamics of steady, one-dimensional premixed H{sub 2}/O{sub 2} flames are studied. The emphasis in this Article is the geometry of the phase space of the dynamical system describing the steady flames. It is shown that steady flames are described by trajectories on the stable manifolds of saddle fixed points. The saddle fixed points correspond to equilibrium points of time-dependent chemical-kinetic systems that are adiabatic and isobaric and whose constant enthalpy matches the asymptotic enthalpy of the flames. The dimensions of the stable manifolds match the dimensions of the chemical-kinetic systems under most conditions, although the dynamics on themmore » are different. It is further shown that the stable manifolds have low-dimensional attractive submanifolds near the saddlepoint. An algorithm for generating trajectories over the spatial domain of these flames is proposed, and it is used to study individual trajectories and trajectory ensembles, whose collective behavior suggests that there are low-dimensional submanifolds away from the saddlepoint.« less
  • Protein {alpha}-helices are ubiquitous secondary structural elements, seldom considered to be stable without tertiary contacts. However, amino acid sequences in proteins that are based on alternating repeats of four glutamic acid (E) residues and four positively charged residues, a combination of arginine (R) and lysine (K), have been shown to form stable {alpha}-helices in a few proteins, in the absence of tertiary interactions. Here, we find that this ER/K motif is more prevalent than previously reported, being represented in proteins of diverse function from archaea to humans. By using molecular dynamics (MD) simulations, we characterize a dynamic pattern of side-chainmore » interactions that extends along the backbone of ER/K {alpha}-helices. A simplified model predicts that side-chain interactions alone contribute substantial bending rigidity (0.5 pN/nm) to ER/K {alpha}-helices. Results of small-angle x-ray scattering (SAXS) and single-molecule optical-trap analyses are consistent with the high bending rigidity predicted by our model. Thus, the ER/K {alpha}-helix is an isolated secondary structural element that can efficiently span long distances in proteins, making it a promising tool in designing synthetic proteins. We propose that the significant rigidity of the ER/K {alpha}-helix can help regulate protein function, as a force transducer between protein subdomains.« less
  • Understanding the structure and dynamics of membrane proteins in their native, hydrophobic environment is important to understanding how these proteins function. EPR spectroscopy in combination with site-directed spin labeling (SDSL) can measure dynamics and structure of membrane proteins in their native lipid environment; however, until now the dynamics measured have been qualitative due to limited knowledge of the nitroxide spin label's intramolecular motion in the hydrophobic environment. Although several studies have elucidated the structural origins of EPR line shapes of water-soluble proteins, EPR spectra of nitroxide spin-labeled proteins in detergents or lipids have characteristic differences from their water-soluble counterparts, suggestingmore » significant differences in the underlying molecular motion of the spin label between the two environments. To elucidate these differences, membrane-exposed {alpha}-helical sites of the leucine transporter, LeuT, from Aquifex aeolicus, were investigated using X-ray crystallography, mutational analysis, nitroxide side chain derivatives, and spectral simulations in order to obtain a motional model of the nitroxide. For each crystal structure, the nitroxide ring of a disulfide-linked spin label side chain (R1) is resolved and makes contacts with hydrophobic residues on the protein surface. The spin label at site I204 on LeuT makes a nontraditional hydrogen bond with the ortho-hydrogen on its nearest neighbor F208, whereas the spin label at site F177 makes multiple van der Waals contacts with a hydrophobic pocket formed with an adjacent helix. These results coupled with the spectral effect of mutating the i {+-} 3, 4 residues suggest that the spin label has a greater affinity for its local protein environment in the low dielectric than on a water-soluble protein surface. The simulations of the EPR spectra presented here suggest the spin label oscillates about the terminal bond nearest the ring while maintaining weak contact with the protein surface. Combined, the results provide a starting point for determining a motional model for R1 on membrane proteins, allowing quantification of nitroxide dynamics in the aliphatic environment of detergent and lipids. In addition, initial contributions to a rotamer library of R1 on membrane proteins are provided, which will assist in reliably modeling the R1 conformational space for pulsed dipolar EPR and NMR paramagnetic relaxation enhancement distance determination.« less