Molecular, thermodynamic, and biological aspects of recognition and function in neurophysin-hormone systems: a model system for the analysis of protein-peptide interactions

Contributions of the interdomain loop, amino terminus, and subunit interface to the ligand-facilitated dimerization of neurophysin: crystal structures and mutation studies of bovine neurophysin-I

Current evidence indicates that the ligand-facilitated dimerization of neurophysin is mediated in part by dimerization-induced changes at the hormone binding site of the unliganded state that increase ligand affinity. To elucidate other contributory factors, we investigated the potential role of neurophysin's short interdomain loop (residues 55-59), particularly the effects of loop residue mutation and of deleting amino-terminal residues 1-6, which interact with the loop and adjacent residues 53-54. The neurophysin studied was bovine neurophysin-I, necessitating determination of the crystal structures of des 1-6 bovine neurophysin-I in unliganded and liganded dimeric states, as well as the structure of its liganded Q58V mutant, in which peptide was bound with unexpectedly increased affinity. Increases in dimerization constant associated with selected loop residue mutations and with deletion of residues 1-6, together with structural data, provided evidence that dimerization of unliganded neurophysin-I is constrained by hydrogen bonding of the side chains of Gln58, Ser56, and Gln55 and by amino terminus interactions, loss or alteration of these hydrogen bonds, and probable loss of amino terminus interactions, contributing to the increased dimerization of the liganded state. An additional intersubunit hydrogen bond from residue 81, present only in the liganded state, was demonstrated as the largest single effect of ligand binding directly on the subunit interface. Comparison of bovine neurophysins I and II indicates broadly similar mechanisms for both, with the exception in neurophysin II of the absence of Gln55 side chain hydrogen bonds in the unliganded state and a more firmly established loss of amino terminus interactions in the liganded state. Evidence is presented that loop status modulates dimerization via long-range effects on neurophysin conformation involving neighboring Phe22 as a key intermediary.

EGFP-tagged vasopressin precursor protein sorting into large dense core vesicles and secretion from PC12 cells

1. Hypothalamic magnocellular neurons synthesize, store, and secrete large quantities of the neuropeptides, vasopressin (VP) and oxytocin (OT), which are synthesized as protein precursors also containing proteins called neurophysins. These protein precursors are sorted through the regulated secretory pathway (RSP), packaged into large dense core vesicles LDCVs, and their peptide products are secreted from nerve terminals in the posterior pituitary. 2. It has been hypothesized that this efficient packaging is dependent on the interaction of the peptide with neurophysin in a complex that forms the granule core. To test this, PC12 cells were transfected with vasopressin precursor DNA constructs that either contained or deleted the neurophysin moiety and tagged with enhanced green fluorescent protein (EGFP) as reporters. The intracellular routing and secretion of the EGFP-tagged VP precursor proteins were studied by in differentiated PC12 cells by fluorescence microscopy, electron microscopic immunocytochemistry, and fluorescent imaging techniques. 3. The data showed that only when the neurophysin was present in the VP precursor construct did the fluorescent fusion protein become routed to the RSP and get efficiently packaged into LDCVs and secreted. These data are consistent with the view that routing of the precursor to LDCVs requires the amino acids that encode the intravesicular chaperone, neurophysin.

Structures of an unliganded neurophysin and its vasopressin complex: implications for binding and allosteric mechanisms

The structures of des 1-6 bovine neurophysin-II in the unliganded state and as its complex with lysine vasopressin were determined crystallographically at resolutions of 2.4 A and 2.3 A, respectively. The structure of the protein component of the vasopressin complex was, with some local differences, similar to that determined earlier of the full-length protein complexed with oxytocin, but relatively large differences, probably intrinsic to the hormones, were observed between the structures of bound oxytocin and bound vasopressin at Gln 4. The structure of the unliganded protein is the first structure of an unliganded neurophysin. Comparison with the liganded state indicated significant binding-induced conformational changes that were the largest in the loop region comprising residues 50-58 and in the 7-10 region. A subtle binding-induced tightening of the subunit interface of the dimer also was shown, consistent with a role for interface changes in neurophysin allosteric mechanism, but one that is probably not predominant. Interface changes are suggested to be communicated from the binding site through the strands of beta-sheet that connect these two regions, in part with mediation by Gly 23. Comparison of unliganded and liganded states additionally reveals that the binding site for the hormone alpha-amino group is largely preformed and accessible in the unliganded state, suggesting that it represents the initial site of hormone protein recognition. The potential molecular basis for its thermodynamic contribution to binding is discussed.

Structural basis of neurophysin hormone specificity: Geometry, polarity, and polarizability in aromatic ring interactions

The structural origins of the specificity of the neurophysin hormone-binding site for an aromatic residue in peptide position 2 were explored by analyzing the binding of a series of peptides in the context of the crystal structure of liganded neurophysin. A new modeling method for describing the van der Waals surface of binding sites assisted in the analysis. Particular attention was paid to the unusually large (5 kcal/mol) difference in binding free energy between Phe and Leu in position 2, a value representing more than three times the maximum expected based on hydrophobicity alone, and additionally remarkable since modeling indicated that the Leu side chain was readily accommodated by the binding pocket. Although evidence was obtained of a weak thermodynamic linkage between the binding interactions of the residue 2 side chain and of the peptide alpha-amino group, two factors are considered central. (1) The bound Leu side chain can establish only one-third of the van der Waals contacts available to a Phe side chain. (2) The bound Phe side chain appears to be additionally stabilized relative to Leu by more favorable dipole and induced dipole interactions with nonaromatic polar and sulfur ligands in the binding pocket, as evidenced by examination of its interactions in the pocket, analysis of the detailed energetics of transfer of Phe and Leu side chains from water to other phases, and comparison with thermodynamic and structural data for the binding of residue 1 side chains in this system. While such polar interactions of aromatic rings have been previously observed, the present results suggest their potential for significant thermodynamic contributions to protein structure and ligand recognition.

