Side-chain Dynamics Research Articles

Toward On-Demand Polymorphic Transitions of Organic Crystals via Side Chain and Lattice Dynamics Engineering.

Controlling polymorphism, namely, the occurrence of multiple crystal forms for a given compound, is still an open technological challenge that needs to be addressed for the reliable manufacturing of crystalline functional materials. Here, we devised a series of 13 organic crystals engineered to embody molecular fragments undergoing specific nanoscale motion anticipated to drive cooperative order-disorder phase transitions. By combining polarized optical microscopy coupled with a heating/cooling stage, differential scanning calorimetry, X-ray diffraction, low-frequency Raman spectroscopy, and calculations (density functional theory and molecular dynamics), we proved the occurrence of cooperative transitions in all the crystalline systems, and we demonstrated how both the molecular structure and lattice dynamics play crucial roles in these peculiar solid-to-solid transformations. These results introduce an efficient strategy to design polymorphic molecular crystalline materials endowed with specific molecular-scale lattice and macroscopic dynamics.

19F NMR relaxation of buried tryptophan side chains suggest anisotropic rotational diffusion of the protein RfaH.

The recent application of 19F NMR in the study of biomolecular structure and dynamics has made it a potentially attractive probe to complement traditional 15N/13C labelled probes for backbone and sidechain dynamics, albeit with some complications. The utility of 15N relaxation rates of rigid backbone amide groups to determine the rotational diffusion tensor of proteins is well established. Here we show that the measured 19F relaxation rates of two buried and possibly immobile 19F labelled tryptophan sidechains for the multidomain protein RfaH, in its closed conformation, are in reasonable agreement with the calculated values, only when anisotropic rotational diffusion of the protein is considered. While the sparsity of 19F relaxation data from a limited number of probes precludes the experimental determination of the rotational diffusion tensor here, these results demonstrate the influence of rotational diffusion anisotropy of proteins on 19F NMR relaxation of rigid tryptophan sidechains, while adding to the expanding literature of 19F NMR relaxation data sets in biomolecules.

Dynamics of side chains in poly(quinoxaline-2,3-diyl)s studied via quasielastic neutron scattering.

The side chain dynamics of poly(quinoxaline-2,3-diyl)s (PQXs) are expected to influence their conformation. To investigate these dynamics experimentally, quasielastic neutron scattering (QENS) was performed for PQXs in deuterated tetrahydrofuran (THF-d8) and deuterated 1,1,2-trichloroethane/THF (1,1,2-TCE-d3/THF-d8), in which they formed right-handed and left-handed helical structures, respectively. The mean-square displacement of the PQX side chains in 1,1,2-TCE-d3/THF-d8 was lower than that in THF-d8. Furthermore, QENS complementary studies and molecular dynamics simulations unraveled a coupling between the main-chain and side chain dynamics of PQXs, suggesting the possibility of controlling the main-chain helical chirality through the dynamics of chiral side chains. These insights present a novel strategy for the design of synthetic helical macromolecules with precise chirality control.

Dynamic protein-protein interactions of KCNQ1 and KCNE1 measured by EPR line shape analysis

KCNQ1, also known as Kv7.1, is a voltage gated potassium channel that associates with the KCNE protein family. Mutations in this protein has been found to cause a variety of diseases including Long QT syndrome, a type of cardiac arrhythmia where the QT interval observed on an electrocardiogram is longer than normal. This condition is often aggravated during strenuous exercise and can cause fainting spells or sudden death. KCNE1 is an ancillary protein that interacts with KCNQ1 in the membrane at varying molar ratios. This interaction allows for the flow of potassium ions to be modulated to facilitate repolarization of the heart. The interaction between these two proteins has been studied previously with cysteine crosslinking and electrophysiology. In this study, electron paramagnetic resonance (EPR) spectroscopy line shape analysis in tandem with site directed spin labeling (SDSL) was used to observe changes in side chain dynamics as KCNE1 interacts with KCNQ1. KCNE1 was labeled at different sites that were found to interact with KCNQ1 based on previous literature, along with sites outside of that range as a control. Once labeled KCNE1 was incorporated into vesicles, KCNQ1 (helices S1-S6) was titrated into the vesicles. The line shape differences observed upon addition of KCNQ1 are indicative of an interaction between the two proteins. This method provides a first look at the interactions between KCNE1 and KCNQ1 from a dynamics perspective using the full transmembrane portion of KCNQ1.

