NMR Methods for Characterizing the Basic Side Chains of Proteins: Electrostatic Interactions, Hydrogen Bonds, and Conformational Dynamics

Solid state NMR spectral editing of histidine, arginine and lysine using Hadamard encoding

The NMR signals from protein sidechains are rich in information about intra- and inter-molecular interactions, but their detection can be complicated due to spectral overlap as well as conformational and hydrogen exchange. In this work, we demonstrate a protocol for multi-dimensional solid-state NMR spectral editing of signals from basic sidechains based on Hadamard matrix encoding. The Hadamard method acquires multi-dimensional experiments in such a way that both the backbone and under-sampled sidechain signals can be decoded for unambiguous editing in the 15 N spectral frequency dimension. All multi-dimensional 15 N-edited solid-state NMR experiments can be acquired using this strategy, thereby accelerating the acquisition of spectra spanning broad frequency bandwidth. Application of these methods to the ferritin nanocage, reveals signals from N atoms from His, Arg, Lys and Trp sidechains, as well as their tightly bound, ordered water molecules. The Hadamard approach adds to the arsenal of spectroscopic approaches for protein NMR signal detection.

DiAMoNDBack: Diffusion-Denoising Autoregressive Model for Non-Deterministic Backmapping of Cα Protein Traces

Coarse-grained molecular models of proteins permit access to length and time scales unattainable by all-atom models and the simulation of processes that occur on long time scales, such as aggregation and folding. The reduced resolution realizes computational accelerations, but an atomistic representation can be vital for a complete understanding of mechanistic details. Backmapping is the process of restoring all-atom resolution to coarse-grained molecular models. In this work, we report DiAMoNDBack (Diffusion-denoising Autoregressive Model for Non-Deterministic Backmapping) as an autoregressive denoising diffusion probability model to restore all-atom details to coarse-grained protein representations retaining only Cα coordinates. The autoregressive generation process proceeds from the protein N-terminus to C-terminus in a residue-by-residue fashion conditioned on the Cα trace and previously backmapped backbone and side-chain atoms within the local neighborhood. The local and autoregressive nature of our model makes it transferable between proteins. The stochastic nature of the denoising diffusion process means that the model generates a realistic ensemble of backbone and side-chain all-atom configurations consistent with the coarse-grained Cα trace. We train DiAMoNDBack over 65k+ structures from the Protein Data Bank (PDB) and validate it in applications to a hold-out PDB test set, intrinsically disordered protein structures from the Protein Ensemble Database (PED), molecular dynamics simulations of fast-folding mini-proteins from DE Shaw Research, and coarse-grained simulation data. We achieve state-of-the-art reconstruction performance in terms of correct bond formation, avoidance of side-chain clashes, and the diversity of the generated side-chain configurational states. We make the DiAMoNDBack model publicly available as a free and open-source Python package.

Biomolecular NMR spectroscopy in the era of artificial intelligence

Biomolecular nuclear magnetic resonance (NMR) spectroscopy and artificial intelligence (AI) have a burgeoning synergy. Deep learning-based structural predictors have forever changed structural biology, yet these tools currently face limitations in accurately characterizing protein dynamics, allostery, and conformational heterogeneity. We begin by highlighting the unique abilities of biomolecular NMR spectroscopy to complement AI-based structural predictions toward addressing these knowledge gaps. We then highlight the direct integration of deep learning approaches into biomolecular NMR methods. AI-based tools can dramatically improve the acquisition and analysis of NMR spectra, enhancing the accuracy and reliability of NMR measurements, thus streamlining experimental processes. Additionally, deep learning enables the development of novel types of NMR experiments that were previously unattainable, expanding the scope and potential of biomolecular NMR spectroscopy. Ultimately, a combination of AI and NMR promises to further revolutionize structural biology on several levels, advance our understanding of complex biomolecular systems, and accelerate drug discovery efforts.

NMR detection and conformational dependence of two, three, and four-bond isotope shifts due to deuteration of backbone amides

