Transforming binding affinities from three dimensions to two with application to cadherin clustering

Membrane-bound receptors often form large assemblies resulting from binding to soluble ligands, cell-surface molecules on other cells and extracellular matrix proteins. For example, the association of membrane proteins with proteins on different cells (trans-interactions) can drive the oligomerization of proteins on the same cell (cis-interactions). A central problem in understanding the molecular basis of such phenomena is that equilibrium constants are generally measured in three-dimensional solution and are thus difficult to relate to the two-dimensional environment of a membrane surface. Here we present a theoretical treatment that converts three-dimensional affinities to two dimensions, accounting directly for the structure and dynamics of the membrane-bound molecules. Using a multiscale simulation approach, we apply the theory to explain the formation of ordered, junction-like clusters by classical cadherin adhesion proteins. The approach features atomic-scale molecular dynamics simulations to determine interdomain flexibility, Monte Carlo simulations of multidomain motion and lattice simulations of junction formation. A finding of general relevance is that changes in interdomain motion on trans-binding have a crucial role in driving the lateral, cis-, clustering of adhesion receptors.     

Calculation of electrostatic potentials in an enzyme active site

To be able to calculate the contributions of individual amino acids to the electrostatic field of a protein would be of considerable value in designing proteins of enhanced or altered function and stability. Recent studies on the serine protease subtilisin provide direct measurements of the electrostatic potential in the active site of the enzyme produced by two charged amino acids. We have used these results to test a recently developed method for the calculation of electrostatic interactions between two specific sites on a protein. The extent of agreement between the theoretical and experimental results suggests that the continuum solvent model used in the calculations reproduces the essential features of the interaction.