Dynamics of Polyalkylfluorene Conjugated Polymers: Insights from Neutron Spectroscopy and Molecular Dynamics Simulations.

The dynamics of the conjugated polymers poly(9,9-dioctylfluorene) (PF8) and poly(9,9-didodecylfluorene) (PF12), differing by the length of their side chains, is investigated in the amorphous phase using the temperature-dependent quasielastic neutron scattering (QENS) technique. The neutron spectroscopy measurements are synergistically underpinned by molecular dynamics (MD) simulations. The probe is focused on the picosecond time scale, where the structural dynamics of both PF8 and PF12 would mainly be dominated by the motions of their side chains. The measurements highlighted temperature-induced dynamics, reflected in the broadening of the QENS spectra upon heating. The MD simulations reproduced well the observations; hence, the neutron measurements validate the MD force fields, the adopted amorphous model structures, and the numerical procedure. As the QENS spectra are dominated by the signal from the hydrogens on the backbones and side chains of PF8 and PF12, extensive analysis of the MD simulations allowed the following: (i) tagging these hydrogens, (ii) estimating their contributions to the self-part of the van Hove functions and hence to the QENS spectra, and (iii) determining the activation energies of the different motions involving the tagged hydrogens. PF12 is found to exhibit QENS spectra broader than those of PF8, indicating a more pronounced motion of the didodecyl chains of PF12 as compared to dioctyl chains of PF8. This is in agreement with the outcome of our MD analysis: (i) confirming a lower glass transition temperature of PF12 compared to PF8, (ii) showing PF12 having a lower density than PF8, and (iii) highlighting lower activation energies of the motions of PF12 in comparison with PF8. This study helped to gain insights into the temperature-induced side-chain dynamics of the PF8 and PF12 conjugated polymers, influencing their stability, which could potentially impact, on the practical side, the performance of the associated optoelectronic active layer.

Analytical Framework to Understand the Origins of Methyl Side-Chain Dynamics in Protein Assemblies.

Side-chain motions play an important role in understanding protein structure, dynamics, protein-protein, and protein-ligand interactions. However, our understanding of protein side-chain dynamics is currently limited by the lack of analytical tools. Here, we present a novel analytical framework employing experimental nuclear magnetic resonance (NMR) relaxation measurements at atomic resolution combined with molecular dynamics (MD) simulation to characterize with a high level of detail the methyl side-chain dynamics in insoluble protein assemblies, using amyloid fibrils formed by the prion HET-s. We use MD simulation to interpret experimental results, where rotameric hops, including methyl group rotation and χ1/χ2 rotations, cannot be completely described with a single correlation time but rather sample a broad distribution of correlation times, resulting from continuously changing local structure in the fibril. Backbone motion similarly samples a broad range of correlation times, from ∼100 ps to μs, although resulting from mostly different dynamic processes; nonetheless, we find that the backbone is not fully decoupled from the side-chain motion, where changes in side-chain dynamics influence backbone motion and vice versa. While the complexity of side-chain motion in protein assemblies makes it very challenging to obtain perfect agreement between experiment and simulation, our analytical framework improves the interpretation of experimental dynamics measurements for complex protein assemblies.

Minimal Peptoid Dynamics Inform Self-Assembly Propensity.

