Conformational aspects of beta-trypsin interaction with substrates and pancreatic trypsin inhibitor. I. Conformational properties of residues in the enzyme-active center and the structure of a non-valent enzyme-substrate complex

The theoretical conformational analysis of potential surfaces of Ser-195, His-57, Asp-102 and Gln-192 side chains in the active center of native beta-trypsin has been carried out. The above residues are shown to exist in low-energy conformational states in the free enzyme and in its nonbonded substrate complexes. Interrelations between the flexibility of the residues and their catalytical functions were revealed. Conformational aspects of interaction of trypsin with N-acetyl-L-lysine and N-acetyl-L-alanyl--L-alanyl--L-lysyl--L-alanyl methyl amide and with the Gly-12--Ile-19 BPTI fragment were analysed. The productive binding of the substrate at the nonbonded complex stage is shown to take place exclusively in the lowest energy conformation of the enzyme-substrate complex. Basing on theoretical and experimental evidence, the problems of primary and secondary specificity of trypsin, and potential properties of the native protein to form a productive nonbonded complex and to react at the subsequent stages of the catalytical act are discussed. Conformational changes in the active center and interactions with a substrate are shown to result from stabilizing enzyme-substrate interactions. Trypsin inhibition by BPTI molecule does not take place at the nonbonded complex stage.     

Theoretical conformational analysis of methylamide of N-acetyl-L-lysine]

The spatial structure of the methylamide of N-acetyl-L-lysine has been analysed taking into account non-bonded and electrostatic interactions, torsional energy, bond angles distortion and hydrogen bonding. Conformational capacities of the backbone and mutual dependence of spatial structures of the backbone and the side chain was described by conformational maps obtained by energy minimisation, the dihedral angles and the bond angles of the side chain being varied for every phi, psi point. Every possible combination for phi, psi, x1-x5-angles was used corresponding to the stable form of the backbone and to torsion potential minima of the initial approximations in the calculation of preferred conformations of the molecule. Comparisons are made between stable forms of the methylamide of N-acetyl-L-lysine and Lys residues in proteins with known structure.     

Complexes of glucoamylase with maltoside heteroanalogues: bound ligand conformations by use of transferred NOE experiments and molecular modeling

Transferred nuclear Overhauser effect (trNOE) experiments have been performed to investigate the conformations of the competitive inhibitors, methyl 5'-thio-4-N-alpha-maltoside 3a and methyl 5'-thio-4-S-alpha-maltoside 4 when bound to the catalytic subunit of the enzyme glucoamylase. These NMR data suggest that, although each of the free ligands populates two conformational families, both heteroanalogues are bound by the enzyme in conformations in the area of the global energy minimum. These conformations have been used as initial points for docking into the active site of the enzyme taken from a X-ray crystal structure of the related glucoamylase-D-gluco-dihydroacarbose 2 complex. Minimization of the resulting complexes has yielded structures for the bound complexes. Corroboration of the structures is provided by fast T(1)(rho)-relaxation effects for certain ligand protons as a result of close contacts with protons in the enzyme active site. The results auger well for the combined use of transferred NOE spectroscopy and molecular modeling based on X-ray crystal structures of complexes of suitable congeners for the rapid analysis of ligand-receptor interactions.     

Atropisomeric amides: Achiral ligands with chiral conformations

A new class of achiral ligands with atropisomeric conformations has been coordinated to titanium(IV). The ligands are ortho-hydroxy benzamide derivatives which are deprotonated on reaction with titanium tetraisopropoxide to furnish Ti(L)(2)(O-iPr)(2) complexes (L=ortho-phenoxy benzamide). In these octahedral titanium compounds, the ortho-phenoxy benzamide ligands chelate to titanium, bonding through the phenoxide oxygen and the amide carbonyl oxygen. The benzamide ligands adopt atropisomeric conformations with an angle between the aryl and amide groups of approximately 35 degrees. The ligand precursor, ligand, and titanium complexes have been characterized by X-ray crystallography. Only one diastereomer of each titanium complex was observed in the solid state structures.     

Solution structure and dynamics of cyclic and acyclic cholinergic agonists.

