Interconversion of functional motions between mesophilic and thermophilic adenylate kinases

Dynamic properties are functionally important in many proteins, including the enzyme adenylate kinase (AK), for which the open/closed transition limits the rate of catalytic turnover. Here, we compare our previously published coarse-grained (double-well Gō) simulation of mesophilic AK from E. coli (AKmeso) to simulations of thermophilic AK from Aquifex aeolicus (AKthermo). In AKthermo, as with AKmeso, the LID domain prefers to close before the NMP domain in the presence of ligand, but LID rigid-body flexibility in the open (O) ensemble decreases significantly. Backbone foldedness in O and/or transition state (TS) ensembles increases significantly relative to AKmeso in some interdomain backbone hinges and within LID. In contact space, the TS of AKthermo has fewer contacts at the CORE-LID interface but a stronger contact network surrounding the CORE-NMP interface than the TS of AKmeso. A "heated" simulation of AKthermo at 375K slightly increases LID rigid-body flexibility in accordance with the "corresponding states" hypothesis. Furthermore, while computational mutation of 7 prolines in AKthermo to their AKmeso counterparts produces similar small perturbations, mutation of these sites, especially positions 8 and 155, to glycine is required to achieve LID rigid-body flexibility and hinge flexibilities comparable to AKmeso. Mutating the 7 sites to proline in AKmeso reduces some hinges' flexibilities, especially hinge 2, but does not reduce LID rigid-body flexibility, suggesting that these two types of motion are decoupled in AKmeso. In conclusion, our results suggest that hinge flexibility and global functional motions alike are correlated with but not exclusively determined by the hinge residues. This mutational framework can inform the rational design of functionally important flexibility and allostery in other proteins toward engineering novel biochemical pathways.     

Conformational transitions in adenylate kinase. Allosteric communication reduces misligation

Large conformational changes in the LID and NMP domains of adenylate kinase (AKE) are known to be key to ligand binding and catalysis, yet the order of binding events and domain motion is not well understood. Combining the multiple available structures for AKE with the energy landscape theory for protein folding, a theoretical model was developed for allostery, order of binding events, and efficient catalysis. Coarse-grained models and nonlinear normal mode analysis were used to infer that intrinsic structural fluctuations dominate LID motion, whereas ligand-protein interactions and cracking (local unfolding) are more important during NMP motion. In addition, LID-NMP domain interactions are indispensable for efficient catalysis. LID domain motion precedes NMP domain motion, during both opening and closing. These findings provide a mechanistic explanation for the observed 1:1:1 correspondence between LID domain closure, NMP domain closure, and substrate turnover. This catalytic cycle has likely evolved to reduce misligation, and thus inhibition, of AKE. The separation of allosteric motion into intrinsic structural fluctuations and ligand-induced contributions can be generalized to further our understanding of allosteric transitions in other proteins.     

Local motions in a benchmark of allosteric proteins

Allosteric proteins have been studied extensively in the last 40 years, but so far, no systematic analysis of conformational changes between allosteric structures has been carried out. Here, we compile a set of 51 pairs of known inactive and active allosteric protein structures from the Protein Data Bank. We calculate local conformational differences between the two structures of each protein using simple metrics, such as backbone and side-chain Cartesian displacement, and torsion angle change and rearrangement in residue-residue contacts. Thresholds for each metric arise from distributions of motions in two control sets of pairs of protein structures in the same biochemical state. Statistical analysis of motions in allosteric proteins quantifies the magnitude of allosteric effects and reveals simple structural principles about allostery. For example, allosteric proteins exhibit substantial conformational changes comprising about 20% of the residues. In addition, motions in allosteric proteins show strong bias toward weakly constrained regions such as loops and the protein surface. Correlation functions show that motions communicate through protein structures over distances averaging 10-20 residues in sequence space and 10-20 A in Cartesian space. Comparison of motions in the allosteric set and a set of 21 nonallosteric ligand-binding proteins shows that nonallosteric proteins also exhibit bias of motion toward weakly constrained regions and local correlation of motion. However, allosteric proteins exhibit twice as much percent motion on average as nonallosteric proteins with ligand-induced motion. These observations may guide efforts to design flexibility and allostery into proteins.     

