The tertiary structure and backbone dynamics of human prolactin

Human prolactin is a 199-residue (23 kDa) protein closely related to growth hormone and placental lactogen with properties and functions resembling both a hormone and a cytokine. As a traditional hormone, prolactin is produced by lactotrophic cells in the pituitary and secreted into the bloodstream where it acts distally to regulate reproduction and promote lactation. Pituitary cells store prolactin in secretory granules organized around large prolactin aggregates, which are produced within the trans layer of the Golgi complex. Extrapituitary prolactin is synthesized by a wide variety of cells but is not stored in secretory granules. Extrapituitary prolactin displays immunomodulatory activities and acts as a growth factor for cancers of the breast, prostate and tissues of the female reproductive system. We have determined the tertiary structure of human prolactin using three-dimensional (3D) and four-dimensional (4D) heteronuclear NMR spectroscopy. As expected, prolactin adopts an "up-up-down-down" four-helical bundle topology and resembles other members of the family of hematopoietic cytokines. Prolactin displays three discrete structural differences from growth hormone: (1) a structured N-terminal loop in contact with the first helix, (2) a missing mini-helix in the loop between the first and second helices, and (3) a shorter loop between the second and third helices lacking the perpendicular mini-helix observed in growth hormone. Residues necessary for functional binding to the prolactin receptor are clustered on the prolactin surface in a position similar to growth hormone. The backbone dynamics of prolactin were investigated by analysis of 15N NMR relaxation phenomena and demonstrated a rigid four-helical bundle with relatively mobile interconnecting loops. Comparison of global macromolecular tumbling at 0.1mM and 1.0mM prolactin revealed reversible oligomerization, which was correlated to dynamic light scattering experiments. The existence of a reversible oligomerization reaction in solution provides insight into previous results describing the in vitro and in vivo aggregation properties of human prolactin.     

A comparative study of the backbone dynamics of two closely related lipid binding proteins: Bovine heart fatty acid binding protein and porcine ileal lipid binding protein

The backbone dynamics of bovine heart fatty acid binding protein (H-FABP) and porcine ileal lipid binding protein (ILBP) were studied by 15N NMR relaxation (T1 and T2) and steady state heteronuclear 15N{1H} NOE measurements. The microdynamic parameters characterizing the backbone mobility were determined using the model-free approach. For H-FABP, the non-terminal backbone amide groups display a rather compact protein structure of low flexibility. In contrast, for ILBP an increased number of backbone amide groups display unusually high internal mobility. Furthermore, the data indicate a higher degree of conformational exchange processes in the sec-msec time range for ILBP compared to H-FABP. These results suggest significant differences in the conformational stability for these two structurally highly homologous members of the fatty acid binding protein family.     

Differences in backbone dynamics of two homologous bacterial albumin-binding modules: implications for binding specificity and bacterial adaptation

Proteins G and PAB are bacterial albumin-binding proteins expressed at the surface of group C and G streptococci and Peptostreptococcus magnus, respectively. Repeated albumin-binding domains, known as GA modules, are found in both proteins. The third GA module of protein G from the group G streptococcal strain G148 (G148-GA3) and the second GA module of protein PAB from P.magnus strain ALB8 (ALB8-GA) exhibit 59% sequence identity and both fold to form three-helix bundle structures that are very stable against thermal denaturation. ALB8-GA binds human serum albumin with higher affinity than G148-GA3, but G148-GA3 shows substantially broader albumin-binding specificity than ALB8-GA. The (15)N nuclear magnetic resonance spin relaxation measurements reported here, show that the two GA modules exhibit mobility on the picosecond-nanosecond time scale in directly corresponding regions (loops and termini). Most residues in G148-GA3 were seen to be involved in conformational exchange processes on the microsecond-millisecond time scale, whereas for ALB8-GA such motions were only identified for the beginning of helix 2 and its preceding loop. Furthermore, and more importantly, hydrogen-deuterium exchange and saturation transfer experiments reveal large differences between the two GA modules with respect to motions on the second-hour time scale. The high degree of similarity between the two GA modules with respect to sequence, structure and stability, and the observed differences in dynamics, binding affinity and binding specificity to different albumins, suggest a distinct correlation between dynamics, binding affinity and binding specificity. Finally, it is noteworthy in this context that the module G148-GA3, which has broad albumin-binding specificity, is expressed by group C and G streptococci known to infect all mammalian species, whereas P.magnus with the ALB8-GA module has been isolated only from humans.     