Two classes of nicotinic cholinergic agonists, which vary in flexibility and electronegativity, have been synthesized, and their structural and dynamic properties have been studied with nuclear magnetic resonance (NMR) spectroscopy. Although the compounds are chemically identical except for the presence or absence of one cyclicizing C--C bond, single channel recording and radioligand binding studies have shown that the cyclic compounds are considerably more potent than the acyclic derivatives (McGroddy, K.A., A.A. Carter, M.M. Tubbert, and R.E. Oswald. 1993. Biophys. J. 64:325-338). Using one- and two-dimensional NMR spectroscopy, we have shown that these molecules exist in two distinct stable conformers, which differ in the orientation of the amide bond. The cyclic 1,1-dimethyl-4-trifluoroacetyl-piperazinium iodide and its trifluoromethyl derivative compounds are symmetric, and the two conformers are of equal energy. The acyclic N,N,N,N'-tetramethyl-N'-acetylethylene-diamine iodide (TED) and its trifluoromethyl derivative derivatives, however, populate two energetically unequal solution conformations. Using variable temperature NMR spectroscopy on these molecules and their uncharged precursors, we have characterized the energetics of amide bond isomerization and have distinguished steric and electrostatic contributions to the equilibrium between the two conformers. The more populated TED conformer has the amide methyl group trans to the carbonyl oxygen, and it is stabilized by an electrostatic attraction between the partially negative carbonyl oxygen and the positively charged quaternary amine nitrogen. As discussed in the accompanying paper (McGroddy, K.A., A.A. Carter, M.M. Tubbert, and R.E. Oswald. 1993. Biophys. J. 64:325-338), the differences in the stable solution structures of the TED derivatives and their interconversion kinetics may be of biological significance.     

The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

Vibrational analysis of the structure of gramicidin A. I. Normal mode analysis.

Normal mode frequencies have been calculated for single-stranded beta 4.4 and beta 6.3 and for double-stranded increases decreases beta 5.6, increases decreases beta 7.2, increases increases beta 5.6, and increases increases beta 7.2 helices that are possible models for the structure of gramicidin A. The force field used in the calculations is one that reproduces the frequencies of model polypeptide chain structures to about +/- 5 cm-1, and is therefore expected to provide meaningful distinctions between these conformations. The calculations predict significant differences in the infrared and Raman spectra of these beta-helices, suggesting that they should be identifiable from their spectra (which is shown in the following paper to be the case). The most sensitive region is that of the amide I frequencies, where the predicted patterns of intense infrared mode, infrared splittings, and intense Raman mode provide a characteristic identification of each of the above structures.     

S-adenosylmethionine conformations in solution and in protein complexes: conformational influences of the sulfonium group

S-Adenosylmethionine (AdoMet) and other sulfonium ions play central roles in the metabolism of all organisms. The conformational preferences of AdoMet and two other biologically important sulfonium ions, S-methylmethionine and dimethylsulfonioproprionic acid, have been investigated by NMR and computational studies. Molecular mechanics parameters for the sulfonium center have been developed for the AMBER force field to permit analysis of NMR results and to enable comparison of the relative energies of the different conformations of AdoMet that have been found in crystal structures of complexes with proteins. S-Methylmethionine and S-dimethylsulfonioproprionate adopt a variety of conformations in aqueous solution; a conformation with an electrostatic interaction between the sulfonium sulfur and the carboxylate group is not noticeably favored, in contrast to the preferred conformation found by in vacuo calculations. Nuclear Overhauser effect measurements and computational results for AdoMet indicate a predominantly anti conformation about the glycosidic bond with a variety of conformations about the methionyl C(alpha)-C(beta) and C(beta)-C(gamma) bonds. An AdoMet conformation in which the positively charged sulfonium sulfur is near an electronegative oxygen in the ribose ring is common. Comparisons of NMR results for AdoMet with those for the uncharged S-adenosylhomocysteine and 5'-methylthioadenosine, and the anionic ATP, indicate that the solution conformations are not dictated mainly by molecular charge. In 20 reported structures of AdoMet.protein complexes, both anti and syn glycosidic torsional angles are found. The methionyl group typically adopts an extended conformation in complexes with enzymes that transfer the methyl group from the sulfonium center, but is more folded in complexes with proteins that do not catalyze reactions involving the sulfur and which can use the sulfonium sulfur solely as a binding site. The conformational energies of AdoMet in these crystal structures are comparable to those found for AdoMet in solution. The sulfonium sulfur is in van der Waals contact with a protein heteroatom in the structures of four proteins, which reflects an energetically favorable contact. Interactions of the sulfonium with aromatic rings are rarely observed.     

