Using NMR diffusion data to validate MD models of disordered proteins: Test case of N-terminal tail of histone H4.

MD simulations can provide uniquely detailed models of intrinsically disordered proteins (IDPs). However, these models need careful experimental validation. The coefficient of translational diffusion Dtr, measurable by pulsed field gradient NMR, offers a potentially useful piece of experimental information related to the compactness of the IDP's conformational ensemble. Here, we investigate, both experimentally and via the MD modeling, the translational diffusion of a 25-residue N-terminal fragment from histone H4 (N-H4). We found that the predicted values of Dtr, as obtained from mean-square displacement of the peptide in the MD simulations, are largely determined by the viscosity of the MD water (which has been reinvestigated as a part of our study). Beyond that, our analysis of the diffusion data indicates that MD simulations of N-H4 in the TIP4P-Ew water give rise to an overly compact conformational ensemble for this peptide. In contrast, TIP4P-D and OPC simulations produce the ensembles that are consistent with the experimental Dtr result. These observations are supported by the analyses of the 15N spin relaxation rates. We also tested a number of empirical methods to predict Dtr based on IDP's coordinates extracted from the MD snapshots. In particular, we show that the popular approach involving the program HYDROPRO can produce misleading results. This happens because HYDROPRO is not intended to predict the diffusion properties of highly flexible biopolymers such as IDPs. Likewise, recent empirical schemes that exploit the relationship between the small-angle x-ray scattering-informed conformational ensembles of IDPs and the respective experimental Dtr values also prove to be problematic. In this sense, the first-principle calculations of Dtr from the MD simulations, such as demonstrated in this work, should provide a useful benchmark for future efforts in this area.

DPPM-Bridged Binuclear Pt(II) Pincer Complexes: Chemistry, Structure, and Photophysics in Solution Revisited.

A series of luminescent binuclear ([dppm{Pt(NNC)}2]2+) and mononuclear ([PPh3Pt(NNC)]+) complexes containing pincer ligands were synthesized and characterized. Photophysical characteristics of both types of complexes were studied in dichloromethane solution. In the solid phase, the binuclear compounds adopt a syn configuration where the {Pt(NNC)} fragments are held together due to intramolecular Pt-Pt bonding and π-stacking of the pincer ligand aromatic systems. Analysis of the complexes' molecular structure in solution by multinuclear NMR spectroscopy showed that the stacked intramolecular configuration is retained in fluid media, which is in complete agreement with a considerable red shift of the emission wavelength due to formation of the intramolecular Pt-Pt bond, leading to the transformation of an emissive excited state to 3MMLCT. It was also found that triethylamine quenches the emission of both types of complexes; the mechanism of quenching is a combination of dynamic and static channels of excited-state deactivation. In the case of binuclear complexes, deprotonation of the dppm methylene bridge by triethylamine also contributes to the chromophore quenching. To explain the observed chemistry of binuclear complex interactions with Et3N, a chemical equilibrium scheme was suggested, which was confirmed by quantitative monitoring of the 31P signal variations as a function of triethylamine concentration.

Insights into Ubiquitin Product Release in Hydrolysis Catalyzed by the Bacterial Deubiquitinase SdeA.

We report the co-crystal structure of the (catalytic Cys)-to-Ala mutant of the deubiquitinase domain of the Legionella pneumophila effector SdeA (SdeADUB) with its ubiquitin (Ub) product. Most of the intermolecular interactions are preserved in this product-bound structure compared to that of the previously characterized complex of SdeADUB with the suicide inhibitor ubiquitin vinylmethyl ester (Ub-VME), whose structure models the acyl-enzyme thioester intermediate. Nuclear magnetic resonance (NMR) titration studies show a chemical shift perturbation pattern that suggests that the same interactions also exist in solution. Isothermal titration calorimetry and NMR titration data reveal that the affinity of wild-type (WT) SdeADUB for Ub is significantly lower than that of the Cys-to-Ala mutant. This is potentially due to repulsive interaction between the thiolate ion of the catalytic Cys residue in WT SdeADUB and the carboxylate group of the C-terminal Gly76 residue in Ub. In the context of SdeADUB catalysis, this electrostatic repulsion arises after the hydrolysis of the scissile isopeptide bond in the acyl-enzyme intermediate and the consequent formation of the C-terminal carboxylic group in the Ub fragment. We hypothesize that this electrostatic repulsion may expedite the release of the Ub product by SdeADUB. We note that similar repulsive interactions may also occur in other deubiquitinases and hydrolases of ubiquitin-like protein modifiers and may constitute a fairly general mechanism of product release within this family. This is a potentially important feature for a family of enzymes that form extensive protein-protein interactions during enzyme-substrate engagement.