Solution structure and backbone dynamics of the Trypanosoma cruzi cysteine protease inhibitor chagasin

A Trypanosoma cruzi cysteine protease inhibitor, termed chagasin, is the first characterized member of a new family of tight-binding cysteine protease inhibitors identified in several lower eukaryotes and prokaryotes but not present in mammals. In the protozoan parasite T.cruzi, chagasin plays a role in parasite differentiation and in mammalian host cell invasion, due to its ability to modulate the endogenous activity of cruzipain, a lysosomal-like cysteine protease. In the present work, we determined the solution structure of chagasin and studied its backbone dynamics by NMR techniques. Structured as a single immunoglobulin-like domain in solution, chagasin exerts its inhibitory activity on cruzipain through conserved residues placed in three loops in the same side of the structure. One of these three loops, L4, predicted to be of variable length among chagasin homologues, is flexible in solution as determined by measurements of (15)N relaxation. The biological implications of structural homology between chagasin and other members of the immunoglobulin super-family are discussed.     

Solution studies of the SH2 domain from the fyn tyrosine kinase: secondary structure,backbone dynamics and protein association

The SH2 domain from Fyn tyrosine kinase, corresponding to residues 155–270 of the human enzyme, was expressed as a GST-fusion protein in a pGEX-E. coli system. After thrombin cleavage and removal of GST, the protein was studied by heteronuclear NMR. Two different phosphotyrosyl-peptides were synthesized and added to the SH2 domain. One peptide corresponded to the regulatory C-terminal tail region of Fyn. Sequence-specific assignment of NMR spectra was achieved using a combination of¹H-¹⁵N-correlated 2D HSQC,¹⁵N-edited 3D TOCSY-HMQC, and¹⁵N-edited 3D NOESY-HMQC spectra. By analysis of the -proton chemical shifts and NOE intensities, the positions of secondary structural elements were determined and found to correspond closely to that seen in the crystal structure of the, homologous, Src-SH2 domain.To investigate the internal dynamics of the protein backbone, T₁ and T₂ relaxation parameters were measured on the free protein, as well as on both peptide complexes. Analytical ultracentrifugation and dynamic light scattering were employed to measure the effect of concentration and peptide-binding on self-association. The results suggest that, at NMR-sample concentrations, the free protein is present in at least dimeric form. Phosphopeptide binding and lower concentration significantly, but not completely, shift the equilibrium towards monomers. The possible role of this protein association in the regulation of the Src-family tyrosine kinases is discussed.Abbreviations SH Src homology - GST glutathione-S-transferase - IPTG isopropyl--D-galactopyranoside - DTT dithiothreitol - PMSF phenyl-methyl-sulphonyl-fluoride - TBS 50 mM Tris, 150 mM NaCl, 5 mM DTT, pH 8.0 - MWCO molecular weight cut off - NMR nuclear magnetic resonance - HSQC heteronuclear single-quantum correlation - NOESY nuclear Overhauser effect spectroscopy     

Solution structure and backbone dynamics of human Raf-1 kinase inhibitor protein

Human Raf-1 kinase inhibitor protein (hRKIP) is a small multi-functional protein of 187 residues. It contains a conserved pocket, which binds a wide range of ligands from various small molecules to distinct proteins. To provide a structural basis for the ligand diversity of RKIP, we herein determined the solution structure of hRKIP, and analyzed its structural dynamics. In solution, hRKIP mainly comprises two antiparallel β sheets, two α helices and two 3₁₀ helices. NMR dynamic analysis reveals that the overall structure of hRKIP is rigid, but its C-terminal helix which is close to the ligand-binding site is mobile. In addition, residues around the ligand-binding pocket exhibit significant conformational exchange on the μs–ms timescale. Conformational flexibility may allow the ligand-binding pocket and the C-terminal helix to adopt various conformations to interact with different substrates. This work may shed light on the underlying molecular mechanisms of how hRKIP recognizes and binds diverse substrate ligands.     

Solution structure and backbone dynamics of the N-terminal region of the calcium regulatory domain from soybean calcium-dependent protein kinase alpha

