Raman and infrared spectroscopic study of gramicidin A conformations

Raman scattering and infrared spectroscopic techniques were used to study the vibrational spectrum and conformation of the membrane channel protein gramicidin A in the solid state, in organic solutions and, using Raman scattering only, in a phospholipid environment. The investigation also includes measurements on head- and tail-group-modifled gramicidin A and a potassium thiocyanate-gramicidin A complex. Tentative identification of the molecular vibrations is proposed on the basis of the data on model compounds. The existence of four distinct conformations of the gramicidin A chain is established: conformation I present in the solid state, and CH₃OH and CD₃OD solutions; conformation II present in films cast from CHCl₃ solution; conformation III present in (CH₃)₂SO and (CD₃)₂SO solutions at concentrations below 0.5 m gramicidin A; and conformation IV present in the potassium thiocyanate-gramicidin A complex. The data obtainable on a gramicidin A-phospholipid suspension indicate a gramicidin A conformation in this environment corresponding either to the conformation I or II. The details of the spectra in the amide I region are shown to be consistent with a β-parallel hydrogen-bonded π_LD helix for conformational I, in terms of the polypeptide vibrational calculations of Nevskaya and co-workers. Conformation II is found to be consistent with an antiparallel double-stranded π_LD helix, while conformations III and IV probably have π-helical structures with larger channel diameters. The data on head- and tail-modified gramicidin A molecules indicate that their conformations are only slightly different from that of gramicidin A in conformation I.     

Refinement of the spatial structure of the gramicidin A ion channel]

The spatial structure of the gramicidin A (GA) transmembrane ion-channel was refined on the base of cross-peak volumes measured in NOESY spectra (mixing time tau m = 100 and 200 ms). The refinement methods included the comparison of experimental cross-peak volumes with those calculated for low-energy GA conformations, dynamic averaging of the low-energy conformation set and restrained energy minimization. Accuracy of the spatial structure determination was estimated by the penalty function Fr defined as a root mean square deviation of interproton distances corresponding to the calculated and experimental cross-peak volumes. As the initial conformation we used the right-handed pi 6,3 LD pi 6,3 LD helix established on the base of NMR data regardless of the cross-peak volumes. The conformation is in a good agreement with NOE cross-peak volumes (Fr 0.2 to 0.5 A depending on NOESY spectrum). For a number of NOEs formed by the side chain protons, distances errors were found as much as 0.5-2.0 A. Restrained energy minimization procedure had little further success. However some of these errors were eliminated by the change in torsional angle chi 2 of D-Leu12 and dynamic averaging of the Val7 side chain conformations. Apparently, majority of deviations of the calculated and experimental cross-peak volumes are due to the intramolecular mobility of GA and cannot be eliminated within the framework of rigid globule model. In summary the spatial structure of GA ion-channel can be thought as a set of low-energy conformations, differing by the side chain torsion angles chi 1 Val7 and chi 2 D-Leu4 and D-Leu10 and the orientation of the C-terminal ethanolamine group. Root mean square differences between the atomic coordinates of conformations are in the range of 0.3-0.8 A.     

Conformation of 3'CMP bound to RNase A using TrNOESY

The conditions for accurately determining distance constraints from TrNOESY data on a small ligand (3'CMP) bound to a small protein (RNase A, <14 kDa) are described. For small proteins, normal TrNOESY conditions of 10:1 ligand:protein or greater can lead to inaccurate structures for the ligand-bound conformation due to the contribution of the free ligand to the TrNOESY signals. By using two ligand:protein ratios (2:1 and 5:1), which give the same distance constraints, a conformation of 3'CMP bound to RNase A was determined (glycosidic torsion angle, chi=-166 degrees ; pseudorotational phase angle, 0 degrees < or = P < or =36 degrees ). Ligand-protein NOESY cross peaks were also observed and used to dock 3'CMP into the binding pocket of the apo-protein (7rsa). After energy minimization, the conformation of the 3'CMP:RNase A complex was similar to the X-ray structure (1 rpf) except that a C3'-endo conformation for the ribose ring (rather than C2'-exo conformation) was found in the TrNOESY structure.     

