This Article Full Text Full Text (PDF) All Versions of this Article: biophysj.106.081224v1 91/5/1710    most recent Alert me when this article is cited Alert me if a correction is posted Services Similar articles in this journal Similar articles in PubMed Alert me to new issues of the journal Download to citation manager Citing Articles Citing Articles via Google Scholar Google Scholar Articles by Zhang, J. Articles by Lazaridis, T. Search for Related Content PubMed PubMed Citation Articles by Zhang, J. Articles by Lazaridis, T.

Calculating the Free Energy of Association of Transmembrane Helices

Department of Chemistry, The City College of New York/The City University of New York, New York, New York 10031

Correspondence: Address reprint requests to Themis Lazaridis, Dept. of Chemistry, The City College of New York, Convent Avenue at 138th St., New York, NY 10031. Tel.: 212-650-8364; Fax: 212-650-6107; E-mail: tlazaridis{at}sci.ccny.cuny.edu.

A large number of experimental studies have been devoted to quantifying the interaction between transmembrane (TM) helices in detergent micelles and, more recently, in bilayers. Theoretical calculation of association free energy of TM helices would be useful for predicting the propensity of given sequences to oligomerize and for understanding the difference between association in micelles and in bilayers. In this article, the theoretical foundation for calculating the standard association free energy of TM helices is laid out and is applied to glycophorin A in both micelles and bilayers. The standard association free energy is decomposed into the effective energy, translational, rotational, and conformational entropy terms. The effective energy of association is obtained by molecular dynamics simulations in an implicit membrane model. The translational and rotational entropy of association is calculated from the probability distribution of the translational and rotational degrees of freedom obtained from the molecular dynamics simulations. The side-chain conformational entropy of association is estimated from the probability distribution obtained by rigid rotation of all side-chain dihedral angles. The calculated standard association free energy of glycophorin A in N-dodecylphosphocholine micelles is in good agreement with the experimental value. The translational entropy cost is larger, whereas the rotational entropy cost is smaller in bilayers than in micelles. The standard association free energy in 1,2-dimyristoyl-sn-glycero-3-phosphocholine bilayers is calculated to be 1.3 kcal/mol more favorable than in N-dodecylphosphocholine micelles, consistent with available experimental data.