The traditional approaches for treating oligopeptides and proteins in molecular mechanics and molecular dynamics assume that the charges are invariant to conformational changes. These charges are fixed at the beginning and are kept constant throughout the whole simulation. Attempts to use these traditional approaches to explain some experimental dependences of oligopeptide conformation on pH have not been very successful. Our goal was focused on a specific aspect of the folding processes: It attempts to understand how the enviroment affects the conformational preference of a short polypeptide with the sequence: Ac-ETGTKAELLAKYEATHK-NHMe, and vice versa; i.e. how the conformation affects the pK's of ionizable groups. In our work a new approach is presented. It considers explicitly the coupling between the conformation of the molecule and the ionization equilibria at a given pH value. Calculations of the solvation free energy and free energy of ionization of this 17-residue (for which the available NMR and CD experimental data indicate that the conformations containing a right-handed alpha-helix segment at low pH are energetically more favorable) were carried out using a fast multigrid boundary element method (MBE). The MBE method was integrated into ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm to compute the coupling between the ionization state and the conformation of the molecule. The results of our calculations using the present method agree quite well with experiments, in contrast with previous applications with standard techniques (using pre-assigned charges at each pH). In particular we have been able to show how the coupling to the conformation leads to different degree of ionization of a given type of residue, for example glutamic acid, at different positions in the amino acid sequence, at a given pH. The results of this study provide a sound basis to discuss the origin of the stability of oligopeptide conformations and its dependence on enviromental conditions. We expect to continue with the development of this approach and further testing of a series of oligopeptides of biological interest.

Agency
National Institute of Health (NIH)
Institute
National Center for Research Resources (NCRR)
Type
Biotechnology Resource Grants (P41)
Project #
5P41RR004293-07
Application #
6282178
Study Section
Project Start
1997-12-01
Project End
1998-11-30
Budget Start
1997-10-01
Budget End
1998-09-30
Support Year
7
Fiscal Year
1998
Total Cost
Indirect Cost
Name
Cornell University
Department
Type
DUNS #
City
Ithaca
State
NY
Country
United States
Zip Code
14850
Chiang, Chi-Tung; Shores, Kevin S; Freindorf, Marek et al. (2008) Size-restricted proton transfer within toluene-methanol cluster ions. J Phys Chem A 112:11559-65
Kazmierkiewicz, Rajmund; Liwo, Adam; Scheraga, Harold A (2003) Addition of side chains to a known backbone with defined side-chain centroids. Biophys Chem 100:261-80
Kazmierkiewicz, Rajmund; Liwo, Adam; Scheraga, Harold A (2002) Energy-based reconstruction of a protein backbone from its alpha-carbon trace by a Monte-Carlo method. J Comput Chem 23:715-23
Liwo, Adam; Arlukowicz, Piotr; Czaplewski, Cezary et al. (2002) A method for optimizing potential-energy functions by a hierarchical design of the potential-energy landscape: application to the UNRES force field. Proc Natl Acad Sci U S A 99:1937-42
Scheraga, Harold A; Pillardy, Jaroslaw; Liwo, Adam et al. (2002) Evolution of physics-based methodology for exploring the conformational energy landscape of proteins. J Comput Chem 23:28-34
Scheraga, Harold A; Vila, Jorge A; Ripoll, Daniel R (2002) Helix-coil transitions re-visited. Biophys Chem 101-102:255-65
Pillardy, J; Arnautova, Y A; Czaplewski, C et al. (2001) Conformation-family Monte Carlo: a new method for crystal structure prediction. Proc Natl Acad Sci U S A 98:12351-6
Vila, J A; Ripoll, D R; Scheraga, H A (2001) Influence of lysine content and pH on the stability of alanine-based copolypeptides. Biopolymers 58:235-46
Pillardy, J; Czaplewski, C; Liwo, A et al. (2001) Recent improvements in prediction of protein structure by global optimization of a potential energy function. Proc Natl Acad Sci U S A 98:2329-33
Czaplewski, C; Rodziewicz-Motowidlo, S; Liwo, A et al. (2000) Molecular simulation study of cooperativity in hydrophobic association. Protein Sci 9:1235-45

Showing the most recent 10 out of 20 publications