Solution structure of a llama single-domain antibody with hydrophobic residues typical of the VH/VL interface

The three-dimensional structure of a llama single-domain antibody BrucD4-4 was established by use of solution NMR spectroscopy. BrucD4-4 has Val, Gly, Leu, and Trp residues at positions 37, 44, 45, and 47, which are considered to be a hallmark to distinguish llama VH from V(H)H fragments at the germline level. In contrast to the murine and human VHs, BrucD4-4 has sufficient solubility, is monomeric in solution, and displays high-quality NMR spectra characteristic of well-structured proteins. Amide proton/deuterium exchange and the (15)N relaxation data showed that BrucD4-4 has a classic protein structure with a well-packed core and comparatively mobile surface loops. The three-dimensional architecture of BrucD4-4 is analogous to that of VHs from murine and human F(v)s and camelid V(H)Hs with two pleated beta-sheets formed by four and five beta-strands. A canonical and undistorted beta-barrel exposes a number of hydrophobic residues into the solvent on the surface of the three-dimensional structure. The eight-residue H3 loop folds over the side chain of Val37 similarly to that in llama V(H)Hs; however, this interaction may be transient due to the H3 conformational flexibility. Overall, the surface characteristics of BrucD4-4 with respect to hydrophobicity appear to lie between the human VH domain from Fv Pot and the llama V(H)H fragment HC-V, which may explain its enhanced solubility allowing NMR structural analysis.     

Influence of Structural Symmetry on Protein Dynamics

Structural symmetry in homooligomeric proteins has intrigued many researchers over the past several decades. However, the implication of protein symmetry is still not well understood. In this study, we performed molecular dynamics (MD) simulations of two forms of trp RNA binding attenuation protein (TRAP), the wild-type 11-mer and an engineered 12-mer, having two different levels of circular symmetry. The results of the simulations showed that the inter-subunit fluctuations in the 11-mer TRAP were significantly smaller than the fluctuations in the 12-mer TRAP while the internal fluctuations were larger in the 11-mer than in the 12-mer. These differences in thermal fluctuations were interpreted by normal mode analysis and group theory. For the 12-mer TRAP, the wave nodes of the normal modes existed at the flexible interface between the subunits, while the 11-mer TRAP had its nodes within the subunits. The principal components derived from the MD simulations showed similar mode structures. These results demonstrated that the structural symmetry was an important determinant of protein dynamics in circularly symmetric homooligomeric proteins.     

Stability and fluctuations of amide hydrogen bonds in a bacterial cytochrome c: a molecular dynamics study

Molecular dynamics (MD) simulations on a bacterial cytochrome c were performed to investigate the lifetime and fluctuations of backbone hydrogen bonds and to correlate these data with protection factors for hydrogen exchange measured by NMR spectroscopy (Bartalesi et al. in Biochemistry, 42:10923–10930, 2003). The MD simulations provide a consistent pattern in that long lifetimes of hydrogen bonds go along with small amplitude fluctuations. In agreement with experiments, differences in stability were found with a rather flexible N-terminal segment as compared with a more rigid C-terminal part. Protection factors of backbone hydrogen exchange correlate strongly with the number of contacts but also with hydrogen-bond occupancy, hydrogen-bond survival times, as well as the inverse of fluctuations of backbone atoms and hydrogen-bond lengths derived from MD simulation data. We observed a conformational transition in the C-terminal loop, and significant motion in the N-terminal loop, which can be interpreted as being the structural units involved in the onset of the protein unfolding process in agreement with experimental evidence on mitochondrial cytochrome c. Electronic Supplementary Material Supplementary material is available for this article at and is accessible for authorized users. Gernot Kieseritzky and Giulia Morra both contributed equally to this work.     

Coarse-grained simulations of conformational dynamics of proteins: Application to apomyoglobin

A coarse-grained dynamic Monte Carlo method is proposed for investigating the conformational dynamics of proteins. Each residue is represented by two interaction sites, one at the α-carbon, and the other on the amino acid sidechain. Geometry and energy parameters extracted from databank structures are used. The method is applied to the crystal structure of apomyoglobin (apo-Mb). Equilibrium and dynamic properties of apo-Mb are characterized within computation times one order of magnitude shorter than conventional molecular dynamics (MD) simulations. The calculated rms fluctuations in α-carbons are in good agreement with crystallographic temperature factors. Regions exhibiting enhanced conformational mobilities are identified. Among the loops connecting the eight helices A to H, the loop CD undergoes the fastest motions, leading to partial unwinding of helix D. Helix G is the most stable helix on the basis of the kinetic stability of dihedral angles, followed by the respective helices A, E, H, and B. These results, in agreement with H/D exchange and two-dimensional NMR experiments, as well as with MD simulations, lend support to the use of the proposed approach as an efficient, yet physically plausible, means of characterizing protein conformational dynamics. Proteins 31:271–281, 1998. © 1998 Wiley-Liss, Inc.     

