Orienting domains in proteins using dipolar couplings measured by liquid-state NMR: differences in solution and crystal forms of maltodextrin binding protein loaded with beta-cyclodextrin

Protein function is often regulated by conformational changes that occur in response to ligand binding or covalent modification such as phosphorylation. In many multidomain proteins these conformational changes involve reorientation of domains within the protein. Although X-ray crystallography can be used to determine the relative orientation of domains, the crystal-state conformation can reflect the effect of crystal packing forces and therefore may differ from the physiologically relevant form existing in solution. Here we demonstrate that the solution-state conformation of a multidomain protein can be obtained from its X-ray structure using an extensive set of dipolar couplings measured by triple-resonance multidimensional NMR spectroscopy in weakly aligning solvent. The solution-state conformation of the 370-residue maltodextrin-binding protein (MBP) loaded with beta-cyclodextrin has been determined on the basis of one-bond (15)N-H(N), (15)N-(13)C', (13)C(alpha)-(13)C', two-bond (13)C'-H(N), and three-bond (13)C(alpha)-H(N) dipolar couplings measured for 280, 262, 276, 262, and 276 residues, respectively. This conformation was generated by applying hinge rotations to various X-ray structures of MBP seeking to minimize the difference between the experimentally measured and calculated dipolar couplings. Consistent structures have been derived in this manner starting from four different crystal forms of MBP. The analysis has revealed substantial differences between the resulting solution-state conformation and its crystal-state counterpart (Protein Data Bank accession code 1DMB) with the solution structure characterized by an 11(+/-1) degrees domain closure. We have demonstrated that the precision achieved in these analyses is most likely limited by small uncertainties in the intradomain structure of the protein (ca 5 degrees uncertainty in orientation of internuclear vectors within domains). In addition, potential effects of interdomain motion have been considered using a number of different models and it was found that the structures derived on the basis of dipolar couplings accurately represent the effective average conformation of the protein.     

On the Bioactive Conformation of a Small Peptide and its Set of Thermodynamically Accessible Conformations

Abstract

In order to investigate the relationship between the bioactive conformation of a peptide and its set of thermodynamically accessible structures in solution, the conformational profile of the tetrapeptide Ac-Pro-Ala-Pro-Tyr-OH was characterized by computational methods. Search of the conformational space was performed within the molecular mechanics framework using the AMBER4.0 force field with an effective dielectric constant of 80. Unique structures of the peptide were compared with its bioactive conformation for the protein Streptomyces griseus Protease A, as taken from the crystal structure of the enzyme-peptide complex. The results show that the bound conformation is close to one of the unique conformations characterized in the conformational search of the isolated peptide. Moreover, the lowest energy minimum characterized in the conformational search exhibits large deviations when compared to the bound conformation of the crystal structure.     

The conformational free-energy map for solvated neocarrabiose

A Ramachandran map of the conformational potential of mean force (pmf) for neocarrabiose in water was obtained using molecular dynamics (MD) simulations with umbrella sampling. The potential energy map calculated in a previous study for this molecule in vacuum exhibited a global minimum located at (phi = 81 degrees, psi = -141 degrees). However, the global minimum on the new pmf map in aqueous solution is located in an area centered around (phi = 175 degrees, psi = 180 degrees), indicating a considerable solvent shift. This new global minimum-energy solution conformation was found to correspond to the experimental value obtained from NMR-NOE measurements, and is also consistent with the experimental crystal structure for neocarrabiose and the fiber diffraction conformation for iota-carrageenan. The global minimum of the solution pmf and its local topology were found to be approximately reproduced by quick vacuum conformational energy mapping using several approximations that mimic solvation effects by de-emphasizing intramolecular hydrogen bonding.     

Solution NMR conformation of glycosaminoglycans

Nuclear magnetic resonance (NMR) spectroscopy has been giving a pivotal contribution to the progress of glycomics, mostly by elucidating the structural, dynamical, conformational and intermolecular binding aspects of carbohydrates. Particularly in the field of conformation, NOE resonances, scalar couplings, residual dipolar couplings, and chemical shift anisotropy offsets have been the principal NMR parameters utilized. Molecular dynamics calculations restrained by NMR-data input are usually employed in conjunction to generate glycosidic bond dihedral angles. Glycosaminoglycans (GAGs) are a special class of sulfated polysaccharides extensively studied worldwide. Besides regulating innumerous physiological processes, these glycans are also widely explored in the global market as either clinical or nutraceutical agents. The conformational aspects of GAGs are key regulators to the quality of interactions with the functional proteins involved in biological events. This report discusses the solution conformation of each GAG type analyzed by one or more of the above-mentioned methods.     

