Incorporating receptor flexibility in the molecular design of protein interfaces

Liwei Li, Shide Liang, Meaghan M. Pilcher, Samy Meroueh

Research output: Contribution to journalArticle

12 Citations (Scopus)

Abstract

The success of antibody-based pharmaceuticals has led to a resurgence in interest in computational structure-based design. Most efforts to date have been on the redesign of existing interfaces. These efforts have mostly neglected the inherent flexibility of the receptor that is critical for binding. In this work, we extend on a previous study to perform a series of designs of protein binding interfaces by incorporating receptor flexibility using an ensemble of conformers collected from explicit-solvent molecular dynamics (MD) simulations. All designer complexes are subjected to 30 ns of MD in explicit solvent to assess for stability for a total of 480 ns of dynamics. This is followed by end-point free energy calculations whereby intermole- cular potential energy, polar and non-polar solvation energy and entropy of ligand and receptor are subtracted from that of the complex and averaged over 320 snapshots collected from each of the 30 ns MD simulations. Our initial effort consisted of redesigning the interface of three well-studied complexes, namely barnase-barstar, lysozyme-antibody D1.3 and trypsin-BPTI. The design was performed with flexible backbone approach. MD simulations revealed that all three complexes remained stable. Interestingly, the redesigned trypsin-BPTI complex was significantly more favorable than the native complex. This was attributed to the favorable electrostatics and entropy that complemented the already favorable non-polar component. Another aspect of this work consisted of grafting the surface of three proteins, namely tenascin, CheY and MBP1 to bind to barnase, trypsin and lysozyme. The process was initially performed using fixed backbone, and more than 300 ns of the explicit-solvent MD simulation revealed some of the complexes to dissociate over the course of the trajectories, whereas others remained stable. Free energy calculations confirmed that the non-polar component of the free energy as computed by summing the van der Waals energy and the non-polar solvation energy was a strong predictor of stability. Four complexes (two stable and two unstable) were selected, and redesigned using multiple conformers collected from the MD simulation. The resulting designer systems were then immersed in explicit solvent and 30 ns of MD was carried out on each. Interestingly, those complexes that were initially stable remained stable, whereas one of the unstable complexes became stable following redesign with flexible backbone. Free energy calculations showed significant improvements in the affinity for most complexes, revealing that the use of multiple conformers in protein design may significantly enhance such efforts.

Original languageEnglish
Pages (from-to)575-586
Number of pages12
JournalProtein Engineering, Design and Selection
Volume22
Issue number9
DOIs
StatePublished - Sep 2009

Fingerprint

Molecular Dynamics Simulation
Molecular dynamics
Proteins
Free energy
Trypsin
Computer simulation
Solvation
Entropy
Muramidase
Antibodies
Enzymes
Tenascin
Potential energy
Static Electricity
Protein Binding
Drug products
Electrostatics
Membrane Proteins
Ligands
Trajectories

Keywords

  • Flexible backbone design
  • Free energy calculations
  • Molecular dynamics simulations
  • Protein engineering
  • Protein-protein interaction

ASJC Scopus subject areas

  • Biochemistry
  • Biotechnology
  • Bioengineering
  • Molecular Biology

Cite this

Incorporating receptor flexibility in the molecular design of protein interfaces. / Li, Liwei; Liang, Shide; Pilcher, Meaghan M.; Meroueh, Samy.

In: Protein Engineering, Design and Selection, Vol. 22, No. 9, 09.2009, p. 575-586.

Research output: Contribution to journalArticle

