Dynamics of the aromatic amino acid residues in the globular conformation of the basic pancreatic trypsin inhibitor (BPTI). II. Semi-empirical energy calculations

Bioactives from medicinal herb against bedaquiline resistant tuberculosis: removing the dark clouds from the horizon

Tuberculosis is a contagious bacterial ailment that primarily affects the lungs and is brought on by the bacterium Mycobacterium tuberculosis (MTB). An antimycobacterial medication called bedaquiline (BQ) is specified to treat multidrug-resistant tuberculosis (MDR-TB). Despite its contemporary use in clinical practice, the mutations (D32 A/G/N/V/P) constrain the potential of BQ by causing transitions in the structural conformation of the atpE subunit-c after binding. In this study, we have taken the benzylisoquinoline alkaloids from thalictrum foliolosum due to its antimicrobial activity reported in prior literature. We used an efficient and optimized structure-based strategy to examine the wild type (WT) and mutated protein upon molecule binding. Our results emphasize the drastic decline in BQ binding affinity of mutant and WT atpE subunit-c complexes compared to thalirugidine (top hit) from thalictrum foliolosum. The decrease in BQ binding free energy is due to electrostatic energy because nearly every atom in a macromolecule harbors a partial charge, and molecules taking part in molecular recognition will interact electrostatically. Similarly, the high potential mean force of thalirugidine than BQ in WT and mutant complexes demonstrated the remarkable ability to eradicate mycobacteria efficiently. Furthermore, the Alamar blue cell viability and ATP determination assay were performed to validate the computational outcomes in search of novel antimycobacterial. Upon closer examination of the ATP determination assay, it became apparent that both BQ and thalirugidine showed similar reductions in ATP levels at their respective MICs, presenting a potential common mechanism of action.

NMR Studies of Aromatic Ring Flips to Probe Conformational Fluctuations in Proteins

Aromatic residues form a significant part of the protein core, where they make tight interactions with multiple surrounding side chains. Despite the dense packing of internal side chains, the aromatic rings of phenylalanine and tyrosine residues undergo 180° rotations, or flips, which are mediated by transient and large-scale "breathing" motions that generate sufficient void volume around the aromatic ring. Forty years after the seminal work by Wagner and Wüthrich, NMR studies of aromatic ring flips are now undergoing a renaissance as a powerful means of probing fundamental dynamic properties of proteins. Recent developments of improved NMR methods and isotope labeling schemes have enabled a number of advances in addressing the mechanisms and energetics of aromatic ring flips. The nature of the transition states associated with ring flips can be described by thermodynamic activation parameters, including the activation enthalpy, activation entropy, activation volume, and also the isothermal volume compressibility of activation. Consequently, it is of great interest to study how ring flip rate constants and activation parameters might vary with protein structure and external conditions like temperature and pressure. The field is beginning to gather such data for aromatic residues in a variety of environments, ranging from surface exposed to buried. In the future, the combination of solution and solid-state NMR spectroscopy together with molecular dynamics simulations and other computational approaches is likely to provide detailed information about the coupled dynamics of aromatic rings and neighboring residues. In this Perspective, we highlight recent developments and provide an outlook toward the future.

OCRE Domains of Splicing Factors RBM5 and RBM10: Tyrosine Ring-Flip Frequencies Determined by Integrated Use of 1 H NMR Spectroscopy and Molecular Dynamics Simulations

