A potential smoothing algorithm accurately predicts transmembrane helix packing

Potential smoothing, a deterministic analog of stochastic simulated annealing, is a powerful paradigm for the solution of conformational search problems that require extensive sampling, and should be a useful tool in computational approaches to structure prediction and refinement. A novel potential smoothing and search (PSS) algorithm has been developed and applied to predict the packing of transmembrane helices. The highlight of this method is the efficient manner in which it circumvents the combinatorial explosion associated with the large number of minima on multidimensional potential energy surfaces in order to converge to the global energy minimum. Here we show how our potential smoothing and search method succeeds in finding the global minimum energy structure for the glycophorin A (GpA) transmembrane helix dimer by optimizing interhelical van der Waals interactions over rigid and semi–rigid helices. Structures obtained from our ab initio predictions are in close agreement with recent experimental data.

[1]  Ron Elber.  Recent Developments In Theoretical Studies Of Proteins , 1996 .

[2]  Jay W. Ponder,et al.  Analysis and Application of Potential Energy Smoothing and Search Methods for Global Optimization , 1998 .

[3]  John E. Straub,et al.  OPTIMIZATION TECHNIQUES WITH APPLICATIONS TO PROTEINS , 1996 .

[4]  D. Langosch,et al.  The dimerization motif of the glycophorin A transmembrane segment in membranes: Importance of glycine residues , 1998, Protein science : a publication of the Protein Society.

[5]  C. Chothia,et al.  Helix to helix packing in proteins. , 1981, Journal of molecular biology.

[6]  H. Scheraga,et al.  On the multiple-minima problem in the conformational analysis of molecules: deformation of the potential energy hypersurface by the diffusion equation method , 1989 .

[7]  Shugo Nakamura,et al.  CONFORMATIONAL ENERGY MINIMIZATION USING A TWO-STAGE METHOD , 1995 .

[8]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[9]  D. Engelman,et al.  Improved prediction for the structure of the dimeric transmembrane domain of glycophorin A obtained through global searching , 1996, Proteins.

[10]  D. Engelman,et al.  Specificity and promiscuity in membrane helix interactions , 1994, Quarterly Reviews of Biophysics.

[11]  Harold A. Scheraga,et al.  Some approaches to the multiple-minima problem in protein folding , 1995, Global Minimization of Nonconvex Energy Functions: Molecular Conformation and Protein Folding.

[12]  James H. Prestegard,et al.  A Transmembrane Helix Dimer: Structure and Implications , 1997, Science.

[13]  H. Scheraga,et al.  Application of the diffusion equation method for global optimization to oligopeptides , 1992 .

[14]  J. Straub,et al.  Global energy minimum searches using an approximate solution of the imaginary time Schroedinger equation , 1993 .

[15]  P. Kraulis.  A program to produce both detailed and schematic plots of protein structures , 1991 .

[16]  D. Engelman,et al.  Membrane protein folding and oligomerization: the two-stage model. , 1990, Biochemistry.

[17]  Roland L. Dunbrack,et al.  Backbone-dependent rotamer library for proteins. Application to side-chain prediction. , 1993, Journal of molecular biology.

[18]  C. Deber,et al.  Guidelines for membrane protein engineering derived from de novo designed model peptides. , 1998, Biopolymers.