Solution X-ray scattering evidence for agonist- and antagonist-induced modulation of cleft closure in a glutamate receptor ligand-binding domain

Agonist-induced conformational changes in the ligand-binding domains (LBD) of glutamate receptor ion channels provide the driving force for molecular rearrangements that mediate channel opening and subsequent desensitization. The resulting regulated transmembrane ion fluxes form the basis for most excitatory neuronal signaling in the brain. Crystallographic analysis of the GluR2 LBD core has revealed a ligand-binding cleft located between two lobes. Channel antagonists stabilize an open cleft, whereas agonists stabilize a closed cleft. The crystal structure of the apo form is similar to the antagonist-bound, open state. To understand the conformational behavior of the LBD in the absence of crystal lattice constraints, and thus better to appreciate the thermodynamic constraints on ligand binding, we have undertaken a solution x-ray scattering study using two different constructs encoding either the core or an extended LBD. In agreement with the GluR2 crystal structures, the LBD is more compact in the presence of agonist than it is in the presence of antagonist. However, the time-averaged conformation of the ligand-free core in solution is intermediate between the open, antagonist-bound state and the closed, agonist-bound state, suggesting a conformational equilibrium. Addition of peptide moieties that connect the core domain to the other functional domains in each channel subunit appears to constrain the conformational equilibrium in favor of the open state.     

Solution structure of diastereomeric thymidine 3'-O(N,N-diisopropyl-O-beta-cyanoethyl-phosphoramidothioate)s by NMR spectroscopy

The conformational properties of a diastereomeric nucleoside-phosphorothioate pair is reported as reflected by different NMR parameters. Configurational assessment is based partly on the different NOE (nuclear overhauser enhancement) effects of the individual isomers and on the trends observed in other NMR parameters. Vicinal carbon-phosphorus couplings reflect the predominance of the epsilon(-) conformation for the Sp isomer and the epsilon(t) conformation for the Rp isomer. The effects of solvent and temperature on these vicinal couplings are described and the results are interpreted in terms of conformational shift towards the preferred epsilon(t) conformation.     

The utility of residual dipolar couplings in detecting motion in carbohydrates: application to sucrose

The solution structure and dynamics of sucrose are examined using a combination of NMR residual dipolar coupling and molecular mechanics force fields. It is found that the alignment tensors of the individual rings are different, and that fitting 35 measured residual dipolar couplings to structures with specific phi, psi values indicates the presence of three major conformations: phi, psi=(120 degrees ,270 degrees), (45 degrees, 300 degrees) and (90 degrees ,180 degrees). Furthermore, fitting two structures simultaneously to the 35 residual dipolar couplings results in a substantial improvement in the fits. The existence of multiple conformations having similar stabilities is a strong indication of motion, due to the interconversion among these states. Results from four molecular mechanics force fields are in general agreement with the experimental results. However, there are major disagreements between force fields. Because fits of residual dipolar couplings to structures are dependent on the force field used to calculate the structures, multiple force fields were used to interpret NMR data. It is demonstrated that the pucker of the fructofuranosyl ring affects the calculated potential energy surface, and the fit to the residual dipolar couplings data. Previously published 13C nuclear relaxation results suggesting that sucrose is rigid are not inconsistent with the present results when motional timescales are considered.     

CDR-H3 loop ensemble in solution – conformational selection upon antibody binding

ABSTRACT

We analyzed pairs of protein-binding, peptide-binding and hapten-binding antibodies crystallized as complex and in the absence of the antigen with and without conformational differences upon binding in the complementarity-determining region (CDR)-H3 loop. Here, we introduce a molecular dynamics-based approach to capture a diverse conformational ensemble of the CDR-H3 loop in solution. The results clearly indicate that the inherently flexible CDR-H3 loop indeed needs to be characterized as a conformational ensemble. The conformational changes of the CDR-H3 loop in all antibodies investigated follow the paradigm of conformation selection, because we observe the experimentally determined binding competent conformation without the presence of the antigen within the ensemble of pre-existing conformational states in solution before binding. We also demonstrate for several examples that the conformation observed in the antibody crystal structure without antigen present is actually selected to bind the carboxyterminal tail region of the antigen-binding fragment (Fab). Thus, special care must be taken when characterizing antibody CDR-H3 loops by Fab X-ray structures, and the possibility that pre-existing conformations are present should always be considered.     

What is the average conformation of bacteriophage T4 lysozyme in solution? A domain orientation study using dipolar couplings measured by solution NMR