Molecular dynamics simulations of peptides and proteins with a continuum electrostatic model based on screened Coulomb potentials

A continuum electrostatics approach for molecular dynamics (MD) simulations of macromolecules is presented and analyzed for its performance on a peptide and a globular protein. The approach incorporates the screened Coulomb potential (SCP) continuum model of electrostatics, which was reported earlier. The model was validated in a broad set of tests some of which were based on Monte Carlo simulations that included single amino acids, peptides, and proteins. The implementation for large-scale MD simulations presented in this article is based on a pairwise potential that makes the electrostatic model suitable for fast analytical calculation of forces. To assess the suitability of the approach, a preliminary validation is conducted, which consists of (i) a 3-ns MD simulation of the immunoglobulin-binding domain of streptococcal protein G, a 56-residue globular protein and (ii) a 3-ns simulation of Dynorphin, a biological peptide of 17 amino acids. In both cases, the results are compared with those obtained from MD simulations using explicit water (EW) molecules in an all-atom representation. The initial structure of Dynorphin was assumed to be an alpha-helix between residues 1 and 9 as suggested from NMR measurements in micelles. The results obtained in the MD simulations show that the helical structure collapses early in the simulation, a behavior observed in the EW simulation and consistent with spectroscopic data that suggest that the peptide may adopt mainly an extended conformation in water. The dynamics of protein G calculated with the SCP implicit solvent model (SCP-ISM) reveals a stable structure that conserves all the elements of secondary structure throughout the entire simulation time. The average structures calculated from the trajectories with the implicit and explicit solvent models had a cRMSD of 1.1 A, whereas each average structure had a cRMSD of about 0.8A with respect to the X-ray structure. The main conformational differences of the average structures with respect to the crystal structure occur in the loop involving residues 8-14. Despite the overall similarity of the simulated dynamics with EW and SCP models, fluctuations of side-chains are larger when the implicit solvent is used, especially in solvent exposed side-chains. The MD simulation of Dynorphin was extended to 40 ns to study its behavior in an aqueous environment. This long simulation showed that the peptide has a tendency to form an alpha-helical structure in water, but the stabilization free energy is too weak, resulting in frequent interconversions between random and helical conformations during the simulation time. The results reported here suggest that the SCP implicit solvent model is adequate to describe electrostatic effects in MD simulation of both peptides and proteins using the same set of parameters. It is suggested that the present approach could form the basis for the development of a reliable and general continuum approach for use in molecular biology, and directions are outlined for attaining this long-term goal.     

Molecular dynamics study of water penetration in staphylococcal nuclease

The ionization properties of Lys and Glu residues buried in the hydrophobic core of staphylococcal nuclease (SN) suggest that the interior of this protein behaves as a highly polarizable medium with an apparent dielectric constant near 10. This has been rationalized previously in terms of localized conformational relaxation concomitant with the ionization of the internal residue, and with contributions by internal water molecules. Paradoxically, the crystal structure of the SN V66E variant shows internal water molecules and the structure of the V66K variant does not. To assess the structural and dynamical character of interior water molecules in SN, a series of 10-ns-long molecular dynamics (MD) simulations was performed with wild-type SN, and with the V66E and V66K variants with Glu66 and Lys66 in the neutral form. Internal water molecules were identified based on their coordination state and characterized in terms of their residence times, average location, dipole moment fluctuations, hydrogen bonding interactions, and interaction energies. The locations of the water molecules that have residence times of several nanoseconds and display small mean-square displacements agree well with the locations of crystallographically observed water molecules. Additional, relatively disordered water molecules that are not observed crystallographically were found in internal hydrophobic locations. All of the interior water molecules that were analyzed in detail displayed a distribution of interaction energies with higher mean value and narrower width than a bulk water molecule. This underscores the importance of protein dynamics for hydration of the protein interior. Further analysis of the MD trajectories revealed that the fluctuations in the protein structure (especially the loop elements) can strongly influence protein hydration by changing the patterns or strengths of hydrogen bonding interactions between water molecules and the protein. To investigate the dynamical response of the protein to burial of charged groups in the protein interior, MD simulations were performed with Glu66 and Lys66 in the charged state. Overall, the MD simulations suggest that a conformational change rather than internal water molecules is the dominant determinant of the high apparent polarizability of the protein interior.     