The 55-residue OCRE domains of the splicing factors RBM5 and RBM10 contain 15 tyrosines in compact, globular folds. At 25 °C, all 15 tyrosines show symmetric ¹ H NMR spectra, with averaged signals for the pairs of δ- and ϵ-ring hydrogens. At 4 °C, two tyrosines were identified as showing ¹ H NMR line-broadening due to lowered frequency of the ring-flipping. For the other 13 tyrosine rings, it was not evident, from the ¹ H NMR data alone, whether they were either all flipping at high frequencies, or whether slowed flipping went undetected due to small chemical-shift differences between pairs of exchanging ring hydrogen atoms. Here, we integrate ¹ H NMR spectroscopy and molecular dynamics (MD) simulations to determine the tyrosine ring-flip frequencies. In the RBM10-OCRE domain, we found that, for 11 of the 15 tyrosines, these frequencies are in the range 2.0×10⁶ to 1.3×10⁸  s^-1 , and we established an upper limit of <1.0×10⁶  s^-1 for the remaining four residues. The experimental data and the MD simulation are mutually supportive, and their combined use extends the analysis of aromatic ring-flip events beyond the limitations of routine ¹ H NMR line-shape analysis into the nanosecond frequency range.

Protein crystal lattices are dynamic assemblies: the role of conformational entropy in the protein condensed phase

Until recently, the occurrence of conformational entropy in protein crystal contacts was considered to be a very unlikely event. A study based on the most accurately refined protein structures demonstrated that side-chain conformational entropy and static disorder might be common in protein crystal lattices. The present investigation uses structures refined using ensemble refinement to show that although paradoxical, conformational entropy is likely to be the major factor in the emergence and integrity of the protein condensed phase. This study reveals that the role of shape entropy and local entropic forces expands beyond the onset of crystallization. For the first time, the complete pattern of intermolecular interactions by protein atoms in crystal lattices is presented, which shows that van der Waals interactions dominate in crystal formation.

Dynamics of Hydrophobic Core Phenylalanine Residues Probed by Solid-State Deuteron NMR

We conducted a detailed investigation of the dynamics of two phenylalanine side chains in the hydrophobic core of the villin headpiece subdomain protein (HP36) in the hydrated powder state over the 298-80 K temperature range. Our main tools were static deuteron NMR measurements of longitudinal relaxation and line shapes supplemented with computational modeling. The temperature dependence of the relaxation times reveals the presence of two main mechanisms that can be attributed to the ring-flips, dominating at high temperatures, and small-angle fluctuations, dominating at low temperatures. The relaxation is nonexponential at all temperatures with the extent of nonexponentiality increasing from higher to lower temperatures. This behavior suggests a distribution of conformers with unique values of activation energies. The central values of the activation energies for the ring-flipping motions are among the smallest reported for aromatic residues in peptides and proteins and point to a very mobile hydrophobic core. The analysis of the widths of the distributions, in combination with the earlier results on the dynamics of flanking methyl groups (Vugmeyster et al. J. Phys. Chem. B 2013, 117, 6129-6137), suggests that the hydrophobic core undergoes slow concerted fluctuations. There is a pronounced effect of dehydration on the ring-flipping motions, which shifts the distribution toward more rigid conformers. The crossover temperature between the regions of dominance of the small-angle fluctuations and ring-flips shifts from 195 K in the hydrated protein to 278 K in the dry one. This result points to the role of solvent in softening the core and highlights aromatic residues as markers of the protein dynamical transitions.

Solvents, interfaces and protein structure

Mean packing densities in protein interiors are comparable to those of most organic solids but the variations between small regions may be substantial. Packing defects may be related to allowed structural fluctuations. Molecular surface areas can be correlated with free energies of transfer between different solvents. The proportionality factor will depend, in general, on the nature of the solute and both the solvents. The changes in solvent-protein interfacial area on chain folding are large and the implied changes in free energy from this solvent-squeezing effect are correspondingly large. The strong tendency to minimize surface area is reflected in the globular shape of most protein molecules or domains in larger structures. The formation of isolated units of secondary structure from an extended chain represents about one half of the eventual total area change. The tendencies of amino acids to form beta-sheets correlate well with the rank-ordered list based on non-polar area change for each residue type. The calculated area changes for helix and sheet formation are not identical in rank order. The rank-ordered list for alpha-helix formation correlates satisfactorily with the probability list prepared from actual structures if glutamic acid and tyrosine are removed. What special characteristics unrelated to surface area these two amino acids might have is not clear. Tertiary structure formation from preformed secondary structural units can be rank ordered on area change and possible nucleation sites can be identified. A prediction scheme for helix-helix interactions is proposed. The hydrophobic force begins to be felt when two helices are about 0.6 nm (6 A) from their final contact positions. Interfacial surface tension is a logical parameter to relate free energy and solvent contact area, but this macroscopic parameter must be used with great caution. It is suggested that water in the deep grooves, characteristic of the active sites of many enzymes, may have a substantially higher fugacity than bulk water as indicated, at least qualitatively, by the Kelvin equation based on surface curvature. Such water would be more easily displaced than its plane surface counterpart and could contribute significantly to ligand-binding energy. This factor would be in addition to the usual solvent entropic effects associated with surface area reduction on association.

