Quantitative J correlation: a new approach for measuring homonuclear three-bond J(HNH.alpha.) coupling constants in 15N-enriched proteins

Electronic Control of Polyproline II Helix Stability via the Identity of Acyl Capping Groups: the Pivaloyl Group Particularly Promotes PPII

The type II polyproline helix (PPII) is a fundamental secondary structure of proteins, important in globular proteins, in intrinsically disordered proteins, and at protein-protein interfaces. PPII is stabilized in part by n→π* interactions between consecutive carbonyls, via electron delocalization between an electron-donor carbonyl lone pair (n) and an electron-acceptor carbonyl (π*) on the subsequent residue. We previously demonstrated that changes to the electronic properties of the acyl donor can predictably modulate the strength of n→π* interactions, with data from model compounds, in solution in chloroform, in the solid state, and computationally. Herein, we examined whether the electronic properties of acyl capping groups could modulate the stability of PPII in peptides in water. In X-PPGY-NH2 peptides (X=10 acyl capping groups), the effect of acyl group identity on PPII was quantified by circular dichroism and NMR spectroscopy. Electron-rich acyl groups promoted PPII relative to the standard acetyl (Ac-) group, with the pivaloyl and iso-butyryl groups most significantly increasing PPII. In contrast, acyl derivatives with electron-withdrawing substituents and the formyl group relatively disfavored PPII. Similar results, though lesser in magnitude, were also observed in X-APPGY-NH2 peptides, indicating that the capping group can impact PPII conformation at both proline and non-proline residues. The pivaloyl group was particularly favorable in promoting PPII. The effects of acyl capping groups were further analyzed in X-DfpPGY-NH2 and X-ADfpPGY-NH2 peptides, Dfp=4,4-difluoroproline. Data on these peptides indicated that acyl groups induced order Piv- > Ac- > For-. These results suggest that greater consideration should be given to the identity of acyl capping groups in inducing structure in peptides.

A Joint Experimental and Theoretical Study on the Structural and Spectroscopic Properties of the Piv-Pro-d-Ser-NHMe Peptide

In this paper, we investigate the secondary structure of the Piv-Pro-d-Ser-NHMe peptide by means of nuclear magnetic resonance (NMR) and electronic circular dichroism (ECD) experiments, in conjunction with theoretical simulations based on molecular dynamics and time-dependent density functional theory calculations including polarizable embedding to account for solvent effects. The various experimental and theoretical protocols are assessed and validated, and are shown to provide a consistent description of the turn structure adopted by this peptide in solution. In addition, a simple fitting procedure is proposed to make the simulated and experimental ECD almost perfectly match. This full methodology is finally tested on another small peptide, enlightening its efficiency and robustness.

Requirements for efficient endosomal escape by designed mini-proteins

ZF5.3 is a compact, rationally designed mini-protein that escapes efficiently from the endosomes of multiple cell types. Despite its small size (27 amino acids), ZF5.3 can be isolated intact from the cytosol of treated cells and guides multiple classes of proteins into the cytosol and/or nucleus. In the best cases, delivery efficiencies reach or exceed 50% to establish nuclear or cytosolic concentrations of 500 nM or higher. But other than the requirement for unfoldable cargo and an intact HOPS complex, there is little known about how ZF5.3 traverses the limiting endocytic membrane. Here we delineate the attributes of ZF5.3 that enable efficient endosomal escape. We confirm that ZF5.3 is stable at pH values between 5.5 and 7.5, with no evidence of unfolding even at temperatures as high as 95 °C. The high-resolution NMR structure of ZF5.3 at pH 5.5, also reported here, shows a canonical p zinc-finger fold with the penta-arg motif integrated seamlessly into the C-terminal α-helix. At lower pH, ZF5.3 unfolds cooperatively as judged by both circular dichroism and high-resolution NMR. Unfolding occurs upon protonation of a single Zn(II)-binding His side chain whose pKa corresponds almost exactly to that of the late endosomal lumen. pH-induced unfolding is essential for endosomal escape, as a ZF5.3 analog that remains folded at pH 4.5 fails to efficiently reach the cytosol, despite high overall uptake. Finally, using reconstituted liposomes, we identify a high-affinity interaction of ZF5.3 with a specific lipid-BMP-that is selectively enriched in the inner leaflet of late endosomal membranes. This interaction is 10-fold stronger at low pH than neutral pH, providing a molecular picture for why escape occurs preferentially and in a HOPS-dependent manner from late endosomal compartments. The requirements for programmed endosomal escape identified here should aid and inform the design of proteins, peptidomimetics, and other macromolecules that reach cytosolic or nuclear targets intact and at therapeutically relevant concentrations.

Structural Insights into Antibacterial Payload Release from Gold Nanoparticles Bound to E. coli Peptide Deformylase

The lack of new antibiotics and the rapidly rising number of pathogens resistant to antibiotics pose a serious problem to mankind. In bacteria, the cell membrane provides the first line of defence to antibiotics by preventing them from reaching their molecular target. To overcome this entrance barrier, it has been suggested[1] that small Gold-Nanoparticles (AuNP) could possibly function as drug delivery systems for antibiotic ligands. Using actinonin-based ligands, we provide here proof-of-principle of AuNP functionalisation, the capability to bind and inhibit the target protein of the ligand, and the possibility to selectively release the antimicrobial payload. To this end, we successfully synthesised AuNP coated with thio-functionalised actinonin and a derivative. Interactions between 15N-enriched His-peptide deformylase 1-147 from E. coli (His-ecPDF 1-147) and compound-coated AuNP were investigated via 2D 1H-15N-HSQC NMR spectra proving the direct binding to His-ecPDF 1-147. More importantly by adding dithiothreitol (DTT), we show that the derivative is successfully released from AuNPs while still bound to His-ecPDF 1-147. Our findings indicate that AuNP-conjugated ligands can address and bind intracellular target proteins. The system introduced here presents a new delivery platform for antibiotics and allows for the easy optimisation of ligand coated AuNPs.

Improved Protein Dynamics and Hydration in the Martini3 Coarse-Grain Model

The Martini coarse-grain force-field has emerged as an important framework to probe cellular processes at experimentally relevant time- and length-scales. However, the recently developed version, the Martini3 force-field with the implemented Go̅ model (Martini3Go̅), as well as previous variants of the Martini model have not been benchmarked and rigorously tested for globular proteins. In this study, we consider three globular proteins, ubiquitin, lysozyme, and cofilin, and compare protein dynamics and hydration with observables from experiments and all-atom simulations. We show that the Martini3Go̅ model is able to accurately model the structural and dynamic features of small globular proteins. Overall, the structural integrity of the proteins is maintained, as validated by contact maps, radii of gyration (Rg), and SAXS profiles. The chemical shifts predicted from the ensemble sampled in the simulations are consistent with the experimental data. Further, a good match is observed in the protein-water interaction energetics, and the hydration levels of the residues are similar to atomistic simulations. However, the protein-water interaction dynamics is not accurately represented and appears to depend on the protein structural complexity, residue specificity, and water dynamics. Our work is a step toward testing and assessing the Martini3Go̅ model and provides insights into future efforts to refine Martini models with improved solvation effects and better correspondence to the underlying all-atom systems.

Design, synthesis, conformational analysis, and biological activity of Cα1-to-Cα6 1,4- and 4,1-disubstituted 1H-[1,2,3]triazol-1-yl-bridged oxytocin analogues