Lysozyme from T4 bacteriophage is comprised of two domains that are both involved in binding substrate. Although wild-type lysozyme has been exclusively crystallized in a closed form that is similar to the peptidoglycan-bound conformation, a more open structure is thought to be required for ligand binding. To determine the relative arrangement of domains within T4 lysozyme in the solution state, dipolar couplings were measured in several different dilute liquid crystalline media by solution NMR methods. The dipolar coupling data were analyzed with a domain orientation procedure described previously that utilizes high- resolution X-ray structures. The cleft between the domains is significantly larger in the average solution structure than what is observed in the X-ray structure of the ligand-free form of the protein (approximately 17 degrees closure from solution to X-ray structures). A comparison of the solution domain orientation with X-ray-derived structures in the protein data base shows that the solution structure resembles a crystal structure obtained for the M6I mutant. Dipolar couplings were also measured on the lysozyme mutant T21C/T142C, which was oxidized to form an inter-domain disulfide bond (T4SS). In this case, the inter-domain solution structure was found to be more closed than was observed in the crystal (approximately 11 degrees). Direct refinement of lysozyme crystal structures with the measured dipolar couplings using the program CNS, establishes that this degree of closure can be accommodated whilst maintaining the inter-domain cystine bond. The differences between the average solution conformations obtained using dipolar couplings and the crystal conformations for both forms of lysozyme investigated in this study illustrate the impact that crystal packing interactions can have on the arrangement of domains within proteins and the importance of alternative methods to X-ray crystallography for evaluating inter-domain structure.     

Solid and solution state conformation of 1L-1-O-acetyl-2,3:5,6-di-O-isopropylidene-chiro-inositol

The X-ray crystal structure of 1L-1-O-acetyl-2,3:5,6-di-O-isopropylidene-chiro-inositol is described. The inositol ring deviates considerably from the ideal chair conformation to a flattened chair. A comparison of its conformation in solution with that in solid was made by the use of 1H NMR. This conformational analysis revealed that the title compound adopts similar conformations in solid state and in solution states irrespective of solvent polarity.     

S-adenosylmethionine conformations in solution and in protein complexes: conformational influences of the sulfonium group

S-Adenosylmethionine (AdoMet) and other sulfonium ions play central roles in the metabolism of all organisms. The conformational preferences of AdoMet and two other biologically important sulfonium ions, S-methylmethionine and dimethylsulfonioproprionic acid, have been investigated by NMR and computational studies. Molecular mechanics parameters for the sulfonium center have been developed for the AMBER force field to permit analysis of NMR results and to enable comparison of the relative energies of the different conformations of AdoMet that have been found in crystal structures of complexes with proteins. S-Methylmethionine and S-dimethylsulfonioproprionate adopt a variety of conformations in aqueous solution; a conformation with an electrostatic interaction between the sulfonium sulfur and the carboxylate group is not noticeably favored, in contrast to the preferred conformation found by in vacuo calculations. Nuclear Overhauser effect measurements and computational results for AdoMet indicate a predominantly anti conformation about the glycosidic bond with a variety of conformations about the methionyl C(alpha)-C(beta) and C(beta)-C(gamma) bonds. An AdoMet conformation in which the positively charged sulfonium sulfur is near an electronegative oxygen in the ribose ring is common. Comparisons of NMR results for AdoMet with those for the uncharged S-adenosylhomocysteine and 5'-methylthioadenosine, and the anionic ATP, indicate that the solution conformations are not dictated mainly by molecular charge. In 20 reported structures of AdoMet.protein complexes, both anti and syn glycosidic torsional angles are found. The methionyl group typically adopts an extended conformation in complexes with enzymes that transfer the methyl group from the sulfonium center, but is more folded in complexes with proteins that do not catalyze reactions involving the sulfur and which can use the sulfonium sulfur solely as a binding site. The conformational energies of AdoMet in these crystal structures are comparable to those found for AdoMet in solution. The sulfonium sulfur is in van der Waals contact with a protein heteroatom in the structures of four proteins, which reflects an energetically favorable contact. Interactions of the sulfonium with aromatic rings are rarely observed.     

Dehydration-induced conformational changes of poly-l-lysine as influenced by drying rate and carbohydrates

