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

Computational modeling multiple conformational states of proteins with residual dipolar coupling data

Solution nuclear magnetic resonance spectroscopy provides unique opportunities to study the structure and dynamics of biomolecules in aqueous environments. While spin relaxation methods are well recognized for their ability to probe timescales of motion, residual dipolar couplings (RDCs) provide access to amplitudes and directions of motion, characteristics that are important to the function of these molecules. Although observed in the 1960s, the acquisition and computational analysis of RDCs has gained significant momentum in recent years, and particularly applications to motion in proteins have become more numerous. This trend may well continue as RDCs can easily leverage structures produced by new computational methods (e.g., AlphaFold) to produce functional descriptions. In this report, we provide examples and a summary of the ways that RDCs have been used to confirm the existence of internal dynamics, characterize the type of dynamics, and recover atomic-scale structural ensembles that define the full range of conformational sampling.

Model free analysis of experimental residual dipolar couplings in small organic compounds

Residual dipolar couplings (RDCs) contain information on the relative arrangement and dynamics of internuclear spin vectors in chemical compounds. Classically, RDC data is analyzed by fitting to structure models, while model-free approaches (MFA) directly relate RDCs to the corresponding internuclear vectors. The recently introduced software TITANIA implements the MFA and extracts structure and dynamics parameters directly from RDCs to facilitate de novo structure refinement for small organic compounds. Encouraged by our previous results on simulated data, we herein focus on the prerequisites and challenges faced when using purely experimental data for this approach. These concern mainly the fact that not all couplings are accessible in all media, leading to voids in the RDC matrix and the concomitant effects on the structure refinement. It is shown that RDC data sets obtained experimentally from currently available alignment media and measurement methods are of sufficient quality to allow relative configuration determination even when the relative configuration of the analyte is completely unknown.

TITANIA: Model-Free Interpretation of Residual Dipolar Couplings in the Context of Organic Compounds

Residual dipolar couplings (RDCs) become increasingly important as additional NMR parameters in the structure elucidation of organic compounds but are usually used in fitting procedures to discriminate between (computed) structures that are in accordance with RDCs and others that can be ruled out. Thus, the determination of configurations requires prior structural information. The direct use of RDCs as restraints to construct structures based on RDCs has only recently begun also in organic compounds. No protocol has been published though that uses the vector and dynamics information available in multialignment data sets directly for the joint determination of conformation and configuration of organic compounds. This is proposed in the current study. We show that by employing these data, even a flat or random start structure converges into the correctly configured structure when employing multiple alignment data sets in our iterative procedure. The requirements in terms of the number of RDCs and alignment media necessary are discussed in detail.

An updated classification and mechanistic insights into ligand binding of the substrate-binding proteins

Substrate-binding proteins (SBPs) mediate ligand translocation and have been classified into seven clusters (A-G). Although the substrate specificities of these clusters are known to some extent, their ligand-binding mechanism(s) remain(s) incompletely understood. In this study, the list of SBPs belonging to different clusters was updated (764 SBPs) compared to the previously reported study (504 SBPs). Furthermore, a new cluster referred to as cluster H was identified. Results reveal that SBPs follow different ligand-binding mechanisms. Intriguingly, the majority of the SBPs follow the "one domain movement" rather than the well-known "Venus Fly-trap" mechanism. Moreover, SBPs of a few clusters display subdomain conformational movement rather than the complete movement of the N- and C-terminal domains.

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand-protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.

Altered conformational sampling along an evolutionary trajectory changes the catalytic activity of an enzyme

Several enzymes are known to have evolved from non-catalytic proteins such as solute-binding proteins (SBPs). Although attention has been focused on how a binding site can evolve to become catalytic, an equally important question is: how do the structural dynamics of a binding protein change as it becomes an efficient enzyme? Here we performed a variety of experiments, including propargyl-DO3A-Gd(III) tagging and double electron-electron resonance (DEER) to study the rigid body protein dynamics of reconstructed evolutionary intermediates to determine how the conformational sampling of a protein changes along an evolutionary trajectory linking an arginine SBP to a cyclohexadienyl dehydratase (CDT). We observed that primitive dehydratases predominantly populate catalytically unproductive conformations that are vestiges of their ancestral SBP function. Non-productive conformational states, including a wide-open state, are frozen out of the conformational landscape via remote mutations, eventually leading to extant CDT that exclusively samples catalytically relevant compact states. These results show that remote mutations can reshape the global conformational landscape of an enzyme as a mechanism for increasing catalytic activity.

Conformational Trapping of a β-Glucosides-Binding Protein Unveils the Selective Two-Step Ligand-Binding Mechanism of ABC Importers

Substrate-binding proteins (SBPs), selectively capture ligand(s) and ensure their translocation via its cognate ATP-binding cassette (ABC) import system. SBPs bind their cognate ligand(s) via an induced-fit mechanism known as the "Venus Fly-trap"; however, this mechanism lacks the atomic details of all conformational landscape as the confirmatory evidence(s) in its support. In this study, we delineate the atomic details of an SBP, β-glucosides-binding protein (βGlyBP) from Thermus thermophilus HB8. The protein βGlyBP is multi-specific and binds to different types of β-glucosides varying in their glycosidic linkages viz. β-1,2; β-1,3; β-1,4 and β-1,6 with a degree of polymerization of 2-5 glucosyl units. Structurally, the protein βGlyBP possesses four subdomains (N1, N2, C1 and C2). The unliganded protein βGlyBP remains in an open state, which closes upon binding to sophorose (SOP2), laminari-oligosaccharides (LAMn), cello-oligosaccharides (CELn), and gentiobiose (GEN2). This study reports, for the first time, four different structural states (open-unliganded, partial-open-unliganded, open-liganded and closed-liganded) of the protein βGlyBP, revealing its conformational changes upon ligand binding and suggesting a two-step induced-fit mechanism. Further, results suggest that the conformational changes of N1 and C1 subdomains drive the ligand binding, unlike that of the whole N- and C-terminal domains (NTD and CTD) as known in the "Venus Fly-trap" mechanism. Additionally, profiling of stereo-selection mechanism for α- and β-glucosides reveals that in the ligand-binding site four secondary structural elements (L1, H1, H2 and H3) drive the ligand selection. In summary, results demonstrate that the details of conformational changes and ligand selection are pre-encoded in the SBPs.

