An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis

W. Wang, W. A. Lim, A. Jakalian, J. Wang, J. Wang, R. Luo, C. I. Bayly, P. A. Kollman

Research output: Contribution to journalArticle

114 Citations (Scopus)

Abstract

The Src-homology-3 (SH3) domain of the Caenorhabditis elegans protein Sem-5 binds prolinerich sequences. It is reported that the SH3 domains broadly accept amide N-substituted residues instead of only recognizing prolines on the basis of side chain shape or rigidity. We have studied the interactions between Sem-5 and its ligands using molecular dynamics (MD), free energy calculations, and sequence analysis. Relative binding free energies, estimated by a method called MM/PBSA, between different substitutions at sites -1, 0, and +2 of the peptide are consistent with the experimental data. A new method to calculate atomic partial charges, AM1-BCC method, is also used in the binding free energy calculations for different N-substitutions at site - 1. The results are very similar to those obtained from widely used RESP charges in the AMBER force field. AM1-BCC charges can be calculated more rapidly for any organic molecule than can the RESP charges. Therefore, their use can enable a broader and more efficient application of the MM/PBSA method in drug design. Examination of each component of the free energy leads to the construction of van der Waals interaction energy profiles for each ligand as well as for wild-type and mutant Sem-5 proteins. The profiles and free energy calculations indicate that the van der Waals interactions between the ligands and the receptor determine whether an N- or a Cα-substituted residue is favored at each site. A VC value (defined as a product of the conservation percentage of each residue and its van der Waals interaction energy with the ligand) is used to identify several residues on the receptor that are critical for specificity and binding affinity. This VC value may have a potential use in identifying crucial residues for any ligand-protein or protein-protein system. Mutations at two of those crucial residues, N190 and N206, are examined. One mutation, N190I, is predicted to reduce the selectivity of the N-substituted residue at site -1 of the ligand and is shown to bind similarly with N- and Cα-substituted residues at that site.

Original languageEnglish (US)
Pages (from-to)3986-3994
Number of pages9
JournalJournal of the American Chemical Society
Volume123
Issue number17
DOIs
StatePublished - 2001

Fingerprint

src Homology Domains
Molecular Dynamics Simulation
Free energy
Sequence Analysis
Molecular dynamics
Ligands
Proteins
Substitution reactions
Caenorhabditis elegans Proteins
Mutation
Drug Design
Amides
Proline
Rigidity
Peptides
Conservation
Molecules
Pharmaceutical Preparations

ASJC Scopus subject areas

  • Chemistry(all)

Cite this

An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis. / Wang, W.; Lim, W. A.; Jakalian, A.; Wang, J.; Wang, J.; Luo, R.; Bayly, C. I.; Kollman, P. A.

In: Journal of the American Chemical Society, Vol. 123, No. 17, 2001, p. 3986-3994.

Research output: Contribution to journalArticle

Wang, W. ; Lim, W. A. ; Jakalian, A. ; Wang, J. ; Wang, J. ; Luo, R. ; Bayly, C. I. ; Kollman, P. A. / An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis. In: Journal of the American Chemical Society. 2001 ; Vol. 123, No. 17. pp. 3986-3994.
@article{49d85d0b3baf49609cd0c389aaf74cd2,
title = "An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis",
abstract = "The Src-homology-3 (SH3) domain of the Caenorhabditis elegans protein Sem-5 binds prolinerich sequences. It is reported that the SH3 domains broadly accept amide N-substituted residues instead of only recognizing prolines on the basis of side chain shape or rigidity. We have studied the interactions between Sem-5 and its ligands using molecular dynamics (MD), free energy calculations, and sequence analysis. Relative binding free energies, estimated by a method called MM/PBSA, between different substitutions at sites -1, 0, and +2 of the peptide are consistent with the experimental data. A new method to calculate atomic partial charges, AM1-BCC method, is also used in the binding free energy calculations for different N-substitutions at site - 1. The results are very similar to those obtained from widely used RESP charges in the AMBER force field. AM1-BCC charges can be calculated more rapidly for any organic molecule than can the RESP charges. Therefore, their use can enable a broader and more efficient application of the MM/PBSA method in drug design. Examination of each component of the free energy leads to the construction of van der Waals interaction energy profiles for each ligand as well as for wild-type and mutant Sem-5 proteins. The profiles and free energy calculations indicate that the van der Waals interactions between the ligands and the receptor determine whether an N- or a Cα-substituted residue is favored at each site. A VC value (defined as a product of the conservation percentage of each residue and its van der Waals interaction energy with the ligand) is used to identify several residues on the receptor that are critical for specificity and binding affinity. This VC value may have a potential use in identifying crucial residues for any ligand-protein or protein-protein system. Mutations at two of those crucial residues, N190 and N206, are examined. One mutation, N190I, is predicted to reduce the selectivity of the N-substituted residue at site -1 of the ligand and is shown to bind similarly with N- and Cα-substituted residues at that site.",
author = "W. Wang and Lim, {W. A.} and A. Jakalian and J. Wang and J. Wang and R. Luo and Bayly, {C. I.} and Kollman, {P. A.}",
year = "2001",
doi = "10.1021/ja003164o",
language = "English (US)",
volume = "123",
pages = "3986--3994",
journal = "Journal of the American Chemical Society",
issn = "0002-7863",
publisher = "American Chemical Society",
number = "17",

}

