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

Interaction Pattern of Arg 62 in the A-Pocket of Differentially Disease-Associated HLA-B27 Subtypes Suggests Distinct TCR Binding Modes

The single amino acid replacement Asp116His distinguishes the two subtypes HLA-B*2705 and HLA-B*2709 which are, respectively, associated and non-associated with Ankylosing Spondylitis, an autoimmune chronic inflammatory disease. The reason for this differential association is so far poorly understood and might be related to subtype-specific HLA:peptide conformations as well as to subtype/peptide-dependent dynamical properties on the nanoscale. Here, we combine functional experiments with extensive molecular dynamics simulations to investigate the molecular dynamics and function of the conserved Arg62 of the α1-helix for both B27 subtypes in complex with the self-peptides pVIPR (RRKWRRWHL) and TIS (RRLPIFSRL), and the viral peptides pLMP2 (RRRWRRLTV) and NPflu (SRYWAIRTR). Simulations of HLA:peptide systems suggest that peptide-stabilizing interactions of the Arg62 residue observed in crystal structures are metastable for both B27 subtypes under physiological conditions, rendering this arginine solvent-exposed and, probably, a key residue for TCR interaction more than peptide-binding. This view is supported by functional experiments with conservative (R62K) and non-conservative (R62A) B*2705 and B*2709 mutants that showed an overall reduction in their capability to present peptides to CD8+ T cells. Moreover, major subtype-dependent differences in the peptide recognition suggest distinct TCR binding modes for the B*2705 versus the B*2709 subtype

Coulomb Interactions between Cytoplasmic Electric Fields and Phosphorylated Messenger Proteins Optimize Information Flow in Cells

Normal cell function requires timely and accurate transmission of information from receptors on the cell membrane (CM) to the nucleus. Movement of messenger proteins in the cytoplasm is thought to be dependent on random walk. However, Brownian motion will disperse messenger proteins throughout the cytosol resulting in slow and highly variable transit times. We propose that a critical component of information transfer is an intracellular electric field generated by distribution of charge on the nuclear membrane (NM). While the latter has been demonstrated experimentally for decades, the role of the consequent electric field has been assumed to be minimal due to a Debye length of about 1 nanometer that results from screening by intracellular Cl- and K+. We propose inclusion of these inorganic ions in the Debye-Huckel equation is incorrect because nuclear pores allow transit through the membrane at a rate far faster than the time to thermodynamic equilibrium. In our model, only the charged, mobile messenger proteins contribute to the Debye length.Using this revised model and published data, we estimate the NM possesses a Debye-Huckel length of a few microns and find this is consistent with recent measurement using intracellular nano-voltmeters. We demonstrate the field will accelerate isolated messenger proteins toward the nucleus through Coulomb interactions with negative charges added by phosphorylation. We calculate transit times as short as 0.01 sec. When large numbers of phosphorylated messenger proteins are generated by increasing concentrations of extracellular ligands, we demonstrate they generate a self-screening environment that regionally attenuates the cytoplasmic field, slowing movement but permitting greater cross talk among pathways. Preliminary experimental results with phosphorylated RAF are consistent with model predictions.This work demonstrates that previously unrecognized Coulomb interactions between phosphorylated messenger proteins and intracellular electric fields will optimize information transfer from the CM to the NM in cells

Predicting Important Residues and Interaction Pathways in Proteins Using Gaussian Network Model: Binding and Stability of HLA Proteins

A statistical thermodynamics approach is proposed to determine structurally and functionally important residues in native proteins that are involved in energy exchange with a ligand and other residues along an interaction pathway. The structure-function relationships, ligand binding and allosteric activities of ten structures of HLA Class I proteins of the immune system are studied by the Gaussian Network Model. Five of these models are associated with inflammatory rheumatic disease and the remaining five are properly functioning. In the Gaussian Network Model, the protein structures are modeled as an elastic network where the inter-residue interactions are harmonic. Important residues and the interaction pathways in the proteins are identified by focusing on the largest eigenvalue of the residue interaction matrix. Predicted important residues match those known from previous experimental and clinical work. Graph perturbation is used to determine the response of the important residues along the interaction pathway. Differences in response patterns of the two sets of proteins are identified and their relations to disease are discussed