Physical and chemical properties of lysozyme

The conformations of lysozyme in crystals and in aqueous solution are discussed and it is shown that the basic conformation is similar in the two states. Certain parts of the molecule have mobility. The reactions of lysozyme with protons, metal ions and some organic reagents are examined in the light of the conformations and their dynamics. The reactions considered are mainly those of tyrosyl, tryptophyl and carboxylate residues. The reactivity data are used in a discussion of the energy states of the reacting side-chains. In particular the reactivity of Glu-35 and its interaction with Trp-108 lead to suggestions for some new aspects in the hypothesis for the mechanism of action of lysozyme. In most respects the X-ray crystal diffraction and the nuclear magnetic resonance solution studies are in accord and complementary.

Six years of protein structure determination by NMR spectroscopy: what have we learned?

Nuclear magnetic resonance (NMR) spectroscopy in solution is a second technique, in addition to X-ray diffraction in single crystals, for the determination of three-dimensional protein structures at atomic resolution. Structures of proteins derived by NMR have now been with us for six years, and here I entertain the following question: what information have we gained that would not be available if X-ray crystallography were still the only method for protein structure determination? Answers include that NMR structures are available of proteins that have not been crystallized, that the two techniques afford different insights into internal mobility of proteins, and that one gets different views of protein hydration and hence the molecular surface when using NMR spectroscopy or X-ray diffraction.

Infrequent cavity-forming fluctuations in HPr from Staphylococcus carnosus revealed by pressure- and temperature-dependent tyrosine ring flips

Infrequent structural fluctuations of a globular protein is seldom detected and studied in detail. One tyrosine ring of HPr from Staphylococcus carnosus, an 88-residue phosphocarrier protein with no disulfide bonds, undergoes a very slow ring flip, the pressure and temperature dependence of which is studied in detail using the on-line cell high-pressure nuclear magnetic resonance technique in the pressure range from 3 MPa to 200 MPa and in the temperature range from 257 K to 313 K. The ring of Tyr6 is buried sandwiched between a beta-sheet and alpha-helices (the water-accessible area is less than 0.26 nm2), its hydroxyl proton being involved in an internal hydrogen bond. The ring flip rates 10(1)-10(5) s(-1) were determined from the line shape analysis of H(delta1, delta2) and H(epsilon1,epsilon2) of Tyr6, giving an activation volume DeltaV++ of 0.044 +/- 0.008 nm3 (27 mL mol(-1)), an activation enthalpy DeltaH++ of 89 +/- 10 kJ mol(-1), and an activation entropy DeltaS++ of 16 +/- 2 JK(-1) mol(-1). The DeltaV++) and DeltaH++ values for HPr found previously for Tyr and Phe ring flips of BPTI and cytochrome c fall within the range of DeltaV(double dagger) of 28 to 51 mL mol(-1) and DeltaH++ of 71 to 155 kJ mol(-1). The fairly common DeltaV++ and DeltaH++ values are considered to represent the extra space or cavity required for the ring flip and the extra energy required to create a cavity, respectively, in the core part of a globular protein. Nearly complete cold denaturation was found to take place at 200 MPa and 257 K independently from the ring reorientation process.

Molecular and biological constraints on ligand-binding affinity and specificity

