Investigation of dipolar-mediated water-protein interactions in microcrystalline Crh by solid-state NMR spectroscopy

Water-protein interactions play a major role in protein folding, structure, and function, and solid-state NMR has recently been shown to be a powerful tool for the site-resolved observation of these interactions in solid proteins. In this article we report investigations on possible water-protein dipolar transfer mechanisms in the microcrystalline deuterated protein Crh by a set of solid-state NMR techniques. Double-quantum (DQ) filtered and edited heteronuclear correlation experiments are used to follow direct dipolar water-protein magnetization transfers. Experimental data reveal no evidence for "solid-like" water molecules, indicating that residence times of solvent molecules are shorter than required for DQ creation, typically a few hundred microseconds. An alternative magnetization pathway, intermolecular cross-relaxation via heteronuclear nuclear Overhauser effects (NOEs), is probed by saturation transfer experiments. The significant additional enhancements observed when irradiating at the water frequency can possibly be attributed to direct heteronuclear water-protein NOEs; however, a contribution from relayed magnetization transfer via chemical exchange or proton-proton dipolar mechanisms cannot be excluded.     

3D structure determination of the Crh protein from highly ambiguous solid-state NMR restraints

In a wide variety of proteins, insolubility presents a challenge to structural biology, as X-ray crystallography and liquid-state NMR are unsuitable. Indeed, no general approach is available as of today for studying the three-dimensional structures of membrane proteins and protein fibrils. We here demonstrate, at the example of the microcrystalline model protein Crh, how high-resolution 3D structures can be derived from magic-angle spinning solid-state NMR distance restraints for fully labeled protein samples. First, we show that proton-mediated rare-spin correlation spectra, as well as carbon-13 spin diffusion experiments, provide enough short, medium, and long-range structural restraints to obtain high-resolution structures of this 2 x 10.4 kDa dimeric protein. Nevertheless, the large number of 13C/15N spins present in this protein, combined with solid-state NMR line widths of about 0.5-1 ppm, induces substantial ambiguities in resonance assignments, preventing 3D structure determination by using distance restraints uniquely assigned on the basis of their chemical shifts. In the second part, we thus demonstrate that an automated iterative assignment algorithm implemented in a dedicated solid-state NMR version of the program ARIA permits to resolve the majority of ambiguities and to calculate a de novo 3D structure from highly ambiguous solid-state NMR data, using a unique fully labeled protein sample. We present, using distance restraints obtained through the iterative assignment process, as well as dihedral angle restraints predicted from chemical shifts, the 3D structure of the fully labeled Crh dimer refined at a root-mean-square deviation of 1.33 A.     

Protein 3D structure determination by high-resolution solid-state NMR

Developments in NMR technology, sample preparation, pulse sequence methodology and structure calculation protocols have recently allowed one to progress towards structure determination at high-resolution of proteins by solid-state NMR spectroscopy. We here report solid-state NMR protocols based on magic-angle-spinning experiments, combined with modified structure calculation protocols, for structure determination of uniformly ¹³C, ¹⁵N isotopically labeled proteins. We demonstrate the use of these protocols to obtain high-resolution structures for the example of the microcrystalline Crh protein. The CHHC, DARR and PAR solid-state NMR experiments, as well as the calculation protocols using the program ARIA, are presented.     

Denatured-state ensemble and the early-stage folding of the G29A mutant of the B-domain of protein A

