Conformational analysis of endothelin-1: effects of solvation free energy

Determinants of ligand binding to cAMP-dependent protein kinase

Protein kinases are essential for the regulation of cellular growth and metabolism. Since their dysfunction leads to debilitating diseases, they represent key targets for pharmaceutical research. The rational design of kinase inhibitors requires an understanding of the determinants of ligand binding to these proteins. In the present study, a theoretical model based on continuum electrostatics and a surface-area-dependent nonpolar term is used to calculate binding affinities of balanol derivatives, H-series inhibitors, and ATP analogues toward the catalytic subunit of cAMP-dependent protein kinase (cAPK or protein kinase A). The calculations reproduce most of the experimental trends and provide insight into the driving forces responsible for binding. Nonpolar interactions are found to govern protein-ligand affinity. Hydrogen bonds represent a negligible contribution, because hydrogen bond formation in the complex requires the desolvation of the interacting partners. However, the binding affinity is decreased if hydrogen-bonding groups of the ligand remain unsatisfied in the complex. The disposition of hydrogen-bonding groups in the ligand is therefore crucial for binding specificity. These observations should be valuable guides in the design of potent and specific kinase inhibitors.

Conformational analysis of biologically active truncated linear analogs of endothelin-1 using NMR and molecular modeling

Some linear truncated analogs of endothelin-1 display potent agonistic activity at the ET(B) receptor, especially when the side chain of Trp21 is N-formylated. Then, the three-dimensional arrangements of six structurally reduced linear analogs, three formylated and three nonformylated, have been investigated by high resolution NMR spectroscopy and molecular modeling, in order to pinpoint the conformational features related to the biological activity. Two-dimensional double-quantum-filtered correlation spectroscopy (DQFCOSY), total correlation spectroscopy (TOCSY) and nuclear Overhauser enhancement spectroscopy (NOESY) were recorded and analyzed for each molecule. Interspatial distance constraints were derived from the intensity of the NOESY connectivities. The formation of hydrogen bonding was monitored from the temperature dependence of the NH chemical shifts. Molecular models calculated by means of distance geometry, simulated annealing and energy minimization, using the NMR constraints, strongly suggested a global elongated structure for the formylated analogs exhibiting biological activity, and a folded arrangement for the unformylated derivatives. Homology comparisons allowed the identification of a beta-turn-like folding of the C-terminal segment Asp18-Trp21 as a probable key-factor for activity.

Conformational analysis of the highly potent bradykinin antagonist Hoe-140 by means of two different computational methods

The AMBER 4.0 force field was used to perform the characterization of the conformational profile of the highly potent bradykinin antagonist Hoe-140 (D-Arg0-Arg1-Pro2-Hyp3-Gly4-Thi5-Ser6-D-++ +Tic7-Oic8-Arg9). The structural features of the peptide were assessed using two different computational methods, both capable to provide a good sampling of the low-energy conformations of the molecule. Specifically, the conformational space of the peptide was explored: i) computing molecular dynamics trajectories in cycles of high (900 K) and low (300 K) temperature and ii) using simulated annealing (SA) in an iterative fashion. Analysis of the structures characterized indicates that most of the low-energy conformations of the peptide exhibit a betaII'-turn motif at its C-terminus, in agreement with previous experimental and theoretical studies. On the other hand, about a 50% of the low-energy conformations characterized also exhibit different beta-turn type motifs at the N-terminus, whereas the rest of the conformations can be described as bends. Finally, in order to get new insights into the structural requirements necessary to design more potent and selective antagonists of bradykinin, present results were compared with those previously reported by this laboratory on the conformational preferences of the native nonapeptide and its DPhe7 analog.

Temperature dependence of histidine ionization constants in myoglobin

The standard enthalpy of ionization of six titratable histidines in horse metaquomyoglobin was determined by repeating proton NMR titrations as a function of temperature and using the van't Hoff relationship. It was found that deltaH degrees varies between 16 and 37 kJ mol(-1) in the protein, compared with a value of 29 kJ mol(-1) in free histidine. The standard entropy change was evaluated by combining the enthalpy and free energy changes derived from the pKa values. Although the entropy change could not be precisely and accurately obtained by this method, it could be established that it spans a wide range, from -60 to 0 J K(-1) mol(-1), about the value of -23 J K(-1) mol(-1) for the free histidine. The entropy change was used within the framework of enthalpy-entropy compensation to partition the solvation component from the standard thermodynamic quantities for each of the titrating residues. It was shown that the partitioning of the values in the protein is not readily understood in terms of solvent accessibility or electrostatic interactions. The contribution of solvation effects to the temperature response appeared to be significant only in the case of His-119 and His-48. The standard quantities were also used to explore the energetics of proton binding in the native state at temperatures below the onset of thermal denaturation.

Kinase conformations: a computational study of the effect of ligand binding

Protein function is often controlled by ligand-induced conformational transitions. Yet, in spite of the increasing number of three-dimensional crystal structures of proteins in different conformations, not much is known about the driving forces of these transitions. As an initial step toward exploring the conformational and energetic landscape of protein kinases by computational methods, intramolecular energies and hydration free energies were calculated for different conformations of the catalytic domain of cAMP-dependent protein kinase (cAPK) with a continuum (Poisson) model for the electrostatics. Three protein kinase crystal structures for ternary complexes of cAPK with the peptide inhibitor PKI(5-24) and ATP or AMP-PNP were modeled into idealized intermediate and open conformations. Concordant with experimental observation, we find that the binding of PKI(5-24) is more effective in stabilizing the closed and intermediate forms of cAPK than ATP. PKI(5-24) seems to drive the final closure of the active site cleft from intermediate to closed state because ATP does not distinguish between these two states. Binding of PKI(5-24) and ATP is energetically additive.

Geometric versus topological clustering: an insight into conformation mapping

Clustering molecular conformations into "families" is a common procedure in conformational analysis of molecular systems. An implicit assumption which often underlies this clustering approach is that the resulting geometric families reflect the energetic structure of the system's potential energy surface. In a broader context we address the question whether structural similarity is correlated with energy basins, i.e., whether conformations that belong to the same energy basin are also geometrically similar. 'Topological mapping' and principal coordinate projections are used here to address this question and to assess the quality of the 'family clustering' procedure. Applying the analysis to a small tetrapeptide it was found that the general correlation that exists between energy basins and structural similarity is not absolute. Clusters generated by the geometric 'family clustering' procedure do not always reflect the underlying energy basins. In particular it was found that the 'family tree' that is generated by the 'family clustering' procedure is completely inconsistent with its real topological counterpart, the 'disconnectivity' graph of this system. It is also demonstrated that principal coordinate analysis is a powerful visualization technique which, at least for this system, works better when distances are measured in dihedral angle space rather than in cartesian space.