Interactions between biological macromolecules have characteristic values of affinity and specificity that are set according to the biological function that is served by the interaction in the organism. Here we examine the molecular mechanisms that are used to achieve the required values of affinity and specificity in various biological systems.

Characterization of rates of ring-flipping in trimethoprim in its ternary complexes with Lactobacillus casei dihydrofolate reductase and coenzyme analogues

NMR measurements have been used to investigate rates of ring-flipping and the activation parameters for the trimethoxybenzyl ring of the antibacterial drug trimethoprim (TMP) bound to Lactobacillus casei dihydrofolate reductase (DHFR) for a series of ternary complexes formed with analogues of the coenzyme NADPH. Rates were obtained at several temperatures from line shape analyses ((13)C-edited HSQC (1)H spectra) and transfer of magnetization measurements (zz-HSQC) on complexes containing 3'-O-[(13)C]trimethoprim. Examination of the structures of the complexes indicates that ring-flipping can only be achieved following major conformational changes and transient fluctuations of the protein and coenzyme structure around the trimethoxybenzyl ring. There is no simple correlation between rates of ring-flipping and binding constants. The presence of the coenzyme nicotinamide ring (in either its reduced or its oxidized forms) in the binding site close to the trimethoxybenzyl ring moiety is the major factor reducing the ring-flipping on coenzyme binding. Thus, the ternary complex with NADPH shows the largest reduction in the rate of ring-flipping (11 +/- 3 s(-)(1) at 298 K) as compared with the binary complex (793 +/- 80 s(-)(1) at 298 K). Complexes with NADPH analogues that either have no nicotinamide ring or are known to have their nicotinamide rings removed from the binding site show the smallest reductions. For the DHFR.TMP.NADP(+) complex where there are two conformations present, very different rates of ring-flipping were observed for the two forms. The activation parameters (DeltaH() and DeltaS()) for the ring-flipping in all the complexes are discussed in terms of the protein-ligand interactions and the possible constraints on the pathway through the transition state.

19F-n.m.r. studies of 3',5'-difluoromethotrexate binding to Lactobacillus casei dihydrofolate reductase. Molecular motion and coenzyme-induced conformational changes

19F-n.m.r. spectroscopy was used to study the binding of 3',5'-difluoromethotrexate to dihydrofolate reductase (tetrahydrofolate dehydrogenase) from Lactobacillus casei. The benzoyl ring of the bound difluoromethotrexate was found to 'flip' about its symmetry axis, and the rate (7.3 X 10(3) s-1 at 298 K) and activation parameters for this process were determined by lineshape analysis of the 19F-n.m.r. spectrum at a series of temperatures in the range 273-308 K. The contributions to the barrier for this process are discussed. Addition of NADP+ or NADPH to form the enzyme-difluoromethotrexate-coenzyme ternary complex led to an increase in the rate of benzoyl ring flipping by a factor of 2.6-2.8-fold, and to substantial changes in the 19F-n.m.r. chemical shifts. The possible nature of the coenzyme-induced conformational changes responsible for these effects is discussed.

Structural and functional aspects of domain motions in proteins

Three distinct categories of large-scale flexibility in proteins have been documented by single-crystal X-ray diffraction studies: the relatively free movement of essentially rigid globular domains that are connected by a flexible segment of polypeptide, the reorientation of essentially rigid domains among a few distinct conformations, and the concerted transition of a contiguous region of the surface of a protein from a disordered state to an ordered state. In a number of examples, well-defined functions can be assigned to these large-scale structural changes. The occurrence of such motions in proteins of known structure is reviewed, and the best-studied examples are discussed in detail to allow a critical evaluation of the methods used to identify and study these motions.

Hydrogen exchange and structural dynamics of proteins and nucleic acids

Though the structures presented in crystallographic models of macromolecules appear to possess rock-like solidity, real proteins and nucleic acids are not particularly rigid. Most structural work to date has centred upon the native state of macromolecules, the most probable macromolecular form. But the native state of a molecule is merely its most abundant form, certainly not its only form. Thermodynamics requires that all other possible structural forms, however improbable, must also exist, albeit with representation corresponding to the factor exp( —Gi/RT) for each state of free energyGi(see Moelwyn-Hughes, 1961), and one appreciates that each molecule within a population of molecules will in time explore the vast ensemble ofpossiblestructural states.