Kinetic Modelling of Transport Inhibition by Substrates in ABC Importers

Prokaryotic ATP-binding cassette (ABC) importers require a substrate-binding protein (SBP) for the capture and delivery of the cognate substrate to the transmembrane domain (TMD) of the transporter. Various biochemical compounds have been identified that bind to the SBP but are not transported. The mechanistic basis for the "non-cognate" substrates not being transported differs. Some non-cognate substrates fail to trigger the appropriate conformational change in the SBP, resulting in loss of affinity for the TMD or the inability to allosterically activate transport. In another mechanism, the SBP cannot release the bound non-cognate substrate. Here, we used rate equations to derive the steady-state transport rate of cognate substrates of an ABC importer and investigated how non-cognate substrates influence this rate. We found that under limiting non-cognate substrate concentrations, the transport rate remains unaltered for each of the mechanisms. In contrast, at saturating substrate and SBP concentrations, the effect of the non-cognate substrate depends heavily on the respective mechanism. For instance, the transport rate becomes zero when the non-cognate substrate cannot be released by the SBP. Yet it remains unaffected when substrate release is possible but the SBP cannot dock onto the TMDs. Our work shows how the different mechanisms of substrate inhibition impact the transport kinetics, which is relevant for understanding and manipulating solute fluxes and hence the propagation of cells in nutritionally complex milieus.

On the complementarity of X-ray and NMR data

X-ray crystallography and NMR contain complementary information for the structural characterization of biological macromolecules. X-ray diffraction is primarily sensitive to the overall shape of the molecule, whereas NMR is mostly sensitive to the atomic detail. Their combination can therefore provide a stronger justification for the resulting structure. For their combination we have recently proposed REFMAC-NMR, which relies on primary data from both techniques for joint refinement. This possibility raises the compelling question of how far the complementarity can be extended. In this paper, we describe an integrative approach to the refinement with NMR data of four X-ray structures of hen-egg-white lysozyme, solved at atomic resolution in four different crystal forms, and we demonstrate that the outcome critically depends on the crystal form itself, reflecting the sensitivity of NMR to fine details.

Conformational and dynamic plasticity in substrate-binding proteins underlies selective transport in ABC importers

Substrate-binding proteins (SBPs) are associated with ATP-binding cassette importers and switch from an open to a closed conformation upon substrate binding, providing specificity for transport. We investigated the effect of substrates on the conformational dynamics of six SBPs and the impact on transport. Using single-molecule FRET, we reveal an unrecognized diversity of plasticity in SBPs. We show that a unique closed SBP conformation does not exist for transported substrates. Instead, SBPs sample a range of conformations that activate transport. Certain non-transported ligands leave the structure largely unaltered or trigger a conformation distinct from that of transported substrates. Intriguingly, in some cases, similar SBP conformations are formed by both transported and non-transported ligands. In this case, the inability for transport arises from slow opening of the SBP or the selectivity provided by the translocator. Our results reveal the complex interplay between ligand-SBP interactions, SBP conformational dynamics and substrate transport.

Recent excitements in protein NMR: Large proteins and biologically relevant dynamics

The advent of Transverse Relaxation Optimized SpectroscopY (TROSY) and perdeuteration allowed biomolecular NMR spectroscopists to overcome the size limitation barrier (approx. 20 kDa) in de novo structure determination of proteins. The utility of these techniques was immediately demonstrated on large proteins and protein complexes (e.g. GroELGroES, ClpP protease, Hsp90-p53, 20S proteasome, etc.). Further, recent methodological developments such as Residual Dipolar Couplings and Paramagnetic Relaxation Enhancement allowed accurate measurement of long-range structural restraints. Additionally, Carr-Purcell-Meiboom-Gill (CPMG), rotating frame relaxation experiments (R1(rho)) and saturation transfer experiments (CEST and DEST) created never-before accessibility to the (mu)s-ms timescale dynamic parameters that led to the deeper understanding of biological processes. Meanwhile, the excitement in the field continued with a series of developments in the fast data acquisition methods allowing rapid structural studies on less stable proteins. This review aims to discuss important developments in the field of biomolecular NMR spectroscopy in the recent past, i.e., in the post TROSY era. These developments not only gave access to the structural studies of large protein assemblies, but also revolutionized tools in the arsenal of today's biomolecular NMR and point to a bright future of biomolecular NMR spectroscopy.

Multiple Conformations of Gal3 Protein Drive the Galactose-Induced Allosteric Activation of the GAL Genetic Switch of Saccharomyces cerevisiae

