NMR spectroscopy on domain dynamics in biomacromolecules

Normal mode analysis and comparative study of intrinsic dynamics of alcohol oxidase enzymes from GMC protein family

Glucose-Methanol-Choline (GMC) family enzymes are very important in catalyzing the oxidation of a wide range of structurally diverse substrates. Enzymes that constitute the GMC family, share a common tertiary fold but < 25% sequence identity. Cofactor FAD, FAD binding signature motif, and similar structural scaffold of the active site are common features of oxidoreductase enzymes of the GMC family. Protein functionality mainly depends on protein three-dimensional structures and dynamics. In this study, we used the normal mode analysis method to search the intrinsic dynamics of GMC family enzymes. We have explored the dynamical behavior of enzymes with unique substrate catabolism and active site characteristics from different classes of the GMC family. Analysis of individual enzymes and comparative ensemble analysis of enzymes from different classes has shown conserved dynamic motion at FAD binding sites. The present study revealed that GMC enzymes share a strong dynamic similarity (Bhattacharyya coefficient >90% and root mean squared inner product >52%) despite low sequence identity across the GMC family enzymes. The study predicts that local deformation energy between atoms of the enzyme may be responsible for the catalysis of different substrates. This study may help that intrinsic dynamics can be used to make meaningful classifications of proteins or enzymes from different organisms.Communicated by Ramaswamy H. Sarma.

On the use of residual dipolar couplings in multi-state structure calculation of two-domain proteins

Residual dipolar couplings (RDCs) are powerful nuclear magnetic resonance (NMR) probes for the structure calculation of biomacromolecules. Typically, an alignment tensor that defines the orientation of the entire molecule relative to the magnetic field is determined either before refinement of individual bond vectors or simultaneously with this refinement. For single-domain proteins this approach works well since all bond vectors can be described within the same coordinate frame, which is given by the alignment tensor. However, novel approaches are sought after for systems where no universal alignment tensor can be used. Here, we present an approach that can be applied to two-domain proteins that enables the calculation of multiple states within each domain as well as with respect to the relative positions of the two domains.

Joint neutron/molecular dynamics vibrational spectroscopy reveals softening of HIV-1 protease upon binding of a tight inhibitor

Biomacromolecules are inherently dynamic, and their dynamics are interwoven into function. The fast collective vibrational dynamics in proteins occurs in the low picosecond timescale corresponding to frequencies of ∼5-50 cm^-1. This sub-to-low THz frequency regime covers the low-amplitude collective breathing motions of a whole protein and vibrations of the constituent secondary structure elements, such as α-helices, β-sheets and loops. We have used inelastic neutron scattering experiments in combination with molecular dynamics simulations to demonstrate the vibrational dynamics softening of HIV-1 protease, a target of HIV/AIDS antivirals, upon binding of a tight clinical inhibitor darunavir. Changes in the vibrational density of states of matching structural elements in the two monomers of the homodimeric protein are not identical, indicating asymmetric effects of darunavir on the vibrational dynamics. Three of the 11 major secondary structure elements contribute over 40% to the overall changes in the vibrational density of states upon darunavir binding. Molecular dynamics simulations informed by experiments allowed us to estimate that the altered vibrational dynamics of the protease would contribute -3.6 kcal mol^-1 at 300 K, or 25%, to the free energy of darunavir binding. As HIV-1 protease drug resistance remains a concern, our results open a new avenue to help establish a direct quantitative link between protein vibrational dynamics and drug resistance.

Reconstruction of Coupled Intra- and Interdomain Protein Motion from Nuclear and Electron Magnetic Resonance

