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

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 ves hypothesis and protein misfolding

Proteins function by changing conformation. These conformational changes, which involve the concerted motion of a large number of atoms are classical events but, in many cases, the triggers are quantum mechani- cal events such as chemical reactions. Here the initial quantum states after the chemical reaction are assumed to be vibrational excited states, something that has been designated as the VES hypothesis. While the dynamics under classical force fields fail to explain the relatively lower structural stability of the proteins associated with misfolding diseases, the application of the VES hy- pothesis to two cases can provide a new explanation for this phenomenon. This explanation relies on the transfer of vibrational energy from water molecules to proteins, a process whose viability is also examined

Polyglutamine Induced Misfolding of Huntingtin Exon1 is Modulated by the Flanking Sequences

Polyglutamine (polyQ) expansion in exon1 (XN1) of the huntingtin protein is linked to Huntington's disease. When the number of glutamines exceeds a threshold of approximately 36–40 repeats, XN1 can readily form amyloid aggregates similar to those associated with disease. Many experiments suggest that misfolding of monomeric XN1 plays an important role in the length-dependent aggregation. Elucidating the misfolding of a XN1 monomer can help determine the molecular mechanism of XN1 aggregation and potentially help develop strategies to inhibit XN1 aggregation. The flanking sequences surrounding the polyQ region can play a critical role in determining the structural rearrangement and aggregation mechanism of XN1. Few experiments have studied XN1 in its entirety, with all flanking regions. To obtain structural insights into the misfolding of XN1 toward amyloid aggregation, we perform molecular dynamics simulations on monomeric XN1 with full flanking regions, a variant missing the polyproline regions, which are hypothesized to prevent aggregation, and an isolated polyQ peptide (Qn). For each of these three constructs, we study glutamine repeat lengths of 23, 36, 40 and 47. We find that polyQ peptides have a positive correlation between their probability to form a β-rich misfolded state and their expansion length. We also find that the flanking regions of XN1 affect its probability to^x_page_count=28 form a β-rich state compared to the isolated polyQ. Particularly, the polyproline regions form polyproline type II helices and decrease the probability of the polyQ region to form a β-rich state. Additionally, by lengthening polyQ, the first N-terminal 17 residues are more likely to adopt a β-sheet conformation rather than an α-helix conformation. Therefore, our molecular dynamics study provides a structural insight of XN1 misfolding and elucidates the possible role of the flanking sequences in XN1 aggregation

Controlling spins in adsorbed molecules by a chemical switch

The development of chemical systems with switchable molecular spins could lead to the architecture of materials with controllable magnetic or spintronic properties. Here, we present conclusive evidence that the spin of an organometallic molecule coupled to a ferromagnetic substrate can be switched between magnetic off and on states by a chemical stimulus. This is achieved by nitric oxide (NO) functioning as an axial ligand of cobalt(II)tetraphenylporphyrin (CoTPP) ferromagnetically coupled to nickel thin-film (Ni(001)). On NO addition, the coordination sphere of Co2+ is modified and a NO–CoTPP nitrosyl complex is formed, which corresponds to an off state of the Co spin. Thermal dissociation of NO from the nitrosyl complex restores the on state of the Co spin. The NO-induced reversible off–on switching of surface-adsorbed molecular spins observed here is attributed to a spin trans effect