Validation of X-ray Crystal Structure Ensemble Representations of SARS-CoV-2 Main Protease by Solution NMR Residual Dipolar Couplings

Considerable debate has focused on whether sampling of molecular dynamics trajectories restrained by crystallographic data can be used to develop realistic ensemble models for proteins in their natural, solution state. For the SARS-CoV-2 main protease, Mpro, we evaluated agreement between solution residual dipolar couplings (RDCs) and various recently reported multi-conformer and dynamic-ensemble crystallographic models. Although Phenix-derived ensemble models showed only small improvements in crystallographic Rfree, substantially improved RDC agreement over fits to a conventionally refined 1.2-Å X-ray structure was observed, in particular for residues with above average disorder in the ensemble. For a set of six lower resolution (1.55-2.19 Å) Mpro X-ray ensembles, obtained at temperatures ranging from 100 to 310 K, no significant improvement over conventional two-conformer representations was found. At the residue level, large differences in motions were observed among these ensembles, suggesting high uncertainties in the X-ray derived dynamics. Indeed, combining the six ensembles from the temperature series with the two 1.2-Å X-ray ensembles into a single 381-member "super ensemble" averaged these uncertainties and substantially improved agreement with RDCs. However, all ensembles showed excursions that were too large for the most dynamic fraction of residues. Our results suggest that further improvements to X-ray ensemble refinement are feasible, and that RDCs provide a sensitive benchmark in such endeavors. Remarkably, a weighted ensemble of 350 PDB Mpro X-ray structures provided slightly better cross-validated agreement with RDCs than any individual ensemble refinement, implying that differences in lattice confinement also limit the fit of RDCs to X-ray coordinates.

Predicting the Sequence-Dependent Backbone Dynamics of Intrinsically Disordered Proteins

Dynamics is a crucial link between sequence and function for intrinsically disordered proteins (IDPs). NMR spin relaxation is a powerful technique for characterizing the sequence-dependent backbone dynamics of IDPs. Of particular interest is the 15N transverse relaxation rate (R2), which reports on slower dynamics (10s of ns up to 1 μs and beyond). NMR and molecular dynamics (MD) simulations have shown that local interactions and secondary structure formation slow down backbone dynamics and raise R2. Elevated R2 has been suggested to be indicators of propensities of membrane association, liquid-liquid phase separation, and other functional processes. Here we present a sequence-based method, SeqDYN, for predicting R2 of IDPs. The R2 value of a residue is expressed as the product of contributing factors from all residues, which attenuate with increasing sequence distance from the central residue. The mathematical model has 21 parameters, representing the correlation length (where the attenuation is at 50%) and the amplitudes of the contributing factors of the 20 types of amino acids. Training on a set of 45 IDPs reveals a correlation length of 5.6 residues, aromatic and long branched aliphatic amino acids and Arg as R2 promotors whereas Gly and short polar amino acids as R2 suppressors. The prediction accuracy of SeqDYN is competitive against that of recent MD simulations using IDP-specific force fields. For a structured protein, SeqDYN prediction represents R2 in the unfolded state. SeqDYN is available as a web server at https://zhougroup-uic.github.io/SeqDYNidp/ for rapid R2 prediction.

Structural Diversity and Dynamics of Human Three-Finger Proteins Acting on Nicotinic Acetylcholine Receptors

Ly-6/uPAR or three-finger proteins (TFPs) contain a disulfide-stabilized β-structural core and three protruding loops (fingers). In mammals, TFPs have been found in epithelium and the nervous, endocrine, reproductive, and immune systems. Here, using heteronuclear NMR, we determined the three-dimensional (3D) structure and backbone dynamics of the epithelial secreted protein SLURP-1 and soluble domains of GPI-anchored TFPs from the brain (Lynx2, Lypd6, Lypd6b) acting on nicotinic acetylcholine receptors (nAChRs). Results were compared with the data about human TFPs Lynx1 and SLURP-2 and snake α-neurotoxins WTX and NTII. Two different topologies of the β-structure were revealed: one large antiparallel β-sheet in Lypd6 and Lypd6b, and two β-sheets in other proteins. α-Helical segments were found in the loops I/III of Lynx2, Lypd6, and Lypd6b. Differences in the surface distribution of charged and hydrophobic groups indicated significant differences in a mode of TFPs/nAChR interactions. TFPs showed significant conformational plasticity: the loops were highly mobile at picosecond-nanosecond timescale, while the β-structural regions demonstrated microsecond-millisecond motions. SLURP-1 had the largest plasticity and characterized by the unordered loops II/III and cis-trans isomerization of the Tyr39-Pro40 bond. In conclusion, plasticity could be an important feature of TFPs adapting their structures for optimal interaction with the different conformational states of nAChRs.

