Molecular dynamics study of major urinary protein-pheromone interactions: a structural model for ligand-induced flexibility increase

Recently, two independent (15)N NMR relaxation studies indicated that in contrast to the decreased flexibility expected for induced-fit interactions, the backbone flexibility of major urinary protein isoform I (MUP-I) slightly increased upon complex formation with its natural pheromone 2-sec-butyl-4,5-dihydrothiazol. We have investigated the subtle details of molecular interactions by molecular dynamics simulations in explicit solvent. The calculated order parameters S(2) for a free- and ligand-bound protein supply evidence that mobility in various regions of MUP-I can be directly related to small conformational changes of the free- and complexed protein resulting from modifications of the hydrogen bonding network.     

Temperature-dependent spectral density analysis applied to monitoring backbone dynamics of major urinary protein-I complexed with the pheromone 2-sec-butyl-4,5-dihydrothiazole*

Backbone dynamics of mouse major urinary protein I (MUP-I) was studied by (15)N NMR relaxation. Data were collected at multiple temperatures for a complex of MUP-I with its natural pheromonal ligand, 2- sec -4,5-dihydrothiazole, and for the free protein. The measured relaxation rates were analyzed using the reduced spectral density mapping. Graphical analysis of the spectral density values provided an unbiased qualitative picture of the internal motions. Varying temperature greatly increased the range of analyzed spectral density values and therefore improved reliability of the analysis. Quantitative parameters describing the dynamics on picosecond to nanosecond time scale were obtained using a novel method of simultaneous data fitting at multiple temperatures. Both methods showed that the backbone flexibility on the fast time scale is slightly increased upon pheromone binding, in accordance with the previously reported results. Zero-frequency spectral density values revealed conformational changes on the microsecond to millisecond time scale. Measurements at different temperatures allowed to monitor temperature dependence of the motional parameters.     

Thermodynamic analysis of binding between mouse major urinary protein-I and the pheromone 2-sec-butyl-4,5-dihydrothiazole

The mouse pheromone 2-sec-butyl-4,5-dihydrothiazole (SBT) binds to an occluded, nonpolar cavity in the mouse major urinary protein-I (MUP-I). The thermodynamics of this interaction have been characterized using isothermal titration calorimetry (ITC). MUP-I-SBT binding is accompanied by a large favorable enthalpy change (DeltaH = -11.2 kcal/mol at 25 degrees C), an unfavorable entropy change (-TDeltaS = 2.8 kcal/mol at 25 degrees C), and a negative heat capacity change [DeltaC(p)() = -165 cal/(mol K)]. Thermodynamic analysis of binding between MUP-I and several 2-alkyl-4,5-dihydrothiazole ligands indicated that the alkyl chain contributes more favorably to the enthalpy and less favorably to the entropy of binding than would be expected on the basis of the hydrophobic desolvation of short-chain alcohols. However, solvent transfer experiments indicated that desolvation of SBT is accompanied by a net unfavorable change in enthalpy (DeltaH = +1.0 kcal/mol) and favorable change in entropy (-TDeltaS = -1.8 kcal/mol). These results are discussed in terms of the possible physical origins of the binding thermodynamics, including (1) hydrophobic desolvation of both the protein and the ligand, (2) formation of a buried water-mediated hydrogen bond network between the protein and ligand, (3) formation of strong van der Waals interactions, and (4) changes in the structure, dynamics, and/or hydration of the protein upon binding.     

Structural basis of pheromone binding to mouse major urinary protein (MUP-I)

The mouse major urinary proteins are pheromone-binding proteins that function as carriers of volatile effectors of mouse physiology and behavior. Crystal structures of recombinant mouse major urinary protein-I (MUP-I) complexed with the synthetic pheromones, 2-sec-butyl-4,5-dihydrothiazole and 6-hydroxy-6-methyl-3-heptanone, have been determined at high resolution. The purification of MUP-I from mouse liver and a high-resolution structure of the natural isolate are also reported. These results show the binding of 6-hydroxy-6-methyl-3-heptanone to MUP-I, unambiguously define ligand orientations for two pheromones within the MUP-I binding site, and suggest how different chemical classes of pheromones can be accommodated within the MUP-I beta-barrel.     

The binding cavity of mouse major urinary protein is optimised for a variety of ligand binding modes