The conformation of hydrated and air-dried poly-l-lysine in thin films was studied using Fourier transform IR spectroscopy in the amide-I region. Hydrated poly-l-lysine has a random coil conformation. Upon slow drying of small droplets of the polypeptide solution over a period of several hours, an extended β-sheet conformation is adopted. This conformational transition can be prevented by fast air-drying within 2–3 min. Slow air-drying in the presence of sucrose also preserves the aqueous conformation and results in the formation of a glassy state. Comparison of shifts of the OH band with temperature indicates that sucrose/poly-l-lysine mixtures form a molecularly more densely packed glassy matrix, having a higher glass transition temperature (T_g), than sucrose alone. Whether direct interaction of sugar and polypeptide or glass formation is involved in the stabilization during slow air-drying was studied by drying in the presence of glucose or dextran. Compared with dextran (and sucrose to a lesser extent), glucose gives superior protection. Dried glucose has the lowest T_g and the best interacting properties. We conclude that either immobilization by fast air-drying or sufficient interaction with a protectant through hydrogen bonding (slow drying) plays the leading role in the preservation of the aqueous protein structure.     

Accounting for Ligand Conformational Restriction in Calculations of Protein-Ligand Binding Affinities

The conformation adopted by a ligand on binding to a receptor may differ from its lowest-energy conformation in solution. In addition, the bound ligand is more conformationally restricted, which is associated with a configurational entropy loss. The free energy change due to these effects is often neglected or treated crudely in current models for predicting binding affinity. We present a method for estimating this contribution, based on perturbation theory using the quasi-harmonic model of Karplus and Kushick as a reference system. The consistency of the method is checked for small model systems. Subsequently we use the method, along with an estimate for the enthalpic contribution due to ligand-receptor interactions, to calculate relative binding affinities. The AMBER force field and generalized Born implicit solvent model is used. Binding affinities were estimated for a test set of 233 protein-ligand complexes for which crystal structures and measured binding affinities are available. In most cases, the ligand conformation in the bound state was significantly different from the most favorable conformation in solution. In general, the correlation between measured and calculated ligand binding affinities including the free energy change due to ligand conformational change is comparable to or slightly better than that obtained by using an empirically-trained docking score. Both entropic and enthalpic contributions to this free energy change are significant.     

The closed and compact domain organization of the 70-kDa human cytochrome P450 reductase in its oxidized state as revealed by NMR

The NADPH cytochrome P450 reductase (CPR), a diflavin enzyme, catalyzes the electron transfer (ET) from NADPH to the substrate P450. The crystal structures of mammalian and yeast CPRs show a compact organization for the two domains containing FMN (flavin mononucleotide) and FAD (flavin adenine dinucleotide), with a short interflavin distance consistent with fast ET from the NADPH-reduced FAD to the second flavin FMN. This conformation, referred as "closed", contrasts with the alternative opened or extended domain arrangements recently described for partially reduced or mutant CPR. Internal domain flexibility in this enzyme is indeed necessary to account for the apparently conflicting requirements of having FMN flavin accessible to both the FAD and the substrate P450 at the same interface. However, how interdomain dynamics influence internal and external ETs in CPR is still largely unknown. Here, we used NMR techniques to explore the global, domain-specific and residue-specific structural and dynamic properties of the nucleotide-free human CPR in solution in its oxidized state. Based on the backbone resonance assignment of this 70-kDa protein, we collected residue-specific (15)N relaxation and (1)H-(15)N residual dipolar couplings. Surprisingly and in contrast with previous studies, the analysis of these NMR data revealed that the CPR exists in a unique and predominant conformation that highly resembles the closed conformation observed in the crystalline state. Based on our findings and the previous observations of conformational equilibria of the CPR in partially reduced states, we propose that the large-scale conformational transitions of the CPR during the catalytic cycle are tightly controlled to ensure optimal electron delivery.     

Disaccharide conformational flexibility. I. An adiabatic potential energy map for sucrose

Constrained conformational energy minimizations have been used to calculate an adiabatic (phi, psi) potential energy surface for the disaccharide sucrose. The inclusion of molecular flexibility in the conformational energy analysis of this disaccharide was found to have a significant effect upon the allowed conformational space of the molecule. Three low-energy regions were identified on the adiabatic energy surface, and two of these regions were found to contain two related local minimum-energy conformations, with similar energies, differing only in the directionality of the intra-residue hydrogen bonds of the glucose portion of the molecule. The known crystal structures of seven molecules containing the sucrose moiety all fall within the region of the primary allowed minimum and are consistent with the relaxed energy map, while these crystal conformations could not be rationalized using energy maps for rigid residue geometries. The greater flexibility of the furanoid ring relative to that of the pyranoid ring contributed significantly to the enlargement of the low-energy region on the adiabatic map. However, in spite of the importance of limited flexibility in understanding the conformation and fluctuations of sucrose, this molecule was found to be considerably more rigid that some other disaccharides, such as maltose and cellobiose, in accord with experimental studies.     

Multiple conformations of E. coli Hsp90 in solution: insights into the conformational dynamics of Hsp90