Oxytocin (OT) is a neurohypophyseal peptide hormone containing a disulphide-bridged pseudocyclic conformation. The biomedical use of OT peptides is limited amongst others by disadvantageous pharmacokinetic parameters. To increase the stability of OT by replacing the disulphide bridge with the stable and more rigid [1,2,3]triazol-1-yl moiety, we employed the Cu2+-catalysed side chain-to-side chain azide-alkyne 1,3-cycloaddition. Here we report the design, synthesis, conformational analysis, and in vitro pharmacological activity of a homologous series of Cα1-to-Cα6 side chain-to-side chain [1,2,3]triazol-1-yl-containing OT analogues differing in the length of the bridge, location, and orientation of the linking moiety. Exploiting this macrocyclisation approach, it was possible to generate a systematic series of compounds providing interesting insight into the structure-conformation-function relationship of OT. Most analogues were able to adopt similar conformation to endogenous OT in water, namely, a type I β-turn. This approach may in the future generate stabilised pharmacological peptide tools to advance understanding of OT physiology.

Phosphorylation motif dictates GPCR C-terminal domain conformation and arrestin interaction

Arrestin-dependent G protein-coupled receptor (GPCR) signaling pathway is regulated by the phosphorylation state of GPCR's C-terminal domain, but the molecular bases of arrestin:receptor interaction are to be further illuminated. Here we investigated the impact of phosphorylation on the conformational features of the C-terminal region from three rhodopsin-like GPCRs, the vasopressin V2 receptor (V2R), the growth hormone secretagogue or ghrelin receptor type 1a (GHSR), and the β2-adernergic receptor (β2AR). Using phosphomimetic variants, we identified pre-formed secondary structure elements, or short linear motifs (SLiMs), that undergo specific conformational transitions upon phosphorylation. Of importance, such conformational transitions appear to favor arrestin-2 binding. Hence, our results suggest a model in which the phosphorylation-dependent structuration of the GPCR C-terminal regions would modulate arrestin binding and therefore signaling outcomes in arrestin-dependent pathways.

An Inherent Difference between Serine and Threonine Phosphorylation: Phosphothreonine Strongly Prefers a Highly Ordered, Compact, Cyclic Conformation

Phosphorylation and dephosphorylation of proteins by kinases and phosphatases are central to cellular responses and function. The structural effects of serine and threonine phosphorylation were examined in peptides and in proteins, by circular dichroism, NMR spectroscopy, bioinformatics analysis of the PDB, small-molecule X-ray crystallography, and computational investigations. Phosphorylation of both serine and threonine residues induces substantial conformational restriction in their physiologically more important dianionic forms. Threonine exhibits a particularly strong disorder-to-order transition upon phosphorylation, with dianionic phosphothreonine preferentially adopting a cyclic conformation with restricted ϕ (ϕ ∼ -60°) stabilized by three noncovalent interactions: a strong intraresidue phosphate-amide hydrogen bond, an n → π* interaction between consecutive carbonyls, and an n → σ* interaction between the phosphate Oγ lone pair and the antibonding orbital of C-Hβ that restricts the χ2 side-chain conformation. Proline is unique among the canonical amino acids for its covalent cyclization on the backbone. Phosphothreonine can mimic proline's backbone cyclization via noncovalent interactions. The preferred torsions of dianionic phosphothreonine are ϕ,ψ = polyproline II helix > α-helix (ϕ ∼ -60°); χ1 = g-; χ2 ∼ +115° (eclipsed C-H/O-P bonds). This structural signature is observed in diverse proteins, including in the activation loops of protein kinases and in protein-protein interactions. In total, these results suggest a structural basis for the differential use and evolution of threonine versus serine phosphorylation sites in proteins, with serine phosphorylation typically inducing smaller, rheostat-like changes, versus threonine phosphorylation promoting larger, step function-like switches, in proteins.

A scorpion toxin affecting sodium channels shows double cis-trans isomerism

Scorpion α-toxins (α-NaTx) inhibiting the inactivation of voltage-gated sodium channels (Nav ) are a well-studied family of small proteins. We previously showed that the structure of α-NaTx specificity module responsible for selective Nav binding is governed by an interplay between the nest and niche protein motifs. Here, we report the solution structure of the toxin Lqq4 from the venom of the scorpion Leiurus quinquestriatus. Unexpectedly, we find that this toxin presents an ensemble of long-lived structurally distinct states. We unequivocally assign these states to the alternative configurations (cis-trans isomers) of two peptide bonds: V56-P57 and C17-G18; neither of the cis isomers has been described in α-NaTx so far. We argue that the native conformational space of α-NaTx is wider than assumed previously.

Lycopene destabilizes preformed Aβ fibrils: Mechanistic insights from all-atom molecular dynamics simulation

The therapeutic strategy employing destabilization of the preformed Aβ fibril by various natural compounds, as studied by experimental and computational methods, has been reported significant in curing Alzheimer's disease (AD). However, lycopene (a carotenoid), from terpenes family, needs investigation for its destabilization potential of Aβ fibril. The highest antioxidant potential and ability to cross blood brain barrier makes lycopene a preferred choice as drug lead for treating AD. The current study focuses on investigating the destabilization potential and underpinning mechanism of lycopene on different polymorphic forms of Aβ fibril via Molecular Dynamics (MD) simulation. The key findings highlight binding of lycopene to the outer surface of the chain F of the fibril (2NAO). Herein G9, K16 and V18 residues were found to be involved in van der Waals with the methyl groups of the lycopene. Additionally, Y10 and F20 residues were observed to interact via π-π interactions with CC bonds of the lycopene. The surface mediated binding of lycopene to the fibril is attributed to the large size and structural rigidity of lycopene along with the bulky size of 2NAO and narrow space of fibrillar cavity. The destabilization of the fibril is evident by breakage of inherent H-bonds and hydrophobic interactions in the presence of one lycopene molecule. The lesser β-sheet content explains disorganization of the fibril and bars the higher order aggregation curbing neurotoxicity of the fibril. The higher concentration of the lycopene is not found to be linearly correlated with the extent of destabilization of the fibril. Lycopene is also observed to destabilize the other polymorphic form of Aβ fibril (2BEG), by accessing the fibrillar cavity and lowering the β-sheet content. The destabilization observed by lycopene on two major polymorphs of Aβ fibril explains its potency towards developing an effective therapeutic approach in treating AD.

Structural and molecular insights into tacrine-benzofuran hybrid induced inhibition of amyloid-β peptide aggregation and BACE1 activity

Amyloid-β (Aβ) aggregation and β-amyloid precursor protein cleaving enzyme 1 (BACE1) are the potential therapeutic drug targets for Alzheimer's disease (AD). A recent study highlighted that tacrine-benzofuran hybrid C1 displayed anti-aggregation activity against Aβ42 peptide and inhibit BACE1 activity. However, the inhibition mechanism of C1 against Aβ42 aggregation and BACE1 activity remains unclear. Thus, molecular dynamics (MD) simulations of Aβ42 monomer and BACE1 with and without C1 were performed to inspect the inhibitory mechanism of C1 against Aβ42 aggregation and BACE1 activity. In addition, a ligand-based virtual screening followed by MD simulations was employed to explore potent new small-molecule dual inhibitors of Aβ42 aggregation and BACE1 activity. MD simulations highlighted that C1 promotes the non aggregating helical conformation in Aβ42 and destabilizes D23-K28 salt bridge that plays a vital role in the self-aggregation of Aβ42. C1 displays a favourable binding free energy (-50.7 ± 7.3 kcal/mol) with Aβ42 monomer and preferentially binds to the central hydrophobic core (CHC) residues. MD simulations highlighted that C1 strongly interacted with the BACE1 active site (Asp32 and Asp228) and active pockets. The scrutiny of interatomic distances among key residues of BACE1 highlighted the close flap (non-active) position in BACE1 on the incorporation of C1. The MD simulations explain the observed high inhibitory activity of C1 against Aβ aggregation and BACE1 in the in vitro studies. The ligand-based virtual screening followed by MD simulations identified CHEMBL2019027 (C2) as a promising dual inhibitor of Aβ42 aggregation and BACE1 activity.Communicated by Ramaswamy H. Sarma.

