Proline Peptide Bond Isomerization in Ubiquitin Under Folding and Denaturing Conditions by Pressure-Jump NMR.

Proline isomerization is widely recognized as a kinetic bottleneck in protein folding, amplified for proteins rich in Pro residues. We introduced repeated hydrostatic pressure jumps between native and pressure-denaturing conditions inside an NMR sample cell to study proline isomerization in the pressure-sensitized L50A ubiquitin mutant. Whereas in two unfolded heptapeptides, X-Pro peptide bonds isomerized ca 1.6-fold faster at 1 bar than at 2.5 kbar, for ubiquitin ca eight-fold faster isomerization was observed for Pro-38 and ca two-fold for Pro-19 and Pro-37 relative to rates measured in the pressure-denatured state. Activation energies for isomerization in pressure-denatured ubiquitin were close to literature values of 20 kcal/mole for denatured polypeptides but showed a substantial drop to 12.7 kcal/mole for Pro-38 at atmospheric pressure. For ubiquitin isomers with a cis E18-P19 peptide bond, the 1-bar NMR spectrum showed sharp resonances with near random coil chemical shifts for the C-terminal half of the protein, characteristic of an unfolded chain, while most of the N-terminal residues were invisible due to exchange broadening, pointing to a metastable partially folded state for this previously recognized 'folding nucleus'. For cis-P37 isomers, a drop in pressure resulted in the rapid loss of nearly all unfolded-state NMR resonances, while the recovery of native state intensity revealed a slow component attributed to cis â†’ trans isomerization of P37. This result implies that the NMR-invisible cis-P37 isomer adopts a molten globule state that encompasses the entire length of the ubiquitin chain, suggestive of a structure that mostly resembles the folded state.

Experimental NOE, Chemical Shift, and Proline Isomerization Data Provide Detailed Insights into Amelotin Oligomerization.

Amelotin is an intrinsically disordered protein (IDP) rich in Pro residues and is involved in hydroxyapatite mineralization. It rapidly oligomerizes under physiological conditions of pH and pressure but reverts to its monomeric IDP state at elevated pressure. We identified a 105-residue segment of the protein that becomes ordered upon oligomerization, and we used pressure-jump NMR spectroscopy to measure long-range NOE contacts that exist exclusively in the oligomeric NMR-invisible state. The kinetics of oligomerization and dissociation were probed at the residue-specific level, revealing that the oligomerization process is initiated in the C-terminal half of the segment. Using pressure-jump NMR, the degree of order in the oligomer at the sites of Pro residues was probed by monitoring changes in cis/trans equilibria relative to the IDP state after long-term equilibration under oligomerizing conditions. Whereas most Pro residues revert to trans in the oligomeric state, Pro-49 favors a cis configuration and three Pro residues retain an unchanged cis fraction, pointing to their local lack of order in the oligomeric state. NOE contacts and secondary 13C chemical shifts in the oligomeric state indicate the presence of an 11-residue α-helix, preceded by a small intramolecular antiparallel ß-sheet, with slower formation of long-range intermolecular interactions to N-terminal residues. Although none of the models generated by AlphaFold2 for the amelotin monomer was consistent with experimental data, subunits of a hexamer generated by AlphaFold-Multimer satisfied intramolecular NOE and chemical shift data and may provide a starting point for developing atomic models for the oligomeric state.

A methyl-TROSY based 13C relaxation dispersion NMR experiment for studies of chemical exchange in proteins.

