Conformational analysis of cysteine-containing peptides and prospects of their application to 213Bi chelating in antitumor therapy

The method of conformational analysis was applied to the spatial structures of peptide analogues of phytochelatins and some fragments of metallothioneins: (Cys-Gly)3, (Cys-Gly)3Asp, (Cys-Gly)3Glu, (Cys-betaAla)3, (Cys-gammaGlu)3, and (Cys-Gly-Gly)3. All the possible low-energy conformations of the molecules were revealed and the role of intra- and inter-residual interactions in the formation of their spatial structures was determined. A different tendency of the molecules under study for acceptance of conformations favorable for binding bismuth ions was shown. Low-energy structures providing an optimum binding of bismuth ion were shown to be most frequent for (Cys-betaAla)3 peptide. Among the analogues of peptide fragments of the metallothioneins, lacking in natural peptides, low-energy pentapeptide CCXXC fragments (where X = Gln, Asn, Phe, Tyr, and Gly) were revealed. In the alpha-helical conformations of these pentapeptides, the distance between the sulfur atoms corresponds to that in Bi2S3. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2006, vol. 32, no. 3; see also http://www.maik.ru.     

Comparative studies for amyloid beta degradation: “Neprilysin vs insulysin”, “monomeric vs aggregate”, and “whole Aβ40 vs its peptide fragments”

Amyloid beta (Aβ) proteins are produced from amyloid precursor protein cleaved by β- and γ-secretases, and are the main components of senile plaques pathologically found in Alzheimer's disease (AD) patient brains. Therefore, the relationship between AD and Aβs has been well studied for both therapeutic and diagnostic purposes. Several enzymes have been reported to degrade Aβs in vivo, with neprilysin (NEP) and insulysin (insulin-degrading enzyme, IDE) being the most prominent. In this article, we describe the mass spectrometric characterization of peptide fragments generated using NEP and IDE, and clarify the differences in digestion specificities between these two enzymes for non-aggregated Aβ₄₀, aggregated Aβ₄₀, and Aβ₄₀ peptide fragments, including Aβ₁₆. Our results allowed identification of all the peptide fragments from non-aggregated Aβ₄₀: NEP, 23 peptide fragments consisting of 2–11 amino-acid residues, 17 cleavage sites; IDE, 23 peptide fragments consisting of 6–33 amino-acid residues, 15 cleavage sites. Also, we confirmed that IDE can digest only whole Aβ₄₀, whereas NEP can digest both Aβ₄₀ and partial structures such as Aβ₁₆ and peptide fragments generated by the digestion of Aβ₄₀ by IDE. Furthermore, we confirmed that IDE and NEP are unable to digest aggregated Aβ₄₀.     

Characterizing aqueous solution conformations of a peptide backbone using Raman optical activity computations

Mounting spectroscopic evidence indicates that alanine predominantly adopts extended polyproline II (PPII) conformations in short polypeptides. Here we analyze Raman optical activity (ROA) spectra of N-acetylalanine-N′-methylamide (Ala dipeptide) in H₂O and D₂O using density functional theory on Monte Carlo (MC) sampled geometries to examine the propensity of Ala dipeptide to adopt compact right-handed (α_R) and left-handed (α_L) helical conformations. The computed ROA spectra based on MC-sampled α_R and PPII peptide conformations contain all the key spectral features found in the measured spectra. However, there is no significant similarity between the measured and computed ROA spectra based on the α_L- and β-conformations sampled by the MC methods. This analysis suggests that Ala dipeptide populates the α_R and PPII conformations but no substantial population of α_L- or β-structures, despite sampling α_L- and β-structures in our MC simulations. Thus, ROA spectra combined with the theoretical analysis allow us to determine the dominant populated structures. Including explicit solute-solvent interactions in the theoretical analysis is essential for the success of this approach.     

Spatial structure of BAM-12P dodecapeptide and its analogues

Theoretical conformational analysis was used to study the spatial structure and conformational properties of the bovine adrenal medulla dodecapeptide BAM-12P (Tyr1-Gly2-Gly3-Phe4-Met5-Arg6-Arg7-Val8-Gly9-Arg10-Pro11-Glu12). Twenty-three low-energy conformations of the BAM-12P backbone were shown to represent the spatial structure of the peptide. The inverse structural problem was solved, and synthetic analogues of BAM-12P were proposed, the spatial structures of which correspond to a set of low-energy potentially physiologically active conformations of the natural dodecapeptide. The English version of the paper: Russian Journal of Bioorganic Chemistry, 2005, vol. 31, no. 3; see also http://www.maik.ru.     