Gal3p is an allosteric monomeric protein that activates the GAL genetic switch of Saccharomyces cerevisiae in response to galactose. Expression of constitutive mutant of Gal3p or overexpression of wild-type Gal3p activates the GAL switch in the absence of galactose. These data suggest that Gal3p exists as an ensemble of active and inactive conformations. Structural data have indicated that Gal3p exists in open (inactive) and closed (active) conformations. However, a mutant of Gal3p that predominantly exists in inactive conformation and is yet capable of responding to galactose has not been isolated. To understand the mechanism of allosteric transition, we have isolated a triple mutant of Gal3p with V273I, T404A, and N450D substitutions, which, upon overexpression, fails to activate the GAL switch on its own but activates the switch in response to galactose. Overexpression of Gal3p mutants with single or double mutations in any of the three combinations failed to exhibit the behavior of the triple mutant. Molecular dynamics analysis of the wild-type and the triple mutant along with two previously reported constitutive mutants suggests that the wild-type Gal3p may also exist in super-open conformation. Furthermore, our results suggest that the dynamics of residue F237 situated in the hydrophobic pocket located in the hinge region drives the transition between different conformations. Based on this study, we suggest that conformational selection mechanism is the driving force in the allosteric transition of Gal3p, which may have implications in other signaling pathways involving monomeric proteins.

Conformational Dynamics in FKBP Domains: Relevance to Molecular Signaling and Drug Design

Among the 22 FKBP domains in the human genome, FKBP12.6 and the first FKBP domains (FK1) of FKBP51 and FKBP52 are evolutionarily and structurally most similar to the archetypical FKBP12. As such, the development of inhibitors with selectivity among these four FKBP domains poses a significant challenge for structure-based design. The pleiotropic effects of these FKBP domains in a range of signaling processes such as the regulation of ryanodine receptor calcium channels by FKBP12 and FKBP12.6 and steroid receptor regulation by the FK1 domains of FKBP51 and FKBP52 amply justify the efforts to develop selective therapies. In contrast to their close structural similarities, these four FKBP domains exhibit a substantial diversity in their conformational flexibility. A number of distinct conformational transitions have been characterized for FKBP12 spanning timeframes from 20 s to 10 ns and in each case these dynamics have been shown to markedly differ from the conformational behavior for one or more of the other three FKBP domains. Protein flexibilitybased inhibitor design could draw upon the transitions that are significantly populated in only one of the targeted proteins. Both the similarities and differences among these four proteins valuably inform the understanding of how dynamical effects propagate across the FKBP domains as well as potentially how such intramolecular transitions might couple to the larger scale transitions that are central to the signaling complexes in which these FKBP domains function.

Unraveling the Conformational Landscape of Ligand Binding to Glucose/Galactose-Binding Protein by Paramagnetic NMR and MD Simulations

Protein dynamics related to function can nowadays be structurally well characterized (i.e., instances obtained by high resolution structures), but they are still ill-defined energetically, and the energy landscapes are only accessible computationally. This is the case for glucose-galactose binding protein (GGBP), where the crystal structures of the apo and holo states provide structural information for the domain rearrangement upon ligand binding, while the time scale and the energetic determinants for such concerted dynamics have been so far elusive. Here, we use GGBP as a paradigm to define a functional conformational landscape, both structurally and energetically, by using an innovative combination of paramagnetic NMR experiments and MD simulations. Anisotropic NMR parameters induced by self-alignment of paramagnetic metal ions was used to characterize the ensemble of conformations adopted by the protein in solution while the rate of interconversion between conformations was elucidated by long molecular dynamics simulation on two states of GGBP, the closed-liganded (holo_cl) and open-unloaded (apo_op) states. Our results demonstrate that, in its apo state, the protein coexists between open-like (68%) and closed-like (32%) conformations, with an exchange rate around 25 ns. Despite such conformational heterogeneity, the presence of the ligand is the ultimate driving force to unbalance the equilibrium toward the holo_cl form, in a mechanism largely governed by a conformational selection mechanism.

How to tackle protein structural data from solution and solid state: An integrated approach

Long-range NMR restraints, such as diamagnetic residual dipolar couplings and paramagnetic data, can be used to determine 3D structures of macromolecules. They are also used to monitor, and potentially to improve, the accuracy of a macromolecular structure in solution by validating or "correcting" a crystal model. Since crystal structures suffer from crystal packing forces they may not be accurate models for the macromolecular structures in solution. However, the presence of real differences should be tested for by simultaneous refinement of the structure using both crystal and solution NMR data. To achieve this, the program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic and paramagnetic NMR data and/or diamagnetic residual dipolar couplings. Inconsistencies between crystal structures and solution NMR data, if any, may be due either to structural rearrangements occurring on passing from the solution to solid state, or to a greater degree of conformational heterogeneity in solution with respect to the crystal. In the case of multidomain proteins, paramagnetic restraints can provide the correct mutual orientations and positions of domains in solution, as well as information on the conformational variability experienced by the macromolecule.

Improved Accuracy from Joint X-ray and NMR Refinement of a Protein-RNA Complex Structure

Integrated experimental approaches play an increasingly important role in structural biology, taking advantage of the complementary information provided by different techniques. In particular, the combination of NMR data with X-ray diffraction patterns may provide accurate and precise information about local conformations not available from average-resolution X-ray structures alone. Here, we refined the structure of a ternary protein-protein-RNA complex comprising three domains, Sxl and Unr, bound to a single-stranded region derived in the msl2 mRNA. The joint X-ray and NMR refinement reveals that-despite the poor quality of the fit found for the original structural model-the NMR data can be largely accommodated within the uncertainty in the atom positioning (structural noise) from the primary X-ray data and that the overall domain arrangements and binding interfaces are preserved on passing from the crystalline state to the solution. The refinement highlights local conformational differences, which provide additional information on specific features of the structure. For example, conformational dynamics and heterogeneity observed at the interface between the CSD1 and the Sxl protein components in the ternary complex are revealed by the combination of NMR and crystallographic data. The joint refinement protocol offers unique opportunities to detect structural differences arising from various experimental conditions and reveals static or dynamic differences in the conformation of the biomolecule between the solution and the crystals.

Investigating protein conformational energy landscapes and atomic resolution dynamics from NMR dipolar couplings: a review