Ca(2+)-dependent protein kinases (CDPKs) are vital Ca(2+)-signaling proteins in plants and protists which have both a kinase domain and a self-contained calcium regulatory calmodulin-like domain (CLD). Despite being very similar to CaM (>40% identity) and sharing the same fold, recent biochemical and structural evidence suggests that the behavior of CLD is distinct from its namesake, calmodulin. In this study, NMR spectroscopy is employed to examine the structure and backbone dynamics of a 168 amino acid Ca(2+)-saturated construct of the CLD (NtH-CLD) in which almost the entire C-terminal domain is exchange broadened and not visible in the NMR spectra. Structural characterization of the N-terminal domain indicates that the first Ca(2+)-binding loop is significantly more open than in a recently reported structure of the CLD complexed with a putative intramolecular binding region (JD) in the CDPK. Backbone dynamics suggest that parts of the third helix exhibit unusually high mobility, and significant exchange, consistent with previous findings that this helix interacts with the C-terminal domain. Dynamics data also show that the "tether" region, consisting of the first 11 amino acids of CLD, is highly mobile and these residues exhibit distinctive beta-type secondary structure, which may help to position the JD and CLD. Finally, the unusual global dynamic behavior of the protein is rationalized on the basis of possible interdomain rearrangements and the highly variable environments of the C- and N-terminal domains.     

Solution structure and backbone dynamics of the defunct domain of calcium vector protein.

CaVP (calcium vector protein) is a Ca(2+) sensor of the EF-hand protein family which is highly abundant in the muscle of Amphioxus. Its three-dimensional structure is not known, but according to the sequence analysis, the protein is composed of two domains, each containing a pair of EF-hand motifs. We determined recently the solution structure of the C-terminal domain (Trp81-Ser161) and characterized the large conformational and dynamic changes induced by Ca(2+) binding. In contrast, the N-terminal domain (Ala1-Asp86) has lost the capacity to bind the metal ion due to critical mutations and insertions in the two calcium loops. In this paper, we report the solution structure of the N-terminal domain and its backbone dynamics based on NMR spectroscopy, nuclear relaxation, and molecular modeling. The well-resolved three-dimensional structure is typical of a pair of EF-hand motifs, joined together by a short antiparallel beta-sheet. The tertiary arrangement of the two EF-hands results in a closed-type conformation, with near-antiparallel alpha-helices, similar to other EF-hand pairs in the absence of calcium ions. To characterize the internal dynamics of the protein, we measured the (15)N nuclear relaxation rates and the heteronuclear NOE effect in (15)N-labeled N-CaVP at a magnetic field of 11.74 T and 298 K. The domain is mainly monomeric in solution and undergoes an isotropic Brownian rotational diffusion with a correlation time of 7.1 ns, in good agreement with the fluorescence anisotropy decay measurements. Data analysis using a model-free procedure showed that the amide backbone groups in the alpha-helices and beta-strands undergo highly restricted movements on a picosecond to nanosecond time scale. The amide groups in Ca(2+) binding loops and in the linker fragment also display rapid fluctuations with slightly increased amplitudes.     

NMR studies of the backbone flexibility and structure of human growth hormone: a comparison of high and low pH conformations

(15)N NMR relaxation parameters and amide (1)H/(2)H-exchange rates have been used to characterize the structural flexibility of human growth hormone (rhGH) at neutral and acidic pH. Our results show that the rigidity of the molecule is strongly affected by the solution conditions. At pH 7.0 the backbone dynamics parameters of rhGH are uniform along the polypeptide chain and their values are similar to those of other folded proteins. In contrast, at pH 2.7 the overall backbone flexibility increases substantially compared to neutral pH and the average order parameter approaches the lower limit expected for a folded protein. However, a significant variation of the backbone dynamics through the molecule indicates that under acidic conditions the mobility of the residues becomes more dependent on their location within the secondary structure units. In particular, the order parameters of certain loop regions decrease dramatically and become comparable to those found in unfolded proteins. Furthermore, the HN-exchange rates at low pH reveal that the residues most protected from exchange are clustered at one end of the helical bundle, forming a stable nucleus. We suggest that this nucleus maintains the overall fold of the protein under destabilizing conditions. We therefore conclude that the acid state of rhGH consists of a structurally conserved, but dynamically more flexible helical core surrounded by an aura of highly mobile, unstructured loops. However, in spite of its prominent flexibility the acid state of rhGH cannot be considered a "molten globule" state because of its high stability. It appears from our work that under certain conditions, a protein can tolerate a considerable increase in flexibility of its backbone, along with an increased penetration of water into its core, while still maintaining a stable folded conformation.     