Crystal structure of human adenylate kinase 4 (L171P) suggests the role of hinge region in protein domain motion

It is well known that motion of LID and NMP-binding (NMP_bind) domains in adenylate kinase (AK) is important in ligand binding and catalysis. However, the nature of such domain motions is poorly characterized. One of the critical hinge regions is hinge IV, which connects the CORE and LID domains. In addition, the hinge IV contains a strictly conserved residue, L171, in the AK family. To investigate the role of hinge IV, crystal structure of human adenylate kinase 4 (AK4) L171P mutant was determined. This mutation dramatically changes the orientation of the LID domain, which could be described as a novel twisted-and-closed conformation in contrast to the open and closed conformations in other AKs. This mutant provides a new example of domain motions in AK family.     

Global Transitions of Proteins Explored by a Multiscale Hybrid Methodology: Application to Adenylate Kinase

Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.     

Global Transitions of Proteins Explored by a Multiscale Hybrid Methodology: Application to Adenylate Kinase

Efficient and accurate mapping of transition pathways is a challenging problem in allosteric proteins. We propose here a to our knowledge new methodology called collective molecular dynamics (coMD). coMD takes advantage of the collective modes of motions encoded by the fold, simultaneously evaluating the interactions and energetics via a full-atomic MD simulation protocol. The basic approach is to deform the structure collectively along the modes predicted by the anisotropic network model, upon selecting them via a Monte Carlo/Metropolis algorithm from among the complete pool of all accessible modes. Application to adenylate kinase, an allosteric enzyme composed of three domains, CORE, LID, and NMP, shows that both open-to-closed and closed-to-open transitions are readily sampled by coMD, with large-scale motions of the LID dominating. An energy-barrier crossing occurs during the NMP movements. The energy barrier originates from a switch between the salt bridges K136-D118 at the LID-CORE interface and K57-E170 and D33-R156 at the CORE-NMP and LID-NMP interfaces, respectively. Despite its simplicity and computing efficiency, coMD yields ensembles of transition pathways in close accord with detailed full atomic simulations, lending support to its utility as a multiscale hybrid method for efficiently exploring the allosteric transitions of multidomain or multimeric proteins.     

Crystal structure of unligated guanylate kinase from yeast reveals GMP-induced conformational changes

The crystal structure of guanylate kinase (GK) from yeast (Saccharomyces cerevisiae) with a non-acetylated N terminus has been determined in its unligated form (apo-GK) as well as in complex with GMP (GK.GMP). The structure of apo-GK was solved with multiwavelength anomalous diffraction data and refined to an R-factor of 0.164 (R(free)=0.199) at 2.3 A resolution. The structure of GK.GMP was determined using the crystal structure of GK with an acetylated N terminus as the search model and refined to an R-factor of 0.156 (R(free)=0.245) at 1.9 A. GK belongs to the family of nucleoside monophosphate (NMP) kinases and catalyzes the reversible phosphoryl transfer from ATP to GMP. Like other NMP kinases, GK consists of three dynamic domains: the CORE, LID, and NMP-binding domains. Dramatic movements of the GMP-binding domain and smaller but significant movements of the LID domain have been revealed by comparing the structures of apo-GK and GK.GMP. apo-GK has a much more open conformation than the GK.GMP complex. Systematic analysis of the domain movements using the program DynDom shows that the large movements of the GMP-binding domain involve a rotation around an effective hinge axis approximately parallel with helix 3, which connects the GMP-binding and CORE domains. The C-terminal portion of helix 3, which connects to the CORE domain, has strikingly higher temperature factors in GK.GMP than in apo-GK, indicating that these residues become more mobile upon GMP binding. The results suggest that helix 3 plays an important role in domain movement. Unlike the GMP-binding domain, which moves toward the active center of the enzyme upon GMP binding, the LID domain moves away from the active center and makes the presumed ATP-binding site more open. Therefore, the LID domain movement may facilitate the binding of MgATP. The structure of the recombinant GK.GMP complex superimposes very well with that of the native GK.GMP complex, indicating that N-terminal acetylation does not have significant impact on the three-dimensional structure of GK.     