Molecular modelling of glycans: three-dimensional structure and protein fraction interaction. The example of rabbit sero-transferrin

On the basis of experimental data and of computer calculations using the Tripos 5.3 force field in order to examine the three-dimensional structures which are sterically feasible and the conformations which are energetically the most favourable, we have designed a program of molecular modelling of biantennary glycans of the N-acetyllactosaminic type (complex type). We demonstrate that, in absence of any interaction with the protein, a high number of glycan conformations exists which can be classified into five basic conformations, four of which have already been described. In fact, in addition to the Y-, T-, "bird" and "broken-wing" conformations, a "back-folded wing" conformation is energetically feasible. In contrast, the glycan linked to the protein is immobilized into only one conformation: the "broken-wing" conformation. Forming a bridge between the two lobes of the peptide chain, it probably contributes to the maintenance of the protein in a biologically active conformation.     

DNA packaging in bacteriophage: is twist important?

We study the packaging of DNA into a bacteriophage capsid using computer simulation, specifically focusing on the potential impact of twist on the final packaged conformation. We perform two dynamic simulations of packaging a polymer chain into a spherical confinement: one where the chain end is rotated as it is fed, and one where the chain is fed without end rotation. The final packaged conformation exhibits distinct differences in these two cases: the packaged conformation from feeding with rotation exhibits a spool-like character that is consistent with experimental and previous theoretical work, whereas feeding without rotation results in a folded conformation inconsistent with a spool conformation. The chain segment density shows a layered structure, which is more pronounced for packaging with rotation. However, in both cases, the conformation is marked by frequent jumps of the polymer chain from layer to layer, potentially influencing the ability to disentangle during subsequent ejection. Ejection simulations with and without Brownian forces show that Brownian forces are necessary to achieve complete ejection of the polymer chain in the absence of external forces.     

On the dependence of molecular conformation on the type of solvent environment: a molecular dynamics study of cyclosporin A

The dependence of the conformation of cyclosporin A (CPA), a cyclic undecapeptide with potent immunosuppressive activity, on the type of solvent environment is examined using the computer simulation method of molecular dynamics (MD). Conformational and dynamic properties of CPA in aqueous solution are obtained from MD simulations of a CPA molecule dissolved in a box with water molecules. Corresponding properties of CPA in apolar solution are obtained from MD simulations of CPA in a box with carbontetrachloride. The results of these simulations in H2O and in CCl4 are compared to each other and to those of previous simulations of crystalline CPA and of an isolated CPA molecule. The conformation of the backbone of the cyclic polypeptide is basically independent of the type of solvent. In aqueous solution the beta-pleated sheet is slightly weaker and the gamma-turn is a bit less pronounced than in apolar solution. Side chains may adopt different conformations in different solvents. In apolar solution the hydrophobic side chain of the MeBmt residue is in an extended conformation with its hydroxyl group hydrogen bonded to the backbone carbonyl group. In aqueous solution this hydrophobic side chain folds over the core of the molecule and the mentioned hydrogen bond is broken in favor of hydrogen bonding to water molecules. The conformation obtained from the MD simulation in CCl4 nicely agrees with experimental atom-atom distance data as obtained from nmr experiments in chloroform. In aqueous solution the relaxation of atomic motion tends to be slower than in apolar solution.     

Nuclear Overhauser effect and cross-relaxation rate determinations of dihedral and transannular interproton distances in the decapeptide tyrocidine A.