Structural and dynamic insights into Mn4Ca cluster-depleted Photosystem II

In the first steps of natural oxygenic photosynthesis, sunlight is used to oxidize water molecules to protons, electrons and molecular oxygen. This reaction takes place on the Mn4Ca cluster located in the reaction centre of Photosystem II (PSII), where the cluster is assembled and continuously repaired through a process known as photoactivation. Understanding the molecular details of such a process has important implications in different fields, in particular inspiring synthesis and repair strategies for artificial photosynthesis devices. In this regard, a detailed structural and dynamic characterization of Photosystem II lacking a Mn4Ca cluster, namelyapoPSII, is a prerequisite for the full comprehension of the photoactivation. Recently, the structure of theapoPSII was resolved at 2.55 Å resolution [Zhanget al.,eLife, 2017,6, e26933], suggesting a pre-organized structure of the protein cavity hosting the cluster. Anyway, the question of whether these findings are a feature of the method used remains open. Here, by means of classical Molecular Dynamics simulations, we characterized the structural and dynamic features of theapoPSII for different protonation states of the cluster cavity. Albeit an overall conformational stability common to all investigated systems, we found significant deviations in the conformation of the side chains of the active site with respect to the X-ray positions. Our findings suggest that not all residues acting as Mn ligands are pre-organized prior to the Mn4Ca formation and previous local conformational changes are required in order to bind the first Mn ion in the high-affinity binding site

Mechanism of Oxygen Evolution and Mn4CaO5 Cluster Restoration in the Natural Water-Oxidizing Catalyst

Water oxidation occurring in the first steps of natural oxygenic photosynthesis is catalyzed by the pigment/protein complex Photosystem II. This process takes place on the Mn4Ca cluster located in the core of Photosystem II and proceeds along the five steps (S0-S4) of the so-called Kok-Joliot cycle until the release of molecular oxygen. The catalytic cycle can therefore be started afresh through insertion of a new water molecule. Here, combining quantum mechanics/molecular mechanics simulations and minimum energy path calculations, we characterized on different spin surfaces the events occurring in the last sector of the catalytic cycle from structural, electronic, and thermodynamic points of view. We found that the process of oxygen evolution and water insertion can be described well by a two-step mechanism, with oxygen release being the rate-limiting step of the process. Moreover, our results allow us to identify the upcoming water molecule required to regenerate the initial structure of the Mn4Ca cluster in the S0 state. The insertion of the water molecule was found to be coupled with the transfer of a proton to a neighboring hydroxide ion, thus resulting in the reconstitution of the most widely accepted model of the S0 state

Protein folding pathways revealed by essential dynamics sampling

The characterization of the protein folding process represents one of the major challenges in molecular biology. Here, a method to simulate the folding process of a protein to its native state is reported, the essential dynamics sampling (EDS) method, and is successfully applied to detecting the correct folding pathways of two small proteins, the all-beta SH3 domain of Src tyrosine kinase transforming protein (SH3) and the alpha/beta B1 domain of streptococcal protein G (GB1). The main idea of the method is that a subset of the natural modes of fluctuation in the native state is key in directing the folding process. A biased molecular dynamics simulation is performed, in which the restrained degrees of freedom are chosen among those obtained by a principal component, or essential dynamics, analysis of the positional fluctuations of the C alpha atoms in the native state. Successful folding is obtained if the restraints are applied only to the eigenvectors with lowest eigenvalues, representing the most rigid quasi-constraint motions. If the essential eigenvectors, the ones accounting for most of the variance, are used, folding is not successful. These results clearly show that the eigenvectors with lowest eigenvalues contain the main mechanical information necessary to drive the folding process, while the essential eigenvectors represent the large concerted motions which can occur without folding/unfolding the protein