Peptoids are structural isomers of natural peptides, with side chain attachment at the amide nitrogen, conferring this class of compounds with the ability to access both cis and trans ω torsions as well as an increased diversity of ψ/φ states with respect to peptides. Sampling within these dimensions is controlled through side chain selection, and an expansive set of viable peptoid residues exists. It has been shown recently that "minimal" di- and tripeptoids with aromatic side chains can self-assemble into highly ordered structures, with size and morphological definition varying as a function of sequence pattern (e.g., XFF and FXF, where X = a nonaromatic peptoid monomer). Aromatic groups, such as phenylalanine, are regularly used in the design of minimal peptide assemblers. In recognition of this, and to draw parallels between these compounds classes, we have developed a series of descriptors for intramolecular dynamics of aromatic side chains to discern whether these dynamics, in a preassembly condition, can be related to experimentally observed nanoscale assemblies. To do this, we have built on the atomistic peptoid force field reported by Weiser and Santiso (CGenFF-WS) through the rigorous fitting of partial charges and the collation of Charmm General Force Field (CGenFF) parameters relevant to these systems. Our study finds that the intramolecular dynamics of side chains, for a given sequence, is dependent on the specific combination of backbone ω torsions and that homogeneity of sampling across these states correlates well with the experimentally observed ability to assemble into nanomorphologies with long-range order. Sequence patterning is also shown to affect sampling, in a manner consistent for both tripeptoids and tripeptides. Additionally, sampling similarities between the nanofiber forming tripeptoid, Nf-Nke-Nf in the cc state, and the nanotube forming dipeptide FF, highlight a structural motif which may be relevant to the emergence of extended linear assemblies. To assess these properties, a variety of computational approaches have been employed.

Side-chain dynamics of the α1B -adrenergic receptor determined by NMR via methyl relaxation.

G protein-coupled receptors (GPCRs) are medically important membrane proteins that sample inactive, intermediate, and active conformational states characterized by relatively slow interconversions (~μs-ms). On a faster timescale (~ps-ns), the conformational landscape of GPCRs is governed by the rapid dynamics of amino acid side chains. Such dynamics are essential for protein functions such as ligand recognition and allostery. Unfortunately, technical challenges have almost entirely precluded the study of side-chain dynamics for GPCRs. Here, we investigate the rapid side-chain dynamics of a thermostabilized α1B -adrenergic receptor (α1B -AR) as probed by methyl relaxation. We determined order parameters for Ile, Leu, and Val methyl groups in the presence of inverse agonists that bind orthosterically (prazosin, tamsulosin) or allosterically (conopeptide ρ-TIA). Despite the differences in the ligands, the receptor's overall side-chain dynamics are very similar, including those of the apo form. However, ρ-TIA increases the flexibility of Ile1764×56 and possibly of Ile2145×49 , adjacent to Pro2155×50 of the highly conserved P5×50 I3×40 F6×44 motif crucial for receptor activation, suggesting differences in the mechanisms for orthosteric and allosteric receptor inactivation. Overall, increased Ile side-chain rigidity was found for residues closer to the center of the membrane bilayer, correlating with denser packing and lower protein surface exposure. In contrast to two microbial membrane proteins, in α1B -AR Leu exhibited higher flexibility than Ile side chains on average, correlating with the presence of Leu in less densely packed areas and with higher protein-surface exposure than Ile. Our findings demonstrate the feasibility of studying receptor-wide side-chain dynamics in GPCRs to gain functional insights.

Synthesis of a 13C-methylene-labeled isoleucine precursor as a useful tool for studying protein side-chain interactions and dynamics

In this study, we present the synthesis and incorporation of a metabolic isoleucine precursor compound for selective methylene labeling. The utility of this novel α-ketoacid isotopologue is shown by incorporation into the protein Brd4-BD1, which regulates gene expression by binding to acetylated histones. High quality single quantum 13C−1 H-HSQC were obtained, as well as triple quantum HTQC spectra, which are superior in terms of significantly increased 13C-T2 times. Additionally, large chemical shift perturbations upon ligand binding were observed. Our study thus proves the great sensitivity of this precursor as a reporter for side-chain dynamic studies and for investigations of CH-π interactions in protein-ligand complexes.

Beta Amyloid Oligomers with Higher Cytotoxicity have Higher Sidechain Dynamics.