The atomistic mechanism of conformational transition in adenylate kinase: a TEE-REX molecular dynamics study

We report on an atomistic molecular dynamics simulation of the complete conformational transition of Escherichia coli adenylate kinase (ADK) using the recently developed TEE-REX algorithm. Two phases characterize the transition pathway of ADK, which folds into the domains CORE and LID and the AMP binding domain AMPbd. Starting from the closed conformation, half-opening of the AMPbd precedes a partially correlated opening of the LID and AMPbd, defining the second phase. A highly stable salt bridge D118-K136 at the LID-CORE interface, contributing substantially to the total nonbonded LID-CORE interactions, was identified as a major factor that stabilizes the open conformation. Alternative transition pathways, such as AMPbd opening following LID opening, seem unlikely, as full transition events were not observed along this pathway. The simulation data indicate a high enthalpic penalty, possibly obstructing transitions along this route.     

CRANKITE: A fast polypeptide backbone conformation sampler

Background

CRANKITE is a suite of programs for simulating backbone conformations of polypeptides and proteins. The core of the suite is an efficient Metropolis Monte Carlo sampler of backbone conformations in continuous three-dimensional space in atomic details.

Methods

In contrast to other programs relying on local Metropolis moves in the space of dihedral angles, our sampler utilizes local crankshaft rotations of rigid peptide bonds in Cartesian space.

Results

The sampler allows fast simulation and analysis of secondary structure formation and conformational changes for proteins of average length.     

Ligand Docking to Intermediate and Close-To-Bound Conformers Generated by an Elastic Network Model Based Algorithm for Highly Flexible Proteins

Incorporating receptor flexibility in small ligand-protein docking still poses a challenge for proteins undergoing large conformational changes. In the absence of bound structures, sampling conformers that are accessible by apo state may facilitate docking and drug design studies. For this aim, we developed an unbiased conformational search algorithm, by integrating global modes from elastic network model, clustering and energy minimization with implicit solvation. Our dataset consists of five diverse proteins with apo to complex RMSDs 4.7–15 Å. Applying this iterative algorithm on apo structures, conformers close to the bound-state (RMSD 1.4–3.8 Å), as well as the intermediate states were generated. Dockings to a sequence of conformers consisting of a closed structure and its “parents” up to the apo were performed to compare binding poses on different states of the receptor. For two periplasmic binding proteins and biotin carboxylase that exhibit hinge-type closure of two dynamics domains, the best pose was obtained for the conformer closest to the bound structure (ligand RMSDs 1.5–2 Å). In contrast, the best pose for adenylate kinase corresponded to an intermediate state with partially closed LID domain and open NMP domain, in line with recent studies (ligand RMSD 2.9 Å). The docking of a helical peptide to calmodulin was the most challenging case due to the complexity of its 15 Å transition, for which a two-stage procedure was necessary. The technique was first applied on the extended calmodulin to generate intermediate conformers; then peptide docking and a second generation stage on the complex were performed, which in turn yielded a final peptide RMSD of 2.9 Å. Our algorithm is effective in producing conformational states based on the apo state. This study underlines the importance of such intermediate states for ligand docking to proteins undergoing large transitions.     

Structural and dynamic studies on ligand-free adenylate kinase from Mycobacterium tuberculosis revealed a closed conformation that can be related to the reduced catalytic activity