The Role of Rotational Motion in Diffusion NMR Experiments on Supramolecular Assemblies: Application to Sup35NM Fibrils.

Pulsed-field gradient (PFG) NMR is an important tool for characterization of biomolecules and supramolecular assemblies. However, for micrometer-sized objects, such as amyloid fibrils, these experiments become difficult to interpret because in addition to translational diffusion they are also sensitive to rotational diffusion. We have constructed a mathematical theory describing the outcome of PFG NMR experiments on rod-like fibrils. To test its validity, we have studied the fibrils formed by Sup35NM segment of the prion protein Sup35. The interpretation of the PFG NMR data in this system is fully consistent with the evidence from electron microscopy. Contrary to some previously expressed views, the signals originating from disordered regions in the fibrils can be readily differentiated from the similar signals representing small soluble species (e.g. proteolytic fragments). This paves the way for diffusion-sorted NMR experiments on complex amyloidogenic samples.

What Drives 15N Spin Relaxation in Disordered Proteins? Combined NMR/MD Study of the H4 Histone Tail.

Backbone (15N) NMR relaxation is one of the main sources of information on dynamics of disordered proteins. Yet, we do not know very well what drives 15N relaxation in such systems, i.e., how different forms of motion contribute to the measurable relaxation rates. To address this problem, we have investigated, both experimentally and via molecular dynamics simulations, the dynamics of a 26-residue peptide imitating the N-terminal portion of the histone protein H4. One part of the peptide was found to be fully flexible, whereas the other part features some transient structure (a hairpin stabilized by hydrogen bonds). The following motional modes proved relevant for 15N relaxation. 1) Sub-picosecond librations attenuate relaxation rates according to S2 ∼0.85-0.90. 2) Axial peptide-plane fluctuations along a stretch of the peptide chain contribute to relaxation-active dynamics on a fast timescale (from tens to hundreds of picoseconds). 3) φ/ψ backbone jumps contribute to relaxation-active dynamics on both fast (from tens to hundreds of picoseconds) and slow (from hundreds of picoseconds to a nanosecond) timescales. The major contribution is from polyproline II (PPII) ↔ ß transitions in the Ramachandran space; in the case of glycine residues, the major contribution is from PPII ↔ (ß) ↔ rPPII transitions, in which rPPII is the mirror-image (right-handed) version of the PPII geometry, whereas ß geometry plays the role of an intermediate state. 4) Reorientational motion of certain (sufficiently long-lived) elements of transient structure, i.e., rotational tumbling, contributes to slow relaxation-active dynamics on ∼1-ns timescale (however, it is difficult to isolate this contribution). In conclusion, recent advances in the area of force-field development have made it possible to obtain viable Molecular Dynamics models of protein disorder. After careful validation against the experimental relaxation data, these models can provide a valuable insight into mechanistic origins of spin relaxation in disordered peptides and proteins.

Coordination to Imidazole Ring Switches on Phosphorescence of Platinum Cyclometalated Complexes: The Route to Selective Labeling of Peptides and Proteins via Histidine Residues.