Li, Liwei ; Liang, Shide ; Pilcher, Meaghan M. ; Meroueh, Samy. / Incorporating receptor flexibility in the molecular design of protein interfaces. In: Protein Engineering, Design and Selection. 2009 ; Vol. 22, No. 9. pp. 575-586.
@article{32da64c0e7fe435396c552eff8e37d87,
title = "Incorporating receptor flexibility in the molecular design of protein interfaces",
abstract = "The success of antibody-based pharmaceuticals has led to a resurgence in interest in computational structure-based design. Most efforts to date have been on the redesign of existing interfaces. These efforts have mostly neglected the inherent flexibility of the receptor that is critical for binding. In this work, we extend on a previous study to perform a series of designs of protein binding interfaces by incorporating receptor flexibility using an ensemble of conformers collected from explicit-solvent molecular dynamics (MD) simulations. All designer complexes are subjected to 30 ns of MD in explicit solvent to assess for stability for a total of 480 ns of dynamics. This is followed by end-point free energy calculations whereby intermole- cular potential energy, polar and non-polar solvation energy and entropy of ligand and receptor are subtracted from that of the complex and averaged over 320 snapshots collected from each of the 30 ns MD simulations. Our initial effort consisted of redesigning the interface of three well-studied complexes, namely barnase-barstar, lysozyme-antibody D1.3 and trypsin-BPTI. The design was performed with flexible backbone approach. MD simulations revealed that all three complexes remained stable. Interestingly, the redesigned trypsin-BPTI complex was significantly more favorable than the native complex. This was attributed to the favorable electrostatics and entropy that complemented the already favorable non-polar component. Another aspect of this work consisted of grafting the surface of three proteins, namely tenascin, CheY and MBP1 to bind to barnase, trypsin and lysozyme. The process was initially performed using fixed backbone, and more than 300 ns of the explicit-solvent MD simulation revealed some of the complexes to dissociate over the course of the trajectories, whereas others remained stable. Free energy calculations confirmed that the non-polar component of the free energy as computed by summing the van der Waals energy and the non-polar solvation energy was a strong predictor of stability. Four complexes (two stable and two unstable) were selected, and redesigned using multiple conformers collected from the MD simulation. The resulting designer systems were then immersed in explicit solvent and 30 ns of MD was carried out on each. Interestingly, those complexes that were initially stable remained stable, whereas one of the unstable complexes became stable following redesign with flexible backbone. Free energy calculations showed significant improvements in the affinity for most complexes, revealing that the use of multiple conformers in protein design may significantly enhance such efforts.",
keywords = "Flexible backbone design, Free energy calculations, Molecular dynamics simulations, Protein engineering, Protein-protein interaction",
author = "Liwei Li and Shide Liang and Pilcher, {Meaghan M.} and Samy Meroueh",
year = "2009",
month = "9",
doi = "10.1093/protein/gzp042",
language = "English",
volume = "22",
pages = "575--586",
journal = "Protein Engineering, Design and Selection",
issn = "1741-0126",
publisher = "Oxford University Press",
number = "9",

}

TY - JOUR

T1 - Incorporating receptor flexibility in the molecular design of protein interfaces

AU - Li, Liwei

AU - Liang, Shide

AU - Pilcher, Meaghan M.

AU - Meroueh, Samy

PY - 2009/9

Y1 - 2009/9