Crystal structure of the N domain of Lon protease from Mycobacterium avium complex

Lon protease is evolutionarily conserved in prokaryotes and eukaryotic organelles. The primary function of Lon is to selectively degrade abnormal and certain regulatory proteins to maintain the homeostasis in vivo. Lon mainly consists of three functional domains and the N-terminal domain is required for the substrate selection and recognition. However, the precise contribution of the N-terminal domain remains elusive. Here, we determined the crystal structure of the N-terminal 192-residue construct of Lon protease from Mycobacterium avium complex at 2.4 å resolution，and measured NMR-relaxation parameters of backbones. This structure consists of two subdomains, the β-strand rich N-terminal subdomain and the five-helix bundle of C-terminal subdomain, connected by a flexible linker，and is similar to the overall structure of the N domain of Escherichia coli Lon even though their sequence identity is only 26%. The obtained NMR-relaxation parameters reveal two stabilized loops involved in the structural packing of the compact N domain and a turn structure formation. The performed homology comparison suggests that structural and sequence variations in the N domain may be closely related to the substrate selectivity of Lon variants. Our results provide the structure and dynamics characterization of a new Lon N domain, and will help to define the precise contribution of the Lon N-terminal domain to the substrate recognition.

Structural insight into the unique dsDNA binding topology of the human ORC2 wing helix domain

The origin recognition complex (ORC) is indispensable for the initiation of DNA replication during the cell cycle. The DNA-binding modes of the human ORC winged helix domain (WHD) remain enigmatic, as the dsDNA recognition sites of archaeal and Saccharomyces cerevisiae ORC WHDs are distinct. Here, we solved the high-resolution crystal structure of the human ORC2 WHD, although its complex with dsDNA is difficult to crystallize due to its weak binding affinities. The near-complete NMR backbone assignments and chemical shift perturbations reveal a new dsDNA binding topology in addition to the conserved β-sheet hairpin region, in which residues show higher dynamics. The key interacting residues (R540, K548, and K549) were validated by mutagenesis studies. Our data suggest that the ORC2 WHD recognizes dsDNA sequences through its flexible β-sheet hairpin as an anchor point, while the rest of the protein adopts various orientations in different species. This weak but real interaction module identified by NMR is useful for the structural reconstruction of large biomolecular complexes using cryo-EM. The binding topology and dynamics of ORC2 WHDs were also underpinned by molecular dynamics simulations.

Dynamic Studies on Intrinsically Disordered Regions of Two Paralogous Transcription Factors Reveal Rigid Segments with Important Biological Functions

Long stretches of intrinsically disordered regions (IDRs) are abundantly present in eukaryotic transcription factors. Although their biological significance is well appreciated, the underlying structural and dynamic mechanisms of their function are still not clear. Using solution NMR spectroscopy, we have studied the structural and dynamic features of two paralogous HOX transcription factors, SCR and DFD, from Drosophila. Both proteins have a conserved DNA-binding homeodomain and a long stretch of functionally important IDR. Using NMR dynamics, we determined flexibility of each residue in these proteins. The flexibility of the residues in the disordered region is not uniform. In both proteins, the IDRs have short stretches of consecutive residues with relatively less flexibility, that is, higher rigidity. We show that one such rigid segment is specifically recognized by another co-transcription factor, thus highlighting the importance of these rigid segments in IDR-mediated protein-protein interactions. Using molecular dynamics simulation, we further show that the rigid segments sample less conformations compared to the rest of the residues in the disordered region. The restrained conformational sampling of these rigid residues should lower the loss in conformational entropy during their interactions with binding partners resulting in sequence specific binding. This work provides experimental evidence of a "rigid-segment" model of IDRs, where functionally important rigid segments are connected by highly flexible linkers. Furthermore, a comparative study of IDRs in paralogous proteins reveals that in spite of low-sequence conservation, the rigid and flexible segments are sequentially maintained to preserve related functions and regulations of these proteins.

NMR analysis of the backbone dynamics of the small GTPase Rheb and its interaction with the regulatory protein FKBP38