Tuberculosis is the leading cause of death worldwide from a single infectious disease. Search of new therapeutic tools requires the discovery and biochemical characterization of new potential targets among the bacterial proteins essential for the survival and virulence. Among them are the nucleoside monophosphate kinases, involved in the nucleotide biosynthesis. In this work, we determined the solution structure of adenylate kinase (AK) from Mycobacterium tuberculosis (AKmt), a protein of 181 residues that was found to be essential for bacterial survival. The structure was calculated by a simulated annealing protocol and energy minimization using experimental restraints, collected by nuclear magnetic resonance spectroscopy. The final, well-defined 20 NMR structures show an average root-mean-square deviation of 0.77 A for the backbone atoms in regular secondary structure segments. The protein has a central CORE domain, composed of a five-stranded parallel beta-sheet surrounded by seven alpha-helices, and two peripheral domains, AMPbd and LID. As compared to other crystallographic structures of free form AKs, AKmt is more compact, with the AMP(bd) domain closer to the CORE of the protein. Analysis of the (15)N relaxation data enabled us to obtain the global rotational correlation time (9.19 ns) and the generalized order parameters (S(2)) of amide vectors along the polypeptide sequence. The protein exhibits restricted movements on a picosecond to nanosecond time scale in the secondary structural regions with amplitudes characterized by an average S(2)() value of 0.87. The loops beta1/alpha1, beta2/alpha2, alpha2/alpha3, alpha3/alpha4, alpha4/beta3, beta3/alpha5, alpha6/alpha7 (LID), alpha7/alpha8, and beta5/alpha9 exhibit rapid fluctuations with enhanced amplitudes. These structural and dynamic features of AKmt may be related to its low catalytic activity that is 10-fold lower than in their eukaryote counterparts.     

Crystal structure of the zinc-, cobalt-, and iron-containing adenylate kinase from <Emphasis Type="Italic">Desulfovibrio gigas</Emphasis>: a novel metal-containing adenylate kinase from Gram-negative bacteria

Adenylate kinases (AK) from Gram-negative bacteria are generally devoid of metal ions in their LID domain. However, three metal ions, zinc, cobalt, and iron, have been found in AK from Gram-negative bacteria. Crystal structures of substrate-free AK from Desulfovibrio gigas with three different metal ions (Zn²⁺, Zn-AK; Co²⁺, Co-AK; and Fe²⁺, Fe-AK) bound in its LID domain have been determined by X-ray crystallography to resolutions 1.8, 2.0, and 3.0 Å, respectively. The zinc and iron forms of the enzyme were crystallized in space group I222, whereas the cobalt-form crystals were C2. The presence of the metals was confirmed by calculation of anomalous difference maps and by X-ray fluorescence scans. The work presented here is the first report of a structure of a metal-containing AK from a Gram-negative bacterium. The native enzyme was crystallized, and only zinc was detected in the LID domain. Co-AK and Fe-AK were obtained by overexpressing the protein in Escherichia coli. Zn-AK and Fe-AK crystallized as monomers in the asymmetric unit, whereas Co-AK crystallized as a dimer. Nevertheless, all three crystal structures are very similar to each other, with the same LID domain topology, the only change being the presence of the different metal atoms. In the absence of any substrate, the LID domain of all holoforms of AK was present in a fully open conformational state. Normal mode analysis was performed to predict fluctuations of the LID domain along the catalytic pathway.     

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.     

Structural characterization of the closed conformation of mouse guanylate kinase

Guanylate kinase (GMPK) is a nucleoside monophosphate kinase that catalyzes the reversible phosphoryl transfer from ATP to GMP to yield ADP and GDP. In addition to phosphorylating GMP, antiviral prodrugs such as acyclovir, ganciclovir, and carbovir and anticancer prodrugs such as the thiopurines are dependent on GMPK for their activation. Hence, structural information on mammalian GMPK could play a role in the design of improved antiviral and antineoplastic agents. Here we present the structure of the mouse enzyme in an abortive complex with the nucleotides ADP and GMP, refined at 2.1 A resolution with a final crystallographic R factor of 0.19 (R(free) = 0.23). Guanylate kinase is a member of the nucleoside monophosphate (NMP) kinase family, a family of enzymes that despite having a low primary structure identity share a similar fold, which consists of three structurally distinct regions termed the CORE, LID, and NMP-binding regions. Previous studies on the yeast enzyme have shown that these parts move as rigid bodies upon substrate binding. It has been proposed that consecutive binding of substrates leads to "closing" of the active site bringing the NMP-binding and LID regions closer to each other and to the CORE region. Our structure, which is the first of any guanylate kinase with both substrates bound, supports this hypothesis. It also reveals the binding site of ATP and implicates arginines 44, 137, and 148 (in addition to the invariant P-loop lysine) as candidates for catalyzing the chemical step of the phosphoryl transfer.     

How different are structurally flexible and rigid binding sites? Sequence and structural features discriminating proteins that do and do not undergo conformational change upon ligand binding