Thermal unfolding of a llama antibody fragment: a two-state reversible process

Camelids produce functional "heavy chain" antibodies which are devoid of light chains and CH1 domains [Hamers-Casterman, C., et al. (1993) Nature 363, 446-448]. It has been shown that the variable domains of these heavy chain antibodies (the V(HH) fragments) are functional at or after exposure to high temperatures, in contrast to conventional antibodies [Linden van der, R. H. J., et al. (1999) Biochim. Biophys. Acta 1431, 37-44]. For a detailed understanding of the higher thermostability of these V(HH) fragments, knowledge of their structure and conformational dynamics is required. As a first step toward this goal, we report here the essentially complete (1)H and (15)N NMR backbone resonance assignments of a llama V(HH) antibody fragment, and an extensive analysis of the structure at higher temperatures. The H-D exchange NMR data at 300 K indicate that the framework of the llama V(HH) fragment is highly protected with a DeltaG(ex) of >5.4 kcal/mol, while more flexibility is observed for surface residues, particularly in the loops and the two outer strands (residues 4-7, 10-13, and 58-60) of the beta-sheet. The CD data indicate a reversible, two-state unfolding mechanism with a melting transition at 333 K and a DeltaH(m) of 56 kcal/mol. H-D exchange studies using NMR and ESI-MS show that below 313 K exchange occurs through local unfolding events whereas above 333 K exchange mainly occurs through global unfolding. The lack of a stable core at high temperatures, observed for V(HH) fragments, has also been observed for conventional antibody fragments. The main distinction between the llama V(HH) fragment and conventional antibody fragments is the reversibility of the thermal unfolding process, explaining its retained functionality after exposure to high temperatures.     

Separating instability from aggregation propensity in γS-crystallin variants

Molecular dynamics (MD) simulations, circular dichroism (CD), and dynamic light scattering (DLS) measurements were used to investigate the aggregation propensity of the eye-lens protein γS-crystallin. The wild-type protein was investigated along with the cataract-related G18V variant and the symmetry-related G106V variant. The MD simulations suggest that local sequence differences result in dramatic differences in dynamics and hydration between these two apparently similar point mutations. This finding is supported by the experimental measurements, which show that although both variants appear to be mostly folded at room temperature, both display increased aggregation propensity. Although the disease-related G18V variant is not the most strongly destabilized, it aggregates more readily than either the wild-type or the G106V variant. These results indicate that γS-crystallin provides an excellent model system for investigating the role of dynamics and hydration in aggregation by locally unfolded proteins.     

DROIDS 1.20: A GUI-Based Pipeline for GPU-Accelerated Comparative Protein Dynamics

Traditional informatics in comparative genomics work only with static representations of biomolecules (i.e., sequence and structure), thereby ignoring the molecular dynamics (MD) of proteins that define function in the cell. A comparative approach applied to MD would connect this very short timescale process, defined in femtoseconds, to one of the longest in the universe: molecular evolution measured in millions of years. Here, we leverage advances in graphics-processing-unit-accelerated MD simulation software to develop a comparative method of MD analysis and visualization that can be applied to any two homologous Protein Data Bank structures. Our open-source pipeline, DROIDS (Detecting Relative Outlier Impacts in Dynamic Simulations), works in conjunction with existing molecular modeling software to convert any Linux gaming personal computer into a “comparative computational microscope” for observing the biophysical effects of mutations and other chemical changes in proteins. DROIDS implements structural alignment and Benjamini-Hochberg-corrected Kolmogorov-Smirnov statistics to compare nanosecond-scale atom bond fluctuations on the protein backbone, color mapping the significant differences identified in protein MD with single-amino-acid resolution. DROIDS is simple to use, incorporating graphical user interface control for Amber16 MD simulations, cpptraj analysis, and the final statistical and visual representations in R graphics and UCSF Chimera. We demonstrate that DROIDS can be utilized to visually investigate molecular evolution and disease-related functional changes in MD due to genetic mutation and epigenetic modification. DROIDS can also be used to potentially investigate binding interactions of pharmaceuticals, toxins, or other biomolecules in a functional evolutionary context as well.     