The folding mechanism of the G29A mutant of the B-domain of protein A (BdpA) has been studied by all-atom molecular dynamics simulation using AMBER force field (ff03) and generalized Born continuum solvent model. Started from the extended chain conformation, a total of 16 simulations (400 ns each) at 300 K captured some early folding events of the G29A mutant of BdpA. In one of the 16 trajectories, the G29A mutant folded within 2.8 A (root mean square) of the wild-type NMR structure. We observed that the fast burial of hydrophobic residues was the driving force to bring the distant residues into close proximity. The initiation of the helix I and III occurred during the stage of hydrophobic collapse. The initiation and growth of the helix II was slow. Both the secondary structure formation and the development of the native tertiary contacts suggested a multistage folding process. Clustering analysis indicated that two helix species (helices I and III) could be intermediates. Further analysis revealed that the hydrophobic residues of partially folded helix II formed nativelike hydrophobic contacts with helices I and III that stabilized a nativelike state and delayed the completion of folding of the entire protein. The details of the early folding process were compared with other theoretical and experimental studies. It was found that a nativelike hydrophobic cluster was formed by residues including F(30), I(31), L(34), L(44), L(45), and A(48) that prevented further development of the native structures, and breaking the hydrophobic cluster like this one contributed to the rate-limiting step. This was in complete agreement with the recent kinetic measurements in which mutations of these residues to Gly and Ala substantially increased the folding rates by as much as 60 times. Apparently, destabilization of nonnative states dramatically enhanced the folding rates.     

On the solid state structure of 4-iodobenzoic acid

The solid-state structure of 4-iodobenzoic acid has been confirmed by variable temperature X-ray diffraction, variable temperature solid-state NMR and differential scanning calorimetry. 4-iodobenzoic acid crystallizes in the space group P2(1)/n, and dimerizes in the solid state about a center of inversion. Using extensive X-ray crystallographic data collections, the placement of the carboxylate H atoms from the residual electron density in difference Fourier maps was determined. The position of the electron density associated with the proton is found to vary with temperature in that the population of the disordered sites changes with varying temperature. Determination of the crystal structure between the temperatures of 248 and 198 K was not possible due to a phase transition, an endothermic event occurring at 230.77 K. The phase transition is also indicated by a change in the relaxation time of the ring carbon atoms in the solid-state NMR data. Though the dominating force in the dimeric unit in the solid state is the presence of strong hydrogen bonds, there are also van der Waals forces present between the iodine atoms. In the layered structure, the iodine-iodine distance is within the van der Waals contact radii, an interaction which causes a deformation in the electron density of the iodine atoms.     

14-helical folding in a cyclobutane-containing beta-tetrapeptide

The efficient synthesis of tetrapeptide 5 containing, in alternation, cyclobutane and beta-alanine residues is described. NMR experiments both at low temperature in CDCl(3) and at 298 K in DMSO-d(6) solutions show the contribution of a strong hydrogen bond in the folded major conformation of 5. Temperature coefficients and diffusion times point out a hydrogen bond involving the NH proton from the cyclobutane residue 1 whereas NOEs manifest the high rigidity of the central fragment of the molecule and are compatible with a 14-membered macrocycle. Theoretical calculations predict a most stable folded conformation corresponding to a 14-helix stabilized by a hydrogen bond between NH(10) in the first residue and OC(25) in the third residue. This structure remains unaltered during the molecular dynamics simulation at 298 K in chloroform. All these results provide evidence for a 14-helical folding and reveal the ability of cis-2-aminocyclobutane carboxylic acid residues to promote folded conformations when incorporated into beta-peptides.     

NMR spectra of a microcrystalline protein at 30 kHz MAS

Proteins are not always available in amounts desirable for solid-state magic-angle spinning (MAS) nuclear-magnetic resonance (NMR) spectroscopy. To maximize the signal-to-noise ratio achievable with small samples, the filling factor must be optimized by using small-diameter MAS rotors. These rotors have the added benefit of allowing higher radio frequency field amplitudes during polarization transfer steps and during decoupling periods as well as allowing higher spinning frequencies. We demonstrate the advantages of relatively fast MAS (30 kHz using a 2.5 mm rotor) compared to MAS at 12 kHz for the 10.4 kDa model protein Crh with 93 residues and show that the signal-to-noise ratio in two-dimensional correlation spectra can be significantly improved by taking advantage of optimized pulse sequences available with rapid MAS.     

Modulation of structure and dynamics by disulfide bond formation in unfolded states

