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

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

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

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

Dynamical characterization of two differentially disease associated MHC class i proteins in complex with viral and self-peptides

Major histocompatibility complex (MHC) class I proteins are expressed on the cell surface where they present foreign and self-peptides to effector cells of the immune system. While an understanding of the structural prerequisites for antigen presentation has already been achieved, insight into subtype- or peptide-dependent dynamical characteristics of a peptide-MHC antigen is so far largely obscure. We approached this problem by employing 400-ns molecular dynamics simulations with two human MHC class I subtypes as model systems: the ankylosing spondylitis-associated HLA-B 27:05 and the non-ankylosing spondylitis-associated HLA-B 27:09. Both proteins differ only by a micropolymorphism at the floor of the peptide binding groove (Asp116His). A viral (pLMP2) and three self-peptides (pVIPR, pGR, and TIS) were evaluated. The stability of the binding grooves was found to be both subtype dependent and peptide dependent. A detachment from the C- and/or N-terminal pockets was observed for all peptides except TIS, resulting in a stabilization of the α1-helix in both TIS-displaying subtypes. Estimates of the entropy associated with the bound peptides showed an increased entropy for pLMP2 presented by B 27:05 as compared to B 27:09, in contrast to the self-peptides. Additionally, the flexibility of the α1-helix that is probably important for receptor binding to the B27:peptide epitope is significantly enhanced for B 27:05. These in silico results show that the dynamic properties of peptide-MHC complexes are affected both by the bound peptide and by micropolymorphisms of the heavy chain. Our findings suggest a role for the conformational flexibility of MHC class I molecules in the context of recognition by receptors on effector cells. © 2011 Elsevier Ltd

Molecular dynamics simulations and kinetic measurements provide insights into the structural requirements of substrate size-dependent specificity of oligogalacturonide oxidase 1 (OGOX1)

Oligogalacturonides (OGs) are pectin fragments released from the breakdown of the homogalacturonan during pathogenesis that act as Damage-Associated Molecular Patterns. OG-oxidase 1 (OGOX1) is an Arabidopsis berberine bridge enzyme-like (BBE-l) oligosaccharide oxidase that oxidizes OGs, impairing their elicitor activity and concomitantly releasing H2O2. The OG-oxidizing activity of OGOX1 is markedly pH-dependent, with optimum pH around 10, and is higher towards OGs with a degree of polymerization higher than two. Here, the molecular determinants of OGOX1 responsible for the binding of OGs with different lengths have been investigated through molecular dynamics simulations and enzyme kinetics studies. OGOX1 was simulated in complex with OGs with different degree of polymerization such as di-, tri-, tetra- and penta-galacturonide, in water solution at alkaline pH. Our simulations revealed that, among the four OGOX1/OG combinations, the pentagalacturonide (OG5) showed the best conformation in the active site to be efficiently oxidized by OGOX1. The optimal conformation can be stabilized by salt-bridges formed between the carboxyl groups of OG5 and five positively charged amino acids of OGOX1, highly conserved in all OGOX paralogs. Our results suggest that these interactions limit the mobility of OG5 as well as longer OGs, contributing to maintain the terminal monomer of OGs in the optimal orientation in order to be oxidized by the enzyme. In accordance with these results, the enzyme efficiency (K-cat/K-M) of OGOX1 on OG5 (40.04) was found to be significantly higher than that on OG4 (13.05) and OG3 (0.6)

Evidence for proton shuffling in a thioredoxin-like protein during catalysis

Proteins of the thioredoxin (Trx) superfamily catalyze disulfide-bond formation, reduction and isomerization in substrate proteins both in prokaryotic and in eukaryotic cells. All members of the Trx family with thiol-disulfide oxidoreductase activity contain the characteristic Cys-X-X-Cys motif in their active site. Here, using Poisson-Boltzmann-based protonation-state calculations based on 100-ns molecular dynamics simulations, we investigate the catalytic mechanism of DsbL, the most oxidizing Trx-like protein known to date. We observed several correlated transitions in the protonation states of the buried active-site cysteine and a neighboring lysine coupled to the exposure of the active-site thiolate. These results support the view of an internal proton shuffling mechanism during oxidation crucial for the uptake of two electrons from the substrate protein. Intramolecular disulfide-bond formation is probably steered by the conformational switch facilitating interaction with the active-site thiolate. A consistent catalytic mechanism for DsbL, probably conferrable to other proteins of the same class, is presented. Our results suggest a functional role of hydration entropy of active-site groups