TY - JOUR

T1 - An analysis of the interactions between the Sem-5 SH3 domain and its ligands using molecular dynamics, free energy calculations, and sequence analysis

AU - Wang, W.

AU - Lim, W. A.

AU - Jakalian, A.

AU - Wang, J.

AU - Wang, J.

AU - Luo, R.

AU - Bayly, C. I.

AU - Kollman, P. A.

PY - 2001

Y1 - 2001

N2 - The Src-homology-3 (SH3) domain of the Caenorhabditis elegans protein Sem-5 binds prolinerich sequences. It is reported that the SH3 domains broadly accept amide N-substituted residues instead of only recognizing prolines on the basis of side chain shape or rigidity. We have studied the interactions between Sem-5 and its ligands using molecular dynamics (MD), free energy calculations, and sequence analysis. Relative binding free energies, estimated by a method called MM/PBSA, between different substitutions at sites -1, 0, and +2 of the peptide are consistent with the experimental data. A new method to calculate atomic partial charges, AM1-BCC method, is also used in the binding free energy calculations for different N-substitutions at site - 1. The results are very similar to those obtained from widely used RESP charges in the AMBER force field. AM1-BCC charges can be calculated more rapidly for any organic molecule than can the RESP charges. Therefore, their use can enable a broader and more efficient application of the MM/PBSA method in drug design. Examination of each component of the free energy leads to the construction of van der Waals interaction energy profiles for each ligand as well as for wild-type and mutant Sem-5 proteins. The profiles and free energy calculations indicate that the van der Waals interactions between the ligands and the receptor determine whether an N- or a Cα-substituted residue is favored at each site. A VC value (defined as a product of the conservation percentage of each residue and its van der Waals interaction energy with the ligand) is used to identify several residues on the receptor that are critical for specificity and binding affinity. This VC value may have a potential use in identifying crucial residues for any ligand-protein or protein-protein system. Mutations at two of those crucial residues, N190 and N206, are examined. One mutation, N190I, is predicted to reduce the selectivity of the N-substituted residue at site -1 of the ligand and is shown to bind similarly with N- and Cα-substituted residues at that site.

AB - The Src-homology-3 (SH3) domain of the Caenorhabditis elegans protein Sem-5 binds prolinerich sequences. It is reported that the SH3 domains broadly accept amide N-substituted residues instead of only recognizing prolines on the basis of side chain shape or rigidity. We have studied the interactions between Sem-5 and its ligands using molecular dynamics (MD), free energy calculations, and sequence analysis. Relative binding free energies, estimated by a method called MM/PBSA, between different substitutions at sites -1, 0, and +2 of the peptide are consistent with the experimental data. A new method to calculate atomic partial charges, AM1-BCC method, is also used in the binding free energy calculations for different N-substitutions at site - 1. The results are very similar to those obtained from widely used RESP charges in the AMBER force field. AM1-BCC charges can be calculated more rapidly for any organic molecule than can the RESP charges. Therefore, their use can enable a broader and more efficient application of the MM/PBSA method in drug design. Examination of each component of the free energy leads to the construction of van der Waals interaction energy profiles for each ligand as well as for wild-type and mutant Sem-5 proteins. The profiles and free energy calculations indicate that the van der Waals interactions between the ligands and the receptor determine whether an N- or a Cα-substituted residue is favored at each site. A VC value (defined as a product of the conservation percentage of each residue and its van der Waals interaction energy with the ligand) is used to identify several residues on the receptor that are critical for specificity and binding affinity. This VC value may have a potential use in identifying crucial residues for any ligand-protein or protein-protein system. Mutations at two of those crucial residues, N190 and N206, are examined. One mutation, N190I, is predicted to reduce the selectivity of the N-substituted residue at site -1 of the ligand and is shown to bind similarly with N- and Cα-substituted residues at that site.

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

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

U2 - 10.1021/ja003164o

DO - 10.1021/ja003164o

M3 - Article

C2 - 11457149

AN - SCOPUS:0034789590

VL - 123

SP - 3986

EP - 3994

JO - Journal of the American Chemical Society

JF - Journal of the American Chemical Society

SN - 0002-7863

IS - 17

ER -