The application of 1H NMR chemical shift calculations to diastereotopic groups in proteins.

We have calculated chemical shifts for a range of diastereotopic protons in proteins (i.e. methylene protons, and the methyl groups of valine and leucine residues), using a recently optimised method for chemical shift calculation. The calculations are based on crystal structure coordinates, and have been compared with experimental stereospecific assignments. The results indicate that chemical shifts can be used to suggest stereospecific assignments with about 80% probability of being correct, in cases where both the experimental and the calculated chemical shift differences between a pair of diastereotopic protons are greater than 0.3 ppm. Inaccurate calculations are shown to be caused in most cases by differences between crystal and solution structures. Furthermore, chemical shift calculations based on NMR structures are shown to be capable of acting as a further constraint on structure, by limiting the range of side-chain conformations adopted in structures calculated from NMR data.     

NMR studies of a series of dehydrodermorphins

The third and fifth aromatic residues of dermorphin, a potent mu-opioid peptide, and of its N-terminal fragments, from the pentapeptide to the parent heptapeptide amide, have been systematically substituted with Z-dehydrophenylalanine (delta-Phe) and/or Phe to investigate the conformation-activity relationship. The characterization in DMSO-d6 at 500 MHz indicates that, in this solvent, all peptides adopt essentially random, extended conformations, as a consequence of the strong solvation. The chemical shift of the methyl group of D-Ala is influenced by the precise orientation of the side chain of the third residue in a fashion that can be correlated to the mu potency, consistently with our model of mu-receptor. However, the complexes of the pentapeptides with 18-crown-6-ether, when dissolved in chloroform, adopt ordered, folded conformations, a behavior that closely parallels the CD observations in methanol.     

Experimental attempt to simulate receptor site environment. A 500-MHz 1H nuclear magnetic resonance study of enkephalin amides

The amides of Leu5-enkephalin, Met5-enkephalin, and three analogues, D-Ala2,Leu5-enkephalin, (AcO)Tyr1,Met5-enkephalin, and (AcO)Tyr1,D-Ala2,Met5-enkephalin, have been studied by means of 1H NMR spectroscopy in two different solvent systems: Me2SO-d6 and CDCl3. In the latter solvent the peptides were dissolved as complexes with 18-crown-6-ether, a coronand that binds strongly to the NH3+ groups. The crown ether complexation and the apolar solvent were used to simulate the anionic subsite of the receptor and the hydrophobic environment of the receptor cavity, respectively. The very unusual amide proton chemical shifts and their temperature coefficients suggest the presence of folded conformations in CDCl3 for all peptides, consistent with several models of opioid receptors and with the crystal structure of Leu5-enkephalin. The differences among the proposed cyclic conformations of the five peptides may be correlated, in part, with their different biological activity. All peptides in Me2SO-d6 are characterized by complex mixtures of extended fully solvated conformations.     

Quantum chemistry analysis of the mechanism of action of proteolytic enzymes. I. The state of the cleavable bond of substrates and intermediates

Electronic parameters of amide and ester bonds in some compounds, modelling substrates of proteolytic enzymes, and electronic properties of corresponding tetrahedral compounds, which are intermediates of the hydrolytic reaction, were calculated by the CNDO/2 method. The nature of substituents and the formation of the hydrogen bond by the carbonyl oxygen atom were shown to have no sufficient influence on the charges and bond orders of the amide group. The dramatic dependence of the amide electronic state from the distort degree of its planar structure was found. The resonance stabilization was shown to be absent in the bicyclic beta-lactams. The pK alpha values of the amide nitrogen atom were calculated at various hybridization states in amides.     