Simple and Effective Conformational Sampling Strategy for Intrinsically Disordered Proteins Using the UNRES Web Server

Intrinsically disordered proteins (IDPs) contain more charged amino acids than folded proteins, resulting in a lack of hydrophobic core(s) and a tendency to adopt rapidly interconverting structures rather than well-defined structures. The structural heterogeneity of IDPs, encoded by the amino acid sequence, is closely related to their unique roles in biological pathways, which require them to interact with different binding partners. Recently, Robustelli and co-workers have demonstrated that a balanced all-atom force field can be used to sample heterogeneous structures of disordered proteins ( Proc. Natl. Acad. Sci. U.S.A. 2018, 115, E4758-E4766). However, such a solution requires extensive computational resources, such as Anton supercomputers. Here, we propose a simple and effective solution to sample the conformational space of IDPs using a publicly available web server, namely, the UNited-RESidue (UNRES) web server. Our proposed solution requires no investment in computational resources and no prior knowledge of UNRES. UNRES Replica Exchange Molecular Dynamics (REMD) simulations were carried out on a set of eight disordered proteins at temperatures spanning from 270 to 430 K. Utilizing the latest UNRES force field designed for structured proteins, with proper selections of temperatures, we were able to produce comparable results to all-atom force fields as reported in work done by Robustelli and co-workers. In addition, NMR observables and the radius of gyration calculated from UNRES ensembles were directly compared with the experimental data to further evaluate the accuracy of the UNRES model at all temperatures. Our results suggest that carrying out the UNRES simulations at optimal temperatures using the UNRES web server can be a good alternative to sample heterogeneous structures of IDPs.

Synthesis and conformational preferences of peptides and proteins with cysteine sulfonic acid

Cysteine sulfonic acid (Cys-SO₃H; cysteic acid) is an oxidative post-translational modification of cysteine, resulting from further oxidation from cysteine sulfinic acid (Cys-SO₂H). Cysteine sulfonic acid is considered an irreversible post-translational modification, which serves as a biomarker of oxidative stress that has resulted in oxidative damage to proteins. Cysteine sulfonic acid is anionic, as a sulfonate (Cys-SO₃^-; cysteate), in the ionization state that is almost exclusively present at physiological pH (pK_a ∼ -2). In order to understand protein structural changes that can occur upon oxidation to cysteine sulfonic acid, we analyzed its conformational preferences, using experimental methods, bioinformatics, and DFT-based computational analysis. Cysteine sulfonic acid was incorporated into model peptides for α-helix and polyproline II helix (PPII). Within peptides, oxidation of cysteine to the sulfonic acid proceeds rapidly and efficiently at room temperature in solution with methyltrioxorhenium (MeReO₃) and H₂O₂. Peptides containing cysteine sulfonic acid were also generated on solid phase using trityl-protected cysteine and oxidation with MeReO₃ and H₂O₂. Using methoxybenzyl (Mob)-protected cysteine, solid-phase oxidation with MeReO₃ and H₂O₂ generated the Mob sulfone precursor to Cys-SO₂^- within fully synthesized peptides. These two solid-phase methods allow the synthesis of peptides containing either Cys-SO₃^- or Cys-SO₂^- in a practical manner, with no solution-phase synthesis required. Cys-SO₃^- had low PPII propensity for PPII propagation, despite promoting a relatively compact conformation in ϕ. In contrast, in a PPII initiation model system, Cys-SO₃^- promoted PPII relative to neutral Cys, with PPII initiation similar to Cys thiolate but less than Cys-SO₂^- or Ala. In an α-helix model system, Cys-SO₃^- promoted α-helix near the N-terminus, due to favorable helix dipole interactions and favorable α-helix capping via a sulfonate-amide side chain-main chain hydrogen bond. Across all peptides, the sulfonate side chain was significantly less ordered than that of the sulfinate. Analysis of Cys-SO₃^- in the PDB revealed a very strong propensity for local (i/i or i/i + 1) side chain-main chain sulfonate-amide hydrogen bonds for Cys-SO₃^-, with >80% of Cys-SO₃^- residues exhibiting these interactions. DFT calculations conducted to explore these conformational preferences indicated that side chain-main chain hydrogen bonds of the sulfonate with the intraresidue amide and/or with the i + 1 amide were favorable. However, hydrogen bonds to water or to amides, as well as interactions with oxophilic metals, were weaker for the sulfonate than the sulfinate, due to lower charge density on the oxygens in the sulfonate.

The efficacy of Paxlovid against COVID-19 is the result of the tight molecular docking between Mpro and antiviral drugs (nirmatrelvir and ritonavir)

Purpose

Currently, a number of medications for coronavirus disease 2019 (COVID-19) treatment are tested in clinical trials; however, credible clinical studies are becoming increasingly difficult to come by. Paxlovid is a ritonavir-boosted nirmatrelvir drug that the U.S. Food and Drug Administration (FDA) authorized for the treatment of COVID-19.

This study aimed to demonstrate the interaction of nirmatrelvir and ritonavir on the active site of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main protease (M^pro).

Materials and methods

To locate the optimal docking between M^pro and antiviral drugs, and to conduct dynamic simulations between atoms in the fusion areas, various bioinformatics and mathematical equations were applied.

Results

According to the docking data, nirmatrelvir has a stronger interaction with M^pro than ritonavir, which has more multiple bonds. Molecular docking of antiviral drugs such as Paxlovid has a significant impact on the treatment of COVID-19 virus.

Conclusions

According to this study, Paxlovid may work on new strains, including Omicron, because the M^pro mutation P132H in the Omicron variant has no direct effect on the protein.

How do salt and lipids affect conformational dynamics of Aβ42 monomers in water?

It is well established that amyloid β-protein (Aβ) self-assembly is involved in triggering of Alzheimer's disease. On the other hand, evidence of physiological function of Aβ interacting with lipids has only begun to emerge. Details of Aβ-lipid interactions, which may underlie physiological and pathological activities of Aβ, are not well understood. Here, the effects of salt and 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) lipids on conformational dynamics of Aβ42 monomer in water are examined by all-atom molecular dynamics (MD). We acquired six sets of 250 ns long MD trajectories for each of the three lipid concentrations (0, 27, and 109 mM) in the absence and presence of 150 mM salt. Ten replica trajectories per set are used to enhance sampling of Aβ42 conformational space. We show that salt facilitates long-range tertiary contacts in Aβ42, resulting in more compact Aβ42 conformations. By contrast, addition of lipids results in lipid-concentration dependent Aβ42 unfolding concomitant with enhanced stability of the turn in the A21-A30 region. At the high lipid concentration, salt enables the N-terminal region of Aβ42 to form long-range tertiary contacts and interact with lipids, which results in formation of a parallel β-strand. Aβ42 forms stable lipid-protein complexes whereby the protein is adhered to the lipid cluster rather than embedded into it. We propose that the inability of Aβ42 monomer to get embedded into the lipid cluster may be important for facilitating repair of leaks in the blood-brain barrier without penetrating and damaging cellular membranes.