The following interproton distances are reported for the decapeptide tyrocidine A in solution: (a) r(phi) distances between NH(i) and H alpha (i), (b) r(psi) distances between NH (i + 1) and H alpha (i), (c) r(phi psi) distances between NH(i + 1) and NH(i), (d) NH in equilibrium NH transannular distances, (e) H alpha in equilibrium H alpha transannular distances, (f) r x 1 distances between H alpha and H beta protons, (g) NH(i) in equilibrium H beta (i) distances, (h) NH (i + 1) in equilibrium H beta (i) distances, (i) carboxamide-backbone protons and carboxamide-side chain proton distances, (j) side chain proton-side chain proton distances. The procedures for distance calculations were: NOE ratios and calibration distances, sigma ratios and calibration distances, and correlation times and sigma parameters. The cross-relaxation parameters were obtained from the product, say, of NOE 1 leads to 2 and the monoselective relaxation rate of proton 2; the NOEs were measured by NOE difference spectroscopy. The data are consistent with a type I beta-turn/ type II' beta-turn/ approximately antiparallel beta-pleated sheet conformation of tyrocidine A in solution and the NOEs, cross-relaxation parameters, and interproton distances serve as distinguishing criteria for beta-turn and beta-pleated sheet conformations. It should be borne in mind that measurement of only r phi and r psi distances for a decapeptide only defines the ( phi, psi)-space in terms of 4(10) possible conformations; the distances b-j served to reduce the degeneracy in possible (phi, psi)-space to one tyrocidine A conformation. The latter conformation is consistent with that derived from scalar coupling constants, hydrogen bonding studies, and proton-chromophore distance measurement, and closely resembles the conformation of gramicidin S.     

Combined use of proton-proton Overhauser enhancements and a distance geometry algorithm for determination of polypeptide conformations. Application to micelle-bound glucagon

In a new approach for the determination of polypeptide conformation, experimental data on intramolecular distances between pairs of hydrogen atoms obtained from nuclear Overhauser enhancement studies are used as input for a distance geometry algorithm. The algorithm determines the limits of the conformation space occupied by the polypeptide chain. The experimental data are used in such a way that the real conformation should in all cases be within these limits. Two important features of the method are that the results do not depend critically on the accuracy of the distance measurements by nuclear Overhauser enhancement studies and that internal mobility of the polypeptide conformation is explicitly taken into consideration. The use of this new procedure is illustrated with a structural study of the region 19-27 of glucagon bound to perdeuterated dodecylphosphocholine micelles.     

Structure of the pyruvate dehydrogenase multienzyme complex E1 component from Escherichia coli at 1.85 A resolution

The crystal structure of the recombinant thiamin diphosphate-dependent E1 component from the Escherichia coli pyruvate dehydrogenase multienzyme complex (PDHc) has been determined at a resolution of 1.85 A. The E. coli PDHc E1 component E1p is a homodimeric enzyme and crystallizes with an intact dimer in an asymmetric unit. Each E1p subunit consists of three domains: N-terminal, middle, and C-terminal, with all having alpha/beta folds. The functional dimer contains two catalytic centers located at the interface between subunits. The ThDP cofactors are bound in the "V" conformation in clefts between the two subunits (binding involves the N-terminal and middle domains), and there is a common ThDP binding fold. The cofactors are completely buried, as only the C2 atoms are accessible from solution through the active site clefts. Significant structural differences are observed between individual domains of E1p relative to heterotetrameric multienzyme complex E1 components operating on branched chain substrates. These differences may be responsible for reported alternative E1p binding modes to E2 components within the respective complexes. This paper represents the first structural example of a functional pyruvate dehydrogenase E1p component from any species. It also provides the first representative example for the entire family of homodimeric (alpha2) E1 multienzyme complex components, and should serve as a model for this class of enzymes.     

Modeling the structure of agitoxin in complex with the Shaker K+ channel: a computational approach based on experimental distance restraints extracted from thermodynamic mutant cycles