Proteins composed of multiple domains allow for structural heterogeneity and interdomain dynamics that may be vital for function. Intradomain structures and dynamics can influence interdomain conformations and vice versa. However, no established structure determination method is currently available that can probe the coupling of these motions. The protein Pin1 contains separate regulatory and catalytic domains that sample "extended" and "compact" states, and ligand binding changes this equilibrium. Ligand binding and interdomain distance have been shown to impact the activity of Pin1, suggesting interdomain allostery. In order to characterize the conformational equilibrium of Pin1, we describe a novel method to model the coupling between intra- and interdomain dynamics at atomic resolution using multistate ensembles. The method uses time-averaged nuclear magnetic resonance (NMR) restraints and double electron-electron resonance (DEER) data that resolve distance distributions. While the intradomain calculation is primarily driven by exact nuclear Overhauser enhancements (eNOEs), J couplings, and residual dipolar couplings (RDCs), the relative domain distribution is driven by paramagnetic relaxation enhancement (PREs), RDCs, interdomain NOEs, and DEER. Our data support a 70:30 population of the compact and extended states in apo Pin1. A multistate ensemble describes these conformations simultaneously, with distinct conformational differences located in the interdomain interface stabilizing the compact or extended states. We also describe correlated conformations between the catalytic site and interdomain interface that may explain allostery driven by interdomain contact.

Heterogeneous dynamics in partially disordered proteins

Importance of disordered protein regions is increasingly recognized in biology, but their characterization remains challenging due to the lack of suitable experimental and theoretical methods. NMR experiments can detect multiple timescale dynamics and structural details of disordered protein regions, but their detailed interpretation is often difficult. Here we combine protein backbone 15N spin relaxation data with molecular dynamics (MD) simulations to detect not only heterogeneous dynamics of large partially disordered proteins but also their conformational ensembles. We observed that the rotational dynamics of folded regions in partially disordered proteins is dominated by similar rigid body rotation as in globular proteins, thereby being largely independent of flexible disordered linkers. Disordered regions, on the other hand, exhibit complex rotational motions with multiple timescales below ∼30 ns which are difficult to detect from experimental data alone, but can be captured by MD simulations. Combining MD simulations and backbone 15N spin relaxation data, measured applying segmental isotopic labeling with salt-inducible split intein, we resolved the conformational ensemble and dynamics of partially disordered periplasmic domain of TonB protein from Helicobacter pylori containing 250 residues. To demonstrate the universality of our approach, it was applied also to the partially disordered region of chicken Engrailed 2. Our results pave the way in understanding how TonB transfers energy from inner membrane to the outer membrane receptors in Gram-negative bacteria, as well as the function of other proteins with disordered domains.

Revisiting biomolecular NMR spectroscopy for promoting small-molecule drug discovery

Recently, there has been increasing interest in new modalities such as therapeutic antibodies and gene therapy at a number of pharmaceutical companies. Moreover, in small-molecule drug discovery at such companies, efforts have focused on hard-to-drug targets such as inhibiting protein-protein interactions. Biomolecular NMR spectroscopy has been used in drug discovery in a variety of ways, such as for the reliable detection of binding and providing three-dimensional structural information for structure-based drug design. The advantages of using NMR spectroscopy have been known for decades (Jahnke in J Biomol NMR 39:87-90, (2007); Gossert and Jahnke in Prog Nucl Magn Reson Spectrosc 97:82-125, (2016)). For tackling hard-to-drug targets and increasing the success in discovering drug molecules, in-depth analysis of drug-target protein interactions performed by biophysical methods will be more and more essential. Here, we review the advantages of NMR spectroscopy as a key technology of biophysical methods and also discuss issues such as using cutting-edge NMR spectrometers and increasing the demand of utilizing conformational dynamics information for promoting small-molecule drug discovery.

Observation of the Unbiased Conformers of Putative DNA-Scaffold Ribosugars