A methyl Transverse Relaxation Optimized Spectroscopy (methyl-TROSY) based, multiple quantum (MQ) 13C Carr-Purcell-Meiboom-Gill (CPMG) relaxation dispersion NMR experiment is described. The experiment is derived from the previously developed MQ 13C-1H CPMG scheme (Korzhnev in J Am Chem Soc 126: 3964-73, 2004) supplemented with a CPMG train of refocusing 1H pulses applied with constant frequency and synchronized with the 13C CPMG pulse train. The optimal 1H 'decoupling' scheme that minimizes the amount of fast-relaxing methyl MQ magnetization present during CPMG intervals, makes use of an XY-4 phase cycling of the refocusing composite 1H pulses. For small-to-medium sized proteins, the MQ 13C CPMG experiment has the advantage over its single quantum (SQ) 13C counterpart of significantly reducing intrinsic, exchange-free relaxation rates of methyl coherences. For high molecular weight proteins, the MQ 13C CPMG experiment eliminates complications in the interpretation of MQ 13C-1H CPMG relaxation dispersion profiles arising from contributions to exchange from differences in methyl 1H chemical shifts between ground and excited states. The MQ 13C CPMG experiment is tested on two protein systems: (1) a triple mutant of the Fyn SH3 domain that interconverts slowly on the chemical shift time scale between the major folded state and an excited state folding intermediate; and (2) the 82-kDa enzyme Malate Synthase G (MSG), where chemical exchange at individual Ile δ1 methyl positions occurs on a much faster time-scale.

Co-Evolutionary Fitness Landscapes for Sequence Design.

Efficient and accurate models to predict the fitness of a sequence would be extremely valuable in protein design. We have explored the use of statistical potentials for the coevolutionary fitness landscape, extracted from known protein sequences, in conjunction with Monte Carlo simulations, as a tool for design. As proof of principle, we created a series of predicted high-fitness sequences for three different protein folds, representative of different structural classes: the GA (all-α) and GB (α/ß) binding domains of streptococcal protein G, and an SH3 (all-ß) domain. We found that most of the designed proteins can fold stably to the target structure, and a structure for a representative of each for GA, GB and SH3 was determined. Several of our designed proteins were also able to bind to native ligands, in some cases with higher affinity than wild-type. Thus, a search using a statistical fitness landscape is a remarkably effective tool for finding novel stable protein sequences.

Tilted, Uninterrupted, Monomeric HIV-1 gp41 Transmembrane Helix from Residual Dipolar Couplings.

Cryo-electron microscopy and X-ray crystallography have shown that the pre- and postfusion states of the HIV-1 gp41 viral coat protein, although very different from one another, each adopt C3 symmetric structures. A stable homotrimeric structure for the transmembrane domain (TM) also was modeled and supported by experimental data. For a C3 symmetric structure, alignment in an anisotropic medium must be axially symmetric, with the unique axis of the alignment tensor coinciding with the C3 axis. However, NMR residual dipolar couplings (RDCs) measured under three different alignment conditions were found to be incompatible with C3 symmetry. Subsequent measurements by paramagnetic relaxation enhancement, analytical ultracentrifugation, and DEER EPR, indicate that the transmembrane domain is monomeric. 15N NMR relaxation data and RDCs show that TM is highly ordered and uninterrupted for a total length of 32 residues, extending well into the membrane proximal external region.

Conformational Heterogeneity in the Activation Mechanism of Bax.

Bax is known for its pro-apoptotic role within the mitochondrial pathway of apoptosis. However, the mechanism for transitioning Bax from cytosolic to membrane-bound oligomer remains elusive. Previous nuclear magnetic resonance (NMR) and electron paramagnetic resonance (EPR) studies defined monomeric Bax as conformationally homogeneous. Yet it has recently been proposed that monomeric Bax exists in equilibrium with a minor state that is distinctly different from its NMR structure. Here, we revisited the structural analysis of Bax using methods uniquely suited for unveiling "invisible" states of proteins, namely, NMR paramagnetic relaxation enhancements and EPR double electron-electron resonance (DEER). Additionally we examined the effect of glycerol, the co-solvent of choice in DEER studies, on the structure of Bax using NMR chemical-shift perturbations and residual dipolar couplings. Based on our combined NMR and EPR results, Bax is a conformationally homogeneous protein prior to its activation.

Long Distance Measurements up to 160†Šin the GroEL Tetradecamer Using Q-Band DEER EPR Spectroscopy.