¹⁵N and ¹HN chemical shift data and ¹⁵N relaxation studies have been used to characterise the binding of N-phenyl-naphthylamine (NPN) to mouse major urinary protein (MUP). NPN binds in the β-barrel cavity of MUP, hydrogen bonding to Tyr120 and making extensive non-bonded contacts with hydrophobic side chains. In contrast to the natural pheromone 2-sec-butyl-4,5-dihydrothiazole, NPN binding gives no change to the overall mobility of the protein backbone of MUP. Comparison with 11 different ligands that bind to MUP shows a range of binding modes involving 16 different residues in the β-barrel cavity. These finding justify why MUP is able to adapt to allow for many successful binding partners.     

Increased protein backbone conformational entropy upon hydrophobic ligand binding

For complexes between proteins and very small hydrophobic ligands, hydrophobic effects alone may be insufficient to outweigh the unfavorable entropic terms resulting from bimolecular association. NMR relaxation experiments indicate that the backbone flexibility of mouse major urinary protein increases upon binding the hydrophobic mouse pheromone 2-sec-butyl-4,5-dihydrothiazole. The associated increase in backbone conformational entropy of the protein appears to make a substantial contribution toward stabilization of the protein-pheromone complex. This term is likely comparable in magnitude to other important free energy contributions to binding and may represent a general mechanism to promote binding of very small ligands to macromolecules.     

Assignments of backbone 1H, 13C, and 15N resonances and secondary structure of ribonuclease H from Escherichia coli by heteronuclear three-dimensional NMR spectroscopy.

The assignments of individual magnetic resonances of backbone nuclei of a larger protein, ribonuclease H from Escherichia coli, which consists of 155 amino acid residues and has a molecular mass of 17.6 kDa are presented. To remove the problem of degenerate chemical shifts, which is inevitable in proteins of this size, three-dimensional NMR was applied. The strategy for the sequential assignment was, first, resonance peaks of amides were classified into 15 amino acid types by 1H-15N HMQC experiments with samples in which specific amino acids were labeled with 15N. Second, the amide 1H-15N peaks were connected along the amino acid sequence by tracing intraresidue and sequential NOE cross peaks. In order to obtain unambiguous NOE connectivities, four types of heteronuclear 3D NMR techniques, 1H-15N-1H 3D NOESY-HMQC, 1H-15N-1H 3D TOCSY-HMQC, 13C-1H-1H 3D HMQC-NOESY, and 13C-1H-1H 3D HMQC-TOCSY, were applied to proteins uniformly labeled either with 15N or with 13C. This method gave a systematic way to assign backbone nuclei (N, NH, C alpha H, and C alpha) of larger proteins. Results of the sequential assignments and identification of secondary structure elements that were revealed by NOE cross peaks among backbone protons are reported.     

1H, 13C and 15N backbone assignments of cyclophilin when bound to cyclosporin A (CsA) and preliminary structural characterization of the CsA binding site.

The backbone 1H, 13C and 15N chemical shifts of cyclophilin (CyP) when bound to cyclosporin A (CsA) have been assigned from heteronuclear two- and three-dimensional NMR experiments involving selectively 15N- and uniformly 15N- and 15N,13C-labeled cyclophilin. From an analysis of the 1H and 15N chemical shifts of CyP that change upon binding to CsA and from CyP/CsA NOEs, we have determined the regions of cyclophilin involved in binding to CsA.     

15N-labeled tRNA. Identification of dihydrouridine in Escherichia coli tRNAfMet, tRNALys, and tRNAPhe by 1H-15N two-dimensional NMR

The N3 imino units of dihydrouridine were identified in samples of 15N-labeled Escherichia coli tRNAfMet, tRNALys, and tRNAPhe by 1H-15N two-dimensional NMR. The peaks for dihydrouridine had high field 1H (9.7-9.8 ppm) and 15N (147.8-149.5 ppm) chemical shifts. Assignments were made by 1H-15N chemical shift correlation based on values obtained in model studies with tri-O-benzoyl- and tri-O-acetyldihydrouridine. The rates of exchange of the imino protons with water suggest that the D-loop in tRNAfMet is less stable than the D-loops in tRNALys or tRNAPhe. Closely spaced peaks were observed for the two dihydrouridines in tRNAPhe in a high resolution spectrum.     

Ligand-induced changes in dynamics in the RT loop of the C-terminal SH3 domain of Sem-5 indicate cooperative conformational coupling

