Solution structure and backbone dynamics of component IV Glycera dibranchiata monomeric hemoglobin-CO

Conformational Dynamics of Phytoglobin BvPgb1.2 from Beta vulgaris ssp. vulgaris

Plant hemoglobins, often referred to as phytoglobins, play important roles in abiotic stress tolerance. Several essential small physiological metabolites can be bound to these heme proteins. In addition, phytoglobins can catalyze a range of different oxidative reactions in vivo. These proteins are often oligomeric, but the degree and relevance of subunit interactions are largely unknown. In this study, we delineate which residues are involved in dimer formation of a sugar beet phytoglobin type 1.2 (BvPgb1.2) using NMR relaxation experiments. E. coli cells harboring a phytoglobin expression vector were cultivated in isotope-labeled (²H, ¹³C and ¹⁵N) M9 medium. The triple-labeled protein was purified to homogeneity using two chromatographic steps. Two forms of BvPgb1.2 were examined, the oxy-form and the more stable cyanide-form. Using three-dimensional triple-resonance NMR experiments, sequence-specific assignments for CN-bound BvPgb1.2 were achieved for 137 backbone amide cross-peaks in the ¹H-¹⁵N TROSY spectrum, which amounts to 83% of the total number of 165 expected cross-peaks. A large proportion of the non-assigned residues are located in α-helixes G and H, which are proposed to be involved in protein dimerization. Such knowledge around dimer formation will be instrumental for developing a better understanding of phytoglobins' roles in planta.

PDB2CD visualises dynamics within protein structures

Proteins tend to have defined conformations, a key factor in enabling their function. Atomic resolution structures of proteins are predominantly obtained by either solution nuclear magnetic resonance (NMR) or crystal structure methods. However, when considering a protein whose structure has been determined by both these approaches, on many occasions, the resultant conformations are subtly different, as illustrated by the examples in this study. The solution NMR approach invariably results in a cluster of structures whose conformations satisfy the distance boundaries imposed by the data collected; it might be argued that this is evidence of the dynamics of proteins when in solution. In crystal structures, the proteins are often in an energy minimum state which can result in an increase in the extent of regular secondary structure present relative to the solution state depicted by NMR, because the more dynamic ends of alpha helices and beta strands can become ordered at the lower temperatures. This study examines a novel way to display the differences in conformations within an NMR ensemble and between these and a crystal structure of a protein. Circular dichroism (CD) spectroscopy can be used to characterise protein structures in solution. Using the new bioinformatics tool, PDB2CD, which generates CD spectra from atomic resolution protein structures, the differences between, and possible dynamic range of, conformations adopted by a protein can be visualised.

Insight into the allosteric mechanism of Scapharca dimeric hemoglobin

Allosteric regulation is an essential function of many proteins that control a variety of different processes such as catalysis, signal transduction, and gene regulation. Structural rearrangements have historically been considered the main means of communication between different parts of a protein. Recent studies have highlighted the importance, however, of changes in protein flexibility as an effective way to mediate allosteric communication across a protein. Scapharca dimeric hemoglobin (HbI) is the simplest possible allosteric system, with cooperative ligand binding between two identical subunits. Thermodynamic equilibrium studies of the binding of oxygen to HbI have shown that cooperativity is an entropically driven effect. The change in entropy of the system observed upon ligand binding may arise from changes in the protein, the ligand, or the water of the system. The goal of this study is to determine the contribution of the change in entropy of the protein backbone to HbI cooperative binding. Molecular dynamics simulations and nuclear magnetic resonance relaxation techniques have revealed that the fast internal motions of HbI contribute to the cooperative binding to carbon monoxide in two ways: (1) by contributing favorably to the free energy of the system and (2) by participating in the cooperative mechanism at the HbI subunit interface. The internal dynamics of the weakly cooperative HbI mutant, F97Y, were also investigated with the same methods. The changes in backbone NH dynamics observed for F97Y HbI upon ligand binding are not as large as for the wild type, in agreement with the reduced cooperativity observed for this mutant. The results of this study indicate that interface flexibility and backbone conformational entropy of HbI participate in and are important for the cooperative mechanism of carbon monoxide binding.

Influence of heme post-translational modification and distal ligation on the backbone dynamics of a monomeric hemoglobin