Synthesis and X-ray crystal structures of two transition metal complexes based on functionalised 1,5-anhydro-2-deoxy-D-galactitol and methyl 2-deoxy-alpha-D-galactopyranoside

Two new ligands of transition metal cations based on galactose-derived scaffolds were synthesised: 1,5-anhydro-2-deoxy-3,4,6-tri-O-(2-picolyl)-D-galactitol and methyl 2-deoxy-3,4,6-tri-O-(2-picolyl)-alpha-D-galactopyranoside. These ligands permitted the isolation as single crystals of a Co(II) and a Ni(II) complex, respectively. The structures of both complexes were determined by X-ray crystallography showing a coordination sphere including sugar-bound oxygen atoms. The sugar-derived ligands were found to be in both cases in high energy conformations in the crystal structures of the complexes. These conformations contain an arrangement of sugar-bound oxygen atoms similar to those observed in polyol-metal and carbohydrate-metal complexes.     

Vibrational CD of the amide II band in some model polypeptides and proteins.

The amide II vibrational CD (VCD) spectra of poly (L-glutamic acid) and poly (L-lysine) in various conformational forms and those of several proteins in H2O have been measured. Characteristic VCD patterns have been observed in the amide II region due to helix, beta-sheet, and coil conformations in polypeptides. Based on their x-ray crystal structures, the proteins studied have been assigned to six categories. Proteins in the same category give rise to similar amide II VCD. While the protein conformational type is indicated using the amide II VCD, discrimination between types is less characteristic than with the previously studied amide I' VCD in D2O.     

Experimental and calculated conformational characteristics of the bicyclic heptapeptide phalloidin

Several conformations generated from approximate potential energy calculations are presented for the bicyclic heptapeptide phalloidin which are consistent with the conformation-dependent information obtained from proton nuclear magnetic resonance measurements performed on phalloidin in dimethylsulfoxide solution. In each conformation, the cysteine amide proton is intramolecularly hydrogen bonded, the tryptophan amide is internally buried and the methyl group of the alanine residue preceding tryptophan is shielded by the tryptophan ring. Thus, phalloidin appears to be a relatively rigid molecule in solution.     

A comparative study on the catalytic properties of guanyl-specific ribonucleases.

The kinetic parameters of reactions catalyzed by four guanyl-specific RNases T1, Pb1, Th1 and Sa were studied comparatively using three types of substrates; guanosine-2',3'-cyclophosphates, GpN dinucleoside phosphates and synthetic polyribonucleotides. The kinetic parameters were shown to be similar in spite of considerable differences in primary structures of these RNases, including amino acid residues of the active sites. Therefore, primary structures of guanyl RNases allow for a considerable number of substitutions (both in the 'recognising' and catalytical parts of the active site) without changes in the catalytical parameters.     

Crystal structures of dialkylglycine decarboxylase inhibitor complexes.

The crystal structures of four inhibitor complexes of dialkylglycine decarboxylase are reported. The enzyme does not undergo a domain closure, as does aspartate aminotransferase, upon inhibitor binding. Two active-site conformations have been observed in previous structures that differ in alkali metal ion content, and two active-site conformations have been shown to coexist in solution when a single type of metal ion is present. There is no indication of coexisting conformers in the structures reported here or in the previously reported structures, and the observed conformation is that expected based on the presence of potassium in the enzyme. Thus, although two active-site conformations coexist in solution, a single conformation, corresponding to the more active enzyme, predominates in the crystal. The structure of 1-aminocyclopropane-1-carboxylate bound in the active site shows the aldimine double bond to the pyridoxal phosphate cofactor to be fully out of the plane of the coenzyme ring, whereas the Calpha-CO2(-) bond lies close to it. This provides an explanation for the observed lack of decarboxylation reactivity with this amino acid. The carboxylate groups of both 1-aminocyclopropane-1-carboxylate and 5'-phosphopyridoxyl-2-methylalanine interact with Ser215 and Arg406 as previously proposed. This demonstrates structurally that alternative binding modes, which constitute substrate inhibition, occur in the decarboxylation half-reaction. The structures of d and l-cycloserine bound to the active-site show that the l-isomer is deprotonated at C(alpha), presumably by Lys272, while the d-isomer is not. This difference explains the approximately 3000-fold greater potency of the l versus the d-isomer as a competitive inhibitor of dialkylglycine decarboxylase.     