We report the effects of peptide binding on the (15)N relaxation rates and chemical shifts of the C-SH3 of Sem-5. (15)N spin-lattice relaxation time (T(1)), spin-spin relaxation time (T(2)), and ((1)H)-(15)N NOE were obtained from heteronuclear 2D NMR experiments. These parameters were then analyzed using the Lipari-Szabo model free formalism to obtain parameters that describe the internal motions of the protein. High-order parameters (S(2) > 0.8) are found in elements of regular secondary structure, whereas some residues in the loop regions show relatively low-order parameters, notably the RT loop. Peptide binding is characterized by a significant decrease in the (15)N relaxation in the RT loop. Concomitant with the change in dynamics is a cooperative change in chemical shifts. The agreement between the binding constants calculated from chemical shift differences and that obtained from ITC indicates that the binding of Sem-5 C-SH3 to its putative peptide ligand is coupled to a cooperative conformational change in which a portion of the binding site undergoes a significant reduction in conformational heterogeneity.     

Dynamic conformational equilibria in the physiological function of the Bombyx mori pheromone-binding protein

The Bombyx mori pheromone-binding protein (BmorPBP) undergoes a pH-dependent conformational transition from a form at basic pH, which contains an open cavity suitable for ligand binding (BmorPBP^B), to a form at pH 4.5, where this cavity is occupied by an additional helix (BmorPBP^A). This helix α7 is formed by the C-terminal dodecapeptide 131-142, which is flexibly disordered on the protein surface in BmorPBP^B and in its complex with the pheromone bombykol. Previous work showed that the ligand-binding cavity cannot accommodate both bombykol and helix α7. Here we further investigated mechanistic aspects of the physiologically crucial ejection of the ligand at lower pH values by solution NMR studies of the variant protein BmorPBP(1-128), where the C-terminal helix-forming tetradecapeptide is removed. The NMR structure of the truncated protein at pH 6.5 corresponds closely to BmorPBP^B. At pH 4.5, BmorPBP(1-128) maintains a B-type structure that is in a slow equilibrium, on the NMR chemical shift timescale, with a low-pH conformation for which a discrete set of ¹⁵N-¹H correlation peaks is NMR unobservable. The full NMR spectrum was recovered upon readjusting the pH of the protein solution to 6.5. These data reveal dual roles for the C-terminal tetradecapeptide of BmorPBP in the mechanism of reversible pheromone binding and transport, where it governs dynamic equilibria between two locally different protein conformations at acidic pH and competes with the ligand for binding to the interior cavity.     

Active site contacts in the purine nucleoside phosphorylase--hypoxanthine complex by NMR and ab initio calculations

Hypoxanthine (Hx) with specific (15)N labels has been used to probe hydrogen-bonding interactions with purine nucleoside phosphorylase (PNP) by NMR spectroscopy. Hx binds to human PNP as the N-7H tautomer, and the N-7H (1)H and (15)N chemical shifts are located at 13.9 and 156.5 ppm, respectively, similar to the solution values. In contrast, the (1)H and (15)N chemical shifts of N-1H in the PNP.Hx complex are shifted downfield by 3.5 and 7.5 ppm to 15.9 and 178.8 ppm, respectively, upon binding. Thus, hydrogen bonding at N-1H is stronger than at N-7H in the complex. Ab initio chemical shift calculations on model systems that simulate Hx in solution and bound to PNP are used to interpret the NMR data. The experimental N-7H chemical shift changes are caused by competing effects of two active site contacts. Hydrogen bonding of Glu201 to N-1H causes upfield shifts of the N-7H group, while the local hydrogen bond (C=O to N-7H from Asn243) causes downfield shifts. The observed N-7H chemical shift can be reproduced by a hydrogen bond distance approximately 0.13 A shorter (but within experimental error) of the experimental value found in the X-ray crystal structure of the bovine PNP.Hx complex. The combined use of NMR and ab initio chemical shift computational analysis provides a novel approach to understand enzyme-ligand interactions in PNP, a target for anticancer agents. This approach has the potential to become a high-resolution tool for structural determination.     

15N NMR assignments and chemical shift analysis of uniformly labeled 15N calbindin D9k in the apo, (Cd2+)1 and (Ca2+)2 states.

15N has been uniformly incorporated into the EF-hand Ca(2+)-binding protein calbindin D9k so that heteronuclear experiments can be used to further characterize the structure and dynamics of the apo, (Cd2+)1 and (Ca2+)2 states of the protein. The 15N NMR resonances were assigned by 2D 15N-resolved 1H experiments, which also allowed the identification of a number of sequential and medium-range 1H-1H contacts that are obscured by chemical shift degeneracy in homonuclear experiments. The 15N chemical shifts are analyzed with respect to correlations with protein secondary structure. In addition, the changes in 15N chemical shift found for the apo----(Cd2+)1----(Ca2+)2 binding sequence confirm that the effects on the protein are mainly associated with chelation of the first ion.     