Essential dynamics/factor analysis for the interpretation of molecular dynamics trajectories

Subject of this work is the analysis of molecular dynamics (MD) trajectories of neurophysins I (NPI) and II (NPII) and their complexes with the neurophyseal nonapeptide hormones oxytocin (OT) and vasopresssin (VP), respectively, simulated in water. NPs serve in the neurosecretory granules as carrier proteins for the hormones before their release to the blood. The starting data consisted of two pairs of different trajectories for each of the (NPII/VP)2 and (NPI/OT)2 heterotetramers and two more trajectories for the NPII2 and NPI2 homodimers (six trajectories in total). Using essential dynamics which, to our judgement, is equivalent to factor analysis, we found that only about 10 degrees of freedom per trajectory are necessary and sufficient to describe in full the motions relevant for the function of the protein. This is consistent with these motions to explain about 90% of the total variance of the system. These principal degrees of freedom represent slow anharmonic motional modes, clearly pointing at distinguished mobility of the atoms involved in the protein's functionality.

Structure-function relationships of the vasopressin prohormone domains

1. In this review the structure-function relationships of the different vasopressin prohormone domains are dated and discussed, with special reference to the neurophysin and glycopeptide domains. 2. The primary structures of the currently known neurophysins and glycopeptide sequences are compared and discussed. 3. The hormone-binding and aggregational properties of neurophysin are reviewed and related to a possible function within the regulated secretory pathway. 4. It is proposed, based on the properties reviewed here as well as our own data shown here, that the sorting of the vasopressin prohormone is initiated by hormone binding, which triggers aggregation of the prohormone into the characteristic dense cores of the regulated secretory pathway. 5. This may suggest that prohormone sorting into the regulated secretory pathway is, in general, determined by noncovalent, intramolecular interactions that promote aggregation.

Molecular modeling of the neurophysin I/oxytocin complex

Neurophysins I and II (NPI and NPII) act in the neurosecretory granules as carrier proteins for the neurophyseal hormones oxytocin (OT) and vasopressin (VP), respectively. The NPI/OT functional unit, believed to be an (NPI/OT)2 heterotetramer, was modeled using low-resolution structure information, viz. the C alpha carbon atom coordinates of the homologous NPII/dipeptide complex (file 1BN2 in the Brookhaven Protein Databank) as a template. Its all-atom representation was obtained using standard modeling tools available within the INSIGHT/Biopolymer modules supplied by Biosym Technologies Inc. A conformation of the NPI-bound OT, similar to that recently proposed in a transfer NOE experiment, was docked into the ligand-binding site by a superposition of its Cys1-Tyr2 fragment onto the equivalent portion of the dipeptide in the template. The starting complex for the initial refinements was prepared by two alternative strategies, termed Model I and Model II, each ending with a approximately 100 ps molecular dynamics (MD) simulation in water using the AMBER 4.1 force field. The free homodimer NPI2 was obtained by removal of the two OT subunits from their sites, followed by a similar structure refinement. The use of Model I, consisting of a constrained simulated annealing, resulted in a structure remarkably similar to both the NPII/dipeptide complex and a recently published solid-state structure of the NPII/OT complex. Thus, Model I is recommended as the method of choice for the preparation of the starting all-atom data for MD. The MD simulations indicate that, both in the homodimer and in the heterotetramer, the 3(10)-helices demonstrate an increased mobility relative to the remaining body of the protein. Also, the C-terminal domains in the NPI2 homodimer are more mobile than the N-terminal ones. Finally, a distinct intermonomer interaction is identified, concentrated around its most prominent, although not unique, contribution provided by an H-bond from Ser25 O gamma in one NPI unit to Glu81 O epsilon in the other unit. This interaction is present in the heterotetramer (NPI/OT)2 and absent or weak in the NPI2 homodimer. We speculate that this interaction, along with the increased mobility of the 3(10)-helices and the carboxy domains, may contribute to the allosteric communication between ligand binding and NPI dimerization.

Crystal structure of a bovine neurophysin II dipeptide complex at 2.8 A determined from the single-wavelength anomalous scattering signal of an incorporated iodine atom

The crystal structure of a dipeptide complex of bovine neurophysin II has been solved at 2.8 A resolution solely by using single-wavelength anomalous scattering data from a single iodinated derivative. The asymmetric unit is an elongated tetramer of dimensions 110 x 40 x 30 A, composed of two dimers related by pseudo twofold symmetry. Each monomer consists of two homologous layers, each with four antiparallel beta-strands. The two regions are connected by a helix followed by a long loop. Monomer-monomer contacts involve antiparallel beta-sheet interactions, which form a dimer with two layers of eight beta-strands. One peptide per monomer occupies the principal hormone-binding pocket formed by part of the amino-terminal region and parts of the connecting helix and loop, with binding to protein consistent with conclusions drawn from solution studies. Dimer-dimer contacts involve the Tyr49 region adjacent to this site. A fifth dipeptide, of unknown biological significance, helps to stabilize one of the monomer-monomer interfaces and the tetramer-tetramer network in the crystal.