Yeast tRNAAsp: codon and wobble codon-anticodon interactions. A transferred nuclear Overhauser enhancement study

The conformations of the ribotrinucleoside bisphosphates GpApC and GpApU, the codon and wobble codon for aspartic acid respectively, bound to yeast tRNAAsp in solution, have been examined by means of time-dependent transferred nuclear Overhauser enhancement measurements to determine distances between bound ligand protons. The conformations of the two bound ribotrinucleoside bisphosphates are shown to be very similar with an overall root-mean-square difference in interproton distances of 0.03 nm. The ribose conformations of all the residues are 3'-endo; the glycosidic bond torsion angles of the A and C residues of GpApC and of the A and U residues of GpApU are in the low anti range. These features are typical of an A-RNA type structure. In contrast, the G residue of both GpApC and GpApU exists as a mixture of syn and anti conformations. The overall conformation of the two bound ribotrinucleoside bisphosphates is also similar to A-RNA and the stability of the complexes is enhanced by extensive base-base stacking interactions. In addition, it is shown that the binding of the codon GpApC to tRNAAsp induces self-association into a multicomplex system consisting of four GpApC-tRNAAsp complexes, whereas the wobble codon GpApU fails to induce any observable self-association.     

DFT study of B-like conformations of deoxydinucleoside monophosphates containing Gua and/or Cyt and their complexes with Na+ cation

B-like minimum energy conformations of deoxydinucleoside monophosphate anions (dDMPs) containing Gua and/or Cyt and their Na+ complexes have been studied by the DFT PW91PW91/DZVP method. The optimized geometry of the dDMPs is in close agreement with experimental observations and the obtained minimum energy conformations are consistent with purine-purine, purine-pyrimidine, and pyrimidine-purine arrangements in crystals of B-DNA duplexes. All the studied systems are characterized by pyramidalization of the amino groups, which participate in the formation of unusual hydrogen bond between the carbonyl oxygen of the second base in the dGpdC, dCpdG dDMPs, and their Na+ complexes. In all the obtained structures the bases assume a nearly parallel disposition to each other and this effect is independent on the degree of their spatial superposition. From this it is concluded that the parallel disposition of the bases in the B-like single-stranded conformations is dictated by the sugar-phosphate backbone. Correspondingly, the base-base interactions attain a secondary role in the formation of these spatial structures. The formation of a weak C6-H6...O5' hydrogen bond between cytosine and the phosphate oxygen is reported, in agreement with experimental observations.     

The crystal structures of the complexes between bovine beta-trypsin and ten P1 variants of BPTI

The high-resolution X-ray structures have been determined for ten complexes formed between bovine beta-trypsin and P1 variants (Gly, Asp, Glu, Gln, Thr, Met, Lys, His, Phe, Trp) of bovine pancreatic trypsin inhibitor (BPTI). All the complexes were crystallised from the same conditions. The structures of the P1 variants Asp, Glu, Gln and Thr, are reported here for the first time in complex with any serine proteinase. The resolution of the structures ranged from 1.75 to 2.05 A and the R-factors were about 19-20 %. The association constants of the mutants ranged from 1.5x10(4) to 1.7x10(13) M-1. All the structures could be fitted into well-defined electron density, and all had very similar global conformations. All the P1 mutant side-chains could be accomodated at the primary binding site, but relative to the P1 Lys, there were small local changes within the P1-S1 interaction site. These comprised: (1) changes in the number and dynamics of water molecules inside the pocket; (2) multiple conformations and non-optimal dihedral angles for some of the P1 side-chains, Ser190 and Gln192; and (3) changes in temperature factors of the pocket walls as well as the introduced P1 side-chain. Binding of the cognate P1 Lys is characterised by almost optimal dihedral angles, hydrogen bonding distances and angles, in addition to considerably lower temperature factors. Thus, the trypsin S1 pocket seems to be designed particularly for lysine binding.