Nuclear magnetic resonance spectroscopy is exquisitely sensitive to protein dynamics. In particular inter-nuclear dipolar couplings, that become measurable in solution when the protein is dissolved in a dilute liquid crystalline solution, report on all conformations sampled up to millisecond timescales. As such they provide the opportunity to describe the Boltzmann distribution present in solution at atomic resolution, and thereby to map the conformational energy landscape in unprecedented detail. The development of analytical methods and approaches based on numerical simulation and their application to numerous biologically important systems is presented.

Solution NMR studies of periplasmic binding proteins and their interaction partners

Periplasmic binding proteins (PBPs) are a crucial part of ATP-binding cassette import systems in Gram-negative bacteria. Central to their function is the ability to undergo a large-scale conformational rearrangement from open-unliganded to closed-liganded, which signals the presence of substrate and starts its translocation. Over the years, PBPs have been extensively studied not only owing to their essential role in nutrient uptake but also because they serve as excellent models for both practical applications (e.g., biosensor technology) and basic research (e.g., allosteric mechanisms). Although much of our knowledge at atomic level has been inferred from the detailed, static pictures afforded by crystallographic studies, nuclear magnetic resonance (NMR) has been able to fill certain gaps in such body of work, particularly with regard to dynamic processes. Here, we review NMR studies on PBPs, and their unique insights on conformation, dynamics, energetics, substrate binding, and interactions with related transport proteins. Based on the analysis of recent paramagnetic NMR results, as well as crystallographic and functional observations, we propose a mechanism that could explain the ability of certain PBPs to achieve a closed conformation in absence of ligand while others seem to remain open until ligand-mediated closure.

Coupling between overall rotational diffusion and domain motions in proteins and its effect on dielectric spectra

In this work, we formulate a closed-form solution of the model of a semirigid molecule for the case of fluctuating and reorienting molecular electric dipole moment. We illustrate with numeric calculations the impact of protein domain motions on dielectric spectra using the example of the 128 kDa protein dimer of Enzyme I. We demonstrate that the most drastic effect occurs for situations when the characteristic time of protein domain dynamics is comparable to the time of overall molecular rotational diffusion. We suggest that protein domain motions could be a possible explanation for the high-frequency contribution that accompanies the major relaxation dispersion peak in the dielectric spectra of protein aqueous solutions. We propose that the presented computational methodology could be used for the simultaneous analysis of dielectric spectroscopy and nuclear magnetic resonance data. Proteins 2015; 83:1571-1581. © 2015 Wiley Periodicals, Inc.

The solution structure, binding properties, and dynamics of the bacterial siderophore-binding protein FepB

The periplasmic binding protein (PBP) FepB plays a key role in transporting the catecholate siderophore ferric enterobactin from the outer to the inner membrane in Gram-negative bacteria. The solution structures of the 34-kDa apo- and holo-FepB from Escherichia coli, solved by NMR, represent the first solution structures determined for the type III class of PBPs. Unlike type I and II PBPs, which undergo large "Venus flytrap" conformational changes upon ligand binding, both forms of FepB maintain similar overall folds; however, binding of the ligand is accompanied by significant loop movements. Reverse methyl cross-saturation experiments corroborated chemical shift perturbation results and uniquely defined the binding pocket for gallium enterobactin (GaEnt). NMR relaxation experiments indicated that a flexible loop (residues 225-250) adopted a more rigid and extended conformation upon ligand binding, which positioned residues for optimal interactions with the ligand and the cytoplasmic membrane ABC transporter (FepCD), respectively. In conclusion, this work highlights the pivotal role that structural dynamics plays in ligand binding and transporter interactions in type III PBPs.

Simultaneous use of solution NMR and X-ray data in REFMAC5 for joint refinement/detection of structural differences

