Glutamyl-isoleucine

Thermodynamic and kinetic studies of biomolecular interactions give insight into specificity of molecular recognition processes and advance rational drug design. Binding of phosphotyrosine (pY)-containing peptides to Src- and Grb2-SH2 domains was investigated using a surface plasmon resonance (SPR)-based method. This SPR assay yielded thermodynamic binding constants in solution, and the kinetic information contained in the SPR signal allowed kinetic analysis, which demonstrated distinct ways for pY ligands to interact with the SH2 domains. The results for binding to Src SH2 were consistent with sequestration of water molecules in the interface of the pYEEI peptide/Src SH2 complex. The results for a pYVNV peptide binding to Grb2 SH2 suggested a conformational change for Grb2 SH2 upon binding, which is not observed for Src SH2. Binding of a cyclic construct, allowing the pYVNV sequence in the bound conformation, did not have the expected entropy advantage. The results suggest an alternative binding mode for this construct, with the hydrophobic ring-closing part interacting with the protein. In all cases, except for full-length Grb2 protein, the affinity for the immobilized peptide at the SPR sensor and in solution was identical. This study demonstrates that SPR thermodynamic and kinetic analysis is a useful strategic tool in drug design.

Small molecules can be discovered or engineered to bind tightly to biologically relevant proteins, and these molecules have proven to be powerful tools for both basic research and therapeutic applications. In many cases, detailed biophysical analyses of the intermolecular binding events are essential for improving the activity of the small molecules. These interactions can often be characterized as straightforward bimolecular binding events, and a variety of experimental and analytical techniques have been developed and refined to facilitate these analyses. Several investigators have recently synthesized heterodimeric molecules that are designed to bind simultaneously with two different proteins to form ternary complexes. These heterodimeric molecules often display compelling biological activity; however, they are difficult to characterize. The bimolecular interaction between one protein and the heterodimeric ligand (primary dissociation constant) can be determined by a number of methods. However, the interaction between that protein-ligand complex and the second protein (secondary dissociation constant) is more difficult to measure due to the noncovalent nature of the original protein-ligand complex. Consequently, these heterodimeric compounds are often characterized in terms of their activity, which is an experimentally dependent metric. We have developed a general quantitative mathematical model that can be used to measure both the primary (protein + ligand) and secondary (protein-ligand + protein) dissociation constants for heterodimeric small molecules. These values are largely independent of the experimental technique used and furthermore provide a direct measure of the thermodynamic stability of the ternary complexes that are formed. Fluorescence polarization and this model were used to characterize the heterodimeric molecule, SLFpYEEI, which binds to both FKBP12 and the Fyn SH2 domain, demonstrating that the model is useful for both predictive as well as ex post facto analytical applications.

Although molecular dynamics simulation methods are useful in the modeling of macromolecular systems, they remain computationally expensive, with production work requiring costly high-performance computing (HPC) resources. We review recent innovations in accelerating molecular dynamics on graphics processing units (GPUs), and we describe GPUGRID, a volunteer computing project that uses the GPU resources of nondedicated desktop and workstation computers. In particular, we demonstrate the capability of simulating thousands of all-atom molecular trajectories generated at an average of 20 ns/day each (for systems of approximately 30 000-80 000 atoms). In conjunction with a potential of mean force (PMF) protocol for computing binding free energies, we demonstrate the use of GPUGRID in the computation of accurate binding affinities of the Src SH2 domain/pYEEI ligand complex by reconstructing the PMF over 373 umbrella sampling windows of 55 ns each (20.5 mus of total data). We obtain a standard free energy of binding of -8.7 +/- 0.4 kcal/mol within 0.7 kcal/mol from experimental results. This infrastructure will provide the basis for a robust system for high-throughput accurate binding affinity prediction.