High-Resolution Conformational Analysis of RGDechi-Derived Peptides Based on a Combination of NMR Spectroscopy and MD Simulations

The crucial role of integrin in pathological processes such as tumor progression and metastasis formation has inspired intense efforts to design novel pharmaceutical agents modulating integrin functions in order to provide new tools for potential therapies. In the past decade, we have investigated the biological proprieties of the chimeric peptide RGDechi, containing a cyclic RGD motif linked to an echistatin C-terminal fragment, able to specifically recognize αvβ3 without cross reacting with αvβ5 and αIIbβ3 integrin. Additionally, we have demonstrated using two RGDechi-derived peptides, called RGDechi1-14 and ψRGDechi, that chemical modifications introduced in the C-terminal part of the peptide alter or abolish the binding to the αvβ3 integrin. Here, to shed light on the structural and dynamical determinants involved in the integrin recognition mechanism, we investigate the effects of the chemical modifications by exploring the conformational space sampled by RGDechi1-14 and ψRGDechi using an integrated natural-abundance NMR/MD approach. Our data demonstrate that the flexibility of the RGD-containing cycle is driven by the echistatin C-terminal region of the RGDechi peptide through a coupling mechanism between the N- and C-terminal regions.

Enhanced Nuclear Magnetic Resonance Spectroscopy with Isotropic Mixing as a Pseudodimension

Chemical analysis based on liquid-state nuclear magnetic resonance spectroscopy exploits numerous observables, mainly chemical shifts, relaxation rates, and internuclear coupling constants. Regarding the latter, the efficiencies of internuclear coherence transfers may be encoded in spectral peak intensities. The dependencies of these intensities on the experimental parameter that influences the transfer, for example, mixing time, are an important source of structural information. Yet, they are costly to measure and difficult to analyze. Here, we show that peak intensity build-up curves in two-dimensional total correlation spectroscopy (2D TOCSY) experiments may be quickly measured by employing nonuniform sampling and that their analysis can be effective if supported by quantum mechanical calculations. Thus, such curves can be used to form a new, third pseudodimension of the TOCSY spectrum. Similarly to the other two frequency dimensions, this one also resolves ambiguities and provides characteristic information. We show how the approach supports the analysis of a fragment of protein Tau Repeat-4 domain. Yet, its potential applications are far broader, including the analysis of complex mixtures or other polymers.

Solution cis-Proline Conformation of IPCs Inhibitor Aureobasidin A Elucidated via NMR-Based Conformational Analysis

Aureobasidin A (abA) is a natural depsipeptide that inhibits inositol phosphorylceramide (IPC) synthases with significant broad-spectrum antifungal activity. abA is known to have two distinct conformations in solution corresponding to trans- and cis-proline (Pro) amide bond rotamers. While the trans-Pro conformation has been studied extensively, cis-Pro conformers have remained elusive. Conformational properties of cyclic peptides are known to strongly affect both potency and cell permeability, making a comprehensive characterization of abA conformation highly desirable. Here, we report a high-resolution 3D structure of the cis-Pro conformer of aureobasidin A elucidated for the first time using a recently developed NMR-driven computational approach. This approach utilizes ForceGen's advanced conformational sampling of cyclic peptides augmented by sparse distance and torsion angle constraints derived from NMR data. The obtained 3D conformational structure of cis-Pro abA has been validated using anisotropic residual dipolar coupling measurements. Support for the biological relevance of both the cis-Pro and trans-Pro abA configurations was obtained through molecular similarity experiments, which showed a significant 3D similarity between NMR-restrained abA conformational ensembles and another IPC synthase inhibitor, pleofungin A. Such ligand-based comparisons can further our understanding of the important steric and electrostatic characteristics of abA and can be utilized in the design of future therapeutics.

Conformational dynamics in the disordered region of human CPEB3 linked to memory consolidation

Background

Current understanding of the molecular basis of memory consolidation points to an important function of amyloid formation by neuronal-specific isoforms of the cytoplasmic polyadenylation element binding (CPEB) protein family. In particular, CPEB is thought to promote memory persistence through formation of self-sustaining prion-like amyloid assemblies at synapses, mediated by its intrinsically disordered region (IDR) and leading to permanent physical alterations at the basis of memory persistence. Although the molecular mechanisms by which amyloid formation takes place in CPEB have been described in invertebrates, the way amyloid formation occurs in the human homolog CPEB3 (hCPEB3) remains unclear. Here, we characterize by NMR spectroscopy the atomic level conformation and ps-ms dynamics of the 426-residue IDR of hCPEB3, which has been associated with episodic memory in humans.

Results

We show that the 426-residue N-terminal region of hCPEB3 is a dynamic, intrinsically disordered region (IDR) which lacks stable folded structures. The first 29 residues, M₁QDDLLMDKSKTQPQPQQQQRQQQQPQP₂₉, adopt a helical + disordered motif, and residues 86–93: P₈₃QQPPPP₉₃, and 166–175: P₁₆₆PPPAPAPQP₁₇₅ form polyproline II (PPII) helices. The (VG)₅ repeat motif is completely disordered, and residues 200–250 adopt three partially populated α-helices. Residues 345–355, which comprise the nuclear localization signal (NLS), form a modestly populated α-helix which may mediate STAT5B binding. These findings allow us to suggest a model for nascent hCPEB3 structural transitions at single residue resolution, advancing that amyloid breaker residues, like proline, are a key difference between functional versus pathological amyloids.

Conclusion

Our NMR spectroscopic analysis of hCPEB3 provides insights into the first structural transitions involved in protein–protein and protein-mRNA interactions. The atomic level understanding of these structural transitions involved in hCPEB3 aggregation is a key first step toward understanding memory persistence in humans, as well as sequence features that differentiate beneficial amyloids from pathological ones.

Areas

Biophysics, Structural Biology, Biochemistry & Neurosciences.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01310-6.

Exploring Nearest Neighbor Interactions and Their Influence on the Gibbs Energy Landscape of Unfolded Proteins and Peptides

The Flory isolated pair hypothesis (IPH) is one of the corner stones of the random coil model, which is generally invoked to describe the conformational dynamics of unfolded and intrinsically disordered proteins (IDPs). It stipulates, that individual residues sample the entire sterically allowed space of the Ramachandran plot without exhibiting any correlations with the conformational dynamics of its neighbors. However, multiple lines of computational, bioinformatic and experimental evidence suggest that nearest neighbors have a significant influence on the conformational sampling of amino acid residues. This implies that the conformational entropy of unfolded polypeptides and proteins is much less than one would expect based on the Ramachandran plots of individual residues. A further implication is that the Gibbs energies of residues in unfolded proteins or polypeptides are not additive. This review provides an overview of what is currently known and what has yet to be explored regarding nearest neighbor interactions in unfolded proteins.

Convergent evolution of a parasite-encoded complement control protein-scaffold to mimic binding of mammalian TGF-β to its receptors, TβRI and TβRII