During oxidative folding, the formation of disulfide bonds has profound effects on guiding the protein folding pathway. Until now, comparatively little is known about the changes in the conformational dynamics in folding intermediates of proteins that contain only a subset of their native disulfide bonds. In this comprehensive study, we probe the conformational landscape of non-native states of lysozyme containing a single native disulfide bond utilizing nuclear magnetic resonance (NMR) spectroscopy, small-angle X-ray scattering (SAXS), circular dichroism (CD) data, and modeling approaches. The impact on conformational dynamics varies widely depending on the loop size of the single disulfide variants and deviates significantly from random coil predictions for both NMR and SAXS data. From these experiments, we conclude that the introduction of single disulfides spanning a large portion of the polypeptide chain shifts the structure and dynamics of hydrophobic core residues of the protein so that these regions exhibit levels of order comparable to the native state on the nanosecond time scale.     

The interaction of lipid modified pseudopeptides with lipid membranes

We have studied the structure of two lipopeptides based on the simple dipeptide building block L-Phe-D-Oxd. These peptides have been reported previously to form fiber-like materials. The lipopeptides synthesized here had the structures C(n)(2)H((2n+1))CO-L-Phe-D-Oxd-OBn or C(n)(2)H((2n+1))CO-D-Phe-L-Oxd-OBn with n = 5 or 11. Addition of the N-terminal lipid modification did not cause a major disturbance of the structures these molecules form. The lipid modifications themselves showed highly rigid structures as inferred from solid-state (2)H NMR. The peptide backbone showed (13)C NMR chemical shifts in agreement with β-sheet secondary structure. Addition of a lipid modification to the N-terminus is a common motif in biology to attach proteins to the membrane. Therefore, we also investigated the lipopeptides in the presence of synthetic POPC bilayers. Two different molecular species were detected under these circumstances: (i) lipopeptide monomers that showed chain order parameters similar to those of the host membrane, (ii) lipopeptide aggregates that exhibited very similar structures and dynamics as the crystalline aggregates. Overall, the lipopeptides showed a well defined and rigid secondary structure that is in agreement with fibrillar aggregates previously detected for those peptides without the lipid modification.     

A conventional an 2DCOS infrared approach to the kinetics of protein misfolding

Cell viability depends on the correct folding of the proteins involved in metabolism. Proteins are synthesized on the endoplasmic reticulum and must follow a pathway to a correct, metastable, tridimensional structure. Changes in structure or in environmental conditions can drive an instability of the folding conditions and produce non-active aggregates that in principle are proteolysed by the cellular mechanisms. However, these aggregates can be even more stable than the native proteins, escaping the cellular control. They can be classified as amorphous, if there is not a well-organized structural pattern, or ordered if a repetitive pattern is produced. These ordered structures, known as fibrils, are involved in many diseases. Infrared spectroscopy is a method of choice to study its formation because it is not affected by turbidity or the formation of high molecular weight aggregates. Moreover, in both cases, two bands characteristic of intermolecular β-sheets allow the monitoring of the aggregate formation. In both cases, the appearance of these bands involves a non-reversible path in protein folding. It has been suggested that a difference in the ordered structures involves an increasing in band intensity. This change can be the origin in variations on the 2DCOS maps. The synchronous map gives an overall idea of the process involved. The asynchronous is more informative because reflects the kinetic changes produced. The outcome of both processes, amorphous or ordered is that 2DCOS can provide a further insight to the knowledge of the kinetic processes giving rise to aggregated structures. This outcome could consist on the order in which the different secondary structures are prone to form the aggregates.     

Mpro-C蛋白三维结构域交换的机理: 来自分子模拟的线索

SARS冠状病毒主蛋白酶(M^pro)在病毒的蛋白酶切过程中发挥着重要作用. M^pro的晶体结构显示它存在两种形式的二聚体: 一种是发生三维结构域交换的形式, 另一种是非交换的形式. M^pro的C端结构域(M^pro-C)单独表达时也能形成与全长M^pro类似的三维结构域交换二聚体. 三维结构域交换通常发生在蛋白质的表面, 但M^pro-C 的结构域交换却发生在疏水核心. 在本文中, 我们利用分子动力学模拟及三维结构域交换预测算法研究了M^pro-C 中被高度埋藏的核心螺旋片段发生交换的机理. 我们发现基于结构与基于序列的已有算法都不能正确预言出M^pro-C和M^pro中发生结构域交换的铰链区位置. 分子模拟结果表明M^pro-C中的交换片段在天然态下埋藏得很好, 但在变性单体中则会被释放并暴露在外面. 因此, 在完全或部分解折叠状态下交换片段的打开有助于促进单体间的相互作用及结构域交换二聚体的形成.     