The cyanobacterium Synechococcus sp. PCC 7002 uses a hemoglobin of the truncated lineage (GlbN) in the detoxification of reactive species generated in the assimilation of nitrate. In view of a sensing or enzymatic role, several states of GlbN are of interest with respect to its structure-activity relationship. Nuclear magnetic resonance spectroscopy was applied to compare the structure and backbone dynamics of six GlbN forms differing in their oxidation state [Fe(II) or Fe(III)], distal ligand to the iron (histidine, carbon monoxide, or cyanide), or heme post-translational modification (b heme or covalently attached heme). Structural properties were assessed with pseudocontact shift calculations. (15)N relaxation data were analyzed by reduced spectral density mapping (picosecond to nanosecond motions) and by inspection of elevated R(2) values (microsecond to millisecond motions). On the picosecond to nanosecond time scale, GlbN exhibited little flexibility and was unresponsive to the differences among the various forms. Regions of slightly higher mobility were the CE turn, the EF loop, and the H-H' kink. In contrast, fluctuations on the microsecond to millisecond time scale depended on the form. Cyanide binding to the ferric state did not enhance motions, whereas reduction to the ferrous bis-histidine state resulted in elevated R(2) values for several amides. This response was attributed, at least in part, to a weakening of the distal histidine coordination. Carbon monoxide binding quenched some of these fluctuations. The results emphasized the role of the distal ligand in dictating backbone flexibility and illustrated the multiple ways in which motions are controlled by the hemoglobin fold.

Applications and limitations of in silico models in drug discovery

Drug discovery in the late twentieth and early twenty-first century has witnessed a myriad of changes that were adopted to predict whether a compound is likely to be successful, or conversely enable identification of molecules with liabilities as early as possible. These changes include integration of in silico strategies for lead design and optimization that perform complementary roles to that of the traditional in vitro and in vivo approaches. The in silico models are facilitated by the availability of large datasets associated with high-throughput screening, bioinformatics algorithms to mine and annotate the data from a target perspective, and chemoinformatics methods to integrate chemistry methods into lead design process. This chapter highlights the applications of some of these methods and their limitations. We hope this serves as an introduction to in silico drug discovery.

Structure and dynamics of Mycobacterium tuberculosis truncated hemoglobin N: insights from NMR spectroscopy and molecular dynamics simulations

The potent nitric oxide dioxygenase (NOD) activity (trHbN-Fe²⁺-O₂ + (•)NO → trHbN-Fe³⁺-OH₂ + NO₃⁻) of Mycobacterium tuberculosis truncated hemoglobin N (trHbN) protects aerobic respiration from inhibition by (•)NO. The high activity of trHbN has been attributed in part to the presence of numerous short-lived hydrophobic cavities that allow partition and diffusion of the gaseous substrates (•)NO and O₂ to the active site. We investigated the relation between these cavities and the dynamics of the protein using solution NMR spectroscopy and molecular dynamics (MD). Results from both approaches indicate that the protein is mainly rigid with very limited motions of the backbone N-H bond vectors on the picoseconds-nanoseconds time scale, indicating that substrate diffusion and partition within trHbN may be controlled by side-chains movements. Model-free analysis also revealed the presence of slow motions (microseconds-milliseconds), not observed in MD simulations, for many residues located in helices B and G including the distal heme pocket Tyr33(B10). All currently known crystal structures and molecular dynamics data of truncated hemoglobins with the so-called pre-A N-terminal extension suggest a stable α-helical conformation that extends in solution. Moreover, a recent study attributed a crucial role to the pre-A helix for NOD activity. However, solution NMR data clearly show that in near-physiological conditions these residues do not adopt an α-helical conformation and are significantly disordered and that the helical conformation seen in crystal structures is likely induced by crystal contacts. Although this lack of order for the pre-A does not disagree with an important functional role for these residues, our data show that one should not assume an helical conformation for these residues in any functional interpretation. Moreover, future molecular dynamics simulations should not use an initial α-helical conformation for these residues in order to avoid a bias based on an erroneous initial structure for the N-termini residues. This work constitutes the first study of a truncated hemoglobin dynamics performed by solution heteronuclear relaxation NMR spectroscopy.

Conformational dynamics of the tetrameric hemoglobin molecule as revealed by hydrogen exchange: I. Effects of pH, temperature, and ligand binding