Structural characterisation of the human alpha-lactalbumin molten globule at high temperature

Molten globules are partially folded forms of proteins thought to be general intermediates in protein folding. The 15N-1H HSQC NMR spectrum of the human alpha-lactalbumin (alpha-LA) molten globule at pH 2 and 20 degrees C is characterised by broad lines which make direct study by NMR methods difficult; this broadening arises from conformational fluctuations throughout the protein on a millisecond to microsecond timescale. Here, we find that an increase in temperature to 50 degrees C leads to a dramatic sharpening of peaks in the 15N-1H HSQC spectrum of human alpha-LA at pH 2. Far-UV CD and ANS fluorescence experiments demonstrate that under these conditions human alpha-LA maintains a high degree of helical secondary structure and the exposed hydrophobic surfaces that are characteristic of a molten globule. Analysis of the H(alpha), H(N) and 15N chemical shifts of the human alpha-LA molten globule at 50 degrees C leads to the identification of regions of native-like helix in the alpha-domain and of non-native helical propensity in the beta-domain. The latter may be responsible for the observed overshoot in ellipticity at 222 nm in kinetic refolding experiments.     

Influence of the completeness of chemical shift assignments on NMR structures obtained with automated NOE assignment

Reliable automated NOE assignment and structure calculation on the basis of a largely complete, assigned input chemical shift list and a list of unassigned NOESY cross peaks has recently become feasible for routine NMR protein structure calculation and has been shown to yield results that are equivalent to those of the conventional, manual approach. However, these algorithms rely on the availability of a virtually complete list of the chemical shifts. This paper investigates the influence of incomplete chemical shift assignments on the reliability of NMR structures obtained with automated NOESY cross peak assignment. The program CYANA was used for combined automated NOESY assignment with the CANDID algorithm and structure calculations with torsion angle dynamics at various degrees of completeness of the chemical shift assignment which was simulated by random omission of entries in the experimental ¹H chemical shift lists that had been used for the earlier, conventional structure determinations of two proteins. Sets of structure calculations were performed choosing the omitted chemical shifts randomly among all assigned hydrogen atoms, or among aromatic hydrogen atoms. For comparison, automated NOESY assignment and structure calculations were performed with the complete experimental chemical shift but under random omission of NOESY cross peaks. When heteronuclear-resolved three-dimensional NOESY spectra are available the current CANDID algorithm yields in the absence of up to about 10% of the experimental ¹H chemical shifts reliable NOE assignments and three-dimensional structures that deviate by less than 2 Å from the reference structure obtained using all experimental chemical shift assignments. In contrast, the algorithm can accommodate the omission of up to 50% of the cross peaks in heteronuclear- resolved NOESY spectra without producing structures with a RMSD of more than 2 Å to the reference structure. When only homonuclear NOESY spectra are available, the algorithm is slightly more susceptible to missing data and can tolerate the absence of up to about 7% of the experimental ¹H chemical shifts or of up to 30% of the NOESY peaks.Abbreviations: BmPBP^A – Bombyx mori pheromone binding protein form A; CYANA – combined assignment and dynamics algorithm for NMR applications; NMR – nuclear magnetic resonance; NOE – nuclear Overhauser effect; NOESY – NOE spectroscopy; RMSD – root-mean-square deviation; WmKT – Williopsis mrakii killer toxin     

Advances towards resonance assignments for uniformly—13C, 15N enriched bacteriorhodopsin at 18.8 T in purple membranes

Solid state NMR spectra from uniformly (13)C, (15)N enriched bacteriorhodospin (bR) purified from H. salinarium were acquired at 18.8 T using magic angle spinning methods. Isolated resonances of 2D (13)C-(13)C spectra exhibited 0.50-0.55 ppm line-widths. Several amino acid types could be assigned, and at least 12 out of 15 Ile peaks could be resolved clearly and identified based on their characteristic chemical shifts and connectivities. This study confirms that high resolution solid state NMR spectra can be obtained for a 248 amino acid uniformly labeled membrane protein in its native membrane environment and indicates that site-specific assignments are likely to be feasible with heteronuclear multidimensional spectra.     

Reaction of horseradish peroxidase with azide and some implications for the heme environmental structure. NMR and kinetic studies.