Computational methods are used to determine the three-dimensional structure of the Agitoxin (AgTx2)-Shaker complex. In a first stage, a large number of models of the complex are generated using high temperature molecular dynamics, accounting for side chain flexibility with distance restraints deduced from thermodynamic analysis of double mutant cycles. Four plausible binding mode candidates are found using this procedure. In a second stage, the quality and validity of the resulting complexes is assessed by examining the stability of the binding modes during molecular dynamics simulations with explicit water molecules and by calculating the binding free energies of mutant proteins using a continuum solvent representation and comparing with experimental data. The docking protocol and the continuum solvent model are validated using the Barstar-Barnase and the lysozyme-antibody D1.2 complexes, for which there are high-resolution structures as well as double mutant data. This combination of computational methods permits the identification of two possible structural models of AgTx2 in complex with the Shaker K+ channel, additional structural analysis providing further evidence in favor of a single model. In this final complex, the toxin is bound to the extracellular entrance of the channel along the pore axis via a combination of hydrophobic, hydrogen bonding, and electrostatic interactions. The magnitude of the buried solvent accessible area corresponding to the protein-protein contact is on the order of 1000 A with roughly similar contributions from each of the four subunits. Some side chains of the toxin adopt different conformation than in the experimental solution structure, indicating the importance of an induced-fit upon the formation of the complex. In particular, the side chain of Lys-27, a residue highly conserved among scorpion toxins, points deep into the pore with its positively charge amino group positioned at the outer binding site for K+. Specific site-directed mutagenesis experiments are suggested to verify and confirm the structure of the toxin-channel complex.     

Molecular basis for bacterial class I release factor methylation by PrmC

Class I release factors bind to ribosomes in response to stop codons and trigger peptidyl-tRNA hydrolysis at the P site. Prokaryotic and eukaryotic RFs share one motif: a GGQ tripeptide positioned in a loop at the end of a stem region that interacts with the ribosomal peptidyl transferase center. The glutamine side chain of this motif is specifically methylated in both prokaryotes and eukaryotes. Methylation in E. coli is due to PrmC and results in strong stimulation of peptide chain release. We have solved the crystal structure of the complex between E. coli RF1 and PrmC bound to the methyl donor product AdoHCy. Both the GGQ domain (domain 3) and the central region (domains 2 and 4) of RF1 interact with PrmC. Structural and mutagenic data indicate a compact conformation of RF1 that is unlike its conformation when it is bound to the ribosome but is similar to the crystal structure of the protein alone.     

Ion binding of cyclolinopeptide A: an NMR and CD conformational study

CD and nmr techniques have been used to study, in acetonitrile solution, the ion-complexing capability of cyclolinopeptide A (CLA), a cyclic nonapeptide of sequence cyclo-(Pro-Pro-Phe-Phe-Leu-Ile-Ile-Leu-Val) endowed with remarkable cytoprotective ability in vitro, and the conformation of the Ba(2+)/CLA complex. At room temperature, CLA in acetonitrile shows a proton nmr spectrum characteristic of the coexistence of many different conformers in intermediate exchange. The backbone contains a cis Pro-Pro bond, with all other peptide bonds in the trans conformation. CLA binds Ba2+ more tightly than the other cations studied, namely K+, Na+, Mg2+, and Ca2+; CD data are indicative of the presence of both 1:2 (sandwich) and 1:1 (equimolar) type complexes, depending on the Ba2+ ion concentration, whereas nmr data are consistent with an equimolar form. The relevant conformational features of the equimolar Ba2+/CLA complex are that the backbone contains all trans peptide bonds, a type I 6----3 beta-turn and a 3----1 gamma-turn (or a distorted 3----9 beta-turn). The global shape of the complexed peptide can be described as a bowl, with the concave (polar) side hosting Ba2+ and the convex side predominantly apolar.     

The intercalation of imidazoacridinones into DNA induces conformational changes in their side chain