Transcriptional regulation depends on the appropriate set of positive and negative signals that provide correct gene expression. Studies in eukaryotic gene expression have shown NC2 to be a general repressor of transcription which blocks the interaction between TFIIB and TBP. However, during the last few years NC2 has been found to bind the transcriptionally active promoters and interact with several positive transcription factors. These data suggested a controversial role of NC2 in transcription. Using in vitro transcription on minimal LTR promoter of HIV-1, we show that removal of NC2 from HeLa nuclear extract increases the yield of transcripts as well as unspecific transcription initiation in a template amount dependent manner. Fractions of HeLa nuclear extract containing NC2 can restore basal transcription repression and precise selection of transcription initiation point. This points to a new role for NC2 as a repressor of inaccurate transcription initiation that allows specific transcription to take place.

PhosphoThr peptide binding globally rigidifies much of the FHA domain from Arabidopsis receptor kinase-associated protein phosphatase

A net increase in the backbone rigidity of the kinase-interacting FHA domain (KI-FHA) from the Arabidopsis receptor kinase-associated protein phosphatase (KAPP) accompanies the binding of a phosphoThr peptide from its CLV1 receptor-like kinase partner, according to (15)N NMR relaxation at 11.7 and 14.1 T. All of the loops of free KI-FHA display evidence of nanosecond-scale motions. Many of these same residues have residual dipolar couplings that deviate from structural predictions. Binding of the CLV1 pT868 peptide seems to reduce nanosecond-scale fluctuations of all loops, including half of the residues of recognition loops. Residues important for affinity are found to be rigid, i.e., conserved residues and residues of the subsite for the key pT+3 peptide position. This behavior parallels SH2 and PTB domain recognition of pTyr peptides. PhosphoThr peptide binding increases KI-FHA backbone rigidity (S(2)) of three recognition loops, a loop nearby, seven strands from the beta-sandwich, and a distal loop. Compensating the trend of increased rigidity, binding enhances fast mobility at a few sites in four loops on the periphery of the recognition surface and in two loops on the far side of the beta-sandwich. Line broadening evidence of microsecond- to millisecond-scale fluctuations occurs across the six-stranded beta-sheet and nearby edges of the beta-sandwich; this forms a network connected by packing of interior side chains and H-bonding. A patch of the slowly fluctuating residues coincides with the site of segment-swapped dimerization in crystals of the FHA domain of human Chfr. Phosphopeptide binding introduces microsecond- to millisecond-scale fluctuations to more residues of the long 8/9 recognition loop of KI-FHA. The rigidity of this FHA domain appears to couple as a whole to pThr peptide binding.

The six zinc fingers of metal-responsive element binding transcription factor-1 form stable and quasi-ordered structures with relatively small differences in zinc affinities

Six Cys(2)His(2) zinc fingers (F1-6) comprise the DNA binding domain of metal-responsive element binding transcription factor-1 (MTF-1). F1-6 is necessary for basal and zinc-induced expression of metallothionein genes. Analysis of NMR structural and dynamic data for an F1-6 protein construct demonstrates that each zinc finger adopts a stable betabetaalpha fold in the presence of stoichiometric Zn(II), provided that all cysteine ligands are in a reduced state. Parallel studies of protein constructs spanning the four N-terminal core DNA binding fingers (F1-4) and two C-terminal low DNA affinity fingers (F5-6) reveal similar stable zinc finger structures. In both the F1-6 and F5-6 proteins, the finger 5 cysteines were found to readily oxidize at neutral pH. Detailed spectral density and hydrodynamic analysis of (15)N relaxation data revealed quasi-ordered anisotropic rotational diffusion properties of the six F1-6 zinc fingers that could influence MTF-1 DNA binding function. A more general effect on the rotational diffusion properties of Cys(2)His(2) zinc fingers was also uncovered that is dependent upon the position of each finger within multifinger domains. Analysis of NMR (1)H-(15)N-heteronuclear single quantum coherence spectral peak intensities measured as a function of added Zn(II) in conjunction with Zn(II) binding modeling studies indicated that the Zn(II) affinities of all MTF-1 zinc fingers are within approximately 10-50-fold. These analyses further suggested that metal sensing by MTF-1 in eukaryotic cells involves multiple zinc fingers and occurs over a 100-fold or less range of accessible Zn(II) concentration.

Structural and dynamic consequences of increasing repeats in a MUC1 peptide tumor antigen