A molecular dynamics study of the 41-56 beta-hairpin from B1 domain of protein G.

The structural and dynamical behavior of the 41-56 beta-hairpin from the protein G B1 domain (GB1) has been studied at different temperatures using molecular dynamics (MD) simulations in an aqueous environment. The purpose of these simulations is to establish the stability of this hairpin in view of its possible role as a nucleation site for protein folding. The conformation of the peptide in the crystallographic structure of the protein GB1 (native conformation) was lost in all simulations. The new equilibrium conformations are stable for several nanoseconds at 300K (>10 ns), 350 K (>6.5 ns), and even at 450 K (up to 2.5 ns). The new structures have very similar hairpin-like conformations with properties in agreement with available experimental nuclear Overhauser effect (NOE) data. The stability of the structure in the hydrophobic core region during the simulations is consistent with the experimental data and provides further evidence for the role played by hydrophobic interactions in hairpin structures. Essential dynamics analysis shows that the dynamics of the peptide at different temperatures spans basically the same essential subspace. The main equilibrium motions in this subspace involve large fluctuations of the residues in the turn and ends regions. Of the six interchain hydrogen bonds, the inner four remain stable during the simulations. The space spanned by the first two eigenvectors, as sampled at 450 K, includes almost all of the 47 different hairpin structures found in the database. Finally, analysis of the hydration of the 300 K average conformations shows that the hydration sites observed in the native conformation are still well hydrated in the equilibrium MD ensemble.     

Conformation and aggregation of M13 coat protein studied by molecular dynamics.

Molecular dynamics (MD) simulations are performed on M13 coat protein, a small membrane protein for which both alpha- and beta-structures have been suggested. The simulations are started from initial conformations that are either monomers or dimers of alpha-helices or U-shaped beta-sheets. The lipid bilayer is represented by a hydrophobic potential. The results are analyzed in terms of stability, energy and secondary structure. The U-shaped beta-structure changes from a planar to a twisted form with larger twist for the monomer than the dimer. The beta-sheet is much more flexible than the alpha-helix as monitored by the root mean square (rms) fluctuations of the C alpha atoms. A comparison of the energies after 100 ps MD simulation shows that of the monomers, the alpha-helix has the lowest energy. The energy difference between alpha- and beta-structures decreases from 266 kJ/mol to 148 kJ/mol, when going from monomers to dimers. It is expected that this difference will decrease with higher aggregation numbers.     

Optimization of designed armadillo repeat proteins by molecular dynamics simulations and NMR spectroscopy

A multidisciplinary approach based on molecular dynamics (MD) simulations using homology models, NMR spectroscopy, and a variety of biophysical techniques was used to efficiently improve the thermodynamic stability of armadillo repeat proteins (ArmRPs). ArmRPs can form the basis of modular peptide recognition and the ArmRP version on which synthetic libraries are based must be as stable as possible. The 42-residue internal Arm repeats had been designed previously using a sequence-consensus method. Heteronuclear NMR revealed unfavorable interactions present at neutral but absent at high pH. Two lysines per repeat were involved in repulsive interactions, and stability was increased by mutating both to glutamine. Five point mutations in the capping repeats were suggested by the analysis of positional fluctuations and configurational entropy along multiple MD simulations. The most stabilizing single C-cap mutation Q240L was inferred from explicit solvent MD simulations, in which water penetrated the ArmRP. All mutants were characterized by temperature- and denaturant-unfolding studies and the improved mutants were established as monomeric species with cooperative folding and increased stability against heat and denaturant. Importantly, the mutations tested resulted in a cumulative decrease of flexibility of the folded state in silico and a cumulative increase of thermodynamic stability in vitro. The final construct has a melting temperature of about 85°C, 14.5° higher than the starting sequence. This work indicates that in silico studies in combination with heteronuclear NMR and other biophysical tools may provide a basis for successfully selecting mutations that rapidly improve biophysical properties of the target proteins.     

Combining NMR ensembles and molecular dynamics simulations provides more realistic models of protein structures in solution and leads to better chemical shift prediction