Current distance measurements between spin-labels on multimeric protonated proteins using double electron-electron resonance (DEER) EPR spectroscopy are generally limited to the 15-60 Šrange. Here we show how DEER experiments can be extended to dipolar evolution times of ca. 80 µs, permitting distances up to 170 Što be accessed in multimeric proteins. The method relies on sparse spin-labeling, supplemented by deuteration of protein and solvent, to minimize the deleterious impact of multispin effects and substantially increase the apparent spin-label phase memory relaxation time, complemented by high sensitivity afforded by measurements at Q-band. We demonstrate the approach using the tetradecameric molecular machine GroEL as an example. Two engineered surface-exposed mutants, R268C and E315C, are used to measure pairwise distance distributions with mean values ranging from 20 to 100 Šand from 30 to 160 Å, respectively, both within and between the two heptameric rings of GroEL. The measured distance distributions are consistent with the known crystal structure of apo GroEL. The methodology presented here should significantly expand the use of DEER for the structural characterization of conformational changes in higher order oligomers.

Insights into the Conformation of the Membrane Proximal Regions Critical to the Trimerization of the HIV-1 gp41 Ectodomain Bound to Dodecyl Phosphocholine Micelles.

The transitioning of the ectodomain of gp41 from a pre-hairpin to a six-helix bundle conformation is a crucial aspect of virus-cell fusion. To gain insight into the intermediary steps of the fusion process we have studied the pH and dodecyl phosphocholine (DPC) micelle dependent trimer association of gp41 by systematic deletion analysis of an optimized construct termed 17-172 (residues 528 to 683 of Env) that spans the fusion peptide proximal region (FPPR) to the membrane proximal external region (MPER) of gp41, by sedimentation velocity and double electron-electron resonance (DEER) EPR spectroscopy. Trimerization at pH 7 requires the presence of both the FPPR and MPER regions. However, at pH 4, the protein completely dissociates to monomers. DEER measurements reveal a partial fraying of the C-terminal MPER residues in the 17-172 trimer while the other regions, including the FPPR, remain compact. In accordance, truncating nine C-terminal MPER residues (675-683) in the 17-172 construct does not shift the trimer-monomer equilibrium significantly. Thus, in the context of the gp41 ectodomain spanning residues 17-172, trimerization is clearly dependent on FPPR and MPER regions even when the terminal residues of MPER unravel. The antibody Z13e1, which spans both the 2F5 and 4E10 epitopes in MPER, binds to 17-172 with a Kd of 1 ± 0.12 µM. Accordingly, individual antibodies 2F5 and 4E10 also recognize the 17-172 trimer/DPC complex. We propose that binding of the C-terminal residues of MPER to the surface of the DPC micelles models a correct positioning of the trimeric transmembrane domain anchored in the viral membrane.

The C34 Peptide Fusion Inhibitor Binds to the Six-Helix Bundle Core Domain of HIV-1 gp41 by Displacement of the C-Terminal Helical Repeat Region.

The conformational transition of the core domain of HIV-1 gp41 from a prehairpin intermediate to a six-helix bundle is responsible for virus-cell fusion. Several inhibitors which target the N-heptad repeat helical coiled-coil trimer that is fully accessible in the prehairpin intermediate have been designed. One such inhibitor is the peptide C34 derived from the C-heptad repeat of gp41 that forms the exterior of the six-helix bundle. Here, using a variety of biophysical techniques, including dye tagging, size-exclusion chromatography combined with multiangle light scattering, double electron-electron resonance EPR spectroscopy, and circular dichroism, we investigate the binding of C34 to two six-helix bundle mimetics comprising N- and C-heptad repeats either without (core(SP)) or with (core(S)) a short spacer connecting the two. In the case of core(SP), C34 directly exchanges with the C-heptad repeat. For core(S), up to two molecules of C34 bind the six-helix bundle via displacement of the C-heptad repeat. These results suggest that fusion inhibitors such as C34 can target a continuum of transitioning conformational states from the prehairpin intermediate to the six-helix bundle prior to the occurrence of irreversible fusion of viral and target cell membranes.