In this study, we have shown that substitution of chloride ligand for imidazole (Im) ring in the cyclometalated platinum complex Pt(phpy)(PPh3)Cl (1; phpy, 2-phenylpyridine; PPh3, triphenylphosphine), which is nonemissive in solution, switches on phosphorescence of the resulting compound. Crystallographic and nuclear magnetic resonance (NMR) spectroscopic studies of the substitution product showed that the luminescence ignition is a result of Im coordination to give the [Pt(phpy)(Im)(PPh3)]Cl complex. The other imidazole-containing biomolecules, such as histidine and histidine-containing peptides and proteins, also trigger luminescence of the substitution products. The complex 1 proved to be highly selective toward the imidazole ring coordination that allows site-specific labeling of peptides and proteins with 1 using the route, which is orthogonal to the common bioconjugation schemes via lysine, aspartic and glutamic acids, or cysteine and does not require any preliminary modification of a biomolecule. The utility of this approach was demonstrated on (i) site-specific modification of the ubiquitin, a small protein that contains only one His residue in its sequence, and (ii) preparation of nonaggregated HSA-based Pt phosphorescent probe. The latter particles easily internalize into the live HeLa cells and display a high potential for live-cell phosphorescence lifetime imaging (PLIM) as well as for advanced correlation PLIM and FLIM experiments.

Microsecond time-scale conformational exchange in proteins: using long molecular dynamics trajectory to simulate NMR relaxation dispersion data.

With the advent of ultra-long MD simulations it becomes possible to model microsecond time-scale protein dynamics and, in particular, the exchange broadening effects (R(ex)) as probed by NMR relaxation dispersion measurements. This new approach allows one to identify the exchanging species, including the elusive "excited states". It further helps to map out the exchange network, which is potentially far more complex than the commonly assumed 2- or 3-site schemes. Under fast exchange conditions, this method can be useful for separating the populations of exchanging species from their respective chemical shift differences, thus paving the way for structural analyses. In this study, recent millisecond-long MD trajectory of protein BPTI (Shaw et al. Science 2010, 330, 341) is employed to simulate the time variation of amide (15)N chemical shifts. The results are used to predict the exchange broadening of (15)N lines and, more generally, the outcome of the relaxation dispersion measurements using Carr-Purcell-Meiboom-Gill sequence. The simulated R(ex) effect stems from the fast (~10-100 µs) isomerization of the C14-C38 disulfide bond, in agreement with the prior experimental findings (Grey et al. J. Am. Chem. Soc. 2003, 125, 14324).

Effect of rotation in NMR diffusion experiments on micron-sized particles: A generalized theoretical treatment.