Ras homolog enriched in brain (Rheb) is a small GTPase that regulates mammalian/mechanistic target of rapamycin complex 1 (mTORC1) and, thereby, cell growth and metabolism. Here we show that cycling between the inactive GDP- and the active GTP-bound state modulates the backbone dynamics of a C-terminal truncated form, RhebΔCT, which is suggested to influence its interactions. We further investigated the interactions between RhebΔCT and the proposed Rheb-binding domain of the regulatory protein FKBP38. The observed weak interactions with the GTP-analogue- (GppNHp-) but not the GDP-bound state, appear to accelerate the GDP to GTP exchange, but only very weakly compared to a genuine GEF. Thus, FKBP38 is most likely not a GEF but a Rheb effector that may function in membrane targeting of Rheb.

Effects of Hinge-region Natural Polymorphisms on Human Immunodeficiency Virus-Type 1 Protease Structure, Dynamics, and Drug Pressure Evolution

Multidrug resistance to current Food and Drug Administration-approved HIV-1 protease (PR) inhibitors drives the need to understand the fundamental mechanisms of how drug pressure-selected mutations, which are oftentimes natural polymorphisms, elicit their effect on enzyme function and resistance. Here, the impacts of the hinge-region natural polymorphism at residue 35, glutamate to aspartate (E35D), alone and in conjunction with residue 57, arginine to lysine (R57K), are characterized with the goal of understanding how altered salt bridge interactions between the hinge and flap regions are associated with changes in structure, motional dynamics, conformational sampling, kinetic parameters, and inhibitor affinity. The combined results reveal that the single E35D substitution leads to diminished salt bridge interactions between residues 35 and 57 and gives rise to the stabilization of open-like conformational states with overall increased backbone dynamics. In HIV-1 PR constructs where sites 35 and 57 are both mutated (e.g. E35D and R57K), x-ray structures reveal an altered network of interactions that replace the salt bridge thus stabilizing the structural integrity between the flap and hinge regions. Despite the altered conformational sampling and dynamics when the salt bridge is disrupted, enzyme kinetic parameters and inhibition constants are similar to those obtained for subtype B PR. Results demonstrate that these hinge-region natural polymorphisms, which may arise as drug pressure secondary mutations, alter protein dynamics and the conformational landscape, which are important thermodynamic parameters to consider for development of inhibitors that target for non-subtype B PR.

Multilevel Changes in Protein Dynamics upon Complex Formation of the Calcium-Loaded S100A4 with a Nonmuscle Myosin IIA Tail Fragment

Dysregulation of Ca²⁺ -binding S100 proteins plays important role in various diseases. The asymmetric complex of Ca²⁺ -bound S100A4 with nonmuscle myosin IIA has high stability and highly increased Ca²⁺ affinity. Here we investigated the possible causes of this allosteric effect by NMR spectroscopy. Chemical shift-based secondary-structure analysis did not show substantial changes for the complex. Backbone dynamics revealed slow-timescale local motions in the H1 helices of homodimeric S100A4; these were less pronounced in the complex form and might be accompanied by an increase in dimer stability. Different mobilities in the Ca²⁺ -coordinating EF-hand sites indicate that they communicate by an allosteric mechanism operating through changes in protein dynamics; this must be responsible for the elevated Ca²⁺ affinity. These multilevel changes in protein dynamics as conformational adaptation allow S100A4 fine-tuning of its protein-protein interactions inside the cell during Ca²⁺ signaling.

Quantitative evaluation of positive ϕ angle propensity in flexible regions of proteins from three-bond J couplings