The mouse intestinal helminth Heligmosomoides polygyrus modulates host immune responses by secreting a transforming growth factor (TGF)-β mimic (TGM), to expand the population of Foxp3⁺ T_regs. TGM comprises five complement control protein (CCP)-like domains, designated D1-D5. Though lacking homology to TGF-β, TGM binds directly to the TGF-β receptors TβRI and TβRII and stimulates the differentiation of naïve T-cells into T_regs. However, the molecular determinants of binding are unclear. Here, we used surface plasmon resonance, isothermal calorimetry, NMR spectroscopy, and mutagenesis to investigate how TGM binds the TGF-β receptors. We demonstrate that binding is modular, with D1-D2 binding to TβRI and D3 binding to TβRII. D1-D2 and D3 were further shown to compete with TGF-β(TβRII)₂ and TGF-β for binding to TβRI and TβRII, respectively. The solution structure of TGM-D3 revealed that TGM adopts a CCP-like fold but is also modified to allow the C-terminal strand to diverge, leading to an expansion of the domain and opening potential interaction surfaces. TGM-D3 also incorporates a long structurally ordered hypervariable loop, adding further potential interaction sites. Through NMR shift perturbations and binding studies of TGM-D3 and TβRII variants, TGM-D3 was shown to occupy the same site of TβRII as bound by TGF-β using both a novel interaction surface and the hypervariable loop. These results, together with the identification of other secreted CCP-like proteins with immunomodulatory activity in H. polygyrus, suggest that TGM is part of a larger family of evolutionarily plastic parasite effector molecules that mediate novel interactions with their host.

Structural Insights into the Intrinsically Disordered GPCR C-Terminal Region, Major Actor in Arrestin-GPCR Interaction

Arrestin-dependent pathways are a central component of G protein-coupled receptor (GPCRs) signaling. However, the molecular processes regulating arrestin binding are to be further illuminated, in particular with regard to the structural impact of GPCR C-terminal disordered regions. Here, we used an integrated biophysical strategy to describe the basal conformations of the C-terminal domains of three class A GPCRs, the vasopressin V2 receptor (V2R), the growth hormone secretagogue or ghrelin receptor type 1a (GHSR) and the β2-adernergic receptor (β2AR). By doing so, we revealed the presence of transient secondary structures in these regions that are potentially involved in the interaction with arrestin. These secondary structure elements differ from those described in the literature in interaction with arrestin. This suggests a mechanism where the secondary structure conformational preferences in the C-terminal regions of GPCRs could be a central feature for optimizing arrestins recognition.

Measuring the 3J HNHa coupling by a simple 2D-intra-HNCA IP/AP-E.COSY with simultaneous encoding of 15N chemical shift and 1J HaCa evolution

The ³J coupling values are commonly used in biomolecular NMR to extract structural information. Here we present a novel intra HNCA IP/AP E.COSY pulse sequence that allows to measure ³J H^NH^a coupling constants by a simple and rapid two-dimensional ¹H-¹⁵N correlation experiment where the ¹⁵N frequency is encoded at the same time as the ¹J H^aC^a coupling evolution. The advantage with respect to the conventional 3D HNCA E.COSY pulse sequence is the dimensionality reduction to a simple 2D experiment, which decreases acquisition time and facilitates data analysis. The performance of this new experiment is demonstrated with an ubiquitin sample at 500 MHz.

Characterization of Amyloidogenic Peptide Aggregability in Helical Subspace

Prototypical amyloidogenic peptides amyloid-β (Aβ) and α-synuclein (αS) can undergo helix-helix associations via partially folded helical conformers, which may influence pathological progression to Alzheimer's (AD) and Parkinson's disease (PD), respectively. At the other extreme, stable folded helical conformers have been reported to resist self-assembly and amyloid formation. Experimental characterisation of such disparities in aggregation profiles due to subtle differences in peptide stabilities is precluded by the conformational heterogeneity of helical subspace. The diverse physical models used in molecular simulations allow sampling distinct regions of the phase space and are extensive in capturing the ensemble of rich helical subspace. Robust and powerful computational predictive methods utilizing network theory and free energy mapping can model the origin of helical population shifts in amyloidogenic peptides, which highlight their inherent aggregability. In this chapter, we discuss computational models, methods, design rules, and strategies to identify the driving force behind helical self-assembly and the molecular origin of aggregation resistance in helical intermediates of Aβ42 and αS. By extensive multiscale mapping of intrapeptide interactions, we show that the computational models can capture features that are otherwise imperceptible to experiments. Our models predict that targeting terminal residues may allow modulation and control of initial pathogenic aggregability of amyloidogenic peptides.

Predictive Modeling of Neurotoxic α-Synuclein Polymorphs

Assembly of monomeric α-synuclein (αS) into aggregation-resistant helically folded tetramers and related multimers is a key target for Parkinson's disease (PD). Protein dynamics hampers experimental characterization of the polymorphism of these structures and so computational modeling and simulation is providing a complementary approach to obtain high-resolution structural information on the assembly of αS and interactions with biological surfaces. These computational techniques are particularly valuable for intrinsically disordered proteins (IDPs) and short-lived peptide and protein assemblies with as yet undetermined 3D structures. Experimental observables such as NMR J-coupling constants and chemical shifts can be predicted directly from simulation data, and compared with available experimental data to generate the most physically realistic atomic-resolution structure. For appropriately validated and benchmarked computational models, macroscopic aggregation properties can be related to the calculated thermodynamic properties at an atomic level. In this chapter, we describe a useful protocol for designing helical αS multimers, especially tetramers, and scanning the peptide-membrane interface for cell-bound αS tetramers. These computationally modeled structures are validated by comparison with the range of available known experimental parameters at time of writing in early 2020, and used to generate predictive design rules to motivate and guide experiments.

Site-Specific Stabilization and Destabilization of α Helical Peptides upon Phosphorylation and O-GlcNAcylation

Helices (α-helix) are the most common type of secondary structure motif present in proteins. In this study, we have investigated the structural influence of phosphorylation and O-GlcNAcylation, common intracellular post-translational modifications (PTMs), on the α-helical conformation. The simulation studies were performed on the Baldwin model α-helical peptide sequence (Ac-AKAAAAKAAAAKAA-NH₂). The Baldwin sequences were chosen due to the availability of site-specific experimental post-translational data for cross-validation with the simulations. The influence of PTMs was examined across the span of the α-helix, namely, at the N-terminus, position 10 (interior region), and the C-terminus for both serine and threonine residues placed at these positions. Molecular dynamics (MD) simulations revealed that phosphorylation and O-GlcNAcylation at the N-terminus lead to the stabilization of the helical conformation. PTMs in the interior or the C-terminus were found to disrupt helicity, with the disruption being more pronounced for PTMs in the interior region, in accordance with experimental studies. It was found that phosphorylation-derived destabilization was mainly due to the formation of an intraresidue HN-PO₃^2- electrostatic interaction and interactions between the phosphate group and the side chain of adjacent lysine residues (NH₃···PO₃^2-). Hydrophobic and steric clashes were the main causes of destabilization in the case of O-GlcNAcylation. The structural disruptions were found to be more pronounced for PTM at the threonine site when compared to the serine site. The salt-bridge-dependent stability of the α-helix was found to be highly position specific, an i → i + 4 interaction stabilizing the helix, with other placements leading to the destabilization of the helix.

Data Mining of Molecular Simulations Suggest Key Amino Acid Residues for Aggregation, Signaling and Drug Action