Imidazoacridinones (IAs) are a new group of highly active antitumor compounds. The intercalation of the IA molecule into DNA is the preliminary step in the mode of action of these compounds. There are no experimental data about the structure of an intercalation complex formed by imidazoacridinones. Therefore the design of new potentially better compounds of this group should employ the molecular modelling techniques. The results of molecular dynamics simulations performed for four IA analogues are presented. Each of the compounds was studied in two systems: i) in water, and ii) in the intercalation complex with dodecamer duplex d(GCGCGCGCGCGC)2. Significant differences in the conformation of the side chain in the two environments were observed for all studied IAs. These changes were induced by electrostatic as well as van der Waals interactions between the intercalator and DNA. Moreover, the results showed that the geometry of the intercalation complex depends on: i) the chemical constitution of the side chain, and ii) the substituent in position 8 of the ring system.     

Variable ranges of interactions in polypeptide conformations with a method to complement molecular modeling.

Formulations of conformational weights for helix-coil transitions can be extended to substantially more complex situations than are usually pursued. General rules for matrix multiplication that depend parametrically on the interaction ranges and numbers of rotamers of residues are presented. The orders of the matrices of statistical weights can be increased with chain length, so that an individual matrix element can represent any specified single conformation, as needed. By the appropriate choice of interaction ranges and numbers of available conformers, approximations can be introduced in which: (1) an average of the conformations of any chain segment is obtained, (2) specific residue-residue interactions are excluded, or (3) the conformation of a part of the chain is restricted or fixed. The method is appropriate for treating specific interactions in peptides and could be used together with available experimental information to develop models of conformational transitions. As such, the methods represent a class of calculations aimed at more rigorous calculations built around known features of a molecule. The aim is to facilitate calculations that bridge the gap between nonquantitative molecular model building and more rigorous but less directed molecular mechanics calculations. The method can directly include any desired longer range of interactions, if the interaction range is not too long to make impossible the manipulation of the requisite matrices. An outline is presented of an application to treat salt bridges in the C peptide of ribonuclease A.     

Two distinct myosin light chain structures are induced by specific variations within the bound IQ motifs-functional implications

IQ motifs are widespread in nature. Mlc1p is a calmodulin-like myosin light chain that binds to IQ motifs of a class V myosin, Myo2p, and an IQGAP-related protein, Iqg1p, playing a role in polarized growth and cytokinesis in Saccharomyces cerevisiae. The crystal structures of Mlc1p bound to IQ2 and IQ4 of Myo2p differ dramatically. When bound to IQ2, Mlc1p adopts a compact conformation in which both the N- and C-lobes interact with the IQ motif. However, in the complex with IQ4, the N-lobe no longer interacts with the IQ motif, resulting in an extended conformation of Mlc1p. The two light chain structures relate to two distinct subfamilies of IQ motifs, one of which does not interact with the N-lobes of calmodulin-like light chains. The correlation between light chain structure and IQ sequence is demonstrated further by sedimentation velocity analysis of complexes of Mlc1p with IQ motifs from Myo2p and Iqg1p. The resulting 'free' N-lobes of myosin light chains in the extended conformation could mediate the formation of ternary complexes during protein localization and/or partner recruitment.     

A stereochemical model of codon-anticodon interactions, in which bases bind via water bridges

Analysis of available experimental data has allowed us to conclude that upon codon-anticodon interaction there is a functional asymmetry between the third nucleotide of the codon and the first nucleotide of the anticodon. We suggest that this phenomenon would be due to the differential codon and anticodon conformation mobilities, which are somehow restricted. Two postulates are proposed to restrict these mobilities respectively: (1) In the codon-anticodon complex the codon is in the fixed A-RNA conformation, and the rearrangement of the sugar-phosphate backbone, that is necessary for any wobble-pair formation occurs only in the anticodon. (2) This rearrangement does not alter the specific system of contacts between the universal uridine and other nucleotides in the anticodon loop, observed in tRNA crystals. Stereochemical analysis revealed that these restrictions make the construction of any Crick's wobble-pair sterically impossible. To resolve this contradiction we propose that two bases in the third position of the codon-anticodon complex can mate not only according to the Crick's scheme, but also with the help of water bridges. Consequently, besides the standard pairs, the following pairs become permitted: UG, UI, UU, CU, GU and AI (it should be noted that GU and AI are weaker that the other ones). All other pairs appear sterically impossible. These conclusions are in good agreement with the literature date.     