NMR spectroscopic characterization of millisecond protein folding by transverse relaxation dispersion measurements

The cold shock protein CspB adopts its native and functional tertiary structure on the millisecond time scale. We employed transverse relaxation NMR methods, which allow a quantitative measurement of the cooperativity of this fast folding reaction on a residue basis. Thereby, chemical exchange contributions to the transverse relaxation rate (R(2)) were observed for every residue of CspB verifying the potential of this method to identify not only local dynamics but also to characterize global events. Toward this end, the homogeneity of the transition state of folding was probed by comparing Chevron plots (i.e., dependence of the apparent folding rate on the denaturant concentration) determined by stopped-flow fluorescence with Chevron plots of six residues acquired by R(2) dispersion experiments. The coinciding results obtained for probes at different locations in the three-dimensional structure of CspB indicate the ability and significance of transverse relaxation NMR to determine Chevron plots on a residue-by-residue basis providing detailed insights on the nature of the transition state of folding.     

Anesthetic modulation of protein dynamics: insight from an NMR study

Mistic (membrane integrating sequence for translation of integral membrane protein constructs) comprises the four-alpha-helix bundle scaffold found in the transmembrane domains of the Cys-loop receptors that are plausible targets for general anesthetics. Nuclear magnetic resonance (NMR) studies of anesthetic halothane interaction with Mistic in dodecyl phosphocholine (DPC) micelles provide an experimental basis for understanding molecular mechanisms of general anesthesia. Halothane was found to interact directly with Mistic, mostly in the interfacial loop regions. Although the presence of halothane had little effect on Mistic structure, (15)N NMR relaxation dispersion measurements revealed that halothane affected Mistic's motion on the microsecond-millisecond time scale. Halothane shifted the equilibrium of chemical exchange in some residues and made the exchange faster or slower in comparison to the original state in the absence of halothane. The motion on the microsecond-millisecond time scale in several residues disappeared in response to the addition of halothane. Most of the residues experiencing halothane-induced dynamics changes also exhibited profound halothane-induced changes in chemical shift, suggesting that dynamics modification of these residues might result from their direct interaction with halothane molecules. Allosteric modulation by halothane also contributed to dynamics changes, as reflected in residues I52 and Y82 where halothane introduction brought about dynamics changes but not chemical shift changes. The study suggests that inhaled general anesthetics could act on proteins via altering protein motion on the microsecond-millisecond time scale, especially motion in the flexible loops that link different alpha helices. The validation of anesthetic effect on protein dynamics that are potentially correlated with protein functions is a critical step in unraveling the mechanisms of anesthetic action on proteins.     

Continuous matrix-assisted refolding of proteins

A refolding reactor was developed for continuous matrix-assisted refolding of proteins. The reactor was composed of an annular chromatography system and an ultrafiltration system to recycle aggregated proteins produced during the refolding reaction. The feed solution containing the denatured protein was continuously fed to the rotating bed perfused with buffer promoting folding of the protein. As the protein passed through the column, it was separated from chaotropic and reducing agents and the refolding process took place. Native proteins and aggregates could be continuously separated due to different molecular size. The exit stream containing aggregates was collected, concentrated by ultrafiltration and recycled to the feed solution. The high concentrations of chaotropic and reducing agents in the feed solution enabled dissociation of the recycled aggregates and consequently were fed again to the refolding reactor. When the initial feed mixture of denatured protein is used up, only buffer-containing chaotropic agents and recycled aggregates are fully converted to native protein. This process resulted in a stoichiometric conversion from the denatured protein to its correctly folded native state. The system was tested with bovine alpha-lactalbumin as model protein. Superdex 75 PrepGrade was used as size-exclusion medium. The yield of 30% active monomer in the batch process was improved to 41% at a recycling rate of 65%. Assuming that the aggregates can be redissolved and recycled into the feed stream in a quantitative manner, a refolding yield close to 100% is possible. The method can be also applied to other chromatographic principles suited for the separation of aggregates.     