The azide complex of horseradish peroxidase was studied by high resolution 1H and 15N NMR spectroscopy and by the temperature-jump method. The heme peripheral methyl proton peaks and the ligand 15N resonance were resolved to show that binding of azide by horseradish peroxidase occurs only in acidic solution below pH 6.5. It was also found that the chemical exchange rate of azide with the ferric enzyme was much faster on the 1H and 15N NMR time scale. This was further substantiated by kinetics of azide binding by horseradish peroxidase where the chemical exchange rate was confirmed to be in the microseconds range at pH 5.0 and 23 degrees C. This rate is salient in usual ligand exchange reactions in hemoproteins so far reported. pH dependences of the first order association and dissociation rate constants were also studied by the temperature-jump method to suggest a strong linkage of the azide binding with a proton uptake of an amino acid residue on the enzyme. These results were compared with the case of horse metmyoglobin and were interpreted to indicate that a heme-linked ionizable group on the enzyme facilitates the fast entry of the ligand to the coordination site. A histidyl residue is a possible candidate for the ionizable group of the enzyme.     

The identification of attractive volatiles in aged male mouse urine

In many species, older males are often preferred mates because they carry "good" genes that account for their viability. In some animals, including mice, which rely heavily on chemical communication, there is some indication that an animal's age can be determined by its scent. In order to identify the attractants in aged male mouse urine, chemical and behavioral studies were performed. We herein show that aged mice have higher levels of 3,4-dehydro-exo- brevicomin (DB), 2-sec-butyl-4,5-dihydrothiazole (BT), and 2-isopropyl-4,5-dihydrothiazole (IT) and a lower level of 6-hydroxy-6-methyl-3-heptanone relative to adult male mice. We also demonstrate that the attraction of females to the odor of male mouse urine is greater when the urine is from aged males. However, the attraction of aged urine odor was offset by the ultrafiltration of adult and aged mouse urine. When DB, BT, and IT were added to adult urine, the attraction of the urine was enhanced. Our results suggest that inbred aged male mice develop an aging odor that is attractive to female mice in an experimental setting and that this attraction is due to increased mouse pheromone signaling.     

A novel regulatory function of selenocysteine lyase, a unique catalyst to modulate major urinary protein

Mouse selenocysteine lyase (SCL) catalyzes the decomposition of -selenocysteine into -alanine and selenium with pyridoxal 5′-phosphate as a coenzyme. When using SCL as bait in a yeast two-hybrid screening method, major urinary protein I (MUP-I) was identified as a protein that interacts with SCL. This interaction was confirmed with an in vitro binding assay. MUP-I is known as a pheromone-binding protein that accommodates volatile effectors to affect the physiology and behavior of mice. We found that the binding of 2-naphthol to MUP-I was significantly inhibited by SCL, suggesting that SCL regulates the binding capacity of MUP-I.     

Solution studies of staphylococcal nuclease H124L. 2. 1H, 13C, and 15N chemical shift assignments for the unligated enzyme and analysis of chemical shift changes that accompany formation of the nuclease-thymidine 3',5'-bisphosphate-calcium ternary complex.

Accurate 1H, 15N, and 13C chemical shift assignments were determined for staphylococcal nuclease H124L (in the absence of inhibitor or activator ion). Backbone 1H and 15N assignments, obtained by analysis of three-dimensional 1H-15N HMQC-NOESY data [Wang, J., Mooberry, E.S., Walkenhorst, W.F., & Markley, J. L. (1992) Biochemistry (preceding paper in this issue)], were refined and extended by a combination of homo- and heteronuclear two-dimensional NMR experiments. Staphylococcal nuclease H124L samples used in the homonuclear 1H NMR studies were at natural isotopic abundance or labeled randomly with 2H (to an isotope level of 50%); nuclease H124L samples used for heteronuclear NMR experiments were labeled uniformly with 15N (to an isotope level greater than 95%) or uniformly with 13C (to an isotope level of 26%). Additional nuclease H124L samples were labeled selectively by incorporating single 15N- or 13C-labeled amino acids. The chemical shifts of uncomplexed enzyme were then compared with those determined previously for the nuclease H124L.pdTp.Ca2+ ternary complex [Wang, J., LeMaster, D. M., & Markley, J.L. (1990) Biochemistry 29, 88-101; Wang, J., Hinck, A.P., Loh, S. N., & Markley, J.L. (1990) Biochemistry 29, 102-113; Wang, J., Hinck, A.P., Loh, S.N., & Markley, J.L. (1990) Biochemistry 29, 4242-4253]. The results reveal that the binding of pdTp and Ca2+ induces large shifts in the resonances of several amino acid segments. These chemical shift changes are interpreted in terms of changes in backbone torsion angles that accompany the binding of pdTp and Ca2+; changes at the binding site appear to be transmitted to other regions of the molecule through networks of hydrogen bonds.