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.

X-ray Crystallography Module in MD Simulation Program Amber 2023. Refining the Models of Protein Crystals.

The MD simulation package Amber offers an attractive platform to refine crystallographic structures of proteins: (i) state-of-the-art force fields help to regularize protein coordinates and reconstruct the poorly diffracting elements of the structure, such as flexible loops; (ii) MD simulations restrained by the experimental diffraction data provide an effective strategy to optimize structural models of protein crystals, including explicitly modeled interstitial solvent as well as crystal contacts. Here, we present the new crystallography module xray, released as a part of the Amber 2023 package. This module contains functions to calculate and scale structure factors (including the contributions from bulk solvent), evaluate the maximum-likelihood-type crystallographic potential, and compute its derivative forces. The X-ray functionality of Amber no longer relies on external dependencies so that the full advantage of GPU acceleration can be taken. This makes it possible to refine in a short time hundreds of crystal models, including supercell models comprised of multiple unit cells. The new automated Amber-based refinement procedure leads to an appreciable improvement in Rfree (in some cases, by as much as 0.067) as well as MolProbity scores.

Histone H4 Tails in Nucleosomes: a Fuzzy Interaction with DNA.

The interaction of positively charged N-terminal histone tails with nucleosomal DNA plays an important role in chromatin assembly and regulation, modulating their susceptibility to post-translational modifications and recognition by chromatin-binding proteins. Here, we report residue-specific 15 N NMR relaxation rates for histone H4 tails in reconstituted nucleosomes. These data indicate that H4 tails are strongly dynamically disordered, albeit with reduced conformational flexibility compared to a free peptide with the same sequence. Remarkably, the NMR observables were successfully reproduced in a 2-µs MD trajectory of the nucleosome. This is an important step toward resolving an apparent inconsistency where prior simulations were generally at odds with experimental evidence on conformational dynamics of histone tails. Our findings indicate that histone H4 tails engage in a fuzzy interaction with nucleosomal DNA, underpinned by a variable pattern of short-lived salt bridges and hydrogen bonds, which persists at low ionic strength (0-100 mM NaCl).

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.

Molecular Dynamics model of peptide-protein conjugation: case study of covalent complex between Sos1 peptide and N-terminal SH3 domain from Grb2.

We have investigated covalent conjugation of VPPPVPPRRRX' peptide (where X' denotes Nε-chloroacetyl lysine) to N-terminal SH3 domain from adapter protein Grb2. Our experimental results confirmed that the peptide first binds to the SH3 domain noncovalently before establishing a covalent linkage through reaction of X' with the target cysteine residue C32. We have also confirmed that this reaction involves a thiolate-anion form of C32 and follows the SN2 mechanism. For this system, we have developed a new MD-based protocol to model the formation of covalent conjugate. The simulation starts with the known coordinates of the noncovalent complex. When two reactive groups come into contact during the course of the simulation, the reaction is initiated. The reaction is modeled via gradual interpolation between the two sets of force field parameters that are representative of the noncovalent and covalent complexes. The simulation proceeds smoothly, with no appreciable perturbations to temperature, pressure or volume, and results in a high-quality MD model of the covalent complex. The validity of this model is confirmed using the experimental chemical shift data. The new MD-based approach offers a valuable tool to explore the mechanics of protein-peptide conjugation and build accurate models of covalent complexes.

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.

A new structural arrangement in proteins involving lysine NH3+ group and carbonyl.

Screening of the Protein Data Bank led to identification of a recurring structural motif where lysine NH3+ group interacts with backbone carbonyl. This interaction is characterized by linear atom arrangement, with carbonyl O atom positioned on the three-fold symmetry axis of the NH3+ group (angle Cε-Nζ-O close to 180°, distance Nζ-O ca. 2.7-3.0 Å). Typically, this linear arrangement coexists with three regular hydrogen bonds formed by lysine NH3+ group (angle Cε-Nζ-acceptor atom close to 109°, distance Nζ-acceptor atom ca. 2.7-3.0 Å). Our DFT calculations using polarizable continuum environment suggest that this newly identified linear interaction makes an appreciable contribution to protein's energy balance, up to 2 kcal/mol. In the context of protein structure, linear interactions play a role in capping the C-termini of α-helices and 310-helices. Of note, linear interaction involving conserved lysine is consistently found in the P-loop of numerous NTPase domains, where it stabilizes the substrate-binding conformation of the P-loop. Linear interaction NH3+ - carbonyl represents an interesting example of ion-dipole interactions that has so far received little attention compared to ion-ion interactions (salt bridges) and dipole-dipole interactions (hydrogen bonds), but nevertheless represents a distinctive element of protein architecture.

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.

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.

Slow conformational exchange and overall rocking motion in ubiquitin protein crystals.

Proteins perform their functions in solution but their structures are most frequently studied inside crystals. Here we probe how the crystal packing alters microsecond dynamics, using solid-state NMR measurements and multi-microsecond MD simulations of different crystal forms of ubiquitin. In particular, near-rotary-resonance relaxation dispersion (NERRD) experiments probe angular backbone motion, while Bloch-McConnell relaxation dispersion data report on fluctuations of the local electronic environment. These experiments and simulations reveal that the packing of the protein can significantly alter the thermodynamics and kinetics of local conformational exchange. Moreover, we report small-amplitude reorientational motion of protein molecules in the crystal lattice with an ~3-5° amplitude on a tens-of-microseconds time scale in one of the crystals, but not in others. An intriguing possibility arises that overall motion is to some extent coupled to local dynamics. Our study highlights the importance of considering the packing when analyzing dynamics of crystalline proteins.X-ray crystallography is the main method for protein structure determination. Here the authors combine solid-state NMR measurements and molecular dynamics simulations and show that crystal packing alters the thermodynamics and kinetics of local conformational exchange as well as overall rocking motion of protein molecules in the crystal lattice.