Conformational Dynamics of apo‐GlnBP Revealed by Experimental and Computational Analysis

The glutamine binding protein (GlnBP) binds l ‐glutamine and cooperates with its cognate transporters during glutamine uptake. Crystal structure analysis has revealed an open and a closed conformation for apo‐ and holo‐GlnBP, respectively. However, the detailed conformational dynamics have remained unclear. Herein, we combined NMR spectroscopy, MD simulations, and single‐molecule FRET techniques to decipher the conformational dynamics of apo‐GlnBP. The NMR residual dipolar couplings of apo‐GlnBP were in good agreement with a MD‐derived structure ensemble consisting of four metastable states. The open and closed conformations are the two major states. This four‐state model was further validated by smFRET experiments and suggests the conformational selection mechanism in ligand recognition of GlnBP.     

Influence of an Unnatural Amino Acid Side Chain on the Conformational Dynamics of Peptides

In this work, a non‐natural amino acid, H‐propargylglycine‐OH (Pra), is chosen to examine the side‐chain effect on the backbone conformation of small peptides. The conformations of two synthesized Pra‐containing tripeptides, Ac‐Pra‐Pra‐NH₂ (PPTP) and Ac‐Pra‐Ala‐NH₂ (PATP), are examined by infrared (IR) spectroscopy in combination with molecular dynamics (MD) simulations and quantum chemical computations. By analyzing the joint distributions of backbone torsional angles, several significant conformations can be identified for the two tripeptides solvated in D₂O. At room temperature, 44 % of PPTP exists in the α‐α conformation and 33 % of PATP exists in the α‐polyproline‐II conformation. Larger structural inhomogeneity is seen in both cases by MD simulations at elevated temperatures. Thus even a small side chain, such as the propargyl group can significantly alter the peptide backbone conformations. The results suggest that there is no overwhelming conformational propensity of the Pra residue in short peptides. IR spectra simulated in the amide‐I region using two different methods, reasonably reproduce the experimental IR spectra and their temperature dependence.     

Conformational Flexibility in the Transmembrane Protein TSPO

The translocator protein (TSPO) is an integral membrane protein that interacts with a wide variety of endogenous ligands, such as cholesterol and porphyrins, and is also the target for several small molecules with substantial in vivo efficacy. When complexed with the TSPO‐specific radioligand (R)‐PK11195, TSPO folds into a rigid five‐helix bundle. However, little is known about the structure and dynamics of TSPO in the absence of high‐affinity ligands. By means of NMR spectroscopy, we show that TSPO exchanges between multiple conformations in the absence of (R)‐PK11195. Extensive motions on time scales from pico‐ to microseconds occur all along the primary sequence of the protein, leading to a loss of stable tertiary interactions and local unfolding of the helical structure in the vicinity of the ligand‐binding site. The flexible nature of TSPO highlights the importance of conformational plasticity in integral membrane proteins.     

Protein Conformational Flexibility from Structure‐Free Analysis of NMR Dipolar Couplings: Quantitative and Absolute Determination of Backbone Motion in Ubiquitin

A robust procedure for the determination of protein‐backbone motions on time scales of pico‐ to milliseconds directly from residual dipolar couplings has been developed that requires no additional scaling relative to external references. The results for ubiquitin (blue in graph: experimental N? HN order parameters) correspond closely to the amplitude, nature, and distribution of motion found in a 400 ns molecular‐dynamics trajectory of ubiquitin (red).

     

Increasing the Affinity of an O-Antigen Polysaccharide Binding Site in Shigella flexneri Bacteriophage Sf6 Tailspike Protein

Broad and unspecific use of antibiotics accelerates spread of resistances. Sensitive and robust pathogen detection is thus important for a more targeted application. Bacteriophages contain a large repertoire of pathogen-binding proteins. These tailspike proteins (TSP) often bind surface glycans and represent a promising design platform for specific pathogen sensors. We analysed bacteriophage Sf6 TSP that recognizes the O-polysaccharide of dysentery-causing Shigella flexneri to develop variants with increased sensitivity for sensor applications. Ligand polyrhamnose backbone conformations were obtained from 2D ¹H,¹H-trNOESY NMR utilizing methine–methine and methine–methyl correlations. They agreed well with conformations obtained from molecular dynamics (MD), validating the method for further predictions. In a set of mutants, MD predicted ligand flexibilities that were in good correlation with binding strength as confirmed on immobilized S. flexneri O-polysaccharide (PS) with surface plasmon resonance. In silico approaches combined with rapid screening on PS surfaces hence provide valuable strategies for TSP-based pathogen sensor design.     