Characterization of the distribution of internal motions in the basic pancreatic trypsin inhibitor using a large number of internal NMR probes

The experimental observations described in this article indicated that a distribution of many different fluctuations is present in a globular protein. These fluctuations were characterized by observation of many natural internal probes such as the labile peptide protons and the aromatic side chains. The conditions which are necessary to get reactions of the internal probes have been discussed in detail. The structural interpretation of the data was facilitated by the development and the use of new NMR techniques which provided the identification of the resonances of all the labile peptide protons. With NOE measurements a distinction between correlated and uncorrelated exchange events was obtained. This enabled us to elucidate the exchange mechanism over a wide range of p2H and temperature and to classify different subsets of fluctuations with respect to their lifetimes. It was further demonstrated that a change of external conditions such as temperature, p2H or pressure can change the distribution of fluctuations in the protein. The mechanisms responsible for rotation of internal aromatic side chains were also found to change with temperature, and mechanistic aspects of these fluctuations were discussed. This demonstration of a manifold of spatial fluctuations in a small protein provides an impression on the kind of fluctuations which have to be expected for larger proteins. When studying protein reactions one should therefore consider the presence of a large number of different, transiently formed, spatial structures available for the partner in the reaction, which may pick out only that structure which will optimally perform a particular reaction with the highest efficiency.

Dynamical theory of activated processes in globular proteins

A method is described for calculating the reaction rate in globular proteins of activated processes such as ligand binding or enzymatic catalysis. The method is based on the determination of the probability that the system is in the transition state and of the magnitude of the reactive flux for transition-state systems. An "umbrella sampling" simulation procedure is outlined for evaluating the transition-state probability. The reactive flux is obtained from an approach described previously for calculating the dynamics of transition-state trajectories. An application to the rotational isomerization of an aromatic ring in the bovine pancreatic trypsin inhibitor is presented. The results demonstrate the feasibility of calculating rate constants for reactions in proteins and point to the importance of solvent effects for reactions that occur near the protein surface.

Structural interpretation of vicinal proton-proton coupling constants 3JH alpha H beta in the basic pancreatic trypsin inhibitor measured by two-dimensional J-resolved NMR spectroscopy

Two-dimensional J-resolved 1H NMR spectroscopy was used to measure the vicinal spin-spin coupling constants 3JH alpha H beta for numerous, previously individually assigned amino acid residues in the basic pancreatic trypsin inhibitor at various temperatures between 30 and 85 degrees JC. An analysis of this data is proposed which enables one to compare the spatial arrangements of individual amino acid side chains in solution and in single crystals of the protein, and which also provides information on the mobility of the side chains in the solution conformation. As a rule, the amino acid side chains in the interior of the protein were found to be locked into unique spatial orientations, with the mobility restricted to rapid rotational fluctuations about this unique value for the dihedral angle chi 1. In most, but not all, instances the data for the interior amino acids indicate identical average conformations for the amino acid side chains in single crystals and in solution. For residues on the protein surface structural rearrangements between crystal and solution appear to be common, and the mobility in the solution conformation may include rapid averaging between two or several distinct, preferentially populated values of chi 1, analogous to the gauche-trans-gauche isomerization in isolated amino acids.

Correlations between internal mobility and stability of globular proteins

