Theoretical prediction of the binding free energy for mutants of replication protein A

Claudio Carra, Janapriya Saha, Francis A. Cucinotta

Research output: Contribution to journalArticle

5 Citations (Scopus)

Abstract

The replication protein A (RPA) is a heterotrimeric (70, 32, and 14 kDa subunits), single stranded DNA (ssDNA) binding protein required for pivotal functions in the cell metabolism, such as chromosomal replication, prevention of hairpin formation, DNA repair and recombination, and signaling after DNA damage. Studies based on deletions and mutations have identified the high affinity ssDNA binding domains in the 70 kDa subunit of RPA, regions A and B. Individually, the domain A and B have a low affinity for ssDNA, while tandems composed of AA, AB, BB, and BA sequences bind the ssDNA with moderate to high affinity. Single and double point mutations on polar residues in the binding domains leads to a reduction in affinity of RPA for ssDNA, in particular when two hydrophilic residues are involved. In view of these results, we performed a study based on molecular dynamics simulation aimed to reproduce the experimental change in binding free energy, ΔΔG, of RPA70 mutants to further elucidate the nature of the protein-ssDNA interaction. The MM-PB(GB) SA methods implemented in Amber10 and the code FoldX were used to estimate the binding free energy. The theoretical and experimental ΔΔG values correlate better when the results are obtained by MM-PBSA calculated on individual trajectories for each mutant. In these conditions, the correlation coefficient between experimental and theoretical ΔΔG reaches a value of 0.95 despite the overestimation of the energy change by one order of magnitude. The decomposition of the MM-GBSA energy per residue allows us to correlate the change of the affinity with the residue polarity and energy contribution to the binding. The method revealed reliable predictions of the change in the affinity in function of mutations, and can be used to identify new mutants with distinct binding properties.

Original languageEnglish (US)
Pages (from-to)3035-3049
Number of pages15
JournalJournal of Molecular Modeling
Volume18
Issue number7
DOIs
StatePublished - Jul 1 2012

Fingerprint

Replication Protein A
Single-Stranded DNA
Free energy
DNA
deoxyribonucleic acid
free energy
proteins
Proteins
affinity
predictions
mutations
DNA-Binding Proteins
Metabolism
Molecular dynamics
Repair
deletion
Trajectories
Decomposition
metabolism
correlation coefficients

Keywords

  • Amber
  • Binding
  • MM-PBSA
  • Molecular dynamics
  • Mutation
  • Replication protein
  • RPA

ASJC Scopus subject areas

  • Catalysis
  • Computer Science Applications
  • Physical and Theoretical Chemistry
  • Organic Chemistry
  • Computational Theory and Mathematics
  • Inorganic Chemistry

Cite this

Theoretical prediction of the binding free energy for mutants of replication protein A. / Carra, Claudio; Saha, Janapriya; Cucinotta, Francis A.

In: Journal of Molecular Modeling, Vol. 18, No. 7, 01.07.2012, p. 3035-3049.

Research output: Contribution to journalArticle

@article{ded385f3aa3d420d8b0142831f959afa,
title = "Theoretical prediction of the binding free energy for mutants of replication protein A",
abstract = "The replication protein A (RPA) is a heterotrimeric (70, 32, and 14 kDa subunits), single stranded DNA (ssDNA) binding protein required for pivotal functions in the cell metabolism, such as chromosomal replication, prevention of hairpin formation, DNA repair and recombination, and signaling after DNA damage. Studies based on deletions and mutations have identified the high affinity ssDNA binding domains in the 70 kDa subunit of RPA, regions A and B. Individually, the domain A and B have a low affinity for ssDNA, while tandems composed of AA, AB, BB, and BA sequences bind the ssDNA with moderate to high affinity. Single and double point mutations on polar residues in the binding domains leads to a reduction in affinity of RPA for ssDNA, in particular when two hydrophilic residues are involved. In view of these results, we performed a study based on molecular dynamics simulation aimed to reproduce the experimental change in binding free energy, ΔΔG, of RPA70 mutants to further elucidate the nature of the protein-ssDNA interaction. The MM-PB(GB) SA methods implemented in Amber10 and the code FoldX were used to estimate the binding free energy. The theoretical and experimental ΔΔG values correlate better when the results are obtained by MM-PBSA calculated on individual trajectories for each mutant. In these conditions, the correlation coefficient between experimental and theoretical ΔΔG reaches a value of 0.95 despite the overestimation of the energy change by one order of magnitude. The decomposition of the MM-GBSA energy per residue allows us to correlate the change of the affinity with the residue polarity and energy contribution to the binding. The method revealed reliable predictions of the change in the affinity in function of mutations, and can be used to identify new mutants with distinct binding properties.",
keywords = "Amber, Binding, MM-PBSA, Molecular dynamics, Mutation, Replication protein, RPA",
author = "Claudio Carra and Janapriya Saha and Cucinotta, {Francis A.}",
year = "2012",
month = "7",
day = "1",
doi = "10.1007/s00894-011-1313-z",
language = "English (US)",
volume = "18",
pages = "3035--3049",
journal = "Journal of Molecular Modeling",
issn = "1610-2940",
publisher = "Springer Verlag",
number = "7",

}