N2 - The success of antibody-based pharmaceuticals has led to a resurgence in interest in computational structure-based design. Most efforts to date have been on the redesign of existing interfaces. These efforts have mostly neglected the inherent flexibility of the receptor that is critical for binding. In this work, we extend on a previous study to perform a series of designs of protein binding interfaces by incorporating receptor flexibility using an ensemble of conformers collected from explicit-solvent molecular dynamics (MD) simulations. All designer complexes are subjected to 30 ns of MD in explicit solvent to assess for stability for a total of 480 ns of dynamics. This is followed by end-point free energy calculations whereby intermole- cular potential energy, polar and non-polar solvation energy and entropy of ligand and receptor are subtracted from that of the complex and averaged over 320 snapshots collected from each of the 30 ns MD simulations. Our initial effort consisted of redesigning the interface of three well-studied complexes, namely barnase-barstar, lysozyme-antibody D1.3 and trypsin-BPTI. The design was performed with flexible backbone approach. MD simulations revealed that all three complexes remained stable. Interestingly, the redesigned trypsin-BPTI complex was significantly more favorable than the native complex. This was attributed to the favorable electrostatics and entropy that complemented the already favorable non-polar component. Another aspect of this work consisted of grafting the surface of three proteins, namely tenascin, CheY and MBP1 to bind to barnase, trypsin and lysozyme. The process was initially performed using fixed backbone, and more than 300 ns of the explicit-solvent MD simulation revealed some of the complexes to dissociate over the course of the trajectories, whereas others remained stable. Free energy calculations confirmed that the non-polar component of the free energy as computed by summing the van der Waals energy and the non-polar solvation energy was a strong predictor of stability. Four complexes (two stable and two unstable) were selected, and redesigned using multiple conformers collected from the MD simulation. The resulting designer systems were then immersed in explicit solvent and 30 ns of MD was carried out on each. Interestingly, those complexes that were initially stable remained stable, whereas one of the unstable complexes became stable following redesign with flexible backbone. Free energy calculations showed significant improvements in the affinity for most complexes, revealing that the use of multiple conformers in protein design may significantly enhance such efforts.

AB - The success of antibody-based pharmaceuticals has led to a resurgence in interest in computational structure-based design. Most efforts to date have been on the redesign of existing interfaces. These efforts have mostly neglected the inherent flexibility of the receptor that is critical for binding. In this work, we extend on a previous study to perform a series of designs of protein binding interfaces by incorporating receptor flexibility using an ensemble of conformers collected from explicit-solvent molecular dynamics (MD) simulations. All designer complexes are subjected to 30 ns of MD in explicit solvent to assess for stability for a total of 480 ns of dynamics. This is followed by end-point free energy calculations whereby intermole- cular potential energy, polar and non-polar solvation energy and entropy of ligand and receptor are subtracted from that of the complex and averaged over 320 snapshots collected from each of the 30 ns MD simulations. Our initial effort consisted of redesigning the interface of three well-studied complexes, namely barnase-barstar, lysozyme-antibody D1.3 and trypsin-BPTI. The design was performed with flexible backbone approach. MD simulations revealed that all three complexes remained stable. Interestingly, the redesigned trypsin-BPTI complex was significantly more favorable than the native complex. This was attributed to the favorable electrostatics and entropy that complemented the already favorable non-polar component. Another aspect of this work consisted of grafting the surface of three proteins, namely tenascin, CheY and MBP1 to bind to barnase, trypsin and lysozyme. The process was initially performed using fixed backbone, and more than 300 ns of the explicit-solvent MD simulation revealed some of the complexes to dissociate over the course of the trajectories, whereas others remained stable. Free energy calculations confirmed that the non-polar component of the free energy as computed by summing the van der Waals energy and the non-polar solvation energy was a strong predictor of stability. Four complexes (two stable and two unstable) were selected, and redesigned using multiple conformers collected from the MD simulation. The resulting designer systems were then immersed in explicit solvent and 30 ns of MD was carried out on each. Interestingly, those complexes that were initially stable remained stable, whereas one of the unstable complexes became stable following redesign with flexible backbone. Free energy calculations showed significant improvements in the affinity for most complexes, revealing that the use of multiple conformers in protein design may significantly enhance such efforts.

KW - Flexible backbone design

KW - Free energy calculations

KW - Molecular dynamics simulations

KW - Protein engineering

KW - Protein-protein interaction

UR - http://www.scopus.com/inward/record.url?scp=75849162018&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=75849162018&partnerID=8YFLogxK

U2 - 10.1093/protein/gzp042

DO - 10.1093/protein/gzp042

M3 - Article

C2 - 19643976

AN - SCOPUS:75849162018

VL - 22

SP - 575

EP - 586

JO - Protein Engineering, Design and Selection

JF - Protein Engineering, Design and Selection

SN - 1741-0126

IS - 9

ER -