Metadynamics study of mutant human interferon-gamma forms

Human interferon-gamma ([email protected]) is an important antiviral and immunomodulating signaling molecule. The upregulation of its production, however, is related to the etiology of certain autoimmune diseases. In the search for a mechanism for suppressing the [email protected] biological activity, we investigated the possibility to obtain mutant derivatives of the protein, capable to bind to [email protected] cellular receptors, but lacking the ability to trigger the biological response inside the cell. In order to preserve the affinity to the receptor, the introduced mutations should not induce conformational changes in the secondary structure of the mutants. Molecular dynamics simulations were performed to study the secondary structure of 100 randomly chosen [email protected] derivatives with substitutions in amino acids 86-88. The stability of the local structure of all [email protected] forms was investigated by means of metadynamics. It was found that some of the mutated forms preserve the local secondary structure and show similar or higher stability of the mutated helix, compared to the native form. The 12 most promising mutants were suggested for experimental investigation.

[1]  A. Laio,et al.  Metadynamics: a method to simulate rare events and reconstruct the free energy in biophysics, chemistry and material science , 2008 .

[2]  K Schulten,et al.  VMD: visual molecular dynamics. , 1996, Journal of molecular graphics.

[3]  P. Subramaniam,et al.  The COOH-terminal nuclear localization sequence of interferon gamma regulates STAT1 alpha nuclear translocation at an intracellular site. , 2000, Journal of cell science.

[4]  D. Jans,et al.  Trafficking and signaling pathways of nuclear localizing protein ligands and their receptors , 2004, BioEssays : news and reviews in molecular, cellular and developmental biology.

[5]  A. Leach.  Molecular Modelling: Principles and Applications , 1996 .

[6]  D. van der Spoel,et al.  GROMACS: A message-passing parallel molecular dynamics implementation , 1995 .

[7]  Alexander D. MacKerell,et al.  All-atom empirical potential for molecular modeling and dynamics studies of proteins. , 1998, The journal of physical chemistry. B.

[8]  Holger Gohlke,et al.  The Amber biomolecular simulation programs , 2005, J. Comput. Chem..

[9]  W. Windsor,et al.  Crystal structure of a complex between interferon-γ and its soluble high-affinity receptor , 1995, Nature.

[10]  S. Narula,et al.  The carboxyl-terminal region of human interferon gamma is important for biological activity: mutagenic and NMR analysis. , 1991, Protein engineering.

[11]  P. Petkov,et al.  Computer simulations of human interferon gamma mutated forms , 2010 .

[12]  W. Delano.  The PyMOL Molecular Graphics System , 2002 .

[13]  Berend Smit,et al.  Understanding Molecular Simulation , 2001 .

[14]  Giovanni Ciccotti,et al.  Molecular dynamics simulation of rigid molecules , 1986 .

[15]  K. Schroder,et al.  Interferon- : an overview of signals, mechanisms and functions , 2004 .

[16]  A. Karshikoff,et al.  Significance of the Putative Upstream Polybasic Nuclear Localisation Sequence for the Biological Activity of Human Interferon-Gamma , 2009 .

[17]  T. Darden,et al.  Particle mesh Ewald: An N⋅log(N) method for Ewald sums in large systems , 1993 .

[18]  Kremer,et al.  Molecular dynamics simulation for polymers in the presence of a heat bath. , 1986, Physical review. A, General physics.

[19]  Portland Press Ltd.  IUPAC-IUB Joint Commission on Biochemical Nomenclature (JCBN). Nomenclature and symbolism for amino acids and peptides. Recommendations 1983 , 1984 .

[20]  A. Karshikoff,et al.  Human interferon gamma: significance of the C-terminal flexible domain for its biological activity. , 2003, Archives of biochemistry and biophysics.

[21]  Francesco Luigi Gervasio,et al.  Comparing the Efficiency of Biased and Unbiased Molecular Dynamics in Reconstructing the Free Energy Landscape of Met-Enkephalin , 2010 .

[22]  F. Winkler,et al.  Observation of an unexpected third receptor molecule in the crystal structure of human interferon-gamma receptor complex. , 2000, Structure.

[23]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[24]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[25]  F. Jay,et al.  The E1 functional epitope of the human interferon gamma is a nuclear targeting signal-like element. Mapping of the E1 epitope. , 1991, The Journal of biological chemistry.

[26]  S E Ealick,et al.  Three-dimensional structure of recombinant human interferon-gamma. , 1991, Science.

[27]  J. Darnell,et al.  Jak-STAT pathways and transcriptional activation in response to IFNs and other extracellular signaling proteins. , 1994, Science.

[28]  A. Karshikoff,et al.  Role of the C‐terminal chain in human interferonγ stability: An electrostatic study , 2001 .

[29]  W. Kabsch,et al.  Dictionary of protein secondary structure: Pattern recognition of hydrogen‐bonded and geometrical features , 1983, Biopolymers.

[30]  Laxmikant V. Kalé,et al.  Scalable molecular dynamics with NAMD , 2005, J. Comput. Chem..

[31]  M. Paddy,et al.  The Carboxyl Terminus of Interferon-γ Contains a Functional Polybasic Nuclear Localization Sequence* , 1999, The Journal of Biological Chemistry.

[32]  G. N. Ramachandran,et al.  Stereochemistry of polypeptide chain configurations. , 1963, Journal of molecular biology.

[33]  S. Skurkovich,et al.  Anti-interferon-gamma antibodies in the treatment of autoimmune diseases. , 2003, Current opinion in molecular therapeutics.