A Non‐damaging Method to Analyze the Configuration and Dynamics of Nitrotyrosines in Proteins

Often, deregulation of protein activity and turnover by tyrosine nitration drives cells toward pathogenesis. Hence, understanding how the nitration of a protein affects both its function and stability is of outstanding interest. Nowadays, most of the in vitro analyses of nitrated proteins rely on chemical treatment of native proteins with an excess of a chemical reagent. One such reagent, peroxynitrite, stands out for its biological relevance. However, given the excess of the nitrating reagent, the resulting in vitro modification could differ from the physiological nitration. Here, we determine unequivocally the configuration of distinct nitrated‐tyrosine rings in single‐tyrosine mutants of cytochrome c. We aimed to confirm the nitration position by a non‐destructive method. Thus, we have resorted to ¹H‐¹⁵N heteronuclear single quantum coherence(HSQC) spectra to identify the ³J(N? H) correlation between a ¹⁵N‐tagged nitro group and the adjacent aromatic proton. Once the chemical shift of this proton was determined, we compared the ¹H‐¹³C HSQC spectra of untreated and nitrated samples. All tyrosines were nitrated at ε positions, in agreement to previous analysis by indirect techniques. Notably, the various nitrotyrosine residues show a different dynamic behaviour that is consistent with molecular dynamics computations.     

Site‐Resolved Backbone and Side‐Chain Intermediate Dynamics in a Carbohydrate‐Binding Module Protein Studied by Magic‐Angle Spinning NMR Spectroscopy

Magic‐angle spinning solid‐state NMR spectroscopy has been applied to study the dynamics of CBM3b–Cbh9A from Clostridium thermocellum (ctCBM3b), a cellulose binding module protein. This 146‐residue protein has a nine‐stranded β‐sandwich fold, in which 35 % of the residues are in the β‐sheet and the remainder are composed of loops and turns. Dynamically averaged ¹H‐¹³C dipolar coupling order parameters were extracted in a site‐specific manner by using a pseudo‐three‐dimensional constant‐time recoupled separated‐local‐field experiment (dipolar‐chemical shift correlation experiment; DIPSHIFT). The backbone‐Cα and Cβ order parameters indicate that the majority of the protein, including turns, is rigid despite having a high content of loops; this suggests that restricted motions of the turns stabilize the loops and create a rigid structure. Water molecules, located in the crystalline interface between protein units, induce an increased dynamics of the interface residues thereby lubricating crystal water‐mediated contacts, whereas other crystal contacts remain rigid.     

Conformational Space of a Polyphilic Molecule with a Fluorophilic Side Chain Integrated in a DPPC Bilayer

We investigate the conformational space of a polyphilic molecule with hydrophilic, lipophilic and fluorophilic parts inserted as a transmembrane agent into a dipalmitoylphosphatidylcholine bilayer by means of all‐atom molecular dynamics simulations. Special focus is put on the competing structural driving forces arising from the hydrophilic, lipophilic and fluorophilic side chains and the aromatic backbone of the polyphile. We observe a significant difference between the lipophilic and the fluorophilic side chains regarding their intramembrane distribution. While the lipophilic groups remain membrane‐centered, the fluorophilic parts tend to orient toward the phosphate headgroups. This trend is important for understanding the influence of polyphile agents on the properties of phospholipid membranes. From a fundamental point of view, our computed distribution functions of the side chains are related to the interplay of sterical, enthalpic and entropic driving forces. Our findings illustrate the potential of rationally designed membrane additives which can be exploited to tune the properties of phospholipid membranes. © 2017 Wiley Periodicals, Inc.     

Conformational selection of the AGA*IA(M) heparin pentasaccharide when bound to the fibroblast growth factor receptor