We have recently developed an analytical framework to interpret the results of pulsed field gradient (PFG) NMR experiments on solution samples of micron-sized amyloid fibrils [Angew. Chem. Int. Ed. 60 (2021) 15445-15451. https://doi.org/10.1002/anie.202102408]. Here we generalize this result by reporting a rigorous theoretical model of such experiments, implemented in a form of efficient computational scheme. In particular, the new treatment fully accounts for the anisotropy of fibrils' translational diffusion and takes into consideration the finite length of the gradient pulses. The results hold not only for the historic spin-echo sequence, but also for the widely used stimulated echo experiment. We have found that fibrils' rotation can attenuate the echo by a factor comparable with that of translation. However, contrary to some recent claims, the rotational mechanism cannot lead to an apparent fast-diffusion situation.

Structural and dynamic origins of ESR lineshapes in spin-labeled GB1 domain: the insights from spin dynamics simulations based on long MD trajectories.

Site-directed spin labeling (SDSL) ESR is a valuable tool to probe protein systems that are not amenable to characterization by x-ray crystallography, NMR or EM. While general principles that govern the shape of SDSL ESR spectra are known, its precise relationship with protein structure and dynamics is still not fully understood. To address this problem, we designed seven variants of GB1 domain bearing R1 spin label and recorded the corresponding MD trajectories (combined length 180 µs). The MD data were subsequently used to calculate time evolution of the relevant spin density matrix and thus predict the ESR spectra. The simulated spectra proved to be in good agreement with the experiment. Further analysis confirmed that the spectral shape primarily reflects the degree of steric confinement of the R1 tag and, for the well-folded protein such as GB1, offers little information on local backbone dynamics. The rotameric preferences of R1 side chain are determined by the type of the secondary structure at the attachment site. The rotameric jumps involving dihedral angles χ1 and χ2 are sufficiently fast to directly influence the ESR lineshapes. However, the jumps involving multiple dihedral angles tend to occur in (anti)correlated manner, causing smaller-than-expected movements of the R1 proxyl ring. Of interest, ESR spectra of GB1 domain with solvent-exposed spin label can be accurately reproduced by means of Redfield theory. In particular, the asymmetric character of the spectra is attributable to Redfield-type cross-correlations. We envisage that the current MD-based, experimentally validated approach should lead to a more definitive, accurate picture of SDSL ESR experiments.

Simple MD-based model for oxidative folding of peptides and proteins.

Significant strides have been recently made to fold peptides and small proteins in silico using MD simulations. However, facilities are currently lacking to include disulfide bonding in the MD models of protein folding. To address this problem, we have developed a simple empirical protocol to model formation of disulfides, which is perturbation-free, retains the same speed as conventional MD simulations and allows one to control the reaction rate. The new protocol has been tested on 15-aminoacid peptide guanylin containing four cysteine residues; the net simulation time using Amber ff14SB force field was 61 µs. The resulting isomer distribution is in qualitative agreement with experiment, suggesting that oxidative folding of guanylin in vitro occurs under kinetic control. The highly stable conformation of the so-called isomer 2(B) has been obtained for full-length guanylin, which is significantly different from the poorly ordered structure of the truncated peptide PDB ID 1GNB. In addition, we have simulated oxidative folding of guanylin within the 94-aminoacid prohormone proguanylin. The obtained structure is in good agreement with the NMR coordinates 1O8R. The proposed modeling strategy can help to explore certain fundamental aspects of protein folding and is potentially relevant for manufacturing of synthetic peptides and recombinant proteins.

Origin of anti-tumor activity of the cysteine-containing GO peptides and further optimization of their cytotoxic properties.

Antitumor GO peptides have been designed as dimerization inhibitors of prominent oncoprotein mucin 1. In this study we demonstrate that activity of GO peptides is independent of the level of cellular expression of mucin 1. Furthermore, these peptides prove to be broadly cytotoxic, causing cell death also in normal cells such as dermal fibroblasts and endometrial mesenchymal stem cells. To explore molecular mechanism of their cytotoxicity, we have designed and tested a number of new peptide sequences containing the key CxC or CxxC motifs. Of note, these sequences bear no similarity to mucin 1 except that they also contain a pair of proximal cysteines. Several of the new peptides turned out to be significantly more potent than their GO prototypes. The results suggest that cytotoxicity of these peptides stems from their (moderate) activity as disulfide oxidoreductases. It is expected that such peptides, which we have termed DO peptides, are involved in disulfide-dithiol exchange reaction, resulting in formation of adventitious disulfide bridges in cell proteins. In turn, this leads to a partial loss of protein function and rapid onset of apoptosis. We anticipate that coupling DO sequences with tumor-homing transduction domains can create a potentially valuable new class of tumoricidal peptides.

Onset of disorder and protein aggregation due to oxidation-induced intermolecular disulfide bonds: case study of RRM2 domain from TDP-43.

We have investigated the behavior of second RNA-recognition motif (RRM2) of neuropathological protein TDP43 under the effect of oxidative stress as modeled in vitro. Toward this end we have used the specially adapted version of H/D exchange experiment, NMR relaxation and diffusion measurements, dynamic light scattering, controlled proteolysis, gel electrophoresis, site-directed mutagenesis and microsecond MD simulations. Under oxidizing conditions RRM2 forms disulfide-bonded dimers that experience unfolding and then assemble into aggregate particles (APs). These particles are strongly disordered, highly inhomogeneous and susceptible to proteolysis; some of them withstand the dithiothreitol treatment. They can recruit/release monomeric RRM2 through thiol-disulfide exchange reactions. By using a combination of dynamic light scattering and NMR diffusion data we were able to approximate the size distribution function for the APs. The key to the observed aggregation behavior is the diminished ability of disulfide-bonded RRM2 dimers to refold and their increased propensity to misfold, which makes them vulnerable to large thermal fluctuations. The emerging picture provides detailed insight on how oxidative stress can contribute to neurodegenerative disease, with unfolding, aggregation, and proteolytic cleavage as different facets of the process.

Comment on "The ligand polyhedral model approach to the mechanism of complete carbonyl exchange in [Rh4(CO)12] and [Rh6(CO)16]" by Brian F. G. Johnson, Dalton Transactions, 2015, 44, DOI: 10.1039/C4DT03360D.

Experimental results (recent IR, DFT calculations and modern multinuclear NMR measurements on Rh-containing clusters, together with earlier VT multinuclear NMR measurements) show that the use of the Ligand Polyhedral Model (LPM) to provide a general mechanism for ligand fluxionality in Transition Metal Carbonyl Clusters (TMCCs) in solution cannot be sustained; instead there are numerous examples of only partial CO-migration over either part or sometimes the whole of the Rh-polyhedron as well as rhodium and carbonyl polyhedral rearrangements of Rh9- and Rh10-TMCCs containing an interstitial P when, in the high temperature limiting spectra, all the metals and all the carbonyls become equivalent and show time-averaged values of (1)J(Rh-P) and (2)J(P-CO) respectively.

Microsecond time-scale dynamics from relaxation in the rotating frame: experiments using spin lock with alternating phase.

A spin lock comprised of radiofrequency pulses with alternating phase, (x) (-x)(x) (-x) , is proposed as a new technique to probe microsecond time-scale dynamics. A series of R1rho measurements using different pulse duration tp allows one to determine exchange rate, kex, the product p(a)p(b)(Delta omega(ab))2 involving populations of the exchanging species, p(a) and p(b), together with chemical shift difference, (Delta omega(ab)), and the strength of the spin-lock field, B1. The interpretation is based on simple analytical expression for R1rho derived on the basis of Redfield theory. The application of the method is demonstrated for partially deuterated molecule of cyclohexane undergoing chair-to-chair interconversion at -9 degrees C.

Paramagnetic relaxation enhancements in unfolded proteins: theory and application to drkN SH3 domain.

Site-directed spin labeling in combination with paramagnetic relaxation enhancement (PRE) measurements is one of the most promising techniques for studying unfolded proteins. Since the pioneering work of Gillespie and Shortle (J Mol Biol 1997;268:158), PRE data from unfolded proteins have been interpreted using the theory that was originally developed for rotational spin relaxation. At the same time, it can be readily recognized that the relative motion of the paramagnetic tag attached to the peptide chain and the reporter spin such as (1)H(N) is best described as a translation. With this notion in mind, we developed a number of models for the PRE effect in unfolded proteins: (i) mutual diffusion of the two tethered spheres, (ii) mutual diffusion of the two tethered spheres subject to a harmonic potential, (iii) mutual diffusion of the two tethered spheres subject to a simulated mean-force potential (Smoluchowski equation); (iv) explicit-atom molecular dynamics simulation. The new models were used to predict the dependences of the PRE rates on the (1)H(N) residue number and static magnetic field strength; the results are appreciably different from the Gillespie-Shortle model. At the same time, the Gillespie-Shortle approach is expected to be generally adequate if the goal is to reconstruct the distance distributions between (1)H(N) spins and the paramagnetic center (provided that the characteristic correlation time is known with a reasonable accuracy). The theory has been tested by measuring the PRE rates in three spin-labeled mutants of the drkN SH3 domain in 2M guanidinium chloride. Two modifications introduced into the measurement scheme-using a reference compound to calibrate the signals from the two samples (oxidized and reduced) and using peak volumes instead of intensities to determine the PRE rates-lead to a substantial improvement in the quality of data. The PRE data from the denatured drkN SH3 are mostly consistent with the model of moderately expanded random-coil protein, although part of the data point toward a more compact structure (local hydrophobic cluster). At the same time, the radius of gyration reported by Choy et al. (J Mol Biol 2002;316:101) suggests that the protein is highly expanded. This seemingly contradictory evidence can be reconciled if one assumes that denatured drkN SH3 forms a conformational ensemble that is dominated by extended conformations, yet also contains compact (collapsed) species. Such behavior is apparently more complex than predicted by the model of a random-coil protein in good solvent/poor solvent.

Synthesis and structural characterization of two novel heterometallic clusters: [Rh4Pt2(CO)11(dppm)2] and [Ru2Rh2Pt2(CO)12(dppm)2].

Two novel heterometallic octahedral clusters [Rh(4)Pt(2)(CO)(11)(dppm)(2)](1) and [Ru(2)Rh(2)Pt(2)(CO)(12)(dppm)(2)](2) were synthesized by the reaction of [Rh(2)Pt(2)(CO)(6)(dppm)(2)] with [Rh(6)(CO)(14)(NCMe)(2)] and Ru(3)(CO)(12), respectively. Solid state structures of 1 and 2 have been established by a single crystal X-ray diffraction study. Two dppm ligands in 1 are bonded to one platinum and three rhodium atoms, which form an equatorial plane of the Rh(4)Pt(2) octahedron. Two rhodium and two platinum atoms bound to the diphosphine ligands in 2 are nonplanar to give an octahedral C2 symmetric Ru(2)Rh(2)Pt(2)(dppm)2 framework. The (31)P NMR investigation of and (1D, (31)P COSY, (31)P-[(103)Rh] HMQC) and simulation of 1D spectral patterns showed that in both clusters the structures of the M(6)(PP)(2) fragments found in the solid state are maintained in solution.

Tripyrrolylphosphine as a unique bridging ligand in the Rh6CO14(mu2-P(NC4H4)3) cluster: structure, bonding, fluxionality, thermodynamics, and kinetics studies.

Tripyrrolylphosphine reacts with the cluster Rh6CO15(NCMe) to afford the disubstituted Rh6CO14(mu2)-P(NC4H4)3) derivative (2) via the Rh6CO15P(NC4H4)3 intermediate (1) with eta(1)-P coordination. In the solid state, 2 has the phosphine occupying a bridging position where it is bonded to two neighboring Rh atoms in the Rh(6) octahedron through the P-atom and an approximately tetrahedral alpha-carbon atom of one of the pyrrolyl rings. This can be described by the interaction of an electron pair, associated with a negative charge on one of the canonical forms of the NC(4)H(4) ring, with the adjacent Rh center. (1)H NMR spectra show that the solid-state structure is retained in solution, but the phosphine is not rigid, and three distinctive dynamic processes are observed. Each of these represents independent hindered rotation of inequivalent pyrrolyl rings about P-N bonds, the ring involved in the interaction with the Rh(6) skeleton displaying the highest activation barrier with deltaH = 15.8 +/- 0.1 kcal mol(-1) and deltaS = 1.4 +/- 0.3 cal K(-1) mol(-1). The assignment has been confirmed by 1H TOCSY and EXSY experiments, and a mechanism is proposed. The formation of 2 from 1 is reversible in the presence of CO, which is highly unusual for bridged clusters. The kinetics of the forward and reverse reactions have been studied, and the values of DeltaH degrees and DeltaS degrees for formation of 2 (+1.3 +/- 0.5 kcal mol(-1) and -9 +/- 2 cal K(-1) mol (-1), respectively) show that the Rh-C bond in the bridge is comparable in strength with the Rh-CO bond it replaces. The intrinsic entropy of 2 is exceptionally unfavorable, overcoming the favorable entropy caused by CO release, and this allows the reversibility of bridge formation. The reactions proceed via a reactive intermediate that may involve agostic bonding of the ring. The reverse reaction has an exceedingly unfavorable activation entropy that emphasizes the unique nature of 2.