The constitution, configuration, and flexibility of the core sugars of DNA molecules alter their function in diverse roles. Conformational itineraries of the ribofuranosides (fs) have long been known to finely determine rates of processing, yet we also know that, strikingly, semifunctional DNAs containing pyranosides (ps) or other configurations can be created, suggesting sufficient but incompletely understood plasticity. The multiple conformers involved in such processes are necessarily influenced by context and environment: solvent, hosts, ligands. Notably, however, to date the unbiased, "naked" conformers have not been experimentally determined. Here, the inherent conformational biases of DNA scaffold deoxyribosides in unsolvated and solvated forms have now been defined using gas-phase microwave and solution-phase NMR spectroscopies coupled with computational analyses and exploitation of critical differences between natural-abundance isotopologues. Serial determination of precise, individual spectra for conformers of these 25 isotopologues in alpha (α-d) and beta (β-d); pyrano (p) and furano (f) methyl 2-deoxy-d-ribosides gave not only unprecedented atomic-level resolution structures of associated conformers but also their quantitative populations. Together these experiments revealed that typical ²E and ³E conformations of the sugar found in complex DNA structures are not inherently populated. Moreover, while both OH-5' and OH-3' are constrained by intramolecular hydrogen bonding in the unnatural αf scaffold, OH-3' is "born free" in the "naked" lowest lying energy conformer of natural βf. Consequently, upon solvation, unnatural αf is strikingly less perturbable (retaining ²T₁ conformation in vacuo and water) than natural βf. Unnatural αp and βp ribosides also display low conformational perturbability. These first experimental data on inherent, unbiased conformers therefore suggest that it is the background of conformational flexibility of βf that may have led to its emergence out of multiple possibilities as the sugar scaffold for "life's code" and suggest a mechanism by which the resulting freedom of OH-3' (and hence accessibility as a nucleophile) in βf may drive preferential processing and complex structure formation, such as replicative propagation of DNA from 5'-to-3'.

Structural and Mutagenesis Studies Evince the Role of the Extended Protuberant Domain of Ribosomal Protein uL10 in Protein Translation

The lateral stalk of ribosomes constitutes the GTPase-associated center and is responsible for recruiting translation factors to the ribosomes. The eukaryotic stalk contains a P-complex, in which one molecule of uL10 (formerly known as P0) protein binds two copies of P1/P2 heterodimers. Unlike bacterial uL10, eukaryotic uL10 has an extended protuberant (uL10ext) domain inserted into the N-terminal RNA-binding domain. Here, we determined the solution structure of the extended protuberant domain of Bombyx mori uL10 by nuclear magnetic resonance spectroscopy. Comparison of the structures of the B. mori uL10ext domain with eRF1-bound and eEF2-bound ribosomes revealed significant structural rearrangement in a "hinge" region surrounding Phe183, a residue conserved in eukaryotic but not in archaeal uL10. ¹⁵N relaxation analyses showed that residues in the hinge region have significantly large values of transverse relaxation rates. To test the role of the conserved phenylalanine residue, we created a yeast mutant strain expressing an F181A variant of uL10. An in vitro translation assay showed that the alanine substitution increased the level of polyphenylalanine synthesis by ∼33%. Taken together, our results suggest that the hinge motion of the uL10ext domain facilitates the binding of different translation factors to the GTPase-associated center during protein synthesis.

Cross-correlated relaxation rates between protein backbone H-X dipolar interactions

The relaxation interference between dipole-dipole interactions of two separate spin pairs carries structural and dynamics information. In particular, when compared to individual dynamic behavior of those spin pairs, such cross-correlated relaxation (CCR) rates report on the correlation between the spin pairs. We have recently mapped out correlated motion along the backbone of the protein GB3, using CCR rates among and between consecutive H^N-N and H^α-C^α dipole-dipole interactions. Here, we provide a detailed account of the measurement of the four types of CCR rates. All rates were obtained from at least two different pulse sequences, of which the yet unpublished ones are presented. Detailed comparisons between the different methods and corrections for unwanted pathways demonstrate that the averaged CCR rates are highly accurate and precise with errors of 1.5-3% of the entire value ranges.

A critical assessment of methods to recover information from averaged data