NMR isotope shifts occur due to small differences in nuclear shielding when nearby atoms are different isotopes. For molecules dissolved in 1:1 H2O:D2O, the resulting mixture of N-H and N-D isotopes leads to a small splitting of resonances from adjacent nuclei. We used multidimensional NMR to measure isotope shifts for the proteins CUS-3iD and CspA. We observed four-bond 4∆N(ND) isotope shifts in high-resolution 2D 15N-TROSY experiments of the perdeuterated proteins that correlate with the torsional angle psi. Three-bond 3∆C'(ND) isotope shifts detected in H(N)CO spectra correlate with the intraresidue H-O distance, and to a lesser extent with the dihedral angle phi. The conformational dependence of the isotope shifts agree with those previously reported in the literature. Both the 4∆N(ND) and 3∆C'(ND) isotope shifts are sensitive to distances between the atoms giving rise to the isotope shifts and the atoms experiencing the splitting, however, these distances are strongly correlated with backbone dihedral angles making it difficult to resolve distance from stereochemical contributions to the isotope shift. H(NCA)CO spectra were used to measure two-bond 2∆C'(ND) isotope shifts and [D]/[H] fractionation factors. Neither parameter showed significant differences for hydrogen-bonded sites, or changes over a 25° temperature range, suggesting they are not sensitive to hydrogen bonding. Finally, the quartet that arises from the combination of 2∆C'(ND) and 3∆C'(ND) isotope shifts in H(CA)CO spectra was used to measure synchronized hydrogen exchange for the sequence neighbors A315-S316 in the protein CUS-3iD. In many of our experiments we observed minor resonances due to the 10% D2O used for the sample deuterium lock, indicating isotope shifts can be a source of spectral heterogeneity in standard NMR experiments. We suggest that applications of isotope shifts such as conformational analysis and correlated hydrogen exchange could benefit from the larger magnetic fields becoming available.

Slow Rotational Dynamics of Cytosine NH2 Groups in Double-Stranded DNA

Aromatic NH₂ groups are essential as hydrogen-bond donors in secondary structures of DNA and RNA. Although rapid rotations of NH₂ groups of adenine and guanine bases were previously characterized, there has been a lack of quantitative information about slow rotations of cytosine NH₂ groups in Watson-Crick base pairs. In this study, using an NMR method we had recently developed, we determined the kinetic rate constants and energy barriers for cytosine NH₂ rotations in a 15-base-pair DNA duplex. Our data show that the rotational dynamics of cytosine NH₂ groups depend on local environments. Qualitative correlation between the ranges of ¹⁵N chemical shifts and rotational time scales for various NH₂ groups of nucleic acids and proteins illuminates a relationship between the partial double-bond character of the C-N bond and the time scale for NH₂ rotations.

Hindered Rotations of Protein Asparagine/Glutamine Side-Chain NH2 Groups: Impact of Hydrogen Bonding with DNA

Hindered rotation about an sp² C-N bond is known to occur in arginine (Arg), asparagine (Asn), and glutamine (Gln) side chains of proteins. However, very little is known about the rotational dynamics of Asn and Gln side-chain NH₂ groups. Here, using a unique NMR method, we quantitatively characterized the hindered rotations of protein Asn/Gln side-chain NH₂ groups. This NMR method yields simple NH₂-selective spectra that allow for an accurate determination of the kinetic rate constants for the hindered rotations. Through the NMR measurements at different temperatures, we investigated the energy barriers that restrict the C-N bond rotations of protein side-chain NH₂ groups. Through a comparison of the kinetic data for the free and DNA-bound states of the Antp homeodomain, we also examined the impact of hydrogen bonding on the hindered rotations of the side-chain NH₂ groups. Our data suggest that the hydrogen bonding increases the energy barriers by 1-6 kJ/mol.

Minimizing Pervasive Artifacts in 4D Covariance Maps for Protein Side Chain NMR Assignments