Hsp90, an essential eukaryotic chaperone, depends upon its intrinsic ATPase activity for function. Crystal structures of the bacterial Hsp90 homolog, HtpG, and the yeast Hsp90 reveal large domain rearrangements between the nucleotide-free and the nucleotide-bound forms. We used small-angle X-ray scattering and recently developed molecular modeling methods to characterize the solution structure of HtpG and demonstrate how it differs from known Hsp90 conformations. In addition to this HtpG conformation, we demonstrate that under physiologically relevant conditions, multiple conformations coexist in equilibrium. In solution, nucleotide-free HtpG adopts a more extended conformation than observed in the crystal, and upon the addition of AMPPNP, HtpG is in equilibrium between this open state and a closed state that is in good agreement with the yeast AMPPNP crystal structure. These studies provide a unique view of Hsp90 conformational dynamics and provide a model for the role of nucleotide in effecting conformational change.     

The structure of cyclolinopeptide A in chloroform refined by RDC measurements

Three‐dimensional structures of molecules traditionally assigned from nuclear Overhauser effects and vicinal coupling constants are recently complemented by measurements of residual dipolar couplings. Residual dipolar couplings measured in a stretched poly(dimethylsiloxane) gel were used to determine the structure of cyclolinopeptide A in chloroform solution at ?50 °C. After structure refinement, conformational details of main cluster were discussed in relation to crystal and nuclear Overhauser effect derived structures. Copyright © 2014 European Peptide Society and John Wiley & Sons, Ltd.     

CD4 binding partially locks the bridging sheet in gp120 but leaves the beta2/3 strands flexible

The structure of the free form HIV gp120, critical for therapeutic agent development, is unavailable due to its high flexibility. Previous thermodynamic data, structural analysis and simulation results have suggested a large conformational change in the core domain upon CD4 binding. The bridging sheet, which consists of four beta-strands with beta20/21 nestling against the inner/outer domains and beta2/3 facing outward, more exposed to the solvent, was proposed to be unfolded in the native state. In order to test this proposition and to characterize the native conformations, we performed potential mean force (PMF) molecular dynamics (MD) simulations on the CD4-bound crystal structure. We pushed the bridging sheet away from the inner and outer domain to explore the accessible conformational space for the bridging sheet. In addition, we performed conventional MD simulations on structures with the bridging sheet partially unfolded to investigate the stability of the association between the inner and outer domains. Based on the free energy profiles, we find that the whole bridging sheet is unlikely to unfold without other concurrent conformational changes. On the other hand, the partial bridging sheet, beta strands 2/3, can switch its conformation from the folded to the unfolded state. Furthermore, relaxation of conformation with partially unfolded bridging sheet through MD simulations leads to a conformation with beta strands 20/21 quickly re-anchoring against the inner and outer domains. Such a conformation, although lacking some of the hydrophobic interactions present in the CD4-bound structure, displayed high stability as further indicated by other restrained MD simulations. The relevance of this conformation to the free form structure and the pathway for conformational change from the free form to the CD4-bound structure is discussed in detail in light of the available unliganded SIV gp120 crystal structure.     

Solid and solution state conformations of (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-allo-inositol and (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-6-O-methyl-allo-inositol

The synthesis and conformational studies of (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-allo-inositol and (+/-)-3-O-acetyl-1,2:4,5-di-O-isopropylidene-6-O-methyl-allo-inositol are described. Solid state conformations of the title compounds have been studied by solving their X-ray crystal structures. The inositol ring in both the compounds deviate considerably from the ideal chair conformation to flattened chair conformation in the solid state. Their conformations in solution were studied by the use of 1H NMR spectroscopy. These conformational analyses revealed that the title compounds adopt similar conformations in solid and solution states irrespective of the solvent polarity.     

Solution structure of human GAIP (Galpha interacting protein): a regulator of G protein signaling.

The solution structure of the human protein GAIP (Galpha interacting protein), a regulator of G protein signaling, has been determined by NMR techniques. Dipolar couplings of the oriented protein in two different liquid crystal media have been used in the structure calculation. The solution structure of GAIP is compared to the crystal structure of an homologous protein from rat (RGS4) complexed to the alpha-subunit of a G protein. Some of RGS4 residues involved in the Galpha-RGS binding interface have similar orientations in GAIP (free form), indicating that upon binding these residues do not suffer conformational rearrangements, and therefore, their role does not seem to be restricted to Galpha interaction but also to RGS folding and stability. We suggest that other structural differences between the two proteins may be related to the process of binding as well as to a distinct efficiency in their respective GTPase activating function.