Conformational heterogeneity is key to the function of many biomacromolecules, but only a few groups have tried to characterize it until recently. Now, thanks to the increased throughput of experimental data and the increased computational power, the problem of the characterization of protein structural variability has become more and more popular. Several groups have devoted their efforts in trying to create quantitative, reliable and accurate protocols for extracting such information from averaged data. We analyze here different approaches, discussing strengths and weaknesses of each. All approaches can roughly be clustered into two groups: those satisfying the maximum entropy principle and those recovering ensembles composed of a restricted number of molecular conformations. In the first case, the solution focuses on the features that are common to all the infinite solutions satisfying the experimental data; in the second case, the reconstructed ensemble shows the conformational regions where a large probability can be placed. The upper limits for conformational probabilities (MaxOcc) can also be calculated. We also give an overview of the mainstream experimental observables, with considerations on the assumptions underlying their usage.

Unveiling Inherent Degeneracies in Determining Population-Weighted Ensembles of Interdomain Orientational Distributions Using NMR Residual Dipolar Couplings: Application to RNA Helix Junction Helix Motifs

A growing number of studies employ time-averaged experimental data to determine dynamic ensembles of biomolecules. While it is well-known that different ensembles can satisfy experimental data to within error, the extent and nature of these degeneracies, and their impact on the accuracy of the ensemble determination remains poorly understood. Here, we use simulations and a recently introduced metric for assessing ensemble similarity to explore degeneracies in determining ensembles using NMR residual dipolar couplings (RDCs) with specific application to A-form helices in RNA. Various target ensembles were constructed representing different domain-domain orientational distributions that are confined to a topologically restricted (<10%) conformational space. Five independent sets of ensemble averaged RDCs were then computed for each target ensemble and a "sample and select" scheme used to identify degenerate ensembles that satisfy RDCs to within experimental uncertainty. We find that ensembles with different ensemble sizes and that can differ significantly from the target ensemble (by as much as ∑Ω ∼ 0.4 where ∑Ω varies between 0 and 1 for maximum and minimum ensemble similarity, respectively) can satisfy the ensemble averaged RDCs. These deviations increase with the number of unique conformers and breadth of the target distribution, and result in significant uncertainty in determining conformational entropy (as large as 5 kcal/mol at T = 298 K). Nevertheless, the RDC-degenerate ensembles are biased toward populated regions of the target ensemble, and capture other essential features of the distribution, including the shape. Our results identify ensemble size as a major source of uncertainty in determining ensembles and suggest that NMR interactions such as RDCs and spin relaxation, on their own, do not carry the necessary information needed to determine conformational entropy at a useful level of precision. The framework introduced here provides a general approach for exploring degeneracies in ensemble determination for different types of experimental data.

Solution structure of the RNA-binding cold-shock domain of the Chlamydomonas reinhardtii NAB1 protein and insights into RNA recognition

Light-harvesting complex (LHC) proteins are among the most abundant proteins on Earth and play critical roles in photosynthesis, both in light capture and in photoprotective mechanisms. The Chlamydomonas reinhardtii nucleic acid-binding protein 1 (NAB1) is a negative regulator of LHC protein translation. Its N-terminal cold-shock domain (CSD) binds to a 13-nt element [CSD consensus sequence (CSDCS)] found in the mRNA of specific LHC proteins associated with Photosystem II (PSII), an interaction which regulates LHC expression and, consequently, PSII-associated antenna size, structure and function. In the present study, we elucidated the solution structure of the NAB1 CSD as determined by heteronuclear NMR. The CSD adopts a characteristic five-stranded anti parallel β-barrel fold. Upon addition of CSDCS RNA, a large number of NMR chemical shift perturbations were observed, corresponding primarily to surface-exposed residues within the highly conserved β2- and β3-strands in the canonical RNA-binding region, but also to residues on β-strand 5 extending the positive surface patch and the overall RNA-binding site. Additional chemical shift perturbations that accompanied RNA binding involved buried residues, suggesting that transcript recognition is accompanied by conformational change. Our results indicate that NAB1 associates with RNA transcripts through a mechanism involving its CSD that is conserved with mechanisms of sequence-specific nucleic acid recognition employed by ancestrally related bacterial cold-shock proteins (CSPs).

