ABSTRACT
In order to investigate conformational preferences of the 21-residue peptide hormone endothelin-1 (ET-1), an extensive conformational search was carried out in vacuo using a combination of high temperature molecular dynamics/annealing and a Monte Carlo/minimization search in torsion angle space. Fully minimized conformations from the search were grouped into families using a clustering technique based on rms fitting over the Cartesian coordinates of the atoms of the peptide backbone of the ring region. A wide range of local energy mining were identified even though two disulfide bridges (Cys1-Cys15 and Cys3-Cys11) constrain the structure of the peptide. Low energy conformers of ET-1 as a nonionized species in vacuo are stabilized by intramolecular interaction of the ring region (residues 1-15) with the tail (residues 16-21). Strained conformations for individual residues are observed. Conformational similarity to protein loops is established by matching to protein crystal structures. In order to assess the influence of aqueous environment on conformational preference, the electrostatic contribution to the solvation energy was calculated for ET-1 as a fully ionized species (Asp8, Lys9, Glu10, Asp18, N- and C-terminus) using a continuum electrostatics model (DelPhi) for each of the conformers generated in vacuo, and the total solvation free energy was estimated by adding a hydrophobic contribution proportional to solvent accessible surface area. Solvation dramatically alters the relative energetics of ET-1 conformers from that calculated in vacuo. Conformers of ET-1 favored by the electrostatic solvation energy in water include conformers with helical secondary structure in the region of residues 9-15. Perhaps of most importance, it was demonstrated that the contribution to solvation by an individual charge depends not only on its solvent accessibility but on the proximity of other charges, i.e., it is a cooperative effect. This was shown by the calculation of electrostatic solvation energy as a function of conformation with individual charges systematically turned "on" and "off". The cooperative effect of multiple charges on solvation demonstrated in this manner calls into question models that relate solvation energy simply to solvent accessibility by atom or residue alone.
Subject(s)
Endothelins/chemistry , Protein Conformation , Amino Acid Sequence , Calorimetry , Least-Squares Analysis , Models, Molecular , Molecular Sequence Data , Monte Carlo Method , Protein Structure, Secondary , Software , Solubility , ThermodynamicsABSTRACT
The binding of the recombinant proteinase inhibitor eglin c from the leech Hirudo medicinalis to serine (pro)enzymes belonging to the chymotrypsin and subtilisin families has been investigated from the thermodynamic viewpoint, between pH 4.5 and 9.5 and from 10 degrees C to 40 degrees C. The affinity of eglin c for the serine (pro)enzymes considered shows the following trend: Leu-proteinase [the leucine specific serine proteinase from spinach (Spinacia oleracea L.) leaves] greater than human leucocyte elastase congruent to human cathepsin G congruent to subtilisin Carlsberg congruent to bovine alpha-chymotrypsin greater than bovine alpha-chymotrypsinogen A congruent to porcine pancreatic elastase congruent to bovine beta-trypsin. The serine (pro)enzyme-inhibitor complex formation is an entropy-driven process. On increasing the pH from 4.5 to 9.5, the affinity of eglin c for the serine (pro)enzymes considered increases thus reflecting the acid pK shift of the invariant hystidyl catalytic residue from approximately to 6.9 in the free serine proteinases and bovine alpha-chymotrypsinogen A to congruent to 5.1 in the serine (pro)enzyme-inhibitor complexes. Considering the known molecular models, the observed binding behaviour of eglin c was related to the inferred stereochemistry of the serine (pro)enzyme-inhibitor contact regions.