Nuclear magnetic resonance (NMR) is a mainstay of biophysical studies that provides atomic level readouts to formulate molecular mechanisms. Side chains are particularly important to derive mechanisms involving proteins as they carry functional groups, but NMR studies of side chains are often limited by challenges in assigning their signals. Here, we designed a novel computational method that combines spectral derivatives and matrix square-rooting to produce reliable 4D covariance maps from routinely acquired 3D spectra and facilitates side chain resonance assignments. Thus, we generate two 4D maps from 3D-HcccoNH and 3D-HCcH-TOCSY spectra that each help overcome signal overlap or sensitivity losses. These 4D maps feature HC-HSQCs of individual side chains that can be paired to assigned backbone amide resonances of individual aliphatic signals, and both are obtained from a single modified covariance calculation. Further, we present 4D maps produced using conventional triple resonance experiments to easily assign asparagine side chain amide resonances. The 4D covariance maps encapsulate the lengthy manual pattern recognition used in traditional assignment methods and distill the information as correlations that can be easily visualized. We showcase the utility of the 4D covariance maps with a 10 kDa peptidyl carrier protein and a 52 kDa cyclization domain from a nonribosomal peptide synthetase.

Quantifying and visualizing weak interactions between anions and proteins

The molecular properties of proteins are influenced by various ions present in the same solution. While site-specific strong interactions between multivalent metal ions and proteins are well characterized, the behavior of other ions that are only weakly interacting with proteins remains elusive. In the current study, using NMR spectroscopy, we have investigated anion-protein interactions for three proteins that are similar in size but differ in overall charge. Using a unique NMR-based approach, we quantified anions accumulated around the proteins. The determined numbers of anions that are electrostatically attracted to the charged proteins were notably smaller than the overall charge valences and were consistent with predictions from the Poisson-Boltzmann theory. This NMR-based approach also allowed us to measure ionic diffusion and characterize the anions interacting with the positively charged proteins. Our data show that these anions rapidly diffuse while bound to the proteins. Using the same experimental approach, we observed the release of the anions from the protein surface upon the formation of the Antp homeodomain-DNA complex. Using paramagnetic relaxation enhancement (PRE), we visualized the spatial distribution of anions around the free proteins and the Antp homeodomain-DNA complex. The obtained PRE data revealed the localization of anions in the vicinity of the highly positively charged regions of the free Antp homeodomain and provided further evidence of the release of anions from the protein surface upon the protein-DNA association. This study sheds light on the dynamic behavior of anions that electrostatically interact with proteins.

Experimental approaches for investigating ion atmospheres around nucleic acids and proteins

Ionic interactions are crucial to biological functions of DNA, RNA, and proteins. Experimental research on how ions behave around biological macromolecules has lagged behind corresponding theoretical and computational research. In the 21st century, quantitative experimental approaches for investigating ionic interactions of biomolecules have become available and greatly facilitated examinations of theoretical electrostatic models. These approaches utilize anomalous small-angle X-ray scattering, atomic emission spectroscopy, mass spectrometry, or nuclear magnetic resonance (NMR) spectroscopy. We provide an overview on the experimental methodologies that can quantify and characterize ions within the ion atmospheres around nucleic acids, proteins, and their complexes.

Searching for Aquamelt Behavior among Silklike Biomimetics during Fibrillation under Flow