MUC1 mucin is a large transmembrane glycoprotein whose extracelluler domain is composed of repeating units of a 20 amino acid sequence. In the cancer associated state, this protein expression becomes upregulated and underglycosylated. Previous studies, which show an enhanced binding of a 5-repeat over a 1-repeat MUC1 peptide to a panel of anti-MUC1 antibodies, have led us to investigate the structural and dynamic consequences of increasing repeat number. Two MUC1 peptides were studied: a 16mer corresponding to slightly less than one full repeat of the MUC1 tandem repeat sequence (GVTSAPDTRPAPGSTA) and a 40mer corresponding to two full repeats of the MUC1 sequence (VTSAPDTRPAPGSTAPPAHG)2. Isotopically labeled versions of these MUC1 peptides were cloned, expressed, purified, and evaluated structurally and dynamically using 15N- and 13C-edited NMR approaches. The data show that MUC1 structure, dynamics, and antibody binding affinity are invariant with increasing repeat number. In light of these results, we conclude that the enhanced antibody affinity of the 5-repeat over the 1-repeat MUC1 peptide is due to multivalency effects, and not due to the development of higher order structure in the longer length peptides. The implications of these results are discussed within the context of a multiple repeat MUC1 breast cancer vaccine design.

Structure and dynamics of the potato carboxypeptidase inhibitor by 1H and 15N NMR

The solution structure and backbone dynamics of the recombinant potato carboxypeptidase inhibitor (PCI) have been characterized by NMR spectroscopy. The structure, determined on the basis of 497 NOE-derived distance constraints, is much better defined than the one reported in a previous NMR study, with an average pairwise backbone root-mean-square deviation of 0.5 A for the well-defined region of the protein, residues 7-37. Many of the side-chains show now well-defined conformations, both in the hydrophobic core and on the surface of the protein. Overall, the solution structure of free PCI is similar to the one that it shows in the crystal of the complex with carboxypeptidase A. However, some local differences are observed in regions 15-21 and 27-29. In solution, the six N-terminal and the two C-terminal residues are rather flexible, as shown by 15N backbone relaxation measurements. The flexibility of the latter segment may have implications in the binding of the inhibitor by the enzyme. All the remaining residues in the protein are essentially rigid (S2 > 0.8) with the exception of two of them at the end of a short 3/10 helix. Despite the small size of the protein, a number of amide protons are protected from exchange with solvent deuterons. The slowest exchanging protons are those in a small two-strand beta-sheet. The unfolding free energies, as calculated from the exchange rates of these protons, are around 5 kcal/mol. Other protected amide protons are located in the segment 7-12, adjacent to the beta-sheet. Although these residues are not in an extended conformation in PCI, the equivalent residues in structurally homologous proteins form a third strand of the central beta-sheet. The amide protons in the 3/10 helix are only marginally protected, indicating that they exchange by a local unfolding mechanism, which is consistent with the increase in flexibility shown by some of its residues. Backbone alignment-based programs for folding recognition, as opposite to disulfide-bond alignments, reveal new proteins of unrelated sequence and function with a similar structure.

Crystal structures of unligated and CN-ligated Glycera dibranchiata monomer ferric hemoglobin components III and IV

Erythrocytes of the marine annelid, Glycera dibranchiata, contain a mixture of monomeric and polymeric hemoglobins. There are three major monomer hemoglobin components, II, III, IV (also called GMH2, 3, and 4), that have been highly purified and well characterized. We have now crystallized GMH3 and GMH4 and determined their structures to 1.4-1.8 A resolution. The structures were determined for these two monomer hemoglobins in the oxidized (Fe3+, ferric, or met-) forms in both the unligated and cyanide-ligated states. This work differs from two published, refined structures of a Glycera dibranchiata monomer hemoglobin, which has a sequence that is substantially different from any bona fide major monomer hemoglobins (GMH2, 3, or 4). The high-resolution crystal structures (presented here) and the previous NMR structure of CO-ligated GMH4, provide a basis for interpreting structure/function details of the monomer hemoglobins. These details include: (1) the strong correlation between temperature factor and NMR dynamics for respective protein forms; (2) the unique nature of the HisE7Leu primary sequence substitutions in GMH3 and GMH4 and their impact on cyanide ion binding kinetics; (3) the LeuB10Phe difference between GMH3 and GMH4 and its impact on ligand binding; and (4) elucidation of changes in the structural details of the distal and proximal heme pockets upon cyanide binding.