NMR studies of dynamic biomolecular conformational ensembles

Multidimensional heteronuclear NMR approaches can provide nearly complete sequential signal assignments of isotopically enriched biomolecules. The availability of assignments together with measurements of spin relaxation rates, residual spin interactions, J-couplings and chemical shifts provides information at atomic resolution about internal dynamics on timescales ranging from ps to ms, both in solution and in the solid state. However, due to the complexity of biomolecules, it is not possible to extract a unique atomic-resolution description of biomolecular motions even from extensive NMR data when many conformations are sampled on multiple timescales. For this reason, powerful computational approaches are increasingly applied to large NMR data sets to elucidate conformational ensembles sampled by biomolecules. In the past decade, considerable attention has been directed at an important class of biomolecules that function by binding to a wide variety of target molecules. Questions of current interest are: "Does the free biomolecule sample a conformational ensemble that encompasses the conformations found when it binds to various targets; and if so, on what time scale is the ensemble sampled?" This article reviews recent efforts to answer these questions, with a focus on comparing ensembles obtained for the same biomolecules by different investigators. A detailed comparison of results obtained is provided for three biomolecules: ubiquitin, calmodulin and the HIV-1 trans-activation response RNA.

Local isotropic diffusion approximation for coupled internal and overall molecular motions in NMR spin relaxation

The present work demonstrates that NMR spin relaxation rate constants for molecules interconverting between states with different diffusion tensors can be modeled theoretically by combining orientational correlation functions for exchanging spherical molecules with locally isotropic approximations for the diffusion anisotropic tensors. The resulting expressions are validated by comparison with correlation functions obtained by Monte Carlo simulations and are accurate for moderate degrees of diffusion anisotropy typically encountered in investigations of globular proteins. The results are complementary to an elegant, but more complex, formalism that is accurate for all degrees of diffusion anisotropy [Ryabov, Y.; Clore, G. M.; Schwieters, C. D. J. Chem. Phys. 2012, 136, 034108].

Probing the structure of human protein disulfide isomerase by chemical cross-linking combined with mass spectrometry

Protein disulfide-isomerase (PDI) is a four-domain flexible protein that catalyzes the formation of disulfide bonds in the endoplasmic reticulum. Here we have analyzed native PDI purified from human placenta by chemical cross-linking followed by mass spectrometry (CXMS). In addition to PDI the sample contained soluble calnexin and ERp72. Extensive cross-linking was observed within the PDI molecule, both intra- and inter-domain, as well as between the different components in the mixture. The high sensitivity of the analysis in the current experiments, combined with a likely promiscuous interaction pattern of the involved proteins, revealed relatively densely populated cross-link heat maps. The established X-ray structure of the monomeric PDI could be confirmed; however, the dimer as presented in the existing models does not seem to be prevalent in solution as modeling on the observed cross-links revealed new models of dimeric PDI. The observed inter-protein cross-links confirmed the existence of a peptide binding area on calnexin that binds strongly both PDI and ERp72. On the other hand, interaction sites on PDI and ERp72 could not be uniquely identified, indicating a more non-specific interaction pattern.

BIOLOGICAL SIGNIFICANCE

The present work demonstrates the use of chemical cross-linking and mass spectrometry (CXMS) for the determination of a solution structure of natural human PDI and its interaction with the chaperones ERp72 and calnexin. The data shows that the dimeric structure of PDI may be more diverse than indicated by present models. We further observe that the temperature influences the cross-linking pattern of PDI, but this does not influence the overall folding pattern of the molecule.