NMR studies of the interaction of chromomycin A3 with small DNA duplexes I

1H and 31P NMR spectral analysis of a chromomycin/d(ATGCAT)2 complex provides strong evidence for a nonintercalative mode of drug binding. Investigation of the imino proton region of the duplex suggests a protection of one of the two guanine imino protons from fast exchange with the bulk water up to at least 45 degrees C by the drug. Subsequent one-dimensional nuclear Overhauser enhancement experiments place the exchangeable chromomycin chromophoric hydroxyl proton less than 0.45 nm from this guanine imino proton and the chromophore 7-methyl less than 0.45 from the internal thymine 6-proton and/or the guanine 8-proton. 1H two-dimensional NMR reveals that the duplex retains a right-handed B conformation but there are distortions at the TGC region of one chain and large deviations in the chemical shift of protons relative to the uncomplexed duplex in the other chain in the same TGC region. The data suggest that the chromomycin chromophore is oriented such that the hydrophilic side of the ring system is proximal to the helix center in the major groove near the TG region while the aromatic side of the ring is oriented away from the helix but is partially protected from the solvent by the aliphatic chain, which bends back over the two aromatic protons. Changes in the 31P spectrum of the duplex on binding of the drug are different from the effect of either actinomycin or netropsin on nucleic acid fragments.     

The structure and activity of a monomeric interferon-gamma:alpha-chain receptor signaling complex

BACKGROUND: Interferon-gamma (IFN-gamma) is a homodimeric cytokine that exerts its various activities by inducing the aggregation of two different receptors. The alpha chain receptor (IFN-gammaRalpha) is a high affinity receptor that binds to IFN-gamma in a symmetric bivalent manner to form a stable, intermediate 1:2 complex. This intermediate forms a binding template for the subsequent binding of two copies of the second receptor, beta chain receptor (IFN-gammaRbeta), producing the active 1:2:2 signaling complex. RESULTS: A single chain monovalent variant of IFN-gamma (scIFN-gamma) was constructed and complexed to one copy of the extracellular domain (ECD) of IFN-gammaRalpha. The structure of this 1:1 complex was determined and the hormone-receptor interface shown to be characterized by a number of hydrophilic interactions mediated by several highly ordered water networks. The scIFN-gamma interface consists of segments from each of the monomer chains of the homodimer. The principal hydrophobic contact of the receptor involves a tripeptide segment of the receptor having an unusual and high energy conformation. Despite containing only one binding site for IFN-gammaRalpha, the monovalent scIFN-gamma molecule has significant activity in antiviral biological assays. CONCLUSIONS: ScIFN-gamma binds the ECD of IFN-gammaRalpha through a highly hydrated interface with an important set of hormone-receptor contacts mediated through structured waters. Although the interface is highly hydrated, it supports tight binding and has a considerable degree of specificity. The biological activity of scIFN-gamma confirms that the scIFN-gamma:IFN-gammaRalpha complex represents a productive intermediate and that it can effectively recruit the other required component, IFN-gammaRbeta, to signal based on the 1:1:1 complex.     

Benzene-di-N-Substituted Carbamates as Conformationally Constrained Substrate Analogs of Cholesterol Esterase