While chemical shifts are invaluable for obtaining structural information from proteins, they also offer one of the rare ways to obtain information about protein dynamics. A necessary tool in transforming chemical shifts into structural and dynamic information is chemical shift prediction. In our previous work we developed a method for 4D prediction of protein ¹H chemical shifts in which molecular motions, the 4th dimension, were modeled using molecular dynamics (MD) simulations. Although the approach clearly improved the prediction, the X-ray structures and single NMR conformers used in the model cannot be considered fully realistic models of protein in solution. In this work, NMR ensembles (NMRE) were used to expand the conformational space of proteins (e.g. side chains, flexible loops, termini), followed by MD simulations for each conformer to map the local fluctuations. Compared with the non-dynamic model, the NMRE+MD model gave 6–17% lower root-mean-square (RMS) errors for different backbone nuclei. The improved prediction indicates that NMR ensembles with MD simulations can be used to obtain a more realistic picture of protein structures in solutions and moreover underlines the importance of short and long time-scale dynamics for the prediction. The RMS errors of the NMRE+MD model were 0.24, 0.43, 0.98, 1.03, 1.16 and 2.39 ppm for ¹Hα, ¹HN, ¹³Cα, ¹³Cβ, ¹³CO and backbone ¹⁵N chemical shifts, respectively. The model is implemented in the prediction program 4DSPOT, available at .     

Separating Instability from Aggregation Propensity in γS-Crystallin Variants

Molecular dynamics (MD) simulations, circular dichroism (CD), and dynamic light scattering (DLS) measurements were used to investigate the aggregation propensity of the eye-lens protein γS-crystallin. The wild-type protein was investigated along with the cataract-related G18V variant and the symmetry-related G106V variant. The MD simulations suggest that local sequence differences result in dramatic differences in dynamics and hydration between these two apparently similar point mutations. This finding is supported by the experimental measurements, which show that although both variants appear to be mostly folded at room temperature, both display increased aggregation propensity. Although the disease-related G18V variant is not the most strongly destabilized, it aggregates more readily than either the wild-type or the G106V variant. These results indicate that γS-crystallin provides an excellent model system for investigating the role of dynamics and hydration in aggregation by locally unfolded proteins.     

Hybrid Monte Carlo with multidimensional replica exchanges: conformational equilibria of the hypervariable regions of a llama VHH antibody domain

Since the structural repertoire of the hypervariable regions of human antibodies is known to be more restricted than what is implied by sequence variability, a common approach to structural prediction is to use a knowledge-based (KB) method, such as the canonical structure model (C. Chothia and A. M. Lesk, Journal of Molecular Biology, 1987, Vol. 196, pp. 901-917). However, this model is less successful when applied to camelid heavy chain antibodies. In this study, molecular simulations were used to examine the conformational equilibria of the hypervariable regions (H1, H2, and H3) of a llama heavy chain variable domain, for which KB predictions are poor. Simulations were carried out using both conventional molecular dynamics (MD) and hybrid Monte Carlo with multidimensional replica exchanges (HYMREX). The advantage of the latter method is its ability to selectively target parts of the Hamiltonian that can most readily improve sampling. A novel variant of HYMREX was implemented in which, besides the temperature, torsional interactions and the range of nonbonded interactions were varied. To compare the sampling abilities of MD and this HYMREX scheme, simulations were started from a misfolded conformational state. Overall, MD yielded final conformations more similar to the initial state, implying quasi-ergodic sampling. In contrast, HYMREX achieved more ergodic sampling, and the majority of conformations that it sampled agreed well with the known crystal structure. The HYMREX simulation results were used to help identify the chief interactions governing the conformational equilibria and to reexamine the key assumptions underlying the KB predictions. The data show that the H1 region exhibited significant conformational freedom, in support of the hypothesis that main-chain structural variability in this region could play a greater role in antigen binding in camelid antibodies than it does in normal antibodies. Key H1 residues and associated inter-loop interactions are conjectured to account for the poor KB predictions.     

REACH coarse-grained biomolecular simulation: transferability between different protein structural classes

Coarse graining of protein interactions provides a means of simulating large biological systems. The REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-graining method, in which the force constants of a residue-scale elastic network model are calculated from the variance-covariance matrix obtained from atomistic molecular dynamics (MD) simulation, involves direct mapping between scales without the need for iterative optimization. Here, the transferability of the REACH force field is examined between protein molecules of different structural classes. As test cases, myoglobin (all α), plastocyanin (all β), and dihydrofolate reductase (α/β) are taken. The force constants derived are found to be closely similar in all three proteins. An MD version of REACH is presented, and low-temperature coarse-grained (CG) REACH MD simulations of the three proteins are compared with atomistic MD results. The mean-square fluctuations of the atomistic MD are well reproduced by the CGMD. Model functions for the CG interactions, derived by averaging over the three proteins, are also shown to produce fluctuations in good agreement with the atomistic MD. The results indicate that, similarly to the use of atomistic force fields, it is now possible to use a single, generic REACH force field for all protein studies, without having first to derive parameters from atomistic MD simulation for each individual system studied. The REACH method is thus likely to be a reliable way of determining spatiotemporal motion of a variety of proteins without the need for expensive computation of long atomistic MD simulations.     