Continuum secondary structure captures protein flexibility

The DSSP program assigns protein secondary structure to one of eight states. This discrete assignment cannot describe the continuum of thermal fluctuations. Hence, a continuous assignment is proposed. Technically, the continuum results from averaging over ten discrete DSSP assignments with different hydrogen bond thresholds. The final continuous assignment for a single NMR model successfully reflected the structural variations observed between all NMR models in the ensemble. The structural variations between NMR models were verified to correlate with thermal motion; these variations were captured by the continuous assignments. Because the continuous assignment reproduces the structural variation between many NMR models from one single model, functionally important variation can be extracted from a single X-ray structure. Thus, continuous assignments of secondary structure may affect future protein structure analysis, comparison, and prediction.

Tissue inhibitor of metalloproteinases-1 undergoes microsecond to millisecond motions at sites of matrix metalloproteinase-induced fit

The N-terminal, matrix metalloproteinase (MMP)-inhibitory fragment of recombinant, human tissue inhibitor of metalloproteinases (TIMP-1) exhibits varied backbone dynamics and rigidity. Most striking is the presence of chemical exchange in the MMP-binding ridge reported to undergo conformational change upon MMP binding. Conformational exchange fluctuations in microseconds to milliseconds map to the sites of MMP-induced fit at residues Val29 through Leu34 of the AB loop and to the Ala65 and Cys70 "hinges" of the CD loop of TIMP-1. Slow chemical exchange is also present at the type I turn of the EF loop at the base of the MMP-binding ridge. These functional slow motions and other fast internal motions are evident from backbone (15)N spin relaxation at 500 and 750 MHz, whether interpreted by the model-free formalism with axial diffusion anisotropy or by the reduced spectral density approach. The conformational exchange is confirmed by its deviation from the trend between R(2) and the cross-correlation rate eta. The magnetic field-dependence indicates that the chemical exchange broadening in the AB and CD loops is fast on the time-scale of chemical shift differences. The conformational exchange rates for most of these exchanging residues, which can closely approach MMP, appear to be a few thousand to several thousand per second. The slow dynamics of the TIMP-1 AB loop contrast the picosecond to nanosecond dynamics reported in the longer TIMP-2 AB loop.

Determinants of backbone dynamics in native BPTI: cooperative influence of the 14-38 disulfide and the Tyr35 side-chain

15Nitrogen relaxation experiments were used to characterize the backbone dynamics of two modified forms of bovine pancreatic trypsin inhibitor (BPTI). In one form, the disulfide between Cys14 and Cys38 in the wild-type protein was selectively reduced and methylated to generate an analog of the final intermediate in the disulfide-coupled folding pathway. The second form was generated by similarly modifying a mutant protein in which Tyr35 was replaced with Gly (Y35G). For both selectively reduced proteins, the overall conformation of native BPTI was retained, and the relaxation data for these proteins were compared with those obtained previously with the native wild-type and Y35G proteins. Removing the disulfide from either protein had only small effects on the observed longitudinal relaxation rates (R1) or heteronuclear cross relaxation rates (nuclear Overhauser effect), suggesting that the 14-38 disulfide has little influence on the fast (ps to ns) backbone dynamics of either protein. In the wild-type protein, the pattern of residues undergoing slower (micros to ms) internal motions, reflected in unusually large transverse relaxation rates (R2), was also largely unaffected by the removal of this disulfide. It thus appears that the large R2 rates previously observed in native wild-type protein are not a direct consequence of isomerization of the 14-38 disulfide. In contrast with the wild-type protein, reducing the disulfide in Y35G BPTI significantly decreased the number of backbone amides displaying large R2 rates. In addition, the frequencies of the backbone motions in the modified protein, estimated from R2 values measured at multiple refocusing delays, appear to span a wider range than those seen in native Y35G BPTI. Together, these observations suggest that the slow internal motions in Y35G BPTI are more independent in the absence of the 14-38 disulfide and that formation of this bond may lead to a substantial loss of conformational entropy. These effects may account for the previous observation that the Y35G substitution greatly destabilizes the disulfide. The results also demonstrate that the disulfide and the buried side-chain influence the dynamics of the folded protein in a highly cooperative fashion, with the effects of removing either being much greater in the absence of the other.