In this paper, we elucidate a generic mechanism behind strain-induced phase transition in aqueous solutions of silk-inspired biomimetics by atomistic molecular dynamics simulations. We show the results of modeling of homopeptides polyglycine Gly₃₀ and polyalanine Ala₃₀ and a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)₅, i.e., the simplest and yet relevant sequences that could mimic the behavior of natural silk under stress conditions. First, we analyze hydrophobicities of the sequences by calculating the Gibbs free energy of hydration and inspecting the interchain hydrogen bonding and hydration by water. Second, the force-extension profiles are scanned and compared with the results for poly(ethylene oxide), the synthetic polymer for which the aquamelt behavior has been proved recently. Additionally, the conformational transitions of oligopeptides from coiled to extended states are characterized by a generalized order parameter and by the dependence of the solvent-accessible surface area of the chains on applied stretching. Fibrillation itself is surveyed using both the two-dimensional interchain pair correlation function and the SAXS/WAXS patterns for the aggregates formed under stress. These are compared with experimental data found in the literature on fibril structure of silk composite materials doped with oligoalanine peptides. Our results show that tensile stress introduced into aqueous oligopeptide solutions facilitates interchain interactions. The oligopeptides display both a greater resistance to extension as compared to poly(ethylene oxide) and a reduced ability for hydrogen bonding of the stretched chains between oligomers and with water. Fiber formation is proved for all simulated objects, but the most structured one is made of a heteropeptide (Gly-Ala-Gly-Ala-Gly-Ser)₅: For this sequence, we obtain the highest degree of the secondary structure motifs in the fiber. We conclude that this is the most promising candidate among considered sequences to find the aquamelt behavior in further experimental studies.

Racemic phosphorothioate as a tool for NMR investigations of protein-DNA complexes

A major driving force for protein-nucleic acid association is electrostatic interactions via ion pairs of the positively charged basic side chains and negatively charged phosphates. For a better understanding of how proteins scan DNA and recognize particular signatures, it is important to gain atomic-level insight into the behavior of basic side chains at the protein-DNA interfaces. NMR spectroscopy is a powerful tool for investigating the structural, dynamic, and kinetic aspects of protein-DNA interactions. However, resonance assignment of basic side-chain cationic moieties at the molecular interfaces remains to be a major challenge. Here, we propose a fast, robust, and inexpensive approach that greatly facilitates resonance assignment of interfacial moieties and also allows for kinetic measurements of protein translocation between two DNA duplexes. This approach utilizes site-specific incorporation of racemic phosphorothioate at the position of a phosphate that interacts with a protein side chain. This modification retains the electric charge of phosphate and therefore is mild, but causes significant chemical shift perturbations for the proximal protein side chains, which facilitates resonance assignment. Due to the racemic nature of the modification, two different chemical shifts are observed for the species with different diastereomers R_P and S_P of the incorporated phosphorothioate group. Kinetic information on the exchange of the protein molecule between R_P and S_P DNA duplexes can be obtained by ¹⁵N_z exchange spectroscopy. We demonstrate the applications of this approach to the Antennapedia homeodomain-DNA complex and the CREB1 basic leucine-zipper (bZIP)-DNA complex.

Hydrogen-exchange kinetics studied through analysis of self-decoupling of nuclear magnetic resonance

Hydrogen exchange between solute and water molecules occurs across a wide range of timescales. Rapid hydrogen-exchange processes can effectively diminish ¹H-¹⁵N scalar couplings. We demonstrate that the self-decoupling of ¹⁵N nuclear magnetic resonance can allow quantitative investigations of hydrogen exchange on a micro- to millisecond timescale, which is relatively difficult to analyze with other methods. Using a Liouvillian matrix incorporating hydrogen exchange as a mechanism for scalar relaxation, the hydrogen exchange rate can be determined from ¹⁵N NMR line shapes recorded with and without ¹H decoupling. Self-decoupling offers a simple approach to analyze the kinetics of hydrogen exchange in a wide range of timescale.

NMR Observation of Intermolecular Hydrogen Bonds between Protein Tyrosine Side-Chain OH and DNA Phosphate Groups

Hydrogen bonds between protein side-chain hydroxyl (OH) and phosphate groups are one of the most common types of intermolecular hydrogen bonds in protein-DNA/RNA complexes. Using NMR spectroscopy, we identified and characterized the hydrogen bonds between tyrosine side-chain OH and DNA phosphate groups in a protein-DNA complex. These OH groups exhibited relatively slow hydrogen-exchange rates and sizable scalar couplings between hydroxyl 1H and DNA phosphate 31P nuclei across the hydrogen bonds. Information about intermolecular hydrogen bonds facilitates investigations of the DNA/RNA recognition by the protein.