Computational compensatory mutation discovery approach: Predicting a PARP1 variant rescue mutation

The prediction of protein mutations that affect function may be exploited for multiple uses. In the context of disease variants, the prediction of compensatory mutations that reestablish functional phenotypes could aid in the development of genetic therapies. In this work, we present an integrated approach that combines coevolutionary analysis and molecular dynamics (MD) simulations to discover functional compensatory mutations. This approach is employed to investigate possible rescue mutations of a poly(ADP-ribose) polymerase 1 (PARP1) variant, PARP1 V762A, associated with lung cancer and follicular lymphoma. MD simulations show PARP1 V762A exhibits noticeable changes in structural and dynamical behavior compared with wild-type (WT) PARP1. Our integrated approach predicts A755E as a possible compensatory mutation based on coevolutionary information, and molecular simulations indicate that the PARP1 A755E/V762A double mutant exhibits similar structural and dynamical behavior to WT PARP1. Our methodology can be broadly applied to a large number of systems where single-nucleotide polymorphisms have been identified as connected to disease and can shed light on the biophysical effects of such changes as well as provide a way to discover potential mutants that could restore WT-like functionality. This can, in turn, be further utilized in the design of molecular therapeutics that aim to mimic such compensatory effect.     

Molecular dynamics simulations as a tool for improving protein stability

Haloalkane dehalogenase (DhlA) was used as a model protein to explore the possibility to use molecular dynamics (MD) simulations as a tool to identify flexible regions in proteins that can serve as a target for stability enhancement by introduction of a disulfide bond. DhlA consists of two domains: an alpha/beta-hydrolase fold main domain and a cap domain composed of five alpha-helices. MD simulations of DhlA showed high mobility in a helix-loop-helix region in the cap domain, involving residues 184-211. A disulfide cross-link was engineered between residue 201 of this flexible region and residue 16 of the main domain. The mutant enzyme showed substantial changes in both thermal and urea denaturation. The oxidized form of the mutant enzyme showed an increase of the apparent transition temperature from 47.5 to 52.5 degrees C, whereas the T(m,app) of the reduced mutant decreased by more than 8 degrees C compared to the wild-type enzyme. Urea denaturation results showed a similar trend. Measurement of the kinetic stability showed that the introduction of the disulfide bond caused a decrease in activation free energy of unfolding of 0.43 kcal mol(-1) compared to the wild-type enzyme and also indicated that the helix-loop-helix region was involved early in the unfolding process. The results show that MD simulations are capable of identifying mobile protein domains that can successfully be used as a target for stability enhancement by the introduction of a disulfide cross-link.     

NMR analysis shows that a b-type variant of Hydrogenobacter thermophilus cytochrome c552 retains its native structure

Conversion of Hydrogenobacter thermophilus cytochrome c(552) into a b-type cytochrome by mutagenesis of both heme-binding cysteines to alanines significantly reduces the stability of the protein (Tomlinson, E. J., and Ferguson, S. J. (2000) Proc. Natl. Acad. Sci. U. S. A. 97, 5156-5160). To understand the effects of this change on the structure and dynamics of the protein, hetero-nuclear (15)N-edited NMR techniques have been used to characterize this b-type variant. The backbone (15)N, (1)H(N), and (1)H(alpha), and (1)H(beta) resonances of the protein have been assigned. Analysis of (3)J(HN)alpha coupling constants, nuclear Overhauser enhancement intensities, and chemical shift index data demonstrates that the four alpha-helices present in the wild-type protein are retained in the b-type variant. Comparison of the chemical shifts for the b-type and wild-type proteins indicates that the tertiary structures of the two proteins are closely similar. Some subtle differences are, however, observed for residues in the N-terminal region and in the vicinity of the heme-binding pocket. Hydrogen exchange studies show that there are 25 backbone amide protons that exchange very slowly in the b-type variant and confirm that the fluctuations within the b-type protein are of a similar extent to those in the wild-type protein. These data demonstrate the notable retention of the native secondary structure and tertiary fold despite the absence of covalent linkages between the heme group and the protein.