Benzene-1,2-, 1,3-, and 1,4-di-N-substituted carbamates (1-15) are synthesized as the constrained analogs of gauche, eclipsed, and anti conformations, respectively, for the glycerol backbones of triacylglycerol. Carbamates 1-15 are characterized as the pseudo substrate inhibitors of cholesterol esterase. Long chain carbamates are more potent inhibitors than short chain ones. Comparison of different geometries for benzene-di-substituted carbamates, such as benzene-1,2-di-N-octylcarbamate (3) (ortho-3), benzene-1,3-di-N-octylcarbamate (8) (meta-8), and benzene-1,4-di-N-octylcarbamates (13) (para-13), indicates that inhibitory potencies are as followed: meta-8 > para-13 > ortho-3. Therefore, we suggest that the preferable conformation for the C(sn-1)-O/C(sn-2)-O glycerol backbone in the enzyme-triacylgycerol complex is the eclipsed conformation. Meanwhile, kinetic data indicate that among ortho, meta, and para carbamates, meta carbamates most resemble the substrate cholesterol ester.     

Thermodynamic criterion for the conformation of P1 residues of substrates and of inhibitors in complexes with serine proteinases.

Eglin c, turkey ovomucoid third domain, and bovine pancreatic trypsin inhibitor (Kunitz) are all standard mechanism, canonical protein inhibitors of serine proteinases. Each of the three belongs to a different inhibitor family. Therefore, all three have the same canonical conformation in their combining loops but differ in their scaffoldings. Eglin c (Leu45 at P1) binds to chymotrypsin much better than its Ala45 variant (the difference in standard free energy changes on binding is -5.00 kcal/mol). Similarly, turkey ovomucoid third domain (Leu18 at P1) binds to chymotrypsin much better than its Ala18 variant (the difference in standard free energy changes on binding is -4.70 kcal/mol). As these two differences are within the +/-400 cal/mol bandwidth (expected from the experimental error), one can conclude that the system is additive. On the basis that isoenergetic is isostructural, we expect that within both the P1 Ala pair and the P1 Leu pair, the conformation of the inhibitor's P1 side chain and of the enzyme's specificity pocket will be identical. This is confirmed, within the experimental error, by the available X-ray structures of complexes of bovine chymotrypsin Aalpha with eglin c () and with turkey ovomucoid third domain (). A comparison can also be made between the structures of P1 (Lys+)15 of bovine pancreatic trypsin inhibitor (Kunitz) ( and ) and of the P1 (Lys+)18 variant of turkey ovomucoid third domain (), both interacting with chymotrypsin. In this case, the conformation of the side chains is strikingly different. Bovine pancreatic trypsin inhibitor with (Lys+)15 at P1 binds to chymotrypsin more strongly than its Ala15 variant (the difference in standard free energy changes on binding is -1.90 kcal/mol). In contrast, turkey ovomucoid third domain variant with (Lys+)18 at P1 binds to chymotrypsin less strongly than its Ala18 variant (the difference in standard free energies of association is 0.95 kcal/mol). In this case, P1 Lys+ is neither isostructural nor isoenergetic. Thus, a thermodynamic criterion for whether the conformation of a P1 side chain in the complex matches that of an already determined one is at hand. Such a criterion may be useful in reducing the number of required X-ray crystallographic structure determinations. More importantly, the criterion can be applied to situations where direct determination of the structure is extremely difficult. Here, we apply it to determine the conformation of the Lys+ side chain in the transition state complex of a substrate with chymotrypsin. On the basis of kcat/KM measurements, the difference in free energies of activation for Suc-AAPX-pna when X is Lys+ and X is Ala is 1.29 kcal/mol. This is in good agreement with the corresponding difference for turkey ovomucoid third domain variants but in sharp contrast to the bovine pancreatic trypsin inhibitor (Kunitz) data. Therefore, we expect that in the transition state complex of this substrate with chymotrypsin, the P1 Lys+ side chain is deeply inserted into the enzyme's specificity pocket as it is in the (Lys+)18 turkey ovomucoid third domain complex with chymotrypsin.