(3)JHNHα and (3)JC'C' couplings can be readily measured in isotopically enriched proteins and were shown to contain precise information on the backbone torsion angles, ϕ, sampled in disordered regions of proteins. However, quantitative interpretation of these couplings required the population of conformers with positive ϕ angles to be very small. Here, we demonstrate that this restriction can be removed by measurement of (3)JC'Hα values. Even though the functional forms of the (3)JC'Hα and (3)JHNHα Karplus equations are the same, large differences in their coefficients enable accurate determination of the fraction of time that positive ϕ angles are sampled. A four-dimensional triple resonance HACANH[C'] E.COSY experiment is introduced to simultaneously measure (3)JC'Hα and (3)JHNC' in the typically very congested spectra of disordered proteins. High resolution in these spectra is obtained by non-uniform sampling (in the 0.1-0.5% range). Application to the intrinsically disordered protein α-synuclein shows that while most residues have close-to-zero positive ϕ angle populations, up to 16% positive ϕ population is observed for Asn residues. Positive ϕ angle populations determined with the new approach agree closely with consensus values from protein coil libraries and prior analysis of a large set of other NMR parameters. The combination of (3)JHNC' and (3)JC'C' provides information about the amplitude of ϕ angle dynamics.

Changing the topology of protein backbone: the effect of backbone cyclization on the structure and dynamics of a SH3 domain

Understanding of the effects of the backbone cyclization on the structure and dynamics of a protein is essential for using protein topology engineering to alter protein stability and function. Here we have determined, for the first time, the structure and dynamics of the linear and various circular constructs of the N-SH3 domain from protein c-Crk. These constructs differ in the length and amino acid composition of the cyclization region. The backbone cyclization was carried out using intein-mediated intramolecular chemical ligation between the juxtaposed N- and the C-termini. The structure and backbone dynamics studies were performed using solution NMR. Our data suggest that the backbone cyclization has little effect on the overall three-dimensional structure of the SH3 domain: besides the termini, only minor structural changes were found in the proximity of the cyclization region. In contrast to the structure, backbone dynamics are significantly affected by the cyclization. On the subnanosecond time scale, the backbone of all circular constructs on average appears more rigid than that of the linear SH3 domain; this effect is observed over the entire backbone and is not limited to the cyclization site. The backbone mobility of the circular constructs becomes less restricted with increasing length of the circularization loop. In addition, significant conformational exchange motions (on the sub-millisecond time scale) were found in the N-Src loop and in the adjacent β-strands in all circular constructs studied in this work. These effects of backbone cyclization on protein dynamics have potential implications for the stability of the protein fold and for ligand binding.

Complete dissociation of the HIV-1 gp41 ectodomain and membrane proximal regions upon phospholipid binding

The envelope glycoprotein gp41 mediates the process of membrane fusion that enables entry of the HIV-1 virus into the host cell. Strong lipid affinity of the ectodomain suggests that its heptad repeat regions play an active role in destabilizing membranes by directly binding to the lipid bilayers and thereby lowering the free-energy barrier for membrane fusion. In such a model, immediately following the shedding of gp120, the N-heptad and C-heptad helices dissociate and melt into the host cell and viral membranes, respectively, pulling the destabilized membranes into juxtaposition, ready for fusion. Post-fusion, reaching the final 6-helix bundle (6 HB) conformation then involves competition between intermolecular interactions needed for formation of the symmetric 6 HB trimer and the membrane affinity of gp41's ectodomain, including its membrane-proximal regions. Our solution NMR study of the structural and dynamic properties of three constructs containing the ectodomain of gp41 with and without its membrane-proximal regions suggests that these segments do not form inter-helical interactions until the very late steps of the fusion process. Interactions between the polar termini of the heptad regions, which are not associating with the lipid surface, therefore may constitute the main driving force initiating formation of the final post-fusion states. The absence of significant intermolecular ectodomain interactions in the presence of dodecyl phosphocholine highlights the importance of trimerization of gp41's transmembrane helix to prevent complete dissociation of the trimer during the course of fusion.

Dynameomics: data-driven methods and models for utilizing large-scale protein structure repositories for improving fragment-based loop prediction

Protein function is intimately linked to protein structure and dynamics yet experimentally determined structures frequently omit regions within a protein due to indeterminate data, which is often due protein dynamics. We propose that atomistic molecular dynamics simulations provide a diverse sampling of biologically relevant structures for these missing segments (and beyond) to improve structural modeling and structure prediction. Here we make use of the Dynameomics data warehouse, which contains simulations of representatives of essentially all known protein folds. We developed novel computational methods to efficiently identify, rank and retrieve small peptide structures, or fragments, from this database. We also created a novel data model to analyze and compare large repositories of structural data, such as contained within the Protein Data Bank and the Dynameomics data warehouse. Our evaluation compares these structural repositories for improving loop predictions and analyzes the utility of our methods and models. Using a standard set of loop structures, containing 510 loops, 30 for each loop length from 4 to 20 residues, we find that the inclusion of Dynameomics structures in fragment-based methods improves the quality of the loop predictions without being dependent on sequence homology. Depending on loop length, ∼ 25-75% of the best predictions came from the Dynameomics set, resulting in lower main chain root-mean-square deviations for all fragment lengths using the combined fragment library. We also provide specific cases where Dynameomics fragments provide better predictions for NMR loop structures than fragments from crystal structures. Online access to these fragment libraries is available at http://www.dynameomics.org/fragments.

A probe to monitor performance of ¹⁵N longitudinal relaxation experiments for proteins in solution

The magnitude of the ¹⁵N longitudinal relaxation rate typically decreases as magnetic field strength increases in globular proteins in solution. Thus, it is important to test the performance of ¹⁵N longitudinal relaxation experiments at high field strength. Herein, a tool to investigate systematic errors in ¹⁵N longitudinal relaxation rate, R₁, is introduced. The tool, a difference in R₁ values between the two components of the ¹H-coupled ¹⁵N magnetizations, R₁(1)-R₁(2), conveniently detects inefficiencies in cancellation of cross correlation between ¹H-¹⁵N dipolar coupling and ¹⁵N chemical shift anisotropy. Experiments, in varying conditions, and simulations of a two-spin system indicate that insufficient cancellation of the cross correlation is due to (1) ¹H pulse imperfection and (2) ¹H off-resonance effect, and (3) is further amplified by residual ¹⁵N transverse magnetization that is caused by the ¹⁵N off-resonance effect. Results also show that this problem can be easily and practically remedied by discarding the initial decay points when recording ¹⁵N longitudinal relaxation in proteins.

Structural basis for the interaction of human β-defensin 6 and its putative chemokine receptor CCR2 and breast cancer microvesicles

Human β-defensins (hBDs) are believed to function as alarm molecules that stimulate the adaptive immune system when a threat is present. In addition to its antimicrobial activity, defensins present other activities such as chemoattraction of a range of different cell types to the sites of inflammation. We have solved the structure of the hBD6 by NMR spectroscopy that contains a conserved β-defensin domain followed by an extended C-terminus. We use NMR to monitor the interaction of hBD6 with microvesicles shed by breast cancer cell lines and with peptides derived from the extracellular domain of CC chemokine receptor 2 (Nt-CCR2) possessing or not possessing sulfation on Tyr26 and Tyr28. The NMR-derived model of the hBD6/CCR2 complex reveals a contiguous binding surface on hBD6, which comprises amino acid residues of the α-helix and β2-β3 loop. The microvesicle binding surface partially overlaps with the chemokine receptor interface. NMR spin relaxation suggests that free hBD6 and the hBD6/CCR2 complex exhibit microsecond-to-millisecond conformational dynamics encompassing the CCR2 binding site, which might facilitate selection of the molecular configuration optimal for binding. These data offer new insights into the structure-function relation of the hBD6-CCR2 interaction, which is a promising target for the design of novel anticancer agents.

Protein backbone and sidechain torsion angles predicted from NMR chemical shifts using artificial neural networks

A new program, TALOS-N, is introduced for predicting protein backbone torsion angles from NMR chemical shifts. The program relies far more extensively on the use of trained artificial neural networks than its predecessor, TALOS+. Validation on an independent set of proteins indicates that backbone torsion angles can be predicted for a larger, ≥90 % fraction of the residues, with an error rate smaller than ca 3.5 %, using an acceptance criterion that is nearly two-fold tighter than that used previously, and a root mean square difference between predicted and crystallographically observed (ϕ, ψ) torsion angles of ca 12º. TALOS-N also reports sidechain χ(1) rotameric states for about 50 % of the residues, and a consistency with reference structures of 89 %. The program includes a neural network trained to identify secondary structure from residue sequence and chemical shifts.

Backbone dynamics of an atypical orphan response regulator protein, Helicobacter pylori 1043

An atypical orphan response regulator protein, HP1043 (HP-RR) in Helicobacter pylori, is proven to be essential for cell growth and does not require the well known phosphorelay scheme. HP-RR was identified as a symmetric dimer with two functional domains, an N-terminal regulatory domain (HP-RR(r)) and a C-terminal effector domain (HP-RR(e)). HP-RR is a new class of response regulator, as a phosphorylation-independent regulator. Previously, we have presented a detailed three-dimensional structure of HP-RR using NMR spectroscopy and X-ray crystallography. In this study, in order to understand the functional importance of flexibilities in HP-RR(r) and HP-RR(e), T1, T2, heteronuclear NOE experiments have been performed and backbone dynamics of HP-RR(r) and HP-RR(e) were investigated. HP-RR(r) is a symmetric dimer and the interface region, α4-β5-α5 of dimer, showed high rigidity (high S (2) values). Site of rearrangements associated with phosphorylation of HP-RR(r) (Ser(75): R ex = 3.382, Ile(95): R ex = 5.228) showed slow chemical exchanges. HP-RR(e) is composed of three α-helices flanked on two sides by anti-parallel β-sheets. Low order parameters as well as conformational exchanges in the centers of loop regions known as the DNA binding site and transcription site of HP-RR(e) suggested that flexibility of HP-RR(e) is essential for interaction with DNA. In conclusion, backbone dynamics information for HP-RR implies that structural flexibilities in HP-RR(r) are necessary for the phosphorylation site and the dynamic nature of HP-RR(e) is essential for the regulation of interaction between protein and DNA.

Measurement of ¹⁵N relaxation rates in perdeuterated proteins by TROSY-based methods

While extracting dynamics parameters from backbone (15)N relaxation measurements in proteins has become routine over the past two decades, it is increasingly recognized that accurate quantitative analysis can remain limited by the potential presence of systematic errors associated with the measurement of (15)N R(1) and R(2) or R(1ρ) relaxation rates as well as heteronuclear (15)N-{(1)H} NOE values. We show that systematic errors in such measurements can be far larger than the statistical error derived from either the observed signal-to-noise ratio, or from the reproducibility of the measurement. Unless special precautions are taken, the problem of systematic errors is shown to be particularly acute in perdeuterated systems, and even more so when TROSY instead of HSQC elements are used to read out the (15)N magnetization through the NMR-sensitive (1)H nucleus. A discussion of the most common sources of systematic errors is presented, as well as TROSY-based pulse schemes that appear free of systematic errors to the level of <1 %. Application to the small perdeuterated protein GB3, which yields exceptionally high S/N and therefore is an ideal test molecule for detection of systematic errors, yields relaxation rates that show considerably less residue by residue variation than previous measurements. Measured R(2)'/R(1)' ratios fit an axially symmetric diffusion tensor with a Pearson's correlation coefficient of 0.97, comparable to fits obtained for backbone amide RDCs to the Saupe matrix.

Solution structure of the phytotoxic protein PcF: the first characterized member of the Phytophthora PcF toxin family

The PcF protein from Phytophthora cactorum is the first member of the "PcF toxin family" from the plant pathogens Phytophthora spp. It is able to induce withering in tomato and strawberry leaves. The lack of sequence similarity with other proteins hampers the identification of the molecular mechanisms responsible for its toxicity. Here, we show that the six cysteines form a disulphide pattern that is exclusive for PcF and essential for the protein withering activity. The NMR solution structure identifies a novel fold among protein effectors: a helix-loop-helix motif. The presence of a negatively charged surface suggests that it might act as a site of electrostatic interaction. Interestingly, a good fold match with Ole e 6, a plant protein with allergenic activity, highlighted the spatial superimposition of a stretch of identical residues. This finding suggests a possible biological activity based on molecular mimicry.

Phosphorylation-induced changes in backbone dynamics of the dematin headpiece C-terminal domain

Dematin is an actin-binding protein abundant in red blood cells and other tissues. It contains a villin-type 'headpiece' F-actin-binding domain at its extreme C-terminus. The isolated dematin headpiece domain (DHP) undergoes a significant conformational change upon phosphorylation. The mutation of Ser74 to Glu closely mimics the phosphorylation of DHP. We investigated motions in the backbone of DHP and its mutant DHPS74E using several complementary NMR relaxation techniques: laboratory frame (15)N NMR relaxation, which is sensitive primarily to the ps-ns time scale, cross-correlated chemical shift modulation NMR relaxation detecting correlated mus-ms time scale motions of neighboring (13)C' and (15)N nuclei, and cross-correlated relaxation of two (15)N-(1)H dipole-dipole interactions detecting slow motions of backbone NH vectors in successive amino acid residues. The results indicate a reduction in mobility upon the mutation in several regions of the protein. The additional salt bridge formed in DHPS74E that links the N- and C-terminal subdomains is likely to be responsible for these changes.

Protein backbone motions viewed by intraresidue and sequential HN-Halpha residual dipolar couplings

Triple resonance E.COSY-based techniques were used to measure intra-residue and sequential H(N)-H(alpha) residual dipolar couplings (RDCs) for the third IgG-binding domain of protein G (GB3), aligned in Pf1 medium. Measurements closely correlate with values predicted on the basis of an NMR structure, previously determined on the basis of a large number of one-bond backbone RDCs measured in five alignment media. However, in particular the sequential H(N)-H(alpha) RDCs are smaller than predicted for a static structure, suggesting a degree of motion for these internuclear vectors that exceeds that of the backbone amide N-H vectors. Of all experimentally determined GB3 structures available, the best correlation between experimental (1)H-(1)H couplings is observed for a GB3 ensemble, previously derived to generate a realistic picture of the conformational space sampled by GB3 (Clore and Schwieters, J Mol Biol 355:879-886, 2006). However, for both NMR and X-ray-derived structures the (1)H-(1)H couplings are found to be systematically smaller than expected on the basis of alignment tensors derived from (15)N-(1)H amide RDCs, assuming librationally corrected N-H bond lengths of 1.041 A.

Effects of denaturants and substitutions of hydrophobic residues on backbone dynamics of denatured staphylococcal nuclease

Analysis of residual dipolar couplings (RDCs) in the Delta131Delta fragment of staphylococcal nuclease has demonstrated that its ensemble-averaged structure is resistant to perturbations such as high concentrations of urea, low pH, and substitution of hydrophobic residues, suggesting that its residual structure is encoded by local side-chain/backbone interactions. In the present study, the effects of these same perturbations on the backbone dynamics of Delta131Delta were examined through (1)H-(15)N relaxation methods. Unlike the global structure reported by RDCs, the transverse relaxation rates R(2) were quite sensitive to denaturing conditions. At pH 5.2, Delta131Delta exhibits an uneven R(2) profile with several characteristic peaks involving hydrophobic chain segments. Protonation of carboxyl side chains by lowering the pH reduces the values of R(2) along the entire chain, yet these characteristic peaks remain. In contrast, high concentrations of urea or the substitution of 10 hydrophobic residues eliminates these peaks and reduces the R(2) values by a greater amount. The combination of low pH and high urea leads to further decreases in R(2). These denaturant-induced increases in backbone mobility are also reflected in decreases in (15)N NOEs and in relaxation interference parameters, with the former reporting an increase in fast motions and the latter a decrease in slow motions. Comparison between the changes in chain dynamics and the corresponding changes in Stokes radius and the patterns of RDCs suggests that regional variations in backbone dynamics in denatured nuclease arise primarily from local contacts between hydrophobic side chains and local interactions involving charged carboxyl groups.

Backbone dynamics of the human MIA protein studied by (15)N NMR relaxation: implications for extended interactions of SH3 domains

The melanoma inhibitory activity (MIA) protein is a clinically valuable marker in patients with malignant melanoma as enhanced values diagnose metastatic melanoma stages III and IV. Here, we report the backbone dynamics of human MIA studied by (15)N NMR relaxation experiments. The folded core of human MIA is found to be rigid, but several loops connecting beta-sheets, such as the RT-loop for example, display increased mobility on picosecond to nanosecond time scales. One of the most important dynamic features is the pronounced flexibility of the distal loop, comprising residues Asp 68 to Ala 75, where motions on time scales up to milliseconds occur. Further, significant exchange contributions are observed for residues of the canonical binding site of SH3 domains including the RT-loop, the n-Src loop, for the loop comprising residues 13 to 19, which we refer to as the"disulfide loop", in part for the distal loop, and the carboxyl terminus of human MIA. The functional importance of this dynamic behavior is discussed with respect to the biological activity of several point mutations of human MIA. The results of this study suggest that the MIA protein and the recently identified highly homologous fibrocyte-derived protein (FDP)/MIA-like (MIAL) constitute a new family of secreted proteins that adopt an SH3 domain-like fold in solution with expanded ligand interactions.

Characterization of binding-induced changes in dynamics suggests a model for sequence-nonspecific binding of ssDNA by replication protein A

Single-stranded-DNA-binding proteins (SSBs) are required for numerous genetic processes ranging from DNA synthesis to the repair of DNA damage, each of which requires binding with high affinity to ssDNA of variable base composition. To gain insight into the mechanism of sequence-nonspecific binding of ssDNA, NMR chemical shift and (15)N relaxation experiments were performed on an isolated ssDNA-binding domain (RPA70A) from the human SSB replication protein A. The backbone (13)C, (15)N, and (1)H resonances of RPA70A were assigned for the free protein and the d-CTTCA complex. The binding-induced changes in backbone chemical shifts were used to map out the ssDNA-binding site. Comparison to results obtained for the complex with d-C(5) showed that the basic mode of binding is independent of the ssDNA sequence, but that there are differences in the binding surfaces. Amide nitrogen relaxation rates (R(1) and R(2)) and (1)H-(15)N NOE values were measured for RPA70A in the absence and presence of d-CTTCA. Analysis of the data using the Model-Free formalism and spectral density mapping approaches showed that the structural changes in the binding site are accompanied by some significant changes in flexibility of the primary DNA-binding loops on multiple timescales. On the basis of these results and comparisons to related proteins, we propose that the mechanism of sequence-nonspecific binding of ssDNA involves dynamic remodeling of the binding surface.

Characterization of the structure and dynamics of amyloidogenic variants of human lysozyme by NMR spectroscopy

The structures and dynamics of the native states of two mutational variants of human lysozyme, I56T and D67H, both associated with non-neuropathic systemic amyloidosis, have been investigated by NMR spectroscopy. The (1)H and (15)N main-chain amide chemical shifts of the I56T variant are very similar to those of the wild-type protein, but those of the D67H variant are greatly altered for 28 residues in the beta-domain. This finding is consistent with the X-ray crystallographic analysis, which shows that the structure of this variant is significantly altered from that of the wild-type protein in this region. The (1)H-(15)N heteronuclear NOE values show that, with the exception of V121, every residue in the wild-type and I56T proteins is located in tightly packed structures characteristic of the native states of most proteins. In contrast, D67H has a region of substantially increased mobility as shown by a dramatic decrease in heteronuclear NOE values of residues near the site of mutation. Despite this unusual flexibility, the D67H variant has no greater propensity to form amyloid fibrils in vivo or in vitro than has I56T. This finding indicates that it is the increased ability of the variants to access partially folded conformations, rather than intrinsic changes in their native state properties, that is the origin of their amyloidogenicity.

Key interactions in the immunoglobulin-like structure of apo-neocarzinostatin: evidence from nuclear magnetic resonance relaxation data and molecular dynamics simulations

The three-dimensional structure of apo-neocarzinostatin (apo-NCS, MW: ca.11000, antitumoral chromophore carrier protein) is based on a seven-stranded antiparallel beta-sandwich, very similar to the immunoglobulin folding domain. We investigated the backbone dynamics of apo-NCS by (13)C-NMR relaxation measurements and molecular dynamics simulation. Model-free parameters determined from the experimental data are compared with a 1.5-nsec molecular simulation of apo-NCS in aqueous solution. This comparison provides an accurate description of both local and collective movements within the protein. This analysis enabled us to correlate dynamic processes with key interactions of this beta-protein. Local motions that could be relevant for the intermolecular association with the ligand are also described.

Structure of the recombinant full-length hamster prion protein PrP(29-231): the N terminus is highly flexible

The prion diseases seem to be caused by a conformational change of the prion protein (PrP) from the benign cellular form PrPC to the infectious scrapie form PrPSc; thus, detailed information about PrP structure may provide essential insights into the mechanism by which these diseases develop. In this study, the secondary structure of the recombinant Syrian hamster PrP of residues 29-231 [PrP(29-231)] is investigated by multidimensional heteronuclear NMR. Chemical shift index analysis and nuclear Overhauser effect data show that PrP(29-231) contains three helices and possibly one short beta-strand. Most striking is the random-coil nature of chemical shifts for residues 30-124 in the full-length PrP. Although the secondary structure elements are similar to those found in mouse PrP fragment PrP(121-231), the secondary structure boundaries of PrP(29-231) are different from those in mouse PrP(121-231) but similar to those found in the structure of Syrian hamster PrP(90-231). Comparison of resonance assignments of PrP(29-231) and PrP(90-231) indicates that there may be transient interactions between the additional residues and the structured core. Backbone dynamics studies done by using the heteronuclear [1H]-15N nuclear Overhauser effect indicate that almost half of PrP(29-231), residues 29-124, is highly flexible. This plastic region could feature in the conversion of PrPC to PrPSc by template-assisted formation of beta-structure.