TY - JOUR

T1 - Theoretical prediction of the binding free energy for mutants of replication protein A

AU - Carra, Claudio

AU - Saha, Janapriya

AU - Cucinotta, Francis A.

PY - 2012/7/1

Y1 - 2012/7/1

N2 - The replication protein A (RPA) is a heterotrimeric (70, 32, and 14 kDa subunits), single stranded DNA (ssDNA) binding protein required for pivotal functions in the cell metabolism, such as chromosomal replication, prevention of hairpin formation, DNA repair and recombination, and signaling after DNA damage. Studies based on deletions and mutations have identified the high affinity ssDNA binding domains in the 70 kDa subunit of RPA, regions A and B. Individually, the domain A and B have a low affinity for ssDNA, while tandems composed of AA, AB, BB, and BA sequences bind the ssDNA with moderate to high affinity. Single and double point mutations on polar residues in the binding domains leads to a reduction in affinity of RPA for ssDNA, in particular when two hydrophilic residues are involved. In view of these results, we performed a study based on molecular dynamics simulation aimed to reproduce the experimental change in binding free energy, ΔΔG, of RPA70 mutants to further elucidate the nature of the protein-ssDNA interaction. The MM-PB(GB) SA methods implemented in Amber10 and the code FoldX were used to estimate the binding free energy. The theoretical and experimental ΔΔG values correlate better when the results are obtained by MM-PBSA calculated on individual trajectories for each mutant. In these conditions, the correlation coefficient between experimental and theoretical ΔΔG reaches a value of 0.95 despite the overestimation of the energy change by one order of magnitude. The decomposition of the MM-GBSA energy per residue allows us to correlate the change of the affinity with the residue polarity and energy contribution to the binding. The method revealed reliable predictions of the change in the affinity in function of mutations, and can be used to identify new mutants with distinct binding properties.

AB - The replication protein A (RPA) is a heterotrimeric (70, 32, and 14 kDa subunits), single stranded DNA (ssDNA) binding protein required for pivotal functions in the cell metabolism, such as chromosomal replication, prevention of hairpin formation, DNA repair and recombination, and signaling after DNA damage. Studies based on deletions and mutations have identified the high affinity ssDNA binding domains in the 70 kDa subunit of RPA, regions A and B. Individually, the domain A and B have a low affinity for ssDNA, while tandems composed of AA, AB, BB, and BA sequences bind the ssDNA with moderate to high affinity. Single and double point mutations on polar residues in the binding domains leads to a reduction in affinity of RPA for ssDNA, in particular when two hydrophilic residues are involved. In view of these results, we performed a study based on molecular dynamics simulation aimed to reproduce the experimental change in binding free energy, ΔΔG, of RPA70 mutants to further elucidate the nature of the protein-ssDNA interaction. The MM-PB(GB) SA methods implemented in Amber10 and the code FoldX were used to estimate the binding free energy. The theoretical and experimental ΔΔG values correlate better when the results are obtained by MM-PBSA calculated on individual trajectories for each mutant. In these conditions, the correlation coefficient between experimental and theoretical ΔΔG reaches a value of 0.95 despite the overestimation of the energy change by one order of magnitude. The decomposition of the MM-GBSA energy per residue allows us to correlate the change of the affinity with the residue polarity and energy contribution to the binding. The method revealed reliable predictions of the change in the affinity in function of mutations, and can be used to identify new mutants with distinct binding properties.

KW - Amber

KW - Binding

KW - MM-PBSA

KW - Molecular dynamics

KW - Mutation

KW - Replication protein

KW - RPA

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

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

U2 - 10.1007/s00894-011-1313-z

DO - 10.1007/s00894-011-1313-z

M3 - Article

VL - 18

SP - 3035

EP - 3049

JO - Journal of Molecular Modeling

JF - Journal of Molecular Modeling

SN - 1610-2940

IS - 7

ER -