The interaction of the synthetic pentasaccharide AGA*IA(M) (GlcNS,6S-GlcA-GlcNS,3S,6S-IdoA2S-GlcNS,6S-Me) with the extracellular Ig2 domain of the fibroblast growth factor receptor (FGFR2) has been studied by NMR and computational methods. Analysis of the heparin pentasaccharide in the free state and in the complex indicates the existence of a conformational selection process. Although an equilibrium exists between the (1)C(4) and (2)S(0) conformers (ratio 60:40) of the 2-O-sulfo-α-L-iduronate ring (IdoA2S) in the free state, FGFR2 selects only the unique twisted-boat (2)S(0) conformation of this IdoA2S residue. In addition, the protein residues involved in the binding with AGA*IA(M) have also been characterized. The NMR results obtained, from both the ligand and protein perspective, were employed to model the bound conformation of the pentasaccharide by a combined docking and molecular dynamic simulation approach.     

Conformational Flexibility and Dynamics of Two (1→6)‐Linked Disaccharides Related to an Oligosaccharide Epitope Expressed on Malignant Tumour Cells

The conformational flexibility and dynamics of two (1→6)‐linked disaccharides that are related to the action of the glycosyl transferase GnT‐V have been investigated. NMR NOE and T‐ROE spectroscopy experiments, conformation‐dependent coupling constants and molecular dynamics (MD) simulations were used in the analyses. To facilitate these studies, the compounds were synthesised as α‐d‐ [6‐¹³C]‐Manp‐OMe derivatives, which reduced the ¹H NMR spectral overlap and facilitated the determination of two‐ and three‐bond ¹H,¹H, ¹H,¹³C and ¹³C,¹³C‐coupling constants. The population distribution for the glycosidic ω torsion angle in α‐d‐ Manp‐(1→6)‐α‐d‐ Manp‐OMe for gt/gg/tg was equal to 45:50:5, whereas in α‐d‐ Manp‐OMe it was determined to be 56:36:8. The dynamic model that was generated for β‐d‐ GlcpNAc‐(1→6)‐α‐d‐ Manp‐OMe by MD simulations employing the PARM22/SU01 CHARMM‐based force field was in very good agreement with experimental observations. β‐d‐ GlcpNAc‐(1→6)‐α‐d‐ Manp‐OMe is described by an equilibrium of populated states in which the ? torsion angle has the exo‐anomeric conformation, the ψ torsion angle an extended antiperiplanar conformation and the ω torsion angle a distribution of populations predominantly between the gauche–trans and the gauche–gauche conformational states (i.e., gt/gg/tg) is equal to 60:35:5, respectively. The use of site‐specific ¹³C labelling in these disaccharides leads to increased spectral dispersion, thereby making NMR spectroscopy based conformational analysis possible that otherwise might be difficult to attain.     

Characterization of Doubly Ionic Hydrogen Bonds in Protic Ionic Liquids by NMR Deuteron Quadrupole Coupling Constants: Differences to H‐bonds in Amides,Peptides, and Proteins

We present the first deuteron quadrupole coupling constants (DQCCs) for selected protic ionic liquids (PILs) measured by solid‐state NMR spectroscopy. The experimental data are supported by dispersion‐corrected density functional theory (DFT‐D3) calculations and molecular dynamics (MD) simulations. The DQCCs of the N−D bond in the triethylammonium cations are the lowest reported for deuterons in PILs, indicating strong hydrogen bonds between ions. The NMR coupling parameters are compared to those in amides, peptides, and proteins. The DQCCs show characteristic behavior with increasing interaction strength of the counterion and variation of the H‐bond motifs. We report the similar presence of the quadrupolar splitting pattern and the narrow liquid line in the NMR spectra over large temperature ranges, indicating the heterogeneous nature of PILs.     

Effect of the Substituents of the Neighboring Ring in the Conformational Equilibrium of Iduronate in Heparin‐like Trisaccharides