The program REFMAC5 from CCP4 was modified to allow the simultaneous use of X-ray crystallographic data and paramagnetic NMR data (pseudocontact shifts and self-orientation residual dipolar couplings) and/or diamagnetic residual dipolar couplings. Incorporation of these long-range NMR restraints in REFMAC5 can reveal differences between solid-state and solution conformations of molecules or, in their absence, can be used together with X-ray crystallographic data for structural refinement. Since NMR and X-ray data are complementary, when a single structure is consistent with both sets of data and still maintains reasonably `ideal' geometries, the reliability of the derived atomic model is expected to increase. The program was tested on five different proteins: the catalytic domain of matrix metalloproteinase 1, GB3, ubiquitin, free calmodulin and calmodulin complexed with a peptide. In some cases the joint refinement produced a single model consistent with both sets of observations, while in other cases it indicated, outside the experimental uncertainty, the presence of different protein conformations in solution and in the solid state.

SedNMR: a web tool for optimizing sedimentation of macromolecular solutes for SSNMR

We have proposed solid state NMR (SSNMR) of sedimented solutes as a novel approach to sample preparation for biomolecular SSNMR without crystallization or other sample manipulations. The biomolecules are confined by high gravity--obtained by centrifugal forces either directly in a SSNMR rotor or in a ultracentrifugal device--into a hydrated non-crystalline solid suitable for SSNMR investigations. When gravity is removed, the sample reverts to solution and can be treated as any solution NMR sample. We here describe a simple web tool to calculate the relevant parameters for the success of the experiment.

Equivalence between Euler angle conventions for the description of tensorial interactions in liquid NMR: application to different software programs

Long-range orientational restraints derived from alignment or rotational diffusion tensors have greatly contributed to the expansion of applications in biomolecular NMR. The orientation of the principal axis system of these tensors is usually described by the so-called Euler angles. However, no clear consensus has emerged concerning the convention of the associated orthogonal rotations. As a result, the different programs that derive or predict them have adopted different conventions, which make comparison between their results difficult. Moreover, the rotation schemes are seldom completely described. Here, we summarize the different conventions, determine which ones are adopted by commonly used software packages, and establish the formal equivalencies between the different calculated Euler angles.

Measuring dynamic and kinetic information in the previously inaccessible supra-τ(c) window of nanoseconds to microseconds by solution NMR spectroscopy

Nuclear Magnetic Resonance (NMR) spectroscopy is a powerful tool that has enabled experimentalists to characterize molecular dynamics and kinetics spanning a wide range of time-scales from picoseconds to days. This review focuses on addressing the previously inaccessible supra-tc window (defined as τ(c) < supra-τ(c) < 40 μs; in which tc is the overall tumbling time of a molecule) from the perspective of local inter-nuclear vector dynamics extracted from residual dipolar couplings (RDCs) and from the perspective of conformational exchange captured by relaxation dispersion measurements (RD). The goal of the first section is to present a detailed analysis of how to extract protein dynamics encoded in RDCs and how to relate this information to protein functionality within the previously inaccessible supra-τ(c) window. In the second section, the current state of the art for RD is analyzed, as well as the considerable progress toward pushing the sensitivity of RD further into the supra-τ(c) scale by up to a factor of two (motion up to 25 μs). From the data obtained with these techniques and methodology, the importance of the supra-τ(c) scale for protein function and molecular recognition is becoming increasingly clearer as the connection between motion on the supra-τ(c) scale and protein functionality from the experimental side is further strengthened with results from molecular dynamics simulations.

Molecular conformations of crystalline L-cysteine determined with vibrational cross angle measurements

The molecular conformations of crystalline L-cysteine prepared in its orthorhombic form were determined by the vibrational cross angle measurements. Its major dihedral angles of chemical bonds determined by this method are consistent with the results from diffraction experiments. In addition, the relative orientations of the chemical bonds associated with the hydrogen atoms of the NH3(+) group and the thiol group are also determined. The results demonstrate that the vibrational cross angle method based on the multiple-mode approach can potentially become a structural tool for determining molecular conformations. The major challenges for the method to become a general molecular structural tool are discussed, and some approaches to address them are proposed.

Mapping the conformational stability of maltose binding protein at the residue scale using nuclear magnetic resonance hydrogen exchange experiments

Being able to differentiate local fluctuations from global folding-unfolding dynamics of a protein is of major interest for improving our understanding of structure-function determinants. The maltose binding protein (MBP), a protein that belongs to the maltose transport system, has a structure composed of two globular domains separated by a rigid-body "hinge bending". Here we determined, by using hydrogen exchange (HX) nuclear magnetic resonance experiments, the apparent stabilization free energies of 101 residues of MBP bound to β-cyclodextrin (MBP-βCD) under native conditions. We observed that the last helix of MBP (helix α14) has a lower protection factor than the rest of the protein. Further, HX experiments were performed using guanidine hydrochloride under subdenaturing conditions to discriminate between local fluctuations and global unfolding events and to determine the MBP-βCD energy landscape. The results show that helix α4 and a part of helices α5 and α6 are clearly grouped into a subdenaturing folding unit and represent a partially folded intermediate under native conditions. In addition, we observed that amide protons located in the hinge between the two globular domains share similar ΔG(gu)(app) and m values and should unfold simultaneously. These observations provide new points of view for improving our understanding of the thermodynamic stability and the mechanisms that drive folding-unfolding dynamics of proteins.

On the use of ultracentrifugal devices for sedimented solute NMR

We have recently proposed sedimented solute NMR (SedNMR) as a solid-state method to access biomolecules without the need of crystallization or other sample manipulation. The drawback of SedNMR is that samples are intrinsically diluted and this is detrimental for the signal intensity. Ultracentrifugal devices can be used to increase the amount of sample inside the rotor, overcoming the intrinsic sensitivity limitation of the method. We designed two different devices and we here report the directions for using such devices and the relevant equations for determining the parameters for sedimentation.

Exploration of multi-state conformational dynamics and underlying global functional landscape of maltose binding protein

An increasing number of biological machines have been revealed to have more than two macroscopic states. Quantifying the underlying multiple-basin functional landscape is essential for understanding their functions. However, the present models seem to be insufficient to describe such multiple-state systems. To meet this challenge, we have developed a coarse grained triple-basin structure-based model with implicit ligand. Based on our model, the constructed functional landscape is sufficiently sampled by the brute-force molecular dynamics simulation. We explored maltose-binding protein (MBP) which undergoes large-scale domain motion between open, apo-closed (partially closed) and holo-closed (fully closed) states responding to ligand binding. We revealed an underlying mechanism whereby major induced fit and minor population shift pathways co-exist by quantitative flux analysis. We found that the hinge regions play an important role in the functional dynamics as well as that increases in its flexibility promote population shifts. This finding provides a theoretical explanation of the mechanistic discrepancies in PBP protein family. We also found a functional “backtracking” behavior that favors conformational change. We further explored the underlying folding landscape in response to ligand binding. Consistent with earlier experimental findings, the presence of ligand increases the cooperativity and stability of MBP. This work provides the first study to explore the folding dynamics and functional dynamics under the same theoretical framework using our triple-basin functional model.

A central goal of biology is to understand the function of the organism and its constituent parts at each of its scales of complexity. Function at the molecular level is often realized by changes in conformation. Unfortunately, experimental explorations of global motions critical for functional conformational changes are still challenging. In the present work, we developed a coarse grained triple-well structure-based model to explore the underlying functional landscape of maltose-binding protein (MBP). By quantitative flux analysis, we uncover the underlying mechanism by which the major induced fit and minor population shift pathways co-exist. Though we have previously lent credence to the assertion that dynamical equilibrium between open and minor closed conformations exist for all the free PBPs, the generality of this rule is still a matter of open debate. We found that the hinge flexibility is favorable to population shift mechanism. This finding provides a theoretical explanation of the mechanism discrepancies in PBP protein family. We also simulated the folding dynamics using this functional multi-basin model which successfully reproduced earlier protein melting experiment. This represents an exciting opportunity to characterize the interplay between folding and function, which is a long-standing question in the community. The theoretical approach employed in this study is general and can be applied to other systems.

Maximum occurrence analysis of protein conformations for different distributions of paramagnetic metal ions within flexible two-domain proteins

Multidomain proteins are composed of rigid domains connected by (flexible) linkers. Therefore, the domains may experience a large degree of reciprocal reorientation. Pseudocontact shifts and residual dipolar couplings arising from one or more paramagnetic metals successively placed in a single metal binding site in the protein can be used as restraints to assess the degree of mobility of the different domains. They can be used to determine the maximum occurrence (MO) of each possible protein conformation, i.e. the maximum weight that such conformations can have independently of the real structural ensemble, in agreement with the provided restraints. In the case of two-domain proteins, the metal ions can be placed all in the same domain, or distributed between the two domains. It has been demonstrated that the quantity of independent information for the characterization of the system is larger when all metals are bound in the same domain. At the same time, it has been shown that there are practical advantages in placing the metals in different domains. Here, it is shown that distributing the metals between the domains provides a tool for defining a coefficient of compatibility among the restraints obtained from different metals, without a significant decrease of the capability of the MO values to discriminate among conformations with different weights.

Advances in nuclear magnetic resonance for drug discovery

Nuclear Magnetic Resonance (NMR) techniques are widely used in the drug discovery process. The primary feature exploited in these investigations is the large difference in mass between drugs and receptors (usually proteins) and the effect this has on the rotational or translational correlation times for drugs bound to their targets. Many NMR parameters, such as the diffusion coefficient, spin diffusion, nuclear Overhauser enhancement, and transverse and longitudinal relaxation times, are strong functions of either the overall tumbling or translation of molecules in solution. This has led to the development of a wide variety of NMR techniques applicable to the elucidation of protein and nucleic acid structure in solution, the screening of drug candidates for binding to a target of choice, and the study of the conformational changes which occur in a target upon drug binding. High-throughput screening by NMR methods has recently received a boost from the introduction of sophisticated computational techniques for reducing the time needed for the acquisition of the primary NMR data for multidimensional studies.

Domain cooperativity in multidomain proteins: what can we learn from molecular alignment in anisotropic media?

Many proteins have modular design with multiple globular domains connected via flexible linkers. As a simple model of such system, we study a tandem construct consisting of two identical SH3 domains and a variable-length Gly/Ser linker. When the linker is short, this construct represents a dumbbell-shaped molecule with limited amount of domain-domain mobility. Due to its elongated shape, this molecule efficiently aligns in steric alignment media. As the length of the linker increases, the two domains become effectively uncoupled and begin to behave as independent entities. Consequently, their degree of alignment drops, approaching that found in the (near-spherical) isolated SH3 domains. To model the dependence of alignment parameters on the length of the interdomain linker, we have generated in silico a series of conformational ensembles representing SH3 tandems with different linker length. These ensembles were subsequently used as input for alignment prediction software PALES. The predicted alignment tensors were compared with the results of experimental measurements using a series of tandem-SH3 samples in PEG/hexanol alignment media. This comparison broadly confirmed the expected trends. At the same time, it has been found that the isolated SH3 domain aligns much stronger than expected. This finding can be attributed to complex morphology of the PEG/hexanol media and/or to weak site-specific interactions between the protein and the media. In the latter case, there are strong indications that electrostatic interactions may play a role. The fact that PEG/hexanol does not behave as a simple steric media should serve as a caution for studies that use PALES as a quantitative prediction tool (especially for disordered proteins). Further progress in this area depends on our ability to accurately model the anisotropic media and its site-specific interactions with protein molecules. Once this ability is improved, it should be possible to use the alignment parameters as a measure of domain-domain cooperativity, thus identifying the situations where two domains transiently interact with each other or become coupled through a partially structured linker.

Measurement of (1)H-(15)N and (1)H-(13)C residual dipolar couplings in nucleic acids from TROSY intensities

Analogous to the recently introduced ARTSY method for measurement of one-bond (1)H-(15)N residual dipolar couplings (RDCs) in large perdeuterated proteins, we introduce methods for measurement of base (13)C-(1)H and (15)N-(1)H RDCs in protonated nucleic acids. Measurements are based on quantitative analysis of intensities in (1)H-(15)N and (13)C-(1)H TROSY-HSQC spectra, and are illustrated for a 71-nucleotide adenine riboswitch. Results compare favorably with those of conventional frequency-based measurements in terms of completeness and convenience of use. The ARTSY method derives the size of the coupling from the ratio of intensities observed in two TROSY-HSQC spectra recorded with different dephasing delays, thereby minimizing potential resonance overlap problems. Precision of the RDC measurements is limited by the signal-to-noise ratio, S/N, achievable in the 2D TROSY-HSQC reference spectrum, and is approximately given by 30/(S/N) Hz for (15)N-(1)H and 65/(S/N) Hz for (13)C-(1)H. The signal-to-noise ratio of both (1)H-(15)N and (1)H-(13)C spectra greatly benefits when water magnetization during the experiments is not perturbed, such that rapid magnetization transfer from bulk water to the nucleic acid, mediated by rapid amino and hydroxyl hydrogen exchange coupled with (1)H-(1)H NOE transfer, allows for fast repetition of the experiment. RDCs in the mutated helix 1 of the riboswitch are compatible with nucleotide-specifically modeled, idealized A-form geometry and a static orientation relative to the helix 2/3 pair, which differs by ca 6° relative to the X-ray structure of the native riboswitch.

Exploring sparsely populated states of macromolecules by diamagnetic and paramagnetic NMR relaxation

Sparsely populated states of macromolecules, characterized by short lifetimes and high free-energies relative to the predominant ground state, often play a key role in many biological, chemical, and biophysical processes. In this review, we briefly summarize various new developments in NMR spectroscopy that permit these heretofore invisible, sparsely populated states to be detected, characterized, and in some instances visualized. Relaxation dispersion spectroscopy yields detailed kinetic information on processes involving species characterized by distinct chemical shifts with lifetimes in the ∼50 μs-10 ms range and populations as low as 0.5%. In the fast exchange regime (time scale less than ∼250-500 μs), the footprint of sparsely populated states can be observed on paramagnetic relaxation enhancement profiles measured on the resonances of the major species, thereby yielding structural information that is directly related to paramagnetic center-nuclei distances from which it is possible, under suitable circumstances, to compute a structure or ensemble of structures for the minor species. Finally, differential transverse relaxation measurements can be used to detect lifetime broadening effects that directly reflect the unidirectional rates for the conversion of NMR-visible into high-molecular weight NMR-invisible species. Examples of these various approaches are presented.

A structural and functional analysis of type III periplasmic and substrate binding proteins: their role in bacterial siderophore and heme transport

In Escherichia coli the Fhu, Fep and Fec transport systems are involved in the uptake of chelated ferric iron-siderophore complexes, whereas in pathogenic strains heme can also be used as an iron source. An essential step in these pathways is the movement of the ferric-siderophore complex or heme from the outer membrane transporter across the periplasm to the cognate cytoplasmic membrane ATP-dependent transporter. This is accomplished in each case by a dedicated periplasmic binding protein (PBP). Ferric-siderophore binding PBPs belong to the PBP protein superfamily and adopt a bilobal type III structural fold in which the two independently folded amino and carboxy terminal domains are linked together by a single long α-helix of approximately 20 amino acids. Recent structural studies reveal how the PBPs of the Fhu, Fep, Fec and Chu systems are able to bind their corresponding ligands. These complex structures will be discussed and placed in the context of our current understanding of the entire type III family of Gram-negative periplasmic binding proteins and related Gram-positive substrate binding proteins.

Fast approximations of the rotational diffusion tensor and their application to structural assembly of molecular complexes

We present and evaluate a rigid-body, deterministic, molecular docking method, called ELMDOCK, that relies solely on the three-dimensional structure of the individual components and the overall rotational diffusion tensor of the complex, obtained from nuclear spin-relaxation measurements. We also introduce a docking method, called ELMPATIDOCK, derived from ELMDOCK and based on the new concept of combining the shape-related restraints from rotational diffusion with those from residual dipolar couplings, along with ambiguous contact/interface-related restraints obtained from chemical shift perturbations. ELMDOCK and ELMPATIDOCK use two novel approximations of the molecular rotational diffusion tensor that allow computationally efficient docking. We show that these approximations are accurate enough to properly dock the two components of a complex without the need to recompute the diffusion tensor at each iteration step. We analyze the accuracy, robustness, and efficiency of these methods using synthetic relaxation data for a large variety of protein-protein complexes. We also test our method on three protein systems for which the structure of the complex and experimental relaxation data are available, and analyze the effect of flexible unstructured tails on the outcome of docking. Additionally, we describe a method for integrating the new approximation methods into the existing docking approaches that use the rotational diffusion tensor as a restraint. The results show that the proposed docking method is robust against experimental errors in the relaxation data or structural rearrangements upon complex formation and is computationally more efficient than current methods. The developed approximations are accurate enough to be used in structure refinement protocols.

Nuclear magnetic resonance provides a quantitative description of protein conformational flexibility on physiologically important time scales

A complete description of biomolecular activity requires an understanding of the nature and the role of protein conformational dynamics. In recent years, novel nuclear magnetic resonance-based techniques that provide hitherto inaccessible detail concerning biomolecular motions occurring on physiologically important time scales have emerged. Residual dipolar couplings (RDCs) provide precise information about time- and ensemble-averaged structural and dynamic processes with correlation times up to the millisecond and thereby encode key information for understanding biological activity. In this review, we present the application of two very different approaches to the quantitative description of protein motion using RDCs. The first is purely analytical, describing backbone dynamics in terms of diffusive motions of each peptide plane, using extensive statistical analysis to validate the proposed dynamic modes. The second is based on restraint-free accelerated molecular dynamics simulation, providing statistically sampled free energy-weighted ensembles that describe conformational fluctuations occurring on time scales from pico- to milliseconds, at atomic resolution. Remarkably, the results from these two approaches converge closely in terms of distribution and absolute amplitude of motions, suggesting that this kind of combination of analytical and numerical models is now capable of providing a unified description of protein conformational dynamics in solution.

NMR spectroscopy to study the dynamics and interactions of CFTR

The cystic fibrosis transmembrane conductance regulator (CFTR) is a multi-domain membrane chloride channel whose activity is regulated by ATP at two nucleotide-binding domains (NBD1 and NBD2) and by phosphorylation of the regulatory (R) region. The NBDs and the R region have functionally relevant motions that are critical for channel gating. Nuclear magnetic resonance (NMR) spectroscopy is a highly useful technique for obtaining information on the structure and interactions of CFTR and is extremely powerful for probing dynamics. NMR approaches for studying CFTR are reviewed, using our previous NBD1 and the R region results to provide examples. These NMR data are yielding insights into the dynamic properties and interactions that facilitate normal CFTR regulation as well as pathological effects of mutations, including the most common disease mutant, deletion of F508 in NBD1.

Contemporary methods in structure determination of membrane proteins by solution NMR

Integral membrane proteins are vital to life, being responsible for information and material exchange between a cell and its environment. Although high-resolution structural information is needed to understand how these functions are achieved, membrane proteins remain an under-represented subset of the protein structure databank. Solution NMR is increasingly demonstrating its ability to help address this knowledge shortfall, with the development of a diverse array of techniques to counter the challenges presented by membrane proteins. Here we document the advances that are helping to define solution NMR as an effective tool for membrane protein structure determination. Developments introduced over the last decade in the production of isotope-labeled samples, reconstitution of these samples into the growing selection of NMR-compatible membrane-mimetic systems, and the approaches used for the acquisition and application of structural restraints from these complexes are reviewed.

Structural assembly of molecular complexes based on residual dipolar couplings

We present and evaluate a rigid-body molecular docking method, called PATIDOCK, that relies solely on the three-dimensional structure of the individual components and the experimentally derived residual dipolar couplings (RDCs) for the complex. We show that, given an accurate ab initio predictor of the alignment tensor from a protein structure, it is possible to accurately assemble a protein-protein complex by utilizing the RDCs' sensitivity to molecular shape to guide the docking. The proposed docking method is robust against experimental errors in the RDCs and computationally efficient. We analyze the accuracy and efficiency of this method using experimental or synthetic RDC data for several proteins, as well as synthetic data for a large variety of protein-protein complexes. We also test our method on two protein systems for which the structure of the complex and steric-alignment data are available (Lys48-linked diubiquitin and a complex of ubiquitin and a ubiquitin-associated domain) and analyze the effect of flexible unstructured tails on the outcome of docking. The results demonstrate that it is fundamentally possible to assemble a protein-protein complex solely on the basis of experimental RDC data and the prediction of the alignment tensor from 3D structures. Thus, despite the purely angular nature of RDCs, they can be converted into intermolecular distance/translational constraints. Additionally, we show a method for combining RDCs with other experimental data, such as ambiguous constraints from interface mapping, to further improve structure characterization of protein complexes.

Protein alignment using cellulose nanocrystals: practical considerations and range of application

Cellulose nanocrystals (CNCs) form liquid crystals in aqueous solution that confer alignment to macromolecules and permit the measurement of residual dipolar couplings. CNCs possess many attractive features as an alignment medium. They are inexpensive, non-toxic, chemically inert, and robust to denaturants and temperature. Despite these advantages, CNCs are seldom employed as an alignment medium and the range of their applicability has not yet been explored. We have re-examined the use of CNCs in biomolecular NMR by analyzing the effects concentration, ionic strength, and temperature on molecular alignment. Stable alignment was obtained over wide ranges of temperature (10-70 degrees C) and pH (2.5-8.0), which makes CNCs potentially very useful in studies of thermophilic proteins and acid-stabilized molecules. Notably, we find that CNC suspensions are very sensitive to the concentrations of biological buffers, which must be taken into account when they are used in NMR analyses. These results have led us to develop a general procedure for preparing aligned samples with CNCs. Using the SH3 domain from the Fyn tyrosine kinase as a model system, we find that CNCs produce an alignment frame collinear with that of the commonly used Pf1 bacteriophage alignment medium, but of opposite magnitude.

Solution NMR studies of Chlorella virus DNA ligase-adenylate

DNA ligases are essential guardians of genome integrity by virtue of their ability to recognize and seal 3'-OH/5'-phosphate nicks in duplex DNA. The substrate binding and three chemical steps of the ligation pathway are coupled to global and local changes in ligase structure, involving both massive protein domain movements and subtle remodeling of atomic contacts in the active site. Here we applied solution NMR spectroscopy to study the conformational dynamics of the Chlorella virus DNA ligase (ChVLig), a minimized eukaryal ATP-dependent ligase consisting of nucleotidyltransferase, OB, and latch domains. Our analysis of backbone (15)N spin relaxation and (15)N,(1)H residual dipolar couplings of the covalent ChVLig-AMP intermediate revealed conformational sampling on fast (picosecond to nanosecond) and slow timescales (microsecond to millisecond), indicative of interdomain and intradomain flexibility. We identified local and global changes in ChVLig-AMP structure and dynamics induced by phosphate. In particular, the chemical shift perturbations elicited by phosphate were clustered in the peptide motifs that comprise the active site. We hypothesize that phosphate anion mimics some of the conformational transitions that occur when ligase-adenylate interacts with the nick 5'-phosphate.

Constructing atomic-resolution RNA structural ensembles using MD and motionally decoupled NMR RDCs

A broad structural landscape often needs to be characterized in order to fully understand how regulatory RNAs perform their biological functions at the atomic level. We present a protocol for visualizing thermally accessible RNA conformations at atomic-resolution and with timescales extending up to milliseconds. The protocol combines molecular dynamics (MD) simulations with experimental residual dipolar couplings (RDCs) measured in partially aligned (13)C/(15)N isotopically enriched elongated RNA samples. The structural ensembles generated in this manner provide insights into RNA dynamics and its role in functionally important transitions.

Probing mutation-induced structural perturbations by refinement against residual dipolar couplings: application to the U4 spliceosomal RNP complex

Confident interpretation of biochemical experiments performed with mutated proteins relies on verification of the integrity of the mutant structures. We present a simple and rapid refinement protocol for comparing the structures of mutated and wild-type proteins. Our approach involves measurement of residual dipolar couplings, and only requires assignment of the backbone resonances of the mutant species. We demonstrate application of the protocol to a mutant of the 15.5K protein, a core component of the U4 spliceosomal ribonucleoprotein (RNP) complex. Confirmation of the unperturbed structure of the mutated protein prompted re-examination of a previous mutagenesis study and indicated that the interpretation of mutant binding affinities in terms of direct interfacial contacts should be applied with caution.