Role of the Subunits Interactions in the Conformational Transitions in Adult Human Hemoglobin: an Explicit Solvent Molecular Dynamics Study

Hemoglobin exhibits allosteric structural changes upon ligand binding due to the dynamic interactions between the ligand binding sites, the amino acids residues and some other solutes present under physiological conditions. In the present study, the dynamical and quaternary structural changes occurring in two unligated (deoxy-) T structures, and two fully ligated (oxy-) R, R2 structures of adult human hemoglobin were investigated with molecular dynamics. It is shown that, in the sub-microsecond time scale, there is no marked difference in the global dynamics of the amino acids residues in both the oxy- and the deoxy- forms of the individual structures. In addition, the R, R2 are relatively stable and do not present quaternary conformational changes within the time scale of our simulations while the T structure is dynamically more flexible and exhibited the T\rightarrow R quaternary conformational transition, which is propagated by the relative rotation of the residues at the {\alpha}1{\beta}2 and {\alpha}2{\beta}1 interface.Comment: Reprinted (adapted) with permission from J. Phys. Chem. B DOI:10.1021/jp3022908. Copyright (2012) American Chemical Societ

Local Simulation Algorithms for Coulomb Interaction

Long ranged electrostatic interactions are time consuming to calculate in molecular dynamics and Monte-Carlo simulations. We introduce an algorithmic framework for simulating charged particles which modifies the dynamics so as to allow equilibration using a local Hamiltonian. The method introduces an auxiliary field with constrained dynamics so that the equilibrium distribution is determined by the Coulomb interaction. We demonstrate the efficiency of the method by simulating a simple, charged lattice gas.Comment: Last figure changed to improve demonstration of numerical efficienc

Polar zipper sequence in the high-affinity hemoglobin of Ascaris suum: amino acid sequence and structural interpretation.

The extracellular hemoglobin of Ascaris has an extremely high oxygen affinity (P50 = 0.004 mmHg). It consists of eight identical subunits of molecular weight 40,600. Their sequence, determined by protein chemistry, shows two tandemly linked globin-like sequences and an 18-residue C-terminal extension. Two N-linked glycosylation sites contain equal ratios of mannose/glucosamine/fucose of 3:2:1. Electron micrographs suggest that the eight subunits form a polyhedron of point symmetry D4, or 42. The C-terminal extension contains a repeat of the sequence Glu-Glu-His-Lys, which would form a pattern of alternate glutamate and histidine side chains on one side and of glutamate and lysine side chains on the other side of a beta strand. We propose that this represents a polar zipper sequence and that the C-terminal extensions are joined in an eight-stranded beta barrel at the center of the molecule, with histidine and glutamate side chains inside and lysine and glutamate side chains outside the barrel compensating each other's charges. The amino acid sequence of Ascaris hemoglobin fails to explain its high oxygen affinity

Structural biology and phylogenetic estimation

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/62633/1/388527a0.pd

Kinetics of stochastically-gated diffusion-limited reactions and geometry of random walk trajectories

In this paper we study the kinetics of diffusion-limited, pseudo-first-order A + B -> B reactions in situations in which the particles' intrinsic reactivities vary randomly in time. That is, we suppose that the particles are bearing "gates" which interchange randomly and independently of each other between two states - an active state, when the reaction may take place, and a blocked state, when the reaction is completly inhibited. We consider four different models, such that the A particle can be either mobile or immobile, gated or ungated, as well as ungated or gated B particles can be fixed at random positions or move randomly. All models are formulated on a $d$-dimensional regular lattice and we suppose that the mobile species perform independent, homogeneous, discrete-time lattice random walks. The model involving a single, immobile, ungated target A and a concentration of mobile, gated B particles is solved exactly. For the remaining three models we determine exactly, in form of rigorous lower and upper bounds, the large-N asymptotical behavior of the A particle survival probability. We also realize that for all four models studied here such a probalibity can be interpreted as the moment generating function of some functionals of random walk trajectories, such as, e.g., the number of self-intersections, the number of sites visited exactly a given number of times, "residence time" on a random array of lattice sites and etc. Our results thus apply to the asymptotical behavior of the corresponding generating functions which has not been known as yet.Comment: Latex, 45 pages, 5 ps-figures, submitted to PR