The underlying biophysical principle governing the cytotoxicity of the oligomeric aggregates of β-amyloid (Aβ) peptides has long been an enigma. Here we show that the size of Aβ40 oligomers can be actively controlled by incubating the peptides in reverse micelles. Our approach allowed for the first time a detailed comparison of the structures and dynamics of two Aβ40 oligomers of different sizes, viz., 10 and 23 nm, by solid-state NMR. From the chemical shift data, we infer that the conformation and/or the chemical environments of the residues from K16 to K28 are different between the 10-nm and 23-nm oligomers. We find that the 10-nm oligomers are more cytotoxic, and the molecular motion of the sidechain of its charged residue K16 is more dynamic. Interestingly, the residue A21 exhibits unusually high structural rigidity. Our data raise an interesting possibility that the cytotoxicity of Aβ40 oligomers could also be correlated to the motional dynamics of the sidechains.

Comprehensive analysis of relaxation decays from high-resolution relaxometry

Relaxometry consists in measuring relaxation rates over orders of magnitude of magnetic fields to probe motions of complex systems. High-resolution relaxometry (HRR) experiments can be performed on conventional high-field NMR magnets equipped with a sample shuttle. During the experiment, the sample shuttle transfers the sample between the high-field magnetic center and a chosen position in the stray field for relaxation during a variable delay, thus using the stray field as a variable field. As the relaxation delay occurs outside of the probe, HRR experiments cannot rely on the control of cross-relaxation pathways, which is standard in high-field relaxation pulse sequences. Thus, decay rates are not pure relaxation rates, which may impair a reliable description of the dynamics. Previously, we took into account cross-relaxation effects in the analysis of high-resolution relaxometry data by applying a correction factor to relaxometry decay rates in order to estimate relaxation rates. These correction factors were obtained from the iterative simulation of the relaxation decay while the sample lies outside of the probe and a preceding analysis of relaxation rates which relies on the approximation of a priori multi-exponential decays by mono-exponential functions. However, an analysis protocol matching directly experimental and simulated relaxometry decays should be more self consistent and more generally applicable as it can accommodate deviations from mono-exponential decays. Here, we introduce Matching INtensities for the Optimization of Timescales and Amplitudes of motions Under Relaxometry (MINOTAUR), a framework for the analysis of high-resolution relaxometry that takes as input the intensity decays at all fields. This approach uses the full relaxation matrix to calculate intensity decays, allowing complex relaxation pathways to be taken into account. Therefore, it eliminates the need for a correction of decay rates and for fitting multi-exponential decays with mono-exponential functions. The MINOTAUR software is designed as a flexible framework where relaxation matrices and spectral density functions corresponding to various models of motions can be defined on a case-by-case basis. The agreement with our previous analyses of protein side-chain dynamics from carbon-13 relaxation is excellent, while providing a more robust analysis tool. We expect MINOTAUR to become the tool of choice for the analysis of high-resolution relaxometry.

Membrane Lipid Composition Influences the Hydration of Proton Half-Channels in FoF1-ATP Synthase.

The membrane lipid composition plays an important role in the regulation of membrane protein activity. To probe its influence on proton half-channels' structure in FoF1-ATP synthase, we performed molecular dynamics simulations with the bacterial protein complex (PDB ID: 6VWK) embedded in three types of membranes: a model POPC, a lipid bilayer containing 25% (in vivo), and 75% (bacterial stress) of cardiolipin (CL). The structure proved to be stable regardless of the lipid composition. The presence of CL increased the hydration of half-channels. The merging of two water cavities at the inlet half-channel entrance and a long continuous chain of water molecules directly to cAsp61 from the periplasm were observed. Minor conformational changes in half-channels with the addition of CL caused extremely rare direct transitions between aGlu219-aAsp119, aGlu219-aHis245, and aGln252-cAsp61. Deeper penetration of water molecules (W1-W3) also increased the proton transport continuity. Stable spatial positions of significant amino acid (AA) residue aAsn214 were found under all simulation conditions indicate a prevailing influence of AA-AA or AA-W interactions on the side-chain dynamics. These results allowed us to put forward a model of the proton movement in ATP synthases under conditions close to in vivo and to evaluate the importance of membrane composition in simulations.