The recent work is surveyed which leads to the suggestions that the conformation of globular proteins in solution corresponds to a dynamic ensemble of rapidly interconverting spatial structures, that clusters of hydrophobic amino acid side chains have an important role in the architecture of protein molecules, and that mechanistic aspects of protein denaturation can be correlated with internal mobility seen in the native conformation. These conclusions resulted originally from high resolution 1H nuclear magnetic resonance (NMR) studies of aromatic ring mobility, exchange of interior amide protons and thermal denaturation of the basic pancreatic trypsin inhibitor and a group of related proteins. Various new approaches to further characterize proteins in solution have now been taken and preliminary data are presented. These include computer graphics to outline hydrophobic clusters in globular protein structures, high resolution 1H-NMR experiments at variable hydrostatic pressure and 13C-NMR relaxation measurements. At the present early stage of these new investigations it appears that the hydrophobic cluster model for globular proteins is compatible with the data obtained.

Internal dynamics of proteins. Short time and long time motions of aromatic sidechains in PTI

Theoretical approaches to the internal dynamics of proteins are outlined and illustrated by application to the aromatic sidechain motions of tyrosines in the bovine pancreatic trypsin inhibitor. High frequency torsional oscillations are obtained from a molecular dynamics simulation, while the longer time ring rotations are analyzed by use of adiabatic energy minimization and special transition-state trajectory techniques.

Nuclear-magnetic-resonance studies of ferrocytochrome c. pH and temperature dependence

The pH dependence and the temperature dependence of the nuclear magnetic resonance spectrum of horse ferrocytochrome c are described. This protein is very stable; it maintains an ordered structure over the pH range 4 to 12 at 25 degrees C and over the temperature range 4 degrees C to 97 degrees C at pH 7.0. The dynamic characteristics of the conformation of ferrocytochrome c were investigated. Particular emphasis was laid on the aromatic resonances and resonances of methyl groups shifted far upfield. Tyr-48 and Phe-46 were found to be relatively immobile whilst a region of the protein close to Ile-57 was found to be relatively flexible.

Side-chain torsional potentials: effect of dipeptide, protein, and solvent environment

Side-chain torsional potentials in the bovine pancreatic trypsin inhibitor are calculated from empirical energy functions by use of the known X-ray structure of the protein and the rigid-geometry mapping technique. The potentials are analyzed to determine the roles and relative importance of contributions from the dipeptide backbone, the protein, and the crystalline environment of solvent and other protein molecules. The structural characteristics of the side chains determine two major patterns of energy surfaces, E(X1,X2): a gamma-branched pattern and a pattern for longer, straight side chains (Arg, Lys, Glu, and Met). Most of the dipeptide potential curves and surfaces have a local minimum corresponding to the side-chain torsional angles in the X-ray structure. Addition of the protein forces sharpens and/or selects from these minima, providing very good agreement with the experimental conformation for most side chains at the surface or in the core of the protein. Inclusion of the crystalline environment produces still better results, especially for the side chains extending away from the protein. The results are discussed in terms of the details of the interactions due to the surrounding, calculated solvent-accessibility figures and the temperature factors derived from the crystallographic refinement of the pancreatic trypsin inhibitor.

Individual assignments of amide proton resonances in the proton NMR spectrum of the basic pancreatic trypsin inhibitor

Studies of proton-proton nuclear Overhauser effects were used to obtain individual assignments of 17 amide proton resonances in the 360 MHz proton nuclear magnetic resonance spectrum of the basic pancreatic trypsin inhibitor. First, optimizing the conditions for obtaining selective nuclear Overhauser effects in the presence of spin diffusion in macromolecules is discussed. Truncated driven nuclear Overhauser experiments were used to assing the amide proton resonances of the beta-sheet in the inhibitor. It is suggested that these techniques could serve quite generally to obtain individual resonance assignments in beta-sheet secondary structures of proteins. Combination of nuclear Overhauser studies with spin decoupling further resulted in individual assignments of the gamma-methyl resonances of the two isoleucines and numerous Calpha and Cbeta protons.

Ring current effects in the conformation dependent NMR chemical shifts of aliphatic protons in the basic pancreatic trypsin inhibitor

