Protein-protein interactions from linear-scaling first-principles quantum-mechanical calculations

¹ Theory of Condensed Matter Group, Cavendish Laboratory, University of Cambridge - J. J. Thomson Avenue, Cambridge CB3 0HE, UK, EU
² School of Chemistry, University of Southampton - Highfield, Southampton SO17 1BJ, UK, EU
³ MRC Cancer Cell Unit Hutchison/MRC Research Centre - Hills Road, Cambridge CB2 0XZ, UK, EU

^a djc56@cam.ac.uk

Received: 30 March 2010
Accepted: 23 July 2010

Abstract

A modification of the MM-PBSA technique for calculating binding affinities of biomolecular complexes is presented. Classical molecular dynamics is used to explore the motion of the extended interface between two peptides derived from the BRC4 repeat of BRCA2 and the eukaryotic recombinase RAD51. The resulting trajectory is sampled using the linear-scaling density functional theory code, onetep, to determine from first principles, and with high computational efficiency, the relative free energies of binding of the ∼2800 atom receptor-ligand complexes. This new method provides the basis for computational interrogation of protein-protein and protein-ligand interactions within fields ranging from chemical biological studies to small-molecule binding behaviour, with both unprecedented chemical accuracy and affordable computational expense.