Spectroscopically Orthogonal Spin Labels in Structural Biology at Physiological Temperatures.

Electron paramagnetic resonance spectroscopy (EPR) is mostly used in structural biology in conjunction with pulsed dipolar spectroscopy (PDS) methods to monitor interspin distances in biomacromolecules at cryogenic temperatures both in vitro and in cells. In this context, spectroscopically orthogonal spin labels were shown to increase the information content that can be gained per sample. Here, we exploit the characteristic properties of gadolinium and nitroxide spin labels at physiological temperatures to study side chain dynamics via continuous wave (cw) EPR at X band, surface water dynamics via Overhauser dynamic nuclear polarization at X band and short-range distances via cw EPR at high fields. The presented approaches further increase the accessible information content on biomolecules tagged with orthogonal labels providing insights into molecular interactions and dynamic equilibria that are only revealed under physiological conditions.

Determining the Predominant Conformations of Mortiamides A-D in Solution Using NMR Data and Molecular Modeling Tools.

Macrocyclic peptidomimetics have been seriously contributing to our arsenal of drugs to combat diseases. The search for nature's discoveries led us to mortiamides A-D (found in a novel fungus from Northern Canada), which is a family of cyclic peptides that clearly have demonstrated impressive pharmaceutical potential. This prompted us to learn more about their solution-state properties as these are central for binding to target molecules. Here, we secured and isolated mortiamide D, and then acquired high-resolution nuclear magnetic resonance (NMR) data to learn more about its structure and dynamics attributes. Sets of two-dimensional NMR experiments provided atomic-level (through-bond and through-space) data to confirm the primary structure, and NMR-driven molecular dynamics (MD) simulations suggested that more than one predominant three-dimensional (3D) structure exist in solution. Further steps of MD simulations are consistent with the finding that the backbones of mortiamides A-C also have at least two prominent macrocyclic shapes, but the side-chain structures and dynamics differed significantly. Knowledge of these solution properties can be exploited for drug design and discovery.

Cover Image

Protein ScienceVolume 32, Issue 6 Cover ImageFree Access Cover Image First published: 12 May 2023 https://doi.org/10.1002/pro.4652AboutPDF ToolsRequest permissionExport citationAdd to favoritesTrack citation ShareShare Give accessShare full text accessShare full-text accessPlease review our Terms and Conditions of Use and check box below to share full-text version of article.I have read and accept the Wiley Online Library Terms and Conditions of UseShareable LinkUse the link below to share a full-text version of this article with your friends and colleagues. Learn more.Copy URL Graphical Abstract Changes in the side-chain dynamics of two residues in TEM-1 β-lactamase, Tyr 105 and Arg 244, upon binding to FTA, an allosteric inhibitor. Janin plots, akin to Ramachandran plots, compare χ1 vs. χ2 side-chain dihedral angles. Upper panels, FTA binding substantially impacts the dynamics of Tyr 105, which lines the catalytic pocket and is essential for recognizing orthosteric ligands. Lower panels, Arg 244 also undergoes notable changes in side-chain conformational sampling upon FTA binding. For details, see: Erich Hellemann et al., Protein Science, 2023;32(4):e4622. https://doi.org/10.1002/pro.4622 Volume32, Issue6June 2023 RelatedInformation

Regulating polystyrene glass transition temperature by varying the hydration levels of aromatic ring/Li+ interaction.

Polymer properties can be altered via lithium ion doping, whereby adsorbed Li+ binds with H2O within the polymer chain. However, direct spectroscopic evidence of the tightness of Li+/H2O binding in the solid state is limited, and the impact of Li+ on polymer sidechain packing is rarely reported. Here, we investigate a polystyrene/H2O/LiCl system using solid-state NMR, from which we determined a dipolar coupling of 11.4 kHz between adsorbed Li+ and H2O protons. This coupling corroborates a model whereby Li+ interacts with the oxygen atom in H2O via charge affinity, which we believe is the main driving force of Li+ binding. We demonstrated the impact of hydrated Li+ on sidechain packing and dynamics in polystyrene using proton-detected solid-state NMR. Experimental data and density functional theory (DFT) simulations revealed that the addition of Li+ and the increase in the hydration levels of Li+, coupled with aromatic ring binding, change the energy barrier of sidechain packing and dynamics and, consequently, changes the glass transition temperature of polystyrene.