Previously, a highly refined crystal structure and energy refined atomic coordinates were obtained for the basic pancreatic trypsin inhibitor, as well as numerous individual resonance assignments in the 1H NMR spectrum. These data were now used to investigate the contributions from the local ring current fields of the aromatic rings to the overall conformation dependent chemical shifts in this globular protein. A program was written which allowed the consideration of certain aspects of internal mobility of the protein, and the different commonly used ring current equa tions were compared. These studies indicate that ring current shifts are the dominant contribution to the observed conformation dependent chemical shifts of the peripheral aliphatic side chain protons. On the other hand, it appears that ring current shifts do not make dominant contributions to the conformation dependent shifts of the backbone alpha- and amide protons or the aromatic protons in the inhibitor. On the basis of the empirical calibration with the peripheral aliphatic side chain protons, the Johnson-Bovey ring current equation was selected for an analysis of the ring geometries of two prolines in the inhibitor.

Structural interpretation of lanthanide binding to the basic pancreatic trypsin inhibitor by 1H NMR at 360 MHz

The weak binding of lanthanides to the five carboxyl groups of the basic pancreatic trypsin inhibitor (hereafter termed "the inhibitor"), has been investigated in detail using high resolution 1H NMR at 360 MHz. Lanthanides bind to the C-terminus with an apparent binding constant of 30 M-1, and thus competitively inhibit the formation of a salt-bridge between the C-terminus and the N-terminus, Lanthanides bind also to the side chain carboxyl groups of Asp 3, Glu 7, Glu 49 and Asp 50, with binding constants of 10--30 M-1. With the use of lanthanides individual resonance assignments for Phe 4 and Phe 45 were obtained in the 1H NMR spectrum of the inhibitor, and for several spin systems previous identifications were independently confirmed. The present experiments also provide a nice illustration for the use of shift reagents to improve the resolution in 1H NMR spectra of proteins. The exchange broadening for Tyr 35 and Phe 45 over the temperature range 4--72 degrees C could thus be observed for almost all the components of these aromatic spin systems and new details on the dynamic properties were obtained also for other aromatic residues.

A 1H nuclear-magnetic-resonance study of the solution conformation of the isoinhibitor K from Helix pomatia

The isoinhibitor K from Helix pomatia was investigated by high-resolution 1H nuclear magnetic resonance at 360 MHZ. Detailed studies of the labile protons and the resonances of the aromatic residues indicated extensive homologies between the spatial structures of the snail inhibitor, the basic pancreatic trypsin inhibitor (Kunitz) from bovine tissue and the cos colostrum trypsin inhibitor. Comparison of these three homologous inhibitors indicated that the overall flexibility of the globular protein conformation is reduced and its stability with respect to thermal denaturation raised when the content of amino acids with charged side chains is increased. It is suggested that this relation between amino acid composition and stability of the globular solution structure might be valid also for other classes of homologous proteins.

Individual assignments of the methyl resonances in the 1H nuclear magnetic resonance spectrum of the basic pancreatic trypsin inhibitor

In earlier work the resonances of the 20 methyl groups in the basic pancreatic trypsin inhibitor (BPTI) had been identified in the 360-MHz 1H nuclear magnetic resonance (NMR) spectra and most of the methyl lines had from spin-decoupling experiments been assigned to the different types of amino acid residues. The assignments to the different amino acid types were now completed by studies of the saturation transfer between the denatured and the globular forms of the inhibitor and by spin-decoupling experiments in nuclear Overhauser enhancement (NOE) difference spectra. These distinguished between the methyl resonances of Ala and Thr. Furthermore, for most of the methyl resonances, individual assignments to specific residues in the amino acid sequence were obtained from measurements of intramolecular proton-proton NOE's, use of lanthanide NMR shift and relaxation probes, and comparative studies of various chemically modified forms of BPTI. These data provide the basis for individual assignments of the methyl 13C NMR lines in BPTI and for detailed investigations of the relations between the spatial structure of the protein and the chemical shifts of the methyl groups. The methyl groups in BPTI are of particular interest since they are located almost exclusively on the surface of the protein and thus represent potential natural NMR probes for studies of the protein-protein interactions in the complexes formed between BPTI and a variety of proteases.