Proteins are dynamic molecules and often undergo conformational change upon ligand binding. It is widely accepted that flexible loop regions have a critical functional role in enzymes. Lack of consideration of binding site flexibility has led to failures in predicting protein functions and in successfully docking ligands with protein receptors. Here we address the question: which sequence and structural features distinguish the structurally flexible and rigid binding sites? We analyze high-resolution crystal structures of ligand bound (holo) and free (apo) forms of 41 proteins where no conformational change takes place upon ligand binding, 35 examples with moderate conformational change, and 22 cases where a large conformational change has been observed. We find that the number of residue-residue contacts observed per-residue (contact density) does not distinguish flexible and rigid binding sites, suggesting a role for specific interactions and amino acids in modulating the conformational changes. Examination of hydrogen bonding and hydrophobic interactions reveals that cases that do not undergo conformational change have high polar interactions constituting the binding pockets. Intriguingly, the large, aromatic amino acid tryptophan has a high propensity to occur at the binding sites of examples where a large conformational change has been noted. Further, in large conformational change examples, hydrophobic-hydrophobic, aromatic-aromatic, and hydrophobic-polar residue pair interactions are dominant. Further analysis of the Ramachandran dihedral angles (phi, psi) reveals that the residues adopting disallowed conformations are found in both rigid and flexible cases. More importantly, the binding site residues adopting disallowed conformations clustered narrowly into two specific regions of the L-Ala Ramachandran map. Examination of the dihedral angles changes upon ligand binding shows that the magnitude of phi, psi changes are in general minimal, although some large changes particularly between right-handed alpha-helical and extended conformations are seen. Our work further provides an account of conformational changes in the dihedral angles space. The findings reported here are expected to assist in providing a framework for predicting protein-ligand complexes and for template-based prediction of protein function.     

Essential dynamics sampling study of adenylate kinase: comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions

Essential dynamics sampling simulations of the domain conformations of unliganded Escherichia coli adenylate kinase have been performed to determine whether the ligand-induced closed-domain conformation is accessible to the open unliganded enzyme. Adenylate kinase is a three- domain protein with a central CORE domain and twoflanking domains, the LID and the NMPbind domains. The sampling simulations were applied to the CORE and NMPbind domain pair and the CORE and LID domain pair separately. One aim is to compare the results to those of a similar study on the enzyme citrate synthase to determine whether a similar domain-locking mechanism operates in adenylate kinase. Although for adenylate kinase the simulations suggest that the closed-domain conformation of the unliganded enzyme is at a slightly higher free energy than the open for both domain pairs, the results are radically different to those found for citrate synthase. In adenylate kinase the targeted domain conformations could always be achieved, whereas this was not the case in citrate synthase due to an apparent free-energy barrier between the open and closed conformations. Adenylate kinase has been classified as a protein that undergoes closure through a hinge mechanism, whereas citrate synthase has been assigned to the shear mechanism. This was quantified here in terms of the change in the number of interdomain contacting atoms upon closure which showed a considerable increase in adenylate kinase. For citrate synthase this number remained largely the same, suggesting that the domain faces slide over each other during closure. This suggests that shear and hinge mechanisms of domain closure may relate to the existence or absence of an appreciable barrier to closure for the unliganded protein, as the latter can hinge comparatively freely, whereas the former must follow a more constrained path. In general though it appears a bias toward keeping the unliganded enzyme in the open-domain conformation may be a common feature of domain enzymes.     

Comparison of normal mode analyses on a small globular protein in dihedral angle space and Cartesian coordinate space

Normal mode analyses on the protein, bovine pancreatic trypsin inhibitor, in dihedral angle space and Cartesian coordinate space are compared. In Cartesian coordinate space it is found that modes of frequencies lower than 30 cm(-1) contribute 80% of the total mean-square fluctuation and are represented almost completely by motions in the dihedral angles. Bond angle and length fluctuations dominate in modes above 200 cm(-1), but contribute less than 2% to the total mean-square fluctuation. In the low-frequency modes a good correspondence between patterns of atomic displacements was found, but on average the root-mean-square fluctuations of the Cartesian coordinate modes are 13% greater than their dihedral angle counterparts. The main effect of fluctuations in the bond angles and lengths, therefore, is to allow the dihedral angles to become more flexible. As the important subspaces determined from the two methods overlap considerably, dihedral angle space analysis can be applied to proteins too large for Cartesian coordinate space analysis.     