Solution structure and backbone dynamics of the holo form of the frenolicin acyl carrier protein

During polyketide biosynthesis, acyl carrier proteins (ACPs) perform the central role of transferring polyketide intermediates between active sites of polyketide synthase. The 4'-phosphopantetheine prosthetic group of a holo-ACP is a long and flexible arm that can reach into different active sites and provide a terminal sulfhydryl group for the attachment of acyl groups through a thioester linkage. We have determined the solution structure and characterized backbone dynamics of the holo form of the frenolicin acyl carrier protein (fren holo-ACP) by nuclear magnetic resonance (NMR). Unambiguous assignments were made for 433 hydrogen atoms, 333 carbon atoms, and 84 nitrogen atoms, representing a total of 94.6% of the assignable atoms in this protein. From 879 meaningful NOEs and 45 angle constraints, a family of 24 structures has been calculated. The solution structure is composed of three major alpha-helices packed in a bundle with three additional short helices in intervening loops; one of the short helices slowly exchanges between two conformations. Superposition of the major helical regions on the mean structure yields average atomic rmsd values of 0.49 +/- 0.09 and 0.91 +/- 0.08 A for backbone and non-hydrogen atoms, respectively. Although the three-helix bundle fold is conserved among acyl carrier proteins involved in fatty acid synthases and polyketide synthases, a detailed comparison revealed that ACPs from polyketide biosynthetic pathways are more related to each other in tertiary fold than to their homologues from fatty acid biosynthetic pathways. Comparison of the free form of ACPs (NMR structures of fren ACP and the Bacillus subtilis ACP) with the substrate-bound form of ACP (crystal structure of butyryl-ACP from Escherichia coli) suggests that conformational exchange plays a role in substrate binding.     

The dynamic properties of the Hepatitis C Virus E2 envelope protein unraveled by molecular dynamics

Hepatitis C Virus (HCV) is one of the most persistent human viruses. Although effective therapeutic approaches have been recently discovered, their use is limited by the elevated costs. Therefore, the development of alternative/complementary strategies is an urgent need. The E2 glycoprotein, the most immunogenic HCV protein, and its variants represent natural candidates to achieve this goal. Here we report an extensive molecular dynamics (MD) analysis of the intrinsic properties of E2. Our data provide interesting clues on the global and local intrinsic dynamic features of the protein. Present MD data clearly indicate that E2 combines a flexible structure with a network of covalent bonds. Moreover, the analysis of the two most important antigenic regions of the protein provides some interesting insights into their intrinsic structural and dynamic properties. Our data indicate that a fluctuating β-hairpin represents a populated state by the region E2^412?423. Interestingly, the analysis of the epitope E2^427?446 conformation, that undergoes a remarkable rearrangement in the simulation, has significant similarities with the structure that the E2^430?442 fragment adopts in complex with a neutralizing antibody. Present data also suggest that the strict conservation of Gly436 in E2 protein of different HCV genotypes is likely dictated by structural restraints. Moreover, the analysis of the E2^412?423 flexibility provides insights into the mechanisms that some antibodies adopt to anchor Trp437 that is fully buried in E2. Finally, the present investigation suggests that MD simulations should systematically complement crystallographic studies on flexible proteins that are studied in combination with antibodies.     

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.     

Solution structure and dynamics of melanoma inhibitory activity protein