Spatial Structure of BAM-12P Dodecapeptide and Its Analogues

Theoretical conformational analysis was used to study the spatial structure and conformational properties of the bovine adrenal medulla dodecapeptide BAM-12P (Tyr1-Gly2-Gly3-Phe4-Met5-Arg6-Arg7-Val8-Gly9-Arg10-Pro11-Glu12). Twenty-three low-energy conformations of the BAM-12P backbone were shown to represent the spatial structure of the peptide. The inverse structural problem was solved, and synthetic analogues of BAM-12P were proposed, the spatial structures of which correspond to a set of low-energy potentially physiologically active conformations of the natural dodecapeptide.__________Translated from Bioorganicheskaya Khimiya, Vol. 31, No. 3, 2005, pp. 245–250.Original Russian Text Copyright © 2005 by Akhmedov, Tagiev, Hasanov, Makhmudova.     

Conformational features responsible for the binding of cyclic analogues of enkephalin to opioid receptors. III. Probable binding conformations of mu-agonists with phenylalanine in position 3.

Theoretical conformational analysis was carried out for several tetrapeptide analogues of beta-casomorphin and dermorphin containing a Phe residue in position 3. Sets of low-energy backbone structures of the mu-selective peptides [N-Me-Phe3, D-Pro4]-morphiceptin and Tyr-D-Orn-Phe-Asp-NH2 were obtained. These sets of structures were compared for geometrical similarity between themselves and with the low-energy conformations found for the delta-selective peptide Tyr-D-Cys-Phe-D-Pen-OH and nonactive peptide Tyr-Orn-Phe-Asp-NH2. Two pairs of geometrically similar conformations of mu-selective peptides, sharing no similarity with the conformations of peptides showing low affinity to the mu-receptor, were selected as two alternative models of probable mu-receptor-bound backbone conformations. Both models share geometrical similarity with the low-energy structures of the linear mu-selective peptide Tyr-D-Ala-Phe-Phe-NH2. Putative binding conformations of Tyr1 and Phe3 side chains are also discussed.     

DFT conformation and energies of amylose fragments at atomic resolution. Part 1: syn forms of α-maltotetraose

DFT optimization studies of 90 syn α-maltotetraose (DP-4) amylose fragments have been carried out at the B3LYP/6-311++G∗∗ level of theory. The DP-4 fragments studied include V-helix, tightly bent conformations, a boat, and a ¹C₄ conformer. The standard hydroxymethyl rotamers (gg, gt, tg) were examined at different locations in the residue sequence, and their influence on the bridge conformations ?/ψ values and conformer energy is described. Hydroxyl groups were considered to be homodromic, that is, they are either in the all clockwise, ‘c’, or all counterclockwise, ‘r’. Energy differences between conformations are examined in order to assess the stability of the different conformations and to identify the sources of energy that dictate amylose polymer formation. A small nearly cyclic compact structure is of low energy as one would expect when these flexible molecules are studied in vacuo. Many conformations in which the only differences are a single hydroxymethyl variation in the residue sequence show similar energies and bridge conformations, with trends being a result of the hydroxymethyl as well as hydroxyl orientation. In general the ‘c’ structures are of lower energy than the ‘r’ structures, although this is only true for the in vacuo state. The solvent dependence on conformational preference of several low-energy DP-4 structures was investigated via the continuum solvation method COSMO. These results suggest that the ‘r’ structures may be favored for fully solvated molecules.     

Conformational analysis of β-methyl-para-nitrophenylalanine stereoisomers of cyclo[D-Pen2, D-Pen5]enkephalin by NMR spectroscopy and conformational energy calculations