Dependence of distance distributions derived from double electron-electron resonance pulsed EPR spectroscopy on pulse-sequence time.

Pulsed double electron-electron resonance (DEER) provides pairwise P(r) distance distributions in doubly spin labeled proteins. We report that in protonated proteins, P(r) is dependent on the length of the second echo period T owing to local environmental effects on the spin-label phase memory relaxation time Tm . For the protein ABD, this effect results in a 1.4 Šincrease in the P(r) maximum from T=6 to 20 µs. Protein A has a bimodal P(r) distribution, and the relative height of the shorter distance peak at T=10 µs, the shortest value required to obtain a reliable P(r), is reduced by 40 % relative to that found by extrapolation to T=0. Our results indicate that data at a series of T values are essential for quantitative interpretation of DEER to determine the extent of the T dependence and to extrapolate the results to T=0. Complete deuteration (99 %) of the protein was accompanied by a significant increase in Tm and effectively abolished the P(r) dependence on T.

Structure and dynamics of full-length HIV-1 capsid protein in solution.

The HIV-1 capsid protein plays a crucial role in viral infectivity, assembling into a cone that encloses the viral RNA. In the mature virion, the N-terminal domain of the capsid protein forms hexameric and pentameric rings, while C-terminal domain homodimers connect adjacent N-terminal domain rings to one another. Structures of disulfide-linked hexamer and pentamer assemblies, as well as structures of the isolated domains, have been solved previously. The dimer configuration in C-terminal domain constructs differs in solution (residues 144-231) and crystal (residues 146-231) structures by ∼30°, and it has been postulated that the former connects the hexamers while the latter links pentamers to hexamers. Here we study the structure and dynamics of full-length capsid protein in solution, comprising a mixture of monomeric and dimeric forms in dynamic equilibrium, using ensemble simulated annealing driven by experimental NMR residual dipolar couplings and X-ray scattering data. The complexity of the system necessitated the development of a novel computational framework that should be generally applicable to many other challenging systems that currently escape structural characterization by standard application of mainstream techniques of structural biology. We show that the orientation of the C-terminal domains in dimeric full-length capsid and isolated C-terminal domain constructs is the same in solution, and we obtain a quantitative description of the conformational space sampled by the N-terminal domain relative to the C-terminal domain on the nano- to millisecond time scale. The positional distribution of the N-terminal domain relative to the C-terminal domain is large and modulated by the oligomerization state of the C-terminal domain. We also show that a model of the hexamer/pentamer assembly can be readily generated with a single configuration of the C-terminal domain dimer, and that capsid assembly likely proceeds via conformational selection of sparsely populated configurations of the N-terminal domain within the capsid protein dimer.

Measurement of 15N relaxation in the detergent-solubilized tetrameric KcsA potassium channel.

A set of TROSY-HNCO (tHNCO)-based 3D experiments is presented for measuring (15)N relaxation parameters in large, membrane-associated proteins, characterized by slow tumbling times and significant spectral overlap. Measurement of backbone (15)N R (1), R (1rho), (15)N-{(1)H} NOE, and (15)N CSA/dipolar cross correlation is demonstrated and applied to study the dynamic behavior of the homotetrameric KcsA potassium channel in SDS micelles under conditions where this channel is in the closed state. The micelle-encapsulated transmembrane domain, KcsA(TM), exhibits a high degree of order, tumbling as an oblate ellipsoid with a global rotational correlation time, tau(c) = 38 +/- 2.5 ns, at 50 degrees C and a diffusion anisotropy, Dparallel/Dperpendicular = 0.79+/-0.05, corresponding to an aspect ratio a/b >/= 1.4. The N- and C-terminal intracellular segments of KcsA exhibit considerable internal dynamics (S (2) values in the 0.2-0.45 range), but are distinctly more ordered than what has been observed for unstructured random coils. Relaxation behavior in these domains confirms the position of the C-terminal helix, and indicates that in SDS micelles, this amphiphilic helix does not associate into a stable homotetrameric helical bundle. The relaxation data indicate the absence of elevated backbone dynamics on the ps-ns time scale for the 5-residue selectivity filter, which selects K(+) ions to enter the channel.