Alzheimer's disease, the most common form of dementia, currently has no cure. There are only temporary treatments that reduce symptoms and the progression of the disease. Alzheimer's disease is characterized by the prevalence of plaques of aggregated amyloid β (Aβ) peptide. Recent treatments to prevent plaque formation have provided little to relieve disease symptoms. Although there have been numerous molecular simulation studies on the mechanisms of Aβ aggregation, the signaling role has been less studied. In this study, a total of over 38,000 simulated structures, generated from molecular dynamics (MD) simulations, exploring different conformations of the Aβ42 mutants and wild-type peptides were used to examine the relationship between Aβ torsion angles and disease measures. Unique methods characterized the data set and pinpointed residues that were associated in aggregation and others associated with signaling. Machine learning techniques were applied to characterize the molecular simulation data and classify how much each residue influenced the predicted variant of Alzheimer's Disease. Orange3 data mining software provided the ability to use these techniques to generate tables and rank the data. The test and score module coupled with the confusion matrix module analyzed data with calculations of specificity and sensitivity. These methods evaluating frequency and rank allowed us to analyze and predict important residues associated with different phenotypic measures. This research has the potential to help understand which specific residues of Aβ should be targeted for drug development.

The 1H, 15N and 13C resonance assignments of the C-terminal domain of Serpine mRNA binding protein 1 (SERBP1)

SERBP1 is a multifunctional mRNA-binding protein that has been shown to play a regulatory role in a number of biological processes such as thrombosis, DNA damage repair, and the cellular response to nutrient deprivation. Additionally, SERBP1 is upregulated in glioblastoma, leukemia as well as liver, prostrate and ovarian cancers where it has been implicated in metastatic disease and poor patient outcomes. SERBP1 binds target mRNA, stabilizing and regulating the post-translational expression of the transcript. Since SERBP1 lacks canonical RNA-binding motifs such as RRM domains or zinc fingers, its target recognition and binding mechanisms are not well understood. Recent reports suggest that it is capable of recognizing both RNA sequence motifs and structured domains. Here we report the production and purification of the intrinsically disordered C-terminal domain of SERBP1, the assignment of the 1H, 13C, 15N backbone resonances of the protein by solution-state NMR, and secondary structure predictions. We show that the protein is not entirely disordered and identify an α-helix that was stable under the experimental conditions. This work is the first step toward understanding the structural basis underpinning the molecular mechanisms of SERBP1 functions, particularly interactions with mRNA targets.

Resonance assignment and secondary structure of DbpA protein from the European species, Borrelia afzelii

Decorin binding proteins (Dbps) mediate attachment of spirochetes in host organisms during the early stages of Lyme disease infection. Previously, different binding mechanisms of Dbps to glycosaminoglycans have been elucidated for the pathogenic species Borrelia burgdorferi sensu stricto and B. afzelii. We are investigating various European Borrelia spirochetes and their interactions at the atomic level using NMR. We report preparative scale recombinant expression of uniformly stable isotope enriched B. afzelii DbpA in Escherichia coli, its chromatographic purification, and solution NMR assignments of its backbone and sidechain 1H, 13C, and 15N atoms. This data was used to predict secondary structure propensity, which we compared to the North American B. burgdorferi sensu stricto and European B. garinii DbpA for which solution NMR structures had been determined previously. Backbone dynamics of DbpA from B. afzelii were elucidated from spin relaxation and heteronuclear NOE experiments. NMR-based secondary structure analysis together with the backbone dynamics characterization provided a first look into structural differences of B. afzelii DbpA compared to the North American species and will serve as the basis for further investigation of how these changes affect interactions with host components.

NMR Spectroscopy in the Conformational Analysis of Peptides: An Overview

BACKGROUND

NMR spectroscopy is one of the most powerful tools to study the structure and interaction properties of peptides and proteins from a dynamic perspective. Knowing the bioactive conformations of peptides is crucial in the drug discovery field to design more efficient analogue ligands and inhibitors of protein-protein interactions targeting therapeutically relevant systems.

OBJECTIVE

This review provides a toolkit to investigate peptide conformational properties by NMR.

METHODS

Articles cited herein, related to NMR studies of peptides and proteins were mainly searched through PubMed and the web. More recent and old books on NMR spectroscopy written by eminent scientists in the field were consulted as well.

RESULTS

The review is mainly focused on NMR tools to gain the 3D structure of small unlabeled peptides. It is more application-oriented as it is beyond its goal to deliver a profound theoretical background. However, the basic principles of 2D homonuclear and heteronuclear experiments are briefly described. Protocols to obtain isotopically labeled peptides and principal triple resonance experiments needed to study them, are discussed as well.

CONCLUSION

NMR is a leading technique in the study of conformational preferences of small flexible peptides whose structure can be often only described by an ensemble of conformations. Although NMR studies of peptides can be easily and fast performed by canonical protocols established a few decades ago, more recently we have assisted to tremendous improvements of NMR spectroscopy to investigate instead large systems and overcome its molecular weight limit.

Structural and dynamic characterization of the C-terminal tail of ErbB2: Disordered but not random

ErbB2 (or HER2) is a receptor tyrosine kinase overexpressed in some breast cancers and associated with poor prognosis. Treatments targeting the receptor extracellular and kinase domains have greatly improved disease outcome in the last 20 years. In parallel, the structures of these domains have been described, enabling better mechanistic understanding of the receptor function and targeted inhibition. However, the ErbB2 disordered C-terminal cytoplasmic tail (CtErbB2) remains very poorly characterized in terms of structure, dynamics, and detailed functional mechanism. Yet, it is where signal transduction is triggered via phosphorylation of tyrosine residues and carried out via interaction with adaptor proteins. Here, we report the first description, to our knowledge, of the ErbB2 disordered tail at atomic resolution using NMR, complemented by small-angle x-ray scattering. We show that although no part of CtErbB2 has any fully populated secondary or tertiary structure, it contains several transient α-helices and numerous transient polyproline II helices, populated up to 20 and 40%, respectively, and low but significant compaction. The presence of some structural elements suggests, along the lines of the results obtained for EGFR (ErbB1), that they may have a functional role in ErbB2's autoregulation processes. In addition, the transient formation of polyproline II helices is compliant with previously suggested interactions with SH3 domains. All in all, our in-depth structural study opens perspectives in the mechanistic understanding of ErbB2.

Structural Insights into the Interaction of the Intrinsically Disordered Co-activator TIF2 with Retinoic Acid Receptor Heterodimer (RXR/RAR)

Retinoic acid receptors (RARs) and retinoid X receptors (RXRs) form heterodimers that activate target gene transcription by recruiting co-activator complexes in response to ligand binding. The nuclear receptor (NR) co-activator TIF2 mediates this recruitment by interacting with the ligand-binding domain (LBD) of NRs trough the nuclear receptor interaction domain (TIF2_NRID) containing three highly conserved α-helical LxxLL motifs (NR-boxes). The precise binding mode of this domain to RXR/RAR is not clear due to the disordered nature of TIF2. Here we present the structural characterization of TIF2_NRID by integrating several experimental (NMR, SAXS, Far-UV CD, SEC-MALS) and computational data. Collectively, the data are in agreement with a largely disordered protein with partially structured regions, including the NR-boxes and their flanking regions, which are evolutionary conserved. NMR and X-ray crystallographic data on TIF2_NRID in complex with RXR/RAR reveal a multisite binding of the three NR-boxes as well as an active role of their flanking regions in the interaction.

A disulfide-stabilised helical hairpin fold in acrorhagin I: An emerging structural motif in peptide toxins