Structural Properties of Polyglutamine Aggregates Investigated via Molecular Dynamics Simulations

Polyglutamine (polyQ) beta-stranded aggregates constitute the hallmark of Huntington disease. The disease is fully penetrant when Q residues are more than 36-40 ("disease threshold"). Here, based on a molecular dynamics study on polyQ helical structures of different shapes and oligomeric states, we suggest that the stability of the aggregates increases with the number of monomers, while it is rather insensitive to the number of Qs in each monomer. However, the stability of the single monomer does depend on the number of side-chain intramolecular H-bonds, and therefore oil the number of Qs. If such number is lower than that of the disease threshold, the beta-stranded monomers are unstable and hence may aggregate with lower probability, consistently with experimental findings. Our results provide a possible interpretation of the apparent polyQ length dependent-toxicity, and they do not support the so-called "structural threshold hypothesis", which supposes a transition from random coil to a beta-sheet structure only above the disease threshold

Notes on bordered Floer homology

This is a survey of bordered Heegaard Floer homology, an extension of the Heegaard Floer invariant HF-hat to 3-manifolds with boundary. Emphasis is placed on how bordered Heegaard Floer homology can be used for computations.Comment: 73 pages, 29 figures. Based on lectures at the Contact and Symplectic Topology Summer School in Budapest, July 2012. v2: Fixed many small typo

Academic freedom: in justification of a universal ideal

This paper examines the justification for, and benefits of, academic freedom to academics, students, universities and the world at large. The paper surveys the development of the concept of academic freedom within Europe, more especially the impact of the reforms at the University of Berlin instigated by Wilhelm von Humboldt. Following from this, the paper examines the reasons why the various facets of academic freedom are important and why the principle should continue to be supported

The structure of a class 3 nonsymbiotic plant haemoglobin from<i>Arabidopsis thaliana</i>reveals a novel N-terminal helical extension

Plant nonsymbiotic haemoglobins fall into three classes, each with distinct properties but all with largely unresolved physiological functions. Here, the first crystal structure of a class 3 nonsymbiotic plant haemoglobin, that fromArabidopsis thaliana, is reported to 1.77 Å resolution. The protein forms a homodimer, with each monomer containing a two-over-two α-helical domain similar to that observed in bacterial truncated haemoglobins. A novel N-terminal extension comprising two α-helices plays a major role in the dimer interface, which occupies the periphery of the dimer–dimer face, surrounding an open central cavity. The haem pocket contains a proximal histidine ligand and an open sixth iron-coordination site with potential for a ligand, in this structure hydroxide, to form hydrogen bonds to a tyrosine or a tryptophan residue. The haem pocket appears to be unusually open to the external environment, with another cavity spanning the entrance of the two haem pockets. The final 23 residues of the C-terminal domain are disordered in the structure; however, these domains in the functional dimer are adjacent and include the only two cysteine residues in the protein sequence. It is likely that these residues form disulfide bondsin vitroand it is conceivable that this C-terminal region may act in a putative complex with a partner moleculein vivo.</jats:p

Spontaneous Quaternary and Tertiary T-R Transitions of Human Hemoglobin in Molecular Dynamics Simulation

We present molecular dynamics simulations of unliganded human hemoglobin (Hb) A under physiological conditions, starting from the R, R2, and T state. The simulations were carried out with protonated and deprotonated HC3 histidines His(β)146, and they sum up to a total length of 5.6µs. We observe spontaneous and reproducible T→R quaternary transitions of the Hb tetramer and tertiary transitions of the α and β subunits, as detected from principal component projections, from an RMSD measure, and from rigid body rotation analysis. The simulations reveal a marked asymmetry between the α and β subunits. Using the mutual information as correlation measure, we find that the β subunits are substantially more strongly linked to the quaternary transition than the α subunits. In addition, the tertiary populations of the α and β subunits differ substantially, with the β subunits showing a tendency towards R, and the α subunits showing a tendency towards T. Based on the simulation results, we present a transition pathway for coupled quaternary and tertiary transitions between the R and T conformations of Hb