Integrating molecular dynamics and co-evolutionary analysis for reliable target prediction and deregulation of the allosteric inhibition of aspartokinase for amino acid production

Deregulation of allosteric inhibition of enzymes is a challenge for strain engineering and has been achieved so far primarily by random mutation and trial-and-error. In this work, we used aspartokinase, an important allosteric enzyme for industrial amino acids production, to demonstrate a predictive approach that combines protein dynamics and evolution for a rational reengineering of enzyme allostery. Molecular dynamic simulation of aspartokinase III (AK3) from Escherichia coli and statistical coupling analysis of protein sequences of the aspartokinase family allowed to identify a cluster of residues which are correlated during protein motion and coupled during the evolution. This cluster of residues forms an interconnected network mediating the allosteric regulation, including most of the previously reported positions mutated in feedback insensitive AK3 mutants. Beyond these mutation positions, we have successfully constructed another twelve targeted mutations of AK3 desensitized toward lysine inhibition. Six threonine-insensitive mutants of aspartokinase I-homoserine dehydrogenase I (AK1-HD1) were also created based on the predictions. The proposed approach can be widely applied for the deregulation of other allosteric enzymes.     

Conformation and dynamics of an Fab'-bound peptide by isotope-edited NMR spectroscopy.

The dynamics and conformation of the peptide antigen MHKDFLEKIGGL bound to the Fab' fragment of the monoclonal antipeptide antibody B13A2, raised against a peptide from myohemerythrin, have been investigated by isotope-edited NMR techniques. The peptides were labeled with 15N (98%) or 13C (99%) at the backbone of individual amino acid residues. Well-resolved amide proton and nitrogen backbone resonances were obtained and assigned for eight of the 12 residues of this bound peptide. Significant resonance line width and chemical shift differences were observed. The 15N and 1H line width variations are attributed to differential backbone mobilities among the bound peptide residues which are consistent with the previously mapped epitope of this peptide antigen. Local structural information was obtained from isotope-directed NOE studies. The approximate distances associated with the experimental NOEs were estimated on the basis of a theoretical NOE analysis involving the relative integrated intensities of the NOE and source peaks. In this way, the sequential NH-NH NOEs obtained for seven of the Fab'-bound peptide residues were shown to correspond to interproton separations of approximately 3 A or less. Such short distances indicate that the backbone dihedral angles of these residues are in the alpha rather than the beta region of phi,psi conformational space; the peptide most likely adopts a helical conformation from F5 to G11 within the antibody combining site. The significance of these results with respect to the type and extent of conformational information obtainable from studies of high molecular weight systems is discussed.     

Binding of retinol induces changes in rat cellular retinol-binding protein II conformation and backbone dynamics

The structure and backbone dynamics of rat holo cellular retinol-binding protein II (holo-CRBP II) in solution has been determined by multidimensional NMR. The final structure ensemble was based on 3980 distance and 30 dihedral angle restraints, and was calculated using metric matrix distance geometry with pairwise Gaussian metrization followed by simulated annealing. The average RMS deviation of the backbone atoms for the final 25 structures relative to their mean coordinates is 0.85(+/-0.09) A. Comparison of the solution structure of holo-CRBP II with apo-CRBP II indicates that the protein undergoes conformational changes not previously observed in crystalline CRBP II, affecting residues 28-35 of the helix-turn-helix, residues 37-38 of the subsequent linker, as well as the beta-hairpin C-D, E-F and G-H loops. The bound retinol is completely buried inside the binding cavity and oriented as in the crystal structure. The order parameters derived from the (15)N T(1), T(2) and steady-state NOE parameters show that the backbone dynamics of holo-CRBP II is restricted throughout the polypeptide. The T(2) derived apparent backbone exchange rate and amide (1)H exchange rate both indicate that the microsecond to second timescale conformational exchange occurring in the portal region of the apo form has been suppressed in the holo form.