Acrorhagin I (U-AITX-Aeq5a) is a disulfide-rich peptide identified in the aggressive organs (acrorhagi) of the sea anemone Actinia equina. Previous studies (Toxicon 2005, 46:768-74) found that the peptide is toxic in crabs, although the structural and functional properties of acrorhagin I have not been reported. In this work, an Escherichia coli (BL21 strain) expression system was established for the preparation of ¹³C,¹⁵N-labelled acrorhagin I, and the solution structure was determined using NMR spectroscopy. Structurally, acrorhagin I is similar to B-IV toxin from the marine worm Cerebratulus lacteus (PDB id 1VIB), with a well-defined helical hairpin structure stabilised by four intramolecular disulfide bonds. The recombinant peptide was tested in patch-clamp electrophysiology assays against voltage-gated potassium and sodium channels, and in bacterial and fungal growth inhibitory assays and haemolytic assays. Acrorhagin I was not active against any of the ion channels tested and showed no activity in functional assays, indicating that this peptide may possess a different biological function. Metal ion interaction studies using NMR spectroscopy showed that acrorhagin I bound zinc and nickel, suggesting that its function might be modulated by metal ions or that it may be involved in regulating metal ion levels and their transport. The similarity between the structure of acrorhagin I and that of B-IV toxin from a marine worm suggests that this fold may prove to be a recurring motif in disulfide-rich peptides from marine organisms.

Structural transitions in Orb2 prion-like domain relevant for functional aggregation in memory consolidation

The recent structural elucidation of ex vivo Drosophila Orb2 fibrils revealed a novel amyloid formed by interdigitated Gln and His residue side chains belonging to the prion-like domain. However, atomic-level details on the conformational transitions associated with memory consolidation remain unknown. Here, we have characterized the nascent conformation and dynamics of the prion-like domain (PLD) of Orb2A using a nonconventional liquid-state NMR spectroscopy strategy based on 13C detection to afford an essentially complete set of 13Cα, 13Cβ, 1Hα, and backbone 13CO and 15N assignments. At pH 4, where His residues are protonated, the PLD is disordered and flexible, except for a partially populated α-helix spanning residues 55-60, and binds RNA oligos, but not divalent cations. At pH 7, in contrast, His residues are predominantly neutral, and the Q/H segments adopt minor populations of helical structure, show decreased mobility and start to self-associate. At pH 7, the His residues do not bind RNA or Ca2+, but do bind Zn2+, which promotes further association. These findings represent a remarkable case of structural plasticity, based on which an updated model for Orb2A functional amyloidogenesis is suggested.

Structural Characteristics of Monomeric Aβ42 on Fibril in the Early Stage of Secondary Nucleation Process

Amyloid-β (Aβ) aggregates are believed to be one of the main causes of Alzheimer's disease. Aβ peptides form fibrils having cross β-sheet structures mainly through primary nucleation, secondary nucleation, and elongation. In particular, self-catalyzed secondary nucleation is of great interest. Here, we investigate the adsorption of Aβ42 peptides to the Aβ42 fibril to reveal a role of adsorption as a part of secondary nucleation. We performed extensive molecular dynamics simulations based on replica exchange with solute tempering 2 (REST2) to two systems: a monomeric Aβ42 in solution and a complex of an Aβ42 peptide and Aβ42 fibril. Results of our simulations show that the Aβ42 monomer is extended on the fibril. Furthermore, we find that the hairpin structure of the Aβ42 monomer decreases but the helix structure increases by adsorption to the fibril surface. These structural changes are preferable for forming fibril-like aggregates, suggesting that the fibril surface serves as a catalyst in the secondary nucleation process. In addition, the stabilization of the helix structure of the Aβ42 monomer on the fibril indicates that the strategy of a secondary nucleation inhibitor design for Aβ40 can also be used for Aβ42.

Resolving overlapped signals with automated FitNMR analytical peak modeling

Nuclear magnetic resonance (NMR) is a valuable tool for determining the structures of molecules and probing their dynamics. A longstanding problem facing both small-molecule and macromolecular NMR is overlapped signals in crowded spectra. To address this, we have developed a method that extracts peak features by fitting analytically derived models of NMR lineshapes. The approach takes into account the effects of truncation, apodization, and the resulting artifacts, while avoiding systematic errors that have affected other models. Even severely overlapped peaks, beyond the point of coalescence, can be distinguished in both simulated and experimental data. We show that the method can measure unresolved backbone scalar couplings directly from a 2D proton-nitrogen spectrum of a de novo designed mini protein. The algorithm is implemented in the FitNMR open-source R package and can be used to analyze nearly any type of single or multidimensional data from small molecules or biomolecules.

Role of Tyr-39 for the Structural Features of α-Synuclein and for the Interaction with a Strong Modulator of Its Amyloid Assembly

Recent studies suggest that Tyr-39 might play a critical role for both the normal function and the pathological dysfunction of α-synuclein (αS), an intrinsically disordered protein involved in Parkinson’s disease. We perform here a comparative analysis between the structural features of human αS and its Y39A, Y39F, and Y39L variants. By the combined application of site-directed mutagenesis, biophysical techniques, and enhanced sampling molecular simulations, we show that removing aromatic functionality at position 39 of monomeric αS leads to protein variants populating more compact conformations, conserving its disordered nature and secondary structure propensities. Contrasting with the subtle changes induced by mutations on the protein structure, removing aromaticity at position 39 impacts strongly on the interaction of αS with the potent amyloid inhibitor phthalocyanine tetrasulfonate (PcTS). Our findings further support the role of Tyr-39 in forming essential inter and intramolecular contacts that might have important repercussions for the function and the dysfunction of αS.

Protein Motional Details Revealed by Complementary Structural Biology Techniques

Proteins depend on defined molecular plasticity for their functionality. How to comprehensively capture dynamics correctly is of ubiquitous biological importance. Approaches commonly used to probe protein dynamics include model-free elucidation of site-specific motion by NMR relaxation, molecular dynamics (MD)-based approaches, and capturing the substates within a dynamic ensemble by recent eNOE-based multiple-structure approaches. Even though MD is sometimes combined with ensemble-averaged NMR restraints, these approaches have largely been developed and used individually. Owing to the different underlying concepts and practical requirements, it has remained unclear how they compare, and how they cross-validate and complement each other. Here, we extract and compare the differential information contents of MD simulations, NMR relaxation measurements, and eNOE-based multi-state structures for the SH3 domain of chicken α-spectrin. The data show that a validated, consistent, and detailed picture is feasible both for timescales and actual conformational states sampled in the dynamic ensemble. This includes the biologically important side-chain plasticity, for which experimentally cross-validated assessment is a significant challenge.

Rapid Biophysical Characterization and NMR Spectroscopy Structural Analysis of Small Proteins from Bacteria and Archaea

Proteins encoded by small open reading frames (sORFs) have a widespread occurrence in diverse microorganisms and can be of high functional importance. However, due to annotation biases and their technically challenging direct detection, these small proteins have been overlooked for a long time and were only recently rediscovered. The currently rapidly growing number of such proteins requires efficient methods to investigate their structure-function relationship. Herein, a method is presented for fast determination of the conformational properties of small proteins. Their small size makes them perfectly amenable for solution-state NMR spectroscopy. NMR spectroscopy can provide detailed information about their conformational states (folded, partially folded, and unstructured). In the context of the priority program on small proteins funded by the German research foundation (SPP2002), 27 small proteins from 9 different bacterial and archaeal organisms have been investigated. It is found that most of these small proteins are unstructured or partially folded. Bioinformatics tools predict that some of these unstructured proteins can potentially fold upon complex formation. A protocol for fast NMR spectroscopy structure elucidation is described for the small proteins that adopt a persistently folded structure by implementation of new NMR technologies, including automated resonance assignment and nonuniform sampling in combination with targeted acquisition.