ARCHE-NOAH: NMR supersequence with five different CEST experiments for studying protein conformational dynamics.

An NMR NOAH-supersequence is presented consisting of five CEST experiments for studying protein backbone and side-chain dynamics by 15N-CEST, carbonyl-13CO-CEST, aromatic-13Car-CEST, 13Cα-CEST, and methyl-13Cmet-CEST. The new sequence acquires the data for these experiments in a fraction of the time required for the individual experiments, saving over four days of NMR time per sample.

Ligands selectively tune the local and global motions of neurotensin receptor 1 (NTS1)

Nuclear magnetic resonance (NMR) studies have revealed that fast methyl sidechain dynamics can report on entropically-driven allostery. Yet, NMR applications have been largely limited to the super-microsecond motional regimes of G protein-coupled receptors (GPCRs). We use13Cε-methionine chemical shift-based global order parameters to test if ligands affect the fast dynamics of a thermostabilized GPCR, neurotensin receptor 1 (NTS1). We establish that the NTS1 solution ensemble includes substates with lifetimes on several, discrete timescales. The longest-lived states reflect those captured in agonist- and inverse agonist-bound crystal structures, separated by large energy barriers. We observe that the rapid fluctuations of individual methionine residues, superimposed on these long-lived states, respond collectively with the degree of fast, global dynamics correlating with ligand pharmacology. This approach lends confidence to interpreting spectra in terms of local structure and methyl dihedral angle geometry. The results suggest a role for sub-microsecond dynamics and conformational entropy in GPCR ligand discrimination.

Pressure, motion, and conformational entropy in molecular recognition by proteins

The thermodynamics of molecular recognition by proteins is a central determinant of complex biochemistry. For over a half-century, detailed cryogenic structures have provided deep insight into the energetic contributions to ligand binding by proteins. More recently, a dynamical proxy based on NMR-relaxation methods has revealed an unexpected richness in the contributions of conformational entropy to the thermodynamics of ligand binding. Here, we report the pressure dependence of fast internal motion within the ribonuclease barnase and its complex with the protein barstar. In what we believe is a first example, we find that protein dynamics are conserved along the pressure-binding thermodynamic cycle. The femtomolar affinity of the barnase-barstar complex exists despite a penalty by −TΔSconf of +11.7 kJ/mol at ambient pressure. At high pressure, however, the overall change in side-chain dynamics is zero, and binding occurs with no conformational entropy penalty, suggesting an important role of conformational dynamics in the adaptation of protein function to extreme environments. Distinctive clustering of the pressure sensitivity is observed in response to both pressure and binding, indicating the presence of conformational heterogeneity involving less efficiently packed alternative conformation(s). The structural segregation of dynamics observed in barnase is striking and shows how changes in both the magnitude and the sign of regional contributions of conformational entropy to the thermodynamics of protein function are possible.

Aromatic ring flips in differently packed ubiquitin protein crystals from MAS NMR and MD

Probing the dynamics of aromatic side chains provides important insights into the behavior of a protein because flips of aromatic rings in a protein’s hydrophobic core report on breathing motion involving a large part of the protein. Inherently invisible to crystallography, aromatic motions have been primarily studied by solution NMR. The question how packing of proteins in crystals affects ring flips has, thus, remained largely unexplored. Here we apply magic-angle spinning NMR, advanced phenylalanine 1H-13C/2H isotope labeling and MD simulation to a protein in three different crystal packing environments to shed light onto possible impact of packing on ring flips. The flips of the two Phe residues in ubiquitin, both surface exposed, appear remarkably conserved in the different crystal forms, even though the intermolecular packing is quite different: Phe4 flips on a ca. 10–20 ns time scale, and Phe45 are broadened in all crystals, presumably due to µs motion. Our findings suggest that intramolecular influences are more important for ring flips than intermolecular (packing) effects.