Inhibition mode of a bisubstrate inhibitor of KDO8P synthase: a frequency-selective REDOR solid-state and solution NMR characterization

In this report the mode of inhibition of mechanism-based inhibitor (2, K(i) = 0.4 microM) of 3-deoxy-d-manno-2-octulosonate-8-phosphate synthase (KDO8PS), which was designed to mimic the combined key features of its natural substrates arabinose-5-phosphate (A5P) and phoshoenolpyruvate (PEP) into a single molecule, was investigated. Our earlier solid-state NMR observations identified the inhibitor to bind in a way that partly mimics A5P, while the phosphonate moiety of its PEP-mimicking part exhibits no interactions with enzyme residues. This result was apparently in disagreement with the competitive inhibition of 2 against PEP and with the later solved crystal structure of KDO8PS-2 binary complex identifying the interactions of its PEP-mimicking part with the enzyme residues that were not detected by solid-state NMR. To solve this discrepancy, further solid-state REDOR NMR and (31)P solution NMR experiments were applied to a variety of enzyme complexes with the substrates and inhibitor. In particular, a novel frequency-selective REDOR experiment was developed and applied. Integration of the solution and solid-state NMR data clearly demonstrates that under conditions of stoichiometric enzyme-ligand ratio at thermodynamic equilibrium (a) PEP binding is unperturbed by the presence of 2 and (b) both PEP and 2 can bind simultaneously to the synthase, i.e., form a ternary complex with PEP occupying its own subsite and 2 occupying A5P's subsite. The latter observation suggests that under the conditions used in our NMR measurements, the inhibition pattern of 2 against PEP should have a mixed type character. Furthermore, the NMR data directly demonstrate the distinction between the relative binding strength of the two moieties of 2: enzyme interactions with PEP-mimicking moiety are much weaker than those with the A5P moiety. This observation is in agreement with KDO8PS-2 crystal structure showing only remote contacts of the phosphonate due to large structural changes of binding site residues. It is concluded that these phosphonate-enzyme interactions evidenced by both (31)P solution NMR and X-ray are too weak to be preserved under the lyophilization of KDO8PS-2 binary complex and therefore are not evidenced by the solid-state REDOR spectra.     

The protein folding network indicates that the ultrafast folding mutant of villin headpiece subdomain has a deeper folding funnel

Protein folding is a dynamic process with continuous transitions among different conformations. In this work, the dynamics in the protein folding network of villin headpiece subdomain (HP35) has been investigated based on multiple reversible folding trajectories of HP35 and its ultrafast folding mutant where sub-angstrom folding was achieved. The four folding states were clearly separated on the network, validating the classification of the states. Examination of the eight conformers with different formation of the individual helices revealed high plasticity of the three helices in all the four states. A consistent feature between the wild type and mutant protein is the dominant conformer 111 (all three helices formed) in the folded state and conformers 111 and 011 (helices II and III formed) in the major intermediate state, indicating the critical role of helices II and III in the folding mechanism. When compared to the wild type, the folding landscape of the ultrafast folding mutant exhibited a deeper folding funnel towards the folded state. The very beginning of the folding (0-10 ns) was very similar for both protein variants but it soon diverged and displayed different folding pathways. Although going through the major intermediate state is the dominant pathway for both, it was also observed that some folding went through the minor intermediate state for the mutant. The intriguing difference resulting from the mutation at two residues in helix III has been carefully analyzed and discussed in details.     

Self-assembly of chiral trans-cyclobutane-containing β-dipeptides into ordered aggregates