Characterization of phospholipid mixed micelles by translational diffusion.

The concentration dependence of the translational self diffusion rate, D (s), has been measured for a range of micelle and mixed micelle systems. Use of bipolar gradient pulse pairs in the longitudinal eddy current delay experiment minimizes NOE attenuation and is found critical for optimizing sensitivity of the translational diffusion measurement of macromolecules and aggregates. For low volume fractions Phi (Phi\\ le 15% v/v) of the micelles, experimental measurement of the concentration dependence, combined with use of the D (s)= D (o)(1-3.2lambdaPhi) relationship, yields the hydrodynamic volume. For proteins, the hydrodynamic volume, derived from D (s) at infinitely dilute concentration, is found to be about 2.6 times the unhydrated molecular volume. Using the data collected for hen egg white lysozyme as a reference, diffusion data for dihexanoyl phosphatidylcholine (DHPC) micelles indicate approximately 27 molecules per micelle, and a critical micelle concentration of 14 mM. Differences in translational diffusion rates for detergent and long chain phospholipids in mixed micelles are attributed to rapid exchange between free and micelle-bound detergent. This difference permits determination of the free detergent concentration, which, for a high detergent to long chain phospholipid molar ratio, is found to depend strongly on this ratio. The hydrodynamic volume of DHPC/POPC bicelles, loaded with an M2 channel peptide homolog, derived from translational diffusion, predicts a rotational correlation time that slightly exceeds the value obtained from peptide (15)N relaxation data.

Molecular basis of sequence-specific single-stranded DNA recognition by KH domains: solution structure of a complex between hnRNP K KH3 and single-stranded DNA.

To elucidate the basis of sequence-specific single-stranded (ss) DNA recognition by K homology (KH) domains, we have solved the solution structure of a complex between the KH3 domain of the transcriptional regulator heterogeneous nuclear ribonucleoprotein K (hnRNP K) and a 10mer ssDNA. We show that hnRNP K KH3 specifically recognizes a tetrad of sequence 5'd-TCCC. The complex is stabilized by a dense network of methyl-oxygen hydrogen bonds involving the methyl groups of three isoleucine residues and the O2 and N3 atoms of the two central cytosine bases. Comparison with the recently solved structure of a specific protein-ssDNA complex involving the KH3 and KH4 domains of the far upstream element (FUSE) binding protein FBP suggests that the amino acid located five residues N-terminal of the invariant GXXG motif, which is characteristic of all KH domains, plays a crucial role in discrimination of the first two bases of the tetrad.

Structure and dynamics of KH domains from FBP bound to single-stranded DNA.

Gene regulation can be tightly controlled by recognition of DNA deformations that are induced by stress generated during transcription. The KH domains of the FUSE-binding protein (FBP), a regulator of c-myc expression, bind in vivo and in vitro to the single-stranded far-upstream element (FUSE), 1,500 base pairs upstream from the c-myc promoter. FBP bound to FUSE acts through TFIIH at the promoter. Here we report the solution structure of a complex between the KH3 and KH4 domains of FBP and a 29-base single-stranded DNA from FUSE. The KH domains recognize two sites, 9-10 bases in length, separated by 5 bases, with KH4 bound to the 5' site and KH3 to the 3' site. The central portion of each site comprises a tetrad of sequence 5'd-ATTC for KH4 and 5'd-TTTT for KH3. Dynamics measurements show that the two KH domains bind as articulated modules to single-stranded DNA, providing a flexible framework with which to recognize transient, moving targets.