Melanoma inhibitory activity (MIA) is a small secreted protein that is implicated in cartilage cell maintenance and melanoma metastasis. It is representative of a recently discovered family of proteins that contain a Src Homologous 3 (SH3) subdomain. While SH3 domains are normally found in intracellular proteins and mediate protein-protein interactions via recognition of polyproline helices, MIA is single-domain extracellular protein, and it probably binds to a different class of ligands.Here we report the assignments, solution structure, and dynamics of human MIA determined by heteronuclear NMR methods. The structures were calculated in a semi-automated manner without manual assignment of NOE crosspeaks, and have a backbone rmsd of 0.38 Å over the ordered regions of the protein. The structure consists of an SH3-like subdomain with N- and C-terminal extensions of approximately 20 amino acids each that together form a novel fold. The rmsd between the solution structure and our recently reported crystal structure is 0.86 Å over the ordered regions of the backbone, and the main differences are localized to the most dynamic regions of the protein. The similarity between the NMR and crystal structures supports the use of automated NOE assignments and ambiguous restraints to accelerate the calculation of NMR structures.     

Solution structure and dynamics of LuxU from Vibrio harveyi, a phosphotransferase protein involved in bacterial quorum sensing

The marine bacterium Vibrio harveyi controls its bioluminescence by a process known as quorum sensing. In this process, autoinducer molecules are detected by membrane-bound sensor kinase/response regulator proteins (LuxN and LuxQ) that relay a signal via a series of protein phosphorylation reactions to another response regulator protein, LuxO. Phosphorylated LuxO indirectly represses the expression of the proteins responsible for bioluminescence. Integral to this quorum sensing process is the function of the phosphotransferase protein, LuxU. LuxU acts to shuttle the phosphate from the membrane-bound proteins, LuxN and LuxQ, to LuxO. LuxU is a 114 amino acid residue monomeric protein. Solution NMR was used to determine the three-dimensional structure of LuxU. LuxU contains a four-helix bundle topology with the active-site histidine residue (His58) located on alpha-helix C and exposed to solution. The active site represents a cluster of positively charged residues located on an otherwise hydrophobic protein face. NMR spin-relaxation experiments identify a collection of flexible residues localized on the same region of LuxU as His58. The studies described here represent the first structural characterization of an isolated, monomeric bacterial phosphotransferase protein.     

Solution structure and backbone dynamics of the non-receptor protein-tyrosine kinase-6 Src homology 2 domain

Human protein-tyrosine kinase-6 (PTK6, also known as breast tumor kinase (Brk)) is a member of the non-receptor protein-tyrosine kinase family and is expressed in two-thirds of all breast tumors. To understand the structural basis of PTK6 function, we have determined the solution structure and backbone dynamics of the PTK6-Src homology 2 (SH2) domain using multidimensional NMR spectroscopy. The solution structure clearly indicates that the SH2 domain of human PTK6 contains a consensus alpha/beta-fold and a Tyr(P) peptide binding surface, which are common to other SH2 domains. However, two of the alpha-helices (alphaA and alphaB) are located on opposite faces of the central beta-sheet. In addition, the topological arrangement of a central four-stranded antiparallel beta-sheet (strands betaA, betaB, betaC, and betaD) differs from that of other Src family members. Backbone dynamics and Tyr(P) peptide titration experiments revealed that the putative ligand binding sites of the PTK6-SH2 domain undergo distinctive internal motions when compared with other regions of the protein. Surface plasmon resonance analysis showed that the Tyr(P) peptide had a dissociation constant of about 60 microm, which is substantially weaker binding than previously reported for Src family members. The solution structure together with data from the ligand binding mode of PTK6-SH2 provides insight into the molecular basis of the autoinhibitory role of PTK6.     

Investigations into sequence and conformational dependence of backbone entropy,inter-basin dynamics and the Flory isolated-pair hypothesis for peptides

The populations and transitions between Ramachandran basins are studied for combinations of the standard 20 amino acids in monomers, dimers and trimers using an implicit solvent Langevin dynamics algorithm and employing seven commonly used force-fields. Both the basin populations and inter-conversion rates are influenced by the nearest neighbor's conformation and identity, contrary to the Flory isolated-pair hypothesis. This conclusion is robust to the choice of force-field, even though the use of different force-fields produces large variations in the populations and inter-conversion rates between the dominant helical, extended beta, and polyproline II basins. The computed variation of conformational and dynamical properties with different force-fields exceeds the difference between explicit and implicit solvent calculations using the same force-field. For all force-fields, the inter-basin transitions exhibit a directional dependence, with most transitions going through extended beta conformation, even when it is the least populated basin. The implications of these results are discussed in the context of estimates for the backbone entropy of single residues, and for the ability of all-atom simulations to reproduce experimental protein folding data.