Unravelling the effect of N(ε)-(carboxyethyl)lysine on the conformation, dynamics and aggregation propensity of α-synuclein

α-Synuclein (αS) aggregation is a hallmark in several neurodegenerative diseases. Among them, Parkinson's disease is highlighted, characterized by the intraneuronal deposition of Lewy bodies (LBs) which causes the loss of dopaminergic neurons. αS is the main component of LBs and in them, it usually contains post-translational modifications. One of them is the formation of advanced glycation end-products (mainly CEL and MOLD) arising from its reaction with methylglyoxal. Despite its biological relevance, there are no data available proving the effect of glycation on the conformation of αS, nor on its aggregation mechanism. This has been hampered by the formation of a heterogeneous set of compounds that precluded conformational studies. To overcome this issue, we have here produced αS homogeneously glycated with CEL. Its use, together with different biophysical techniques and molecular dynamics simulations, allowed us to study for the first time the effect of glycation on the conformation of a protein. CEL extended the conformation of the N-terminal domain as a result of the loss of transient N-/C-terminal long-range contacts while increasing the heterogeneity of the conformational population. CEL also inhibited the αS aggregation, but it was not able to disassemble preexisting amyloid fibrils, thus proving that CEL found on LBs must be formed in a later event after aggregation.

We study the effect of an advanced glycation end product (N(ε)-(carboxyethyl)lysine), found on the Lewy bodies of people suffering from Parkinson’s disease, on the conformational and aggregation features of alpha-synuclein.

Peptide-based short single β-strand mimics without hydrogen bonding or aggregation

Generation of short peptides with a single β-strand structure in solution is difficult. Herein, we design a new class of single β-strand peptidic mimics that are stable without self-aggregation in protic and non-protic solvents. Introduction of one present β-strand mimic can induce and propagate the β-strand structure for at least a penta-peptide sequence.

Solution structure of the cytoplasmic domain of NhaP2 a K+/H+ antiporter from Vibrio cholera

NhaP2 is a K⁺/H⁺ antiporter from Vibrio cholerae which consists of a transmembrane domain and a cytoplasmic domain of approximately 200 amino acids, both of which are required for cholera infectivity. Here we present the solution structure for a 165 amino acid minimal cytoplasmic domain (P2MIN) form of the protein. The structure reveals a compact N-terminal domain which resembles a Regulator of Conductance of K⁺ channels (RCK) domain connected to a more open C-terminal domain via a flexible 20 amino acid linker. NMR titration experiments showed that the protein binds ATP through its N-terminal domain, which was further supported by waterLOGSY and Saturation Transfer Difference NMR experiments. The two-domain organisation of the protein was confirmed by BIOSAXS, which also revealed that there are no detectable-ATP-induced conformational changes in the protein structure. Finally, in contrast to all known RCK domain structures solved to date, the current work shows that the protein is a monomer.

Structural preferences of cysteine sulfinic acid: The sulfinate engages in multiple local interactions with the peptide backbone

Cysteine sulfinic acid (Cys-SO₂^-) is a non-enzymatic oxidative post-translational modification (PTM) that has been identified in hundreds of proteins. However, the effects of cysteine sulfination are in most cases poorly understood. Cys-SO₂^- is structurally distinctive, with long sulfur-carbon and sulfur-oxygen bonds, and with tetrahedral geometry around sulfur due to its lone pair. Cys-SO₂^- thus has a unique range of potential interactions with the protein backbone which could facilitate protein structural changes. Herein, the structural effects of cysteine oxidation to the sulfinic acid were investigated in model peptides and folded proteins using NMR spectroscopy, circular dichroism, bioinformatics, and computational studies. In the PDB, Cys-SO₂^- shows a greater preference for α-helix than Cys. In addition, Cys-SO₂^- is more commonly found in structures with φ > 0, including in multiple types of β-turn. Sulfinate oxygens engage in hydrogen bonds with adjacent (i or i + 1) amide hydrogens. Over half of sulfinates have at least one hydrogen bond with an adjacent amide, and several structures have hydrogen bonds with both adjacent amides. Alternately, sulfur or either oxygen can act as an electron donor for n→π* interactions with the backbone carbonyl of the same residue, as indicated by frequent S⋯CO or O⋯CO distances below the sums of their van der Waals radii in protein structures. In peptides, Cys-SO₂^- favored α-helical structure at the N-terminus, consistent with helix dipole effects and backbone hydrogen bonds with the sulfinate promoting α-helix. Cys-SO₂^- has only modestly greater polyproline II helix propensity than Cys-SH, likely due to competition from multiple side chain-backbone interactions. Cys-SO₂^- stabilizes the i+1 position of a β-turn relative to Cys-SH. Within proteins, the range of side chain-main chain interactions available to Cys-SO₂^- compared to Cys-SH provides a basis for potential changes in protein structure and function due to cysteine oxidation to the sulfinic acid.

Assessing the One-Bond Cα-H Spin-Spin Coupling Constants in Proteins: Pros and Cons of Different Approaches

In the present work, we explore three different approaches for the computation of the one-bond spin-spin coupling constants (SSCC) 1JCαH in proteins: density functional theory (DFT) calculations, a Karplus-like equation, and Gaussian process regression. The main motivation of this work is to select the best method for fast and accurate computation of the 1JCαH SSCC, for its use in everyday applications in protein structure validation, refinement, and/or determination. Our initial results showed a poor agreement between the DFT-computed and observed 1JCαH SSCC values. Further analysis leads us to the understanding that the model chosen for the DFT computations is inappropriate and that more complex models will require a higher, if not prohibitively, computational cost. Finally, we show that the Karplus-like equation and Gaussian Process regression provide faster and more accurate results than DFT-based calculations.

Conformational Priming of RepA-WH1 for Functional Amyloid Conversion Detected by NMR Spectroscopy

How proteins with a stable globular fold acquire the amyloid state is still largely unknown. RepA, a versatile plasmidic DNA binding protein from Pseudomonas savastanoi, is functional as a transcriptional repressor or as an initiator or inhibitor of DNA replication, the latter via assembly of an amyloidogenic oligomer. Its N-terminal domain (WH1) is responsible for discrimination between these functional abilities by undergoing insufficiently understood structural changes. RepA-WH1 is a stable dimer whose conformational dynamics had not been explored. Here, we have studied it through NMR {¹H}-¹⁵N relaxation and H/D exchange kinetics measurements. The N- and the C-terminal α-helices, and the internal amyloidogenic loop, are partially unfolded in solution. S4-indigo, a small inhibitor of RepA-WH1 amyloidogenesis, binds to and tethers the N-terminal α-helix to a β-hairpin that is involved in dimerization, thus providing evidence for a priming role of fraying ends and dimerization switches in the amyloidogenesis of folded proteins.

Folding dynamics of Aβ42 monomer at pH 4.0–7.5 with and without physiological salt conditions – does the β1 or β2 region fold first?

The seeding region of Aβ42 monomer is jointly affected by the solution acidity, ionic distribution of the salt, and charged residues.