Based on the structure of the regular heparin, we have prepared a smart library of heparin‐like trisaccharides by incorporating some sulfate groups in the sequence α‐D ‐GlcNS‐ (1‐4)‐α‐L ‐Ido2S‐(1‐4)‐α‐D ‐GlcN. According to the 3D structure of heparin, which features one helix turn every four residues, this fragment corresponds to the minimum binding motif. We have performed a complete NMR study and found that the trisaccharides have a similar 3D structure to regular heparin itself, but their spectral properties are such that allow to extract very detailed information about distances and coupling constants as they are isotropic molecules. The characteristic conformational equilibrium of the central iduronate ring has been analyzed combining NMR and molecular dynamics and the populations of the conformers of the central iduronate ring have been calculated. We have found that in those compounds lacking the sulfate group at position 6 of the reducing end glucosamine, the population of ²S₀ of the central iduronate residue is sensitive to the temperature decreasing to 19 % at 278 K. On the contrary, the trisaccharides with 6‐O‐sulfate in the reducing end glucosamine keep the level of population constant with temperature circa 40 % of ²S₀ similar to that observed at room temperature. Another structural feature that has been revealed through this analysis is the larger flexibility of the L ‐IdoAS‐ D ‐GlcN glycosidic linkage, compared with the D ‐GlcNS‐L ‐IdoA. We propose that this is the point where the heparin chain is bended to form structures far from the regular helix known as kink that have been proposed to play an important role in the specificity of the heparin–protein interaction.     

Discovery of Potential Antiviral Compounds against Hendra Virus by Targeting Its Receptor-Binding Protein (G) Using Computational Approaches

Hendra virus (HeV) belongs to the paramyxoviridae family of viruses which is associated with the respiratory distress, neurological illness, and potential fatality of the affected individuals. So far, no competitive approved therapeutic substance is available for HeV. For that reason, the current research work was conducted to propose some novel compounds, by adopting a Computer Aided Drug Discovery approach, which could be used to combat HeV. The G attachment Glycoprotein (Ggp) of HeV was selected to achieve the primary objective of this study, as this protein makes the entry of HeV possible in the host cells. Briefly, a library of 6000 antiviral compounds was screened for potential drug-like properties, followed by the molecular docking of short-listed compounds with the Protein Data Bank (PDB) structure of Ggp. Docked complexes of top two hits, having maximum binding affinities with the active sites of Ggp, were further considered for molecular dynamic simulations of 200 ns to elucidate the results of molecular docking analysis. MD simulations and Molecular Mechanics Energies combined with the Generalized Born and Surface Area (MMGBSA) or Poisson–Boltzmann and Surface Area (MMPBSA) revealed that both docked complexes are stable in nature. Furthermore, the same methodology was used between lead compounds and HeV Ggp in complex with its functional receptor in human, Ephrin-B2. Surprisingly, no major differences were found in the results, which demonstrates that our identified compounds can also perform their action even when the Ggp is attached to the Ephrin-B2 ligand. Therefore, in light of all of these results, we strongly suggest that compounds (S)-5-(benzylcarbamoyl)-1-(2-(4-methyl-2-phenylpiperazin-1-yl)-2-oxoethyl)-6-oxo-3,6-dihydropyridin-1-ium-3-ide and 5-(cyclohexylcarbamoyl)-1-(2-((2-(3-fluorophenyl)-2-methylpropyl)amino)-2-oxoethyl)-6-oxo-3,6-dihydropyridin-1-ium-3-ide could be considered as potential therapeutic agents against HeV; however, further in vitro and in vivo experiments are required to validate this study.     

Conformational Dynamics in the Core of Human Y145Stop Prion Protein Amyloid Probed by Relaxation Dispersion NMR

Microsecond to millisecond timescale backbone dynamics of the amyloid core residues in Y145Stop human prion protein (PrP) fibrils were investigated by using ¹⁵N rotating frame (R_1ρ) relaxation dispersion solid-state nuclear magnetic resonance spectroscopy over a wide range of spin-lock fields. Numerical simulations enabled the experimental relaxation dispersion profiles for most of the fibril core residues to be modelled by using a two-state exchange process with a common exchange rate of 1000 s⁻¹, corresponding to protein backbone motion on the timescale of 1 ms, and an excited-state population of 2 %. We also found that the relaxation dispersion profiles for several amino acids positioned near the edges of the most structured regions of the amyloid core were better modelled by assuming somewhat higher excited-state populations (∼5–15 %) and faster exchange rate constants, corresponding to protein backbone motions on the timescale of ∼100–300 μs. The slow backbone dynamics of the core residues were evaluated in the context of the structural model of human Y145Stop PrP amyloid.