Solution conformations of β-methyl-para-nitrophenylalanine⁴ analogues of the potent δ-opioid peptide cyclo[D-Pen², D-Pen⁵]enkephalin (DPDPE) were studied by combined use of nmr and conformational energy calculations. Nuclear Overhauser effect connectivities and ³J_HNC^α_H coupling constants measured for the (2S, 3S)-, (2S, 3R)-, and (2R, 3R)-stereoisomers of[β-Me-p-NO₂Phe⁴]DPDPE in DMSO were compared with low energy conformers obtained by energy minimization in the Empirical Conformational Energy Program for Peptides #2 force field. The conformers that satisfied all available nmr data were selected as probable solution conformations of these peptides. Side-chain rotamer populations, established using homonuclear (³J_H^α_H^β) and heteronuclear (³J_H^αC^γ) coupling constants and ¹³C chemical shifts, show that the β-methyl substituent eliminates one of the three staggered rotamers of the torsion angle x¹ for each stereoisomer of the β-Me-p-NO₂Phe⁴. Similar solution conformations were suggested for the L-Phe⁴-containing (2S, 3S)- and (2S, 3R)-stereoisomers. Despite some local differences, solution conformations of L- and D-Phe⁴-containing analogues have a common shape of the peptide backbone and allow similar orientations of the main δ-opioid pharmacophores. This type of structure differs from several models of the solution conformations of DPDPE, and from the model of biologically active conformations of DPDPE suggested earlier. The latter model is allowed for the potent (2S, 3S)- and (2S, 3R)-stereoisomers of [β-Me-p-NO₂Phe⁴] DPDPE, but it is forbidden for the less active (2R, 3R)- and (2R, 3S)-stereoisomers. It was concluded that the biologically active stereoisomers of [β-Me-p-No₂Phe⁴] DPDPE in the δ-receptor-bound state may assume a conformation different from their favorable conformations in DMSO. © 1996 John Wiley & Sons, Inc.     

Ramachandran analysis of conserved glycyl residues in homologous proteins of known structure

High conservation of glycyl residues in homologous proteins is fairly frequent. It is commonly understood that glycine tends to be highly conserved either because of its unique Ramachandran angles or to avoid steric clash that would arise with a larger side chain. Using a database of aligned 3D structures of homologous proteins we identified conserved Gly in 288 alignment positions from 85 families. Ninety‐six of these alignment positions correspond to conserved Gly residue with (φ, ψ) values allowed for non‐glycyl residues. Reasons for this observation were investigated by in‐silico mutation of these glycyl residues to Ala. We found in 94% of the cases a short contact exists between the C^β atom of the introduced Ala with the atoms which are often distant in the primary structure. This suggests the lack of space even for a short side chain thereby explaining high conservation of glycyl residues even when they adopt (φ, ψ) values allowed for Ala. In 189 alignment positions, the conserved glycyl residues adopt (φ, ψ) values which are disallowed for Ala. In‐silico mutation of these Gly residues to Ala almost always results in steric hindrance involving C^β atom of Ala as one would expect by comparing Ramachandran maps for Ala and Gly. Rare occurrence of the disallowed glycyl conformations even in ultrahigh resolution protein structures are accompanied by short contacts in the crystal structures and such disallowed conformations are not conserved in the homologues. These observations raise the doubt on the accuracy of such glycyl conformations in proteins.     

Influence of local interactions on protein structure. I. Conformational energy studies of N-acetyl-N′-methylamides of pro-X and X-pro dipeptides

Conformational energy calculations using an Empirical Conformational Energy Program for Peptides (ECEPP) were carried out on the N-acetyl-N′-methylamides of Pro-X, where X = Ala, Asn, Asp, Gly, Leu, Phe, Ser, and Val, and of X-Pro, where X = Ala, Asn, Gly, and Pro. The conformational energy was minimized from starting conformations which included all combinations of low-energy single-residue minima and several standard bend structures. It was found that almost all resulting minima are combinations of low-energy single-residue minima, suggesting that intra residue interactions predominate in determining conformation. The calculations also indicate, however, that inter residue interactions can be important. In addition, librational entropy was found to influence the relative stabilities of some minima. Because of the existence of 10–100 low-energy minima for each dipeptide, the normalized statistical weight of an individual minimum rarely exceeds 0.3, suggesting that these dipeptides have considerable conformational flexibility and exist as statistical ensembles of low-energy structures. The propensity of each dipeptide to form bend conformations was calculated, and the results were compared with available experimental data. It was found that bends are favored in Pro-X dipeptides because ?_Pro is fixed by the pyrrolidine ring in a conformation which is frequently found in bends, but that bends are not favored in X-Pro dipeptides because interactions between the X residue and the pyrrolidine ring restrict the X residue to conformations which are not usually found in bends.     