Two chiral synthetic β-dipeptides have been constructed, one with two trans-cyclobutane residues and the other with one trans and one cis fragment, 1 and 2, respectively, and investigated to get insight into the non-covalent interactions responsible for their self-assembly to form ordered aggregates, as well into parameters such as their morphology and size. Experimental evidence of the formation of these assemblies was provided by spectroscopy, microscopy and X-ray diffraction experiments that suggest the formation of nanoscale helical aggregates. This process involves a conformational change in the molecules of each dipeptide with respect to the preferred conformation of the isolated molecules in solution. A high-resolution NMR spectroscopy study allowed the determination of the dynamics of the gelation process in [D(8)]toluene and the sol-gel transition temperature, which was around 270 K in this solvent at a concentration of 15 mM. NMR spectroscopy experiments also provided some information about conformational changes involved in the sol-gel transition and also suggested a different gel packing for each dipeptide. These observations have been nicely explained by computational studies. The self-assembly of the molecules has been modelled and suggested a head-to-head molecular arrangement for 1 and a head-to-tail arrangement for 2 to give helical structures corresponding to hydrogen-bonded single chains. These chains interact with one another in an antiparallel way to afford bundles, the significant geometry parameters of which fit well to the main peaks observed in wide-angle X-ray diffraction spectra of the aggregates in the solid state.     

Aromatic Cluster Sensor of Protein Folding: Near‐UV Electronic Circular Dichroism Bands Assigned to Fold Compactness

Both far‐ and near‐UV electronic circular dichroism (ECD) spectra have bands sensitive to thermal unfolding of Trp and Tyr residues containing proteins. Beside spectral changes at 222 nm reporting secondary structural variations (far‐UV range), L_b bands (near‐UV range) are applicable as 3D‐fold sensors of protein's core structure. In this study we show that both L_b(Tyr) and L_b(Trp) ECD bands could be used as sensors of fold compactness. ECD is a relative method and thus requires NMR referencing and cross‐validation, also provided here. The ensemble of 204 ECD spectra of Trp‐cage miniproteins is analysed as a training set for “calibrating” Trp?Tyr folded systems of known NMR structure. While in the far‐UV ECD spectra changes are linear as a function of the temperature, near‐UV ECD data indicate a non‐linear and thus, cooperative unfolding mechanism of these proteins. Ensemble of ECD spectra deconvoluted gives both conformational weights and insight to a protein folding?unfolding mechanism. We found that the L_b²⁹³ band is reporting on the 3D‐structure compactness. In addition, the pure near‐UV ECD spectrum of the unfolded state is described here for the first time. Thus, ECD folding information now validated can be applied with confidence in a large thermal window (5≤T≤85 °C) compared to NMR for studying the unfolding of Trp?Tyr residue pairs. In conclusion, folding propensities of important proteins (RNA polymerase II, ubiquitin protein ligase, tryptase‐inhibitor etc.) can now be analysed with higher confidence.     

Magic-angle spinning solid-state NMR spectroscopy of the beta1 immunoglobulin binding domain of protein G (GB1): 15N and 13C chemical shift assignments and conformational analysis

Magic-angle spinning solid-state NMR (SSNMR) studies of the beta1 immunoglobulin binding domain of protein G (GB1) are presented. Chemical shift correlation spectra at 11.7 T (500 MHz 1H frequency) were employed to identify signals specific to each amino acid residue type and to establish backbone connectivities. High sensitivity and resolution facilitated the detection and assignment of every 15N and 13C site, including the N-terminal (M1) 15NH3, the C-terminal (E56) 13C', and side-chain resonances from residues exhibiting fast-limit conformational exchange near room temperature. The assigned spectra lend novel insight into the structure and dynamics of microcrystalline GB1. Secondary isotropic chemical shifts report on conformation, enabling a detailed comparison of the microcrystalline state with the conformation of single crystals and the protein in solution; the consistency of backbone conformation in these three preparations is the best among proteins studied so far. Signal intensities and line widths vary as a function of amino acid position and temperature. High-resolution spectra are observed near room temperature (280 K) and at <180 K, whereas resolution and sensitivity greatly degrade substantially near 210 K; the magnitude of this effect is greatest among the side chains of residues at the intermolecular interface of the microcrystal lattice, which we attribute to intermediate-rate translational diffusion of solvent molecules near the glass transition. These features of GB1 will enable its use as an excellent model protein not only for SSNMR methods development but also for fundamental studies of protein thermodynamics in the solid state.