NMR investigations of protein–carbohydrate interactions: insights into the topology of the bound conformation of a lactose isomer and β-galactosyl xyloses to mistletoe lectin and galectin-1

A hallmark of oligosaccharides is their often limited spatial flexibility, allowing them to access a distinct set of conformers in solution. Viewing each individual or even the complete ensemble of conformations as potential binding partner(s) for lectins in protein–carbohydrate interactions, it is pertinent to address the question on the characteristics of bound state conformation(s) in solution. Also, it is possible that entering the lectin’s binding site distorts the low-energy topology of a glycosidic linkage. As a step to delineate the strategy of ligand selection for galactosides, a common physiological docking point, we have performed a NMR study on two non-homologous lectins showing identical monosaccharide specificity. Thus, the conformation of lactose analogues bound to bovine heart galectin-1 and to mistletoe lectin in solution has been determined by transferred nuclear Overhauser effect measurements. It is demonstrated that the lectins select the syn conformation of lactose and various structural analogues (Galβ(1→4)Xyl, Galβ(1→3)Xyl, Galβ(1→2)Xyl, and Galβ(1→3)Glc) from the ensemble of presented conformations. No evidence for conformational distortion was obtained. Docking of the analogues to the modeled binding sites furnishes explanations, in structural terms, for exclusive recognition of the syn conformer despite the non-homologous design of the binding sites.     

De novo design,synthesis and solution conformational study of two didehydroundecapeptides: effect of nature and number of amino acids interspersed between Phe residues

De novo design of peptides and proteins has recently surfaced as an approach for investigating protein structure and function. This approach vitally tests our knowledge of protein folding and function, while also laying the groundwork for the fabrication of proteins with properties not precedented in nature. The success relies heavily on the ability to design relatively short peptides that can espouse stable secondary structures. To this end, substitution with α,β‐didehydroamino acids, especially α,β‐didehydrophenylalanine (Δ^zPhe), comes in use for spawning well‐defined structural motifs. Introduction of ΔPhe induces β‐bends in small and 3₁₀‐helices in longer peptide sequences. The present work aims to investigate the effect of nature and the number of amino acids interspersed between two ΔPhe residues in two model undecapeptides, Ac‐Gly‐Ala‐ΔPhe‐Ile‐Val‐ΔPhe‐Ile‐Val‐ΔPhe‐Ala‐Gly‐NH₂ (I) and Boc‐Val‐ΔPhe‐Phe‐Ala‐Phe‐ΔPhe‐Phe‐Leu‐Ala‐ΔPhe‐Gly‐OMe (II). Peptide I was synthesized using solid‐phase chemistry and characterized using circular dichroism spectroscopy. Peptide II was synthesized using solution‐phase chemistry and characterized using circular dichroism and nuclear magnetic resonance spectroscopy. Peptide I was designed to examine the effect of incorporating β‐strand‐favoring residues like valine and isoleucine as spacers between two ΔPhe residues on the final conformation of the resulting peptide. Circular dichroism studies on this peptide have shown the existence of a 3₁₀‐helical conformation. Peptide II possesses three amino acids as spacers between ΔPhe residues and has been reported to adopt a mixed 3₁₀/α‐helical conformation using circular dichroism and nuclear magnetic resonance spectroscopy studies. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Comparison by 1H-nmr spectroscopy of the conformation of the 2600 dalton antifreeze glycopeptide of polar cod with that of the high molecular weight antifreeze glycoprotein

The antifreeze glycopeptide (AFGP-8) from polar cod, B. saida, is a 14-amino acid polypeptide having alternating glycotripeptide sequences of Ala-[Gal(β1 → 3)GalNAc(β1 → O)]-Thr-Pro and Ala-[Gal(β1 → 3)GalNAc(β1 → O)]-Thr-Ala, with alanyl residues at amino and carboxy terminals. Conformational studies of AFGP-8 have been carried out by ¹H-nmr and empirical energy calculations to investigate the difference in its antifreeze behavior from that of the more active high-molecular weight AFGP 1-4 of P. borchgrevinki. The ¹H-nmr spectra, including the resonances of the exchangeable amide protons, were assigned by two-dimensional correlated spectroscopy (COSY), one-dimensional difference decoupling, and nuclear Overhauser effect (NOE) measurements. For the four threonyl residues, the amide proton coupling constants and the small coupling constants between H^α and H^β indicate similar conformations, despite significant chemical shift differences. The strong NOE between the α protons and the amide protons of the residue following together with large temperature coefficients of chemical shifts, indicate an extended conformation not consisting of α-helix, turns or bends. Energy computations indicate several low-energy conformations consistent with the observed coupling constants for ?. Among these, a left-handed helical conformation with three repeating residues per turn has been proposed, which is in accordance with the observed NOE between the methyl group of the α-GalNAc and Ala H^βs. While the observed Overhauser effects in the threonyl side chain suggest a certain amount of conformational averaging, the effect involving the acetmido methyl of α-GalNAc and H^βs of Ala indicate that it as is a major conformer. In view of the close similarity between the conformations of AFGP-8 and the more active antifreeze polymer, AFGP 1-4, we propose that the difference in their activities is due to the length of the regular repeating structure with glycosylation at every third amino acid residue, and not due to any fundamental difference in their conformations.     

The conformational properties of α,β‐dehydroamino acids with a C‐terminal ester group

α,β‐Dehydroamino acid esters occur in nature. To investigate their conformational properties, a systematic theoretical analysis was performed on the model molecules Ac‐ΔXaa‐OMe [ΔXaa = ΔAla, (E)‐ΔAbu, (Z)‐ΔAbu, ΔVal] at the B3LYP/6‐311+ + G(d,p) level in the gas phase as well as in chloroform and water solutions with the self‐consistent reaction field‐polarisable continuum model method. The Fourier transform IR spectra in CCl₄ and CHCl₃ have been analysed as well as the analogous solid state conformations drawn from The Cambridge Structural Database. The ΔAla residue has a considerable tendency to adopt planar conformations C5 (?, ψ ≈ ? 180°, 180°) and β2 (?, ψ ≈ ? 180°, 0°), regardless of the environment. The ΔVal residue prefers the conformation β2 (?, ψ ≈ ? 120°, 0°) in a low polar environment, but the conformations α (?, ψ ≈ ? 55°, 35°) and β (?, ψ ≈ ? 55°, 145°) when the polarity increases. The ΔAbu residues reveal intermediate properties, but their conformational dispositions depend on configuration of the side chain of residue: (E)‐ΔAbu is similar to ΔAla, whereas (Z)‐ΔAbu to ΔVal. Results indicate that the low‐energy conformation β2 is the characteristic feature of dehydroamino acid esters. The studied molecules constitute conformational patterns for dehydroamino acid esters with various side chain substituents in either or both Z and E positions. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Studies on the conformation of amino acids. X. Conformations of norvalyl, leucyl and aromatic side groups in a dipeptide unit

Backbone-side group conformations of amino acid residues including one or two δ-carbons in the side group have been investigated. Conformational energies of norvalyl, leucyl, phenylalanyl, tyrosyl, tryptophenyl, and histidinyl side groups in a dipeptide unit have been calculated by using classical energy expressions. The side group conformations about the C^α—C^β and C^β—C^γ bonds are restricted to specific values of the respective rotational angles. Thus, most favourable positions of γ- and δ-atoms of a linear side-chain (norvalyl) are restricted to (γ_I, δ_II) (γ_II, δ_I), (γ_II, δ_II), (γ_III, δ_II), and (γ_III, δ_III), whereas those of the side-chain branching at a sp³ γ-atom (leueyl) are further restricted. It is also shown that there is a definite correlation between the orientations of the two peptide planes and that of the planar group of the aromatic side chain of phenylalanyl type residues. The studies bring out an important fact that while the γ-atoms have definite and characteristic effects on the backbone rotational angles ? and ψ, the δ atoms and beyond have no effects on the preferred ? and ψ values. Thus, the preferred backbone conformations are independent of the preferred side group conformations beyond the γ-atom and vice versa. The observed ?, ψ, χ¹, and χ² values of amino acids, simple peptides, and of the three protein molecules lysozyme, myoglobin, and chymotrypsin have been compared with the theoretical predictions, and the agreement is found to be excellent.     

Studies of conformational isomerism in α-melanocyte stimulating hormone by design of cyclic analogues

Results of energy calculations for α-MSH (α-melanocyte stimulating hormone, Ac-Ser¹-Tyr²-Ser³-Met⁴-Glu⁵-His⁶-Phe⁷-Arg⁸-Trp⁹-Gly¹⁰-Lys¹¹-Pro¹²-Val¹³-NH₂) and [D -Phe⁷]α-MSH were used for design of cyclic peptides with the general aim to stabilize different conformational isomers of the parent compound. The minimal structural modifications of the conformationally flexible Gly¹⁰ residue, as substitutions for L -Ala, D -Ala, or Aib (replacing of hydrogen atoms by methyl groups), were applied to obtain octa- and heptapeptide analogues of α-MSH(4–11) and α-MSH(5–11), which were cyclized by lactam bridges between the side chains in positions 5 and 11. Some of these analogues, namely those with substitutions of the Gly¹⁰ residue with L -Ala or Aib, showed biological activity potencies on frog skin comparable to the potency of the parent tridecapeptide hormone. Additional energy calculations for designed cyclic analogues were used for further refinement of the model for the biologically active conformations of the His-Phe-Arg-Trp “message” sequence within the sequences of α-MSH and [D -Phe⁷]α-MSH. In such conformations the aromatic moieties of the side chains of the His⁶, L/D -Phe⁷, and Trp⁹ residues form a continuous hydrophobic “surface,” presumably interacting with a complementary receptor site. This feature is characteristic for low-energy conformers of active cyclic analogues, but it is absent in the case of inactive analogues. This particular spatial arrangement of functional groups involved in the message sequence is very close for α-MSH and [D -Phe⁷]α-MSH, as well as for biologically active cyclic analogues despite differences of dihedral angle values for corresponding low-energy conformations. © 1998 John Wiley & Sons, Inc. Biopoly 46: 155–167, 1998     

β-Turn conformations in crystal structures of model peptides containing α,α-Di-n-propylglycine and α,α-Di-n-butylglycine

The crystal state conformations of three peptides containing the α,α-dialkylated residues. α,α-di-n-propylglycine (Dpg) and α,α-di-n-butylglycine (Dbg), have been established by x-ray diffraction. Boc-Ala-Dpg-Alu-OMe (I) and Boc-Ala-Dbg-Ala-OMe (III) adopt distorted type II β-turn conformations with Ala (1) and Dpg/Dbg (2) as the corner residues. In both peptides the conformational angles at the Dxg residue (I: ? = 66.2°, ψ = 19.3°; III: ? = 66.5°. ψ = 21.1°) deviate appreciably from ideal values for the i + 2 residue in a type II β-turn. In both peptides the observed (N…O) distances between the Boc CO and Ala (3) NH groups are far too long (1: 3.44 Å: III: 3.63 Å) for an intramolecular 4 → 1 hydrogen bond. Boc-Ala-Dpg-Ata-NHMe (II) crystallizes with two independent molecules in the asymmetric unit. Both molecules HA and HB adopt consecutive β-turn (type III-III in HA and type III-I in IIB) or incipient 3₁₀-helical structures, stabilized by two intramolecular 4 → 1 hydrogen bonds. In all four molecules the bond angle N-C^α-C′ (τ) at the Dxg residues are ≥ 110°. The observation of conformational angles in the helical region of ?,ψ space at these residues is consistent with theoretical predictions. © 1995 John Wiley & Sons, Inc.     

Approach to studying the structure-activity relationship of native peptides. I. Structural organization of the opioid peptide Tyr-Glu-Gly-Phe-Met-Arg-Gly-Leu and its artificial analogs]

Theoretical conformational analysis was carried out for the octapeptide Tyr1-Gly2-Gly3-Phe4-Met5-Arg6-Gly7-Leu8. Possible structure of the opioid peptide under physiological conditions may be described by a set of low-energy conformations belonging to 14 different forms of the backbone. The solution of the "reverse conformational problem" for the opioid peptide enables one to predict the modified amino acid sequences (Ala2, D-Ala2, Ala3, D-Ala3, Ala7, D-Ala7, MeMet5, MeArg6-analogues) which may assume one of the low-energy states of the native hormone. The influence of the solute was not taken into account in our calculations.