Strategy for improving circular dichroism spectra deconvolution accuracy for macrocyclic peptides in drug discovery

Peptide therapeutics have emerged as an appealing modality in the pharmaceutical industry. Understanding peptide conformation in solution remains one of the most critical areas for peptide drug development. Circular dichroism (CD) spectroscopy is a useful technique to study the secondary structure of proteins and peptides, but the current approaches are limited to protein-focused models to predict high-order structures of peptides, and the models were built based on X-ray crystallography instead of solution-based technique, as a result, such models may have poor predictions for peptides. In this study, we present a novel CD deconvolution model to determine peptide conformation in solution. To quantitatively obtain secondary structure information using CD, a calibration model is needed beforehand to establish the relationship between each secondary structure feature and the corresponding CD response. A reference set containing the majority of cyclic peptides with known structures from solution-state NMR spectroscopy was used to build the calibration model for CD deconvolution. Improved prediction accuracy on the secondary structure determination for cyclic peptides was achieved by this model compared to the commercial standard model using commercially available platforms. This new CD deconvolution method is crucial for peptide conformational analysis in solution, and has the potential to greatly accelerate peptide drug candidate optimization in the pharmaceutical drug discovery field.

A story of peptides, lipophilicity and chromatography - back and forth in time

Peptides, as part of the beyond the rule of 5 (bRo5) chemical space, represent a unique class of pharmaceutical compounds. Because of their exceptional position in the chemical space between traditional small molecules (molecular weight (MW) < 500 Da) and large therapeutic proteins (MW > 5000 Da), peptides became promising candidates for targeting challenging binding sites, including even targets traditionally considered as undruggable - e.g. intracellular protein-protein interactions. However, basic knowledge about physicochemical properties that are important for a drug to be membrane permeable is missing but would enhance the drug discovery process of bRo5 molecules. Consequently, there is a demand for quick and simple lipophilicity determination methods for peptides. In comparison to the traditional lipophilicity determination methods via shake flask and in silico prediction, chromatography-based methods could have multiple benefits such as the requirement of low analyte amount, insensitivity to impurities and high throughput. Herein we elucidate the role of peptide lipophilicity and different lipophilicity values. Further, we summarize peptide analysis via common chromatographic techniques, in specific reversed phase liquid chromatography, hydrophilic interaction liquid chromatography and supercritical fluid chromatography and their role in drug discovery and development process.

Quantitative molecular simulations

All-atom simulations can provide molecular-level insights into the dynamics of gas-phase, condensed-phase and surface processes. One important requirement is a sufficiently realistic and detailed description of the underlying intermolecular interactions. The present perspective provides an overview of the present status of quantitative atomistic simulations from colleagues' and our own efforts for gas- and solution-phase processes and for the dynamics on surfaces. Particular attention is paid to direct comparison with experiment. An outlook discusses present challenges and future extensions to bring such dynamics simulations even closer to reality.

Atomistic Simulations for Reactions and Vibrational Spectroscopy in the Era of Machine Learning─Quo Vadis?

Atomistic simulations using accurate energy functions can provide molecular-level insight into functional motions of molecules in the gas and in the condensed phase. This Perspective delineates the present status of the field from the efforts of others and some of our own work and discusses open questions and future prospects. The combination of physics-based long-range representations using multipolar charge distributions and kernel representations for the bonded interactions is shown to provide realistic models for the exploration of the infrared spectroscopy of molecules in solution. For reactions, empirical models connecting dedicated energy functions for the reactant and product states allow statistically meaningful sampling of conformational space whereas machine-learned energy functions are superior in accuracy. The future combination of physics-based models with machine-learning techniques and integration into all-purpose molecular simulation software provides a unique opportunity to bring such dynamics simulations closer to reality.

Multipolar Force Fields for Amide-I Spectroscopy from Conformational Dynamics of the Alanine Trimer

The dynamics and spectroscopy of N-methyl-acetamide (NMA) and trialanine in solution are characterized from molecular dynamics simulations using different energy functions, including a conventional point charge (PC)-based force field, one based on a multipolar (MTP) representation of the electrostatics, and a semiempirical DFT method. For the 1D infrared spectra, the frequency splitting between the two amide-I groups is 10 cm^-1 from the PC, 13 cm^-1 from the MTP, and 47 cm^-1 from self-consistent charge density functional tight-binding (SCC-DFTB) simulations, compared with 25 cm^-1 from experiment. The frequency trajectory required for the frequency fluctuation correlation function (FFCF) is determined from individual normal mode (INM) and full normal mode (FNM) analyses of the amide-I vibrations. The spectroscopy, time-zero magnitude of the FFCF C(t = 0), and the static component Δ₀² from simulations using MTP and analysis based on FNM are all consistent with experiments for (Ala)₃. Contrary to this, for the analysis excluding mode-mode coupling (INM), the FFCF decays to zero too rapidly and for simulations with a PC-based force field, the Δ₀² is too small by a factor of two compared with experiments. Simulations with SCC-DFTB agree better with experiment for these observables than those from PC-based simulations. The conformational ensemble sampled from simulations using PCs is consistent with the literature (including P_II, β, α_R, and α_L), whereas that covered by the MTP-based simulations is dominated by P_II with some contributions from β and α_R. This agrees with and confirms recently reported Bayesian-refined populations based on 1D infrared experiments. FNM analysis together with a MTP representation provides a meaningful model to correctly describe the dynamics of hydrated trialanine.

Low-Frequency Vibronic Mixing Modulates the Excitation Energy Flow in Bacterial Light-Harvesting Complex II

Oscillatory features observed in two-dimensional electronic spectroscopy (2DES) manifest coherent vibrational and electronic dynamics and even the interplay of them. Recently, we developed a 2DES technique utilizing a pair of synchronized repetition-frequency-stabilized lasers, which enables the wide dynamic range measurements of 2DES signals rapidly. Here, we apply this dual-laser 2DES technique to investigate the electronic energy transfer (EET) process in bacterial light-harvesting complex II consisting of B800 and B850 circular aggregates at ambient temperature, and the coherent vibrational wavepakcet associated with the EET between the two aggregates are measured. Examining the principal component analysis of the time-resolved 2DES signal, we found that the EET from B800 to low-lying B850 states is modulated by a low-frequency (156 cm^-1) vibrational mode of the exciton donor (B800). This observation suggests that the donor transition density is modulated by this vibration, which, in turn, modulates the energy transfer dynamics.

Structural basis of SARS-CoV-2 spike protein induced by ACE2

Motivation

The recent emergence of the novel SARS-coronavirus 2 (SARS-CoV-2) and its international spread pose a global health emergency. The spike (S) glycoprotein binds ACE2 and promotes SARS-CoV-2 entry into host cells. The trimeric S protein binds the receptor using the receptor-binding domain (RBD) causing conformational changes in S protein that allow priming by host cell proteases. Unraveling the dynamic structural features used by SARS-CoV-2 for entry might provide insights into viral transmission and reveal novel therapeutic targets. Using structures determined by X-ray crystallography and cryo-EM, we performed structural analysis and atomic comparisons of the different conformational states adopted by the SARS-CoV-2-RBD.

Results

Here, we determined the key structural components induced by the receptor and characterized their intramolecular interactions. We show that κ-helix (polyproline-II) is a predominant structure in the binding interface and in facilitating the conversion to the active form of the S protein. We demonstrate a series of conversions between switch-like κ-helix and β-strand, and conformational variations in a set of short α-helices which affect the hinge region. These conformational changes lead to an alternating pattern in conserved disulfide bond configurations positioned at the hinge, indicating a possible disulfide exchange, an important allosteric switch implicated in viral entry of various viruses, including HIV and murine coronavirus. The structural information presented herein enables to inspect and understand the important dynamic features of SARS-CoV-2-RBD and propose a novel potential therapeutic strategy to block viral entry. Overall, this study provides guidance for the design and optimization of structure-based intervention strategies that target SARS-CoV-2.

Availability

We have implemented the proposed methods in an R package freely available at https://github.com/Grantlab/bio3d

Supplementary information

Supplementary data are available at Bioinformatics online.

Refinement of Peptide Conformational Ensembles by 2D IR Spectroscopy: Application to Ala‒Ala‒Ala

Characterizing ensembles of intrinsically disordered proteins is experimentally challenging because of the ill-conditioned nature of ensemble determination with limited data and the intrinsic fast dynamics of the conformational ensemble. Amide I two-dimensional infrared (2D IR) spectroscopy has picosecond time resolution to freeze structural ensembles as needed for probing disordered-protein ensembles and conformational dynamics. Also, developments in amide I computational spectroscopy now allow a quantitative and direct prediction of amide I spectra based on conformational distributions drawn from molecular dynamics simulations, providing a route to ensemble refinement against experimental spectra. We performed a Bayesian ensemble refinement method on Ala-Ala-Ala against isotope-edited Fourier-transform infrared spectroscopy and 2D IR spectroscopy and tested potential factors affecting the quality of ensemble refinements. We found that isotope-edited 2D IR spectroscopy provides a stringent constraint on Ala-Ala-Ala conformations and returns consistent conformational ensembles with the dominant ppII conformer across varying prior distributions from many molecular dynamics force fields and water models. The dominant factor influencing ensemble refinements is the systematic frequency uncertainty from spectroscopic maps. However, the uncertainty of conformer populations can be significantly reduced by incorporating 2D IR spectra in addition to traditional Fourier-transform infrared spectra. Bayesian ensemble refinement against isotope-edited 2D IR spectroscopy thus provides a route to probe equilibrium-complex protein ensembles and potentially nonequilibrium conformational dynamics.

Molecular dynamics-derived rotamer libraries for d-amino acids within homochiral and heterochiral polypeptides

Computational resources have contributed to the design and engineering of novel proteins by integrating genomic, structural and dynamic aspects of proteins. Non-canonical amino acids, such as d-amino acids, expand the available sequence space for designing and engineering proteins; however, the rotamer libraries for d-amino acids are usually constructed as the mirror images of l-amino acid rotamer libraries, an assumption that has not been tested. To this end, we have performed molecular dynamics (MD) simulations of model host-guest peptide systems containing d-amino acids. Our simulations systematically address the applicability of the mirror image convention as well as the effects of neighboring residue chirality. Rotamer libraries derived from these systems provide realistic rotamer distributions suitable for use in both rational and computational design workflows. Our simulations also address the impact of chirality on the intrinsic conformational preferences of amino acids, providing fundamental insights into the relationship between chirality and biomolecular dynamics. While d-amino acids are rare in naturally occurring proteins, they are used in designed proteins to stabilize a desired conformation, increase bioavailability or confer favorable biochemical and physical attributes. Here, we present d-amino acid rotamer libraries derived from MD simulations of alanine-based host-guest pentapeptides and show how certain residues can deviate from mirror image symmetry. Our simulations directly model d-amino acids as guest residues within the chiral l-Ala and d-Ala pentapeptide series to explicitly incorporate any contributions resulting from the chiralities of neighboring residues.

Quantum Mechanics Study on Hydrophilic and Hydrophobic Interactions in the Trivaline-Water System

With the aim to elucidate hydrophobic effects in the unfolded state of peptides, DFT-M062X computations on the Val₃H⁺· nH₂O ( n up to 22) clusters have been accomplished. As far as the main chain is concerned, four conformers with β-strand and/or polyproline type II conformations, PPII (indicated as β-β, β-PPII, PPII-β, and PPII-PPII), have been found by changing the ϕ and ψ angles. For bare peptide, the side chain (isopropyl) of each residue can independently take on three different orientations with negligible effects on energetics. The great isopropyl spatial separations in β-β and β-PPII conformers allow for the construction of synergic and extensive water-water and water-peptide H-bonding in the minimal hydration Val₃H⁺·22H₂O models without significant steric encumbrance. Conversely, due to the proximity of the isopropyl of the central residue with the other two, some restrictions in the water shell construction around the peptide become evident for the PPII-PPII conformer and the number of energetically accessible structures decreases. This is indicative of correlated motion involving isopropyls and backbone mediated by water molecules, the origin of the nearest neighbor effects. Comparing the thermodynamic data of Ala₃H⁺·22H₂O and Val₃H⁺·22H₂O, what emerges is that both hydration enthalpy and entropy drive the β-strand stability of the latter.

Quantum Mechanics Approach to Hydration Energies and Structures of Alanine and Dialanine

A systematic approach to the phenomena related to hydration of biomolecules is reported at the state of the art of electronic-structure methods. Large-scale CCSD(T), MP4-SDQ, MP2, and DFT(M06-2X) calculations for some hydrated complexes of alanine and dialanine (Ala⋅13 H₂ O, Ala₂ H⁺ ⋅18 H₂ O, and Ala₂ ⋅18 H₂ O) are compared with experimental data and other elaborate modeling to assess the reliability of a simple bottom-up approach. The inclusion of a minimal number of water molecules for microhydration of the polar groups together with the polarizable continuum model is sufficient to reproduce the relative bulk thermodynamic functions of the considered biomolecules. These quantities depend on the adopted electronic-structure method, which should be chosen with great care. Nevertheless, the computationally feasible MP2 and M06-2X functionals with the aug-cc-pVTZ basis set satisfactorily reproduce values derived by high-level CCSD(T) and MP4-SDQ methods, and thus they are suitable for future developments of more elaborate and hence more biochemically significant peptides.

Structure of Penta-Alanine Investigated by Two-Dimensional Infrared Spectroscopy and Molecular Dynamics Simulation

We have studied the structure of (Ala)5, a model unfolded peptide, using a combination of 2D IR spectroscopy and molecular dynamics (MD) simulation. Two different isotopomers, each bis-labeled with (13)C═O and (13)C═(18)O, were strategically designed to shift individual site frequencies and uncouple neighboring amide-I' modes. 2D IR spectra taken under the double-crossed ⟨π/4, -π/4, Y, Z⟩ polarization show that the labeled four-oscillator systems can be approximated by three two-oscillator systems. By utilizing the different polarization dependence of diagonal and cross peaks, we extracted the coupling constants and angles between three pairs of amide-I' transition dipoles through spectral fitting. These parameters were related to the peptide backbone dihedral angles through DFT calculated maps. The derived dihedral angles are all located in the polyproline-II (ppII) region of the Ramachandran plot. These results were compared to the conformations sampled by Hamiltonian replica-exchange MD simulations with three different CHARMM force fields. The C36 force field predicted that ppII is the dominant conformation, consistent with the experimental findings, whereas C22/CMAP predicted similar population for α+, β, and ppII, and the polarizable Drude-2013 predicted dominating β structure. Spectral simulation based on MD representative conformations and structure ensembles demonstrated the need to include multiple 2D spectral features, especially the cross-peak intensity ratio and shape, in structure determination. Using 2D reference spectra defined by the C36 structure ensemble, the best spectral simulation is achieved with nearly 100% ppII population, although the agreement with the experimental cross-peak intensity ratio is still insufficient. The dependence of population determination on the choice of reference structures/spectra and the current limitations on theoretical modeling relating peptide structures to spectral parameters are discussed. Compared with the previous results on alanine based oligopeptides, the dihedral angles of our fitted structure, and the most populated ppII structure from the C36 simulation are in good agreement with those suggesting a major ppII population. Our results provide further support for the importance of ppII conformation in the ensemble of unfolded peptides.

Interplay between Intrinsic Conformational Propensities and Intermolecular Interactions in the Self-Assembly of Short Surfactant-like Peptides Composed of Leucine/Isoleucine

To study how the conformational propensities of individual amino acid residues, primary structures (i.e., adjacent residues and molecular lengths), and intermolecular interactions of peptides affect their self-assembly properties, we report the use of replica exchange molecular dynamics (REMD) to investigate the monomers, dimers, and trimers of a series of short surfactant-like peptides (I3K, L3K, L4K, and L5K). For four-residue peptides X3K (I3K and L3K), the results show that their different aggregation behaviors arise from the different intrinsic conformational propensities of isoleucine and leucine. For LmK peptides (L3K, L4K, and L5K), the molecular length is found to dictate their aggregation via primarily modulating intermolecular interactions. Increasing the number of hydrophobic amino acid residues of LmK peptides enhances their intermolecular H-bonding and promotes the formation of β-strands in dimer and trimer aggregates, overwhelming the intrinsic preference of Leu for helical structures. Thus, the interplay between the conformational propensities of individual amino acid residues for secondary structures and molecular interactions determines the self-assembly properties of the peptides, and the competition between intramolecular and intermolecular H-bonding interactions determines the probability of β-sheet alignment of peptide molecules. These results are validated by comparing simulated and experimental CD spectra of the peptides. This study will aid the design of short peptide amphiphiles and improve the mechanistic understanding of their self-assembly behavior.

Interfacial water at the trialanine hydrophilic surface: a DFT electronic structure and bottom-up investigation

A model describing a network of hydrogen bonded water-trialanine has been developed to estimate hydration effects on various conformers of the peptide.

Intrinsic α-helical and β-sheet conformational preferences: a computational case study of alanine

A fundamental question in protein science is what is the intrinsic propensity for an amino acid to be in an α-helix, β-sheet, or other backbone dihedral angle ( ϕ-ψ) conformation. This question has been hotly debated for many years because including all protein crystal structures from the protein database, increases the probabilities for α-helical structures, while experiments on small peptides observe that β-sheet-like conformations predominate. We perform molecular dynamics (MD) simulations of a hard-sphere model for Ala dipeptide mimetics that includes steric interactions between nonbonded atoms and bond length and angle constraints with the goal of evaluating the role of steric interactions in determining protein backbone conformational preferences. We find four key results. For the hard-sphere MD simulations, we show that (1) β-sheet structures are roughly three and half times more probable than α-helical structures, (2) transitions between α-helix and β-sheet structures only occur when the backbone bond angle τ (NCα C) is greater than 110°, and (3) the probability distribution of τ for Ala conformations in the "bridge" region of ϕ-ψ space is shifted to larger angles compared to other regions. In contrast, (4) the distributions obtained from Amber and CHARMM MD simulations in the bridge regions are broader and have increased τ compared to those for hard sphere simulations and from high-resolution protein crystal structures. Our results emphasize the importance of hard-sphere interactions and local stereochemical constraints that yield strong correlations between ϕ-ψ conformations and τ.

Neighboring residue effects in terminally blocked dipeptides: implications for residual secondary structures in intrinsically unfolded/disordered proteins

For nuclear magnetic resonance (NMR)-based protein structure determinations, the random coil chemical shifts are very important because the secondary and tertiary protein structure predictions become possible by examining deviations of measured chemical shifts from those reference chemical shift values. In addition, neighboring residue effects on chemical shifts and J-coupling constants are crucial in understanding the nature of conformational propensities exhibited by unfolded or intrinsically disordered proteins. We recently reported the 1D NMR results for a complete set of terminally blocked dipeptides (Oh KI, Jung YS, Hwang GS, Cho M. J Biomol NMR 2012;53:25-41), but the NMR resonance assignments were not possible so that the average chemical shifts and J-coupling constants were only considered. In the present work, to thoroughly investigate the neighboring residue effects and random coil chemical shifts we extend the previous studies with 2D NMR, and measured all the (3) J(HNHα) values and H(α) and H(N) chemical shifts of the same set of terminally blocked dipeptides that are free from structural effects like secondary structure, hydrogen-bond, long-range backbone, and side-chain interactions. In particular, the preceding and following residue effects on amino-acid backbone conformational propensities are revealed and directly compared with previous works on either short peptides or empirical chemical shift database.

The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development

The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.

Identification of polyproline II regions derived from the proline-rich nuclear receptor coactivators PNRC and PNRC2: new insights for ERα coactivator interactions

Protein-protein interactions are crucial for signal transductions required for cell differentiation and proliferation. Their modulation is therefore key to the development of therapeutic alternatives, particularly in the context of cancer. According to literature data, the polyproline-rich nuclear receptor coactivators PNRC and PNRC2 interact with estrogen receptor (ERα) through their PxxP SH3-binding motifs. In a search to identify the molecular features governing this interaction, we explored using electronic circular dichroism (ECD) spectroscopy and molecular dynamics (MD) calculations, the capacity of a range of putative biologically active peptides derived from these proteins and containing this PxxP motif(s) to form polyproline II (PPII) domains. An additional more exhaustive structural study on a lead PPII peptide was also performed using 2D nuclear magnetic resonance (NMR) spectroscopy. With the exception of one of all the investigated peptides (PNRC-D), binding assays failed to detect any affinity for Grb2 SH3 domains, suggesting that PPII motifs issued from Grb2 antagonists have a binding mode distinct from those derived from Grb2 agonists. Instead, the peptides revealed a competitive binding ability against a synthetic peptide (ERα17p) with a putative PPII-cognate domain located within a coregulator recruitment region of ERα (AF-2 site). Our work, which constitutes the first structure-related interaction study concerning PNRC and PNRC2, supports not only the existence of PxxP-induced PPII sequences in these coregulators, but also confirms the presence of a PPII recognition site in the AF-2 of the steroid receptor ERα, a region important for transcription regulation.

pH-Independence of trialanine and the effects of termini blocking in short peptides: a combined vibrational, NMR, UVCD, and molecular dynamics study

Several lines of evidence now well establish that unfolded peptides in general, and alanine in specific, have an intrinsic preference for the polyproline II (pPII) conformation. Investigation of local order in the unfolded state is, however, complicated by experimental limitations and the inherent dynamics of the system, which has in some cases yielded inconsistent results from different types of experiments. One method of studying these systems is the use of short model peptides, and specifically short alanine peptides, known for predominantly sampling pPII structure in aqueous solution. Recently, He et al. ( J. Am. Chem. Soc. 2012 , 134 , 1571 - 1576 ) proposed that unblocked tripeptides may not be suitable models for studying conformational propensities in unfolded peptides due to the presence of end effect, that is, electrostatic interactions between investigated amino acid residues and terminal charges. To determine whether changing the protonation states of the N- and C-termini influence the conformational manifold of the central amino acid residue in tripeptides, we have examined the pH-dependence of unblocked trialanine and the conformational preferences of alanine in the alanine dipeptide. To this end, we measured and globally analyzed amide I' band profiles and NMR J-coupling constants. We described conformational distributions as the superposition of two-dimensional Gaussian distributions assignable to specific subspaces of the Ramachandran plot. Results show that the conformational ensemble of trialanine as a whole, and the pPII content (χpPII = 0.84) in particular, remains practically unaffected by changing the protonation state. We found that compared to trialanine, the alanine dipeptide has slightly lower pPII content (χpPII = 0.74) and an ensemble more reminiscent of the unblocked Gly-Ala-Gly model peptide. In addition, a two-state thermodynamic analysis of the conformational sensitive Δε(T) and (3)J(H(N)H(α))(T) data obtained from electronic circular dichroism and H NMR spectra indicate that the free energy landscape of trialanine is similar in all protonation states. MD simulations for the investigated peptides corroborate this notion and show further that the hydration shell around unblocked trialanine is unaffected by the protonation/deprotonation of the C-terminal group. In contrast, the alanine dipeptide shows a reduced water density around the central residue as well as a less ordered hydration shell, which decreases the pPII propensity and reduces the lifetime of sampled conformations.

Conformational distributions of denatured and unstructured proteins are similar to those of 20 × 20 blocked dipeptides

Understanding intrinsic conformational preferences of amino-acids in unfolded proteins is important for elucidating the underlying principles of their stability and re-folding on biological timescales. Here, to investigate the neighbor interaction effects on the conformational propensities of amino-acids, we carried out (1)H NMR experiments for a comprehensive set of blocked dipeptides and measured the scalar coupling constants between alpha protons and amide protons as well as their chemical shifts. Detailed inspection of these NMR properties shows that, irrespective of amino-acid side-chain properties, the distributions of the measured coupling constants and chemical shifts of the dipeptides are comparatively narrow, indicating small variances of their conformation distributions. They are further compared with those of blocked amino-acids (Ac-X-NHMe), oligopeptides (Ac-GGXGG-NH(2)), and native (lysozyme), denatured (lysozyme and outer membrane protein X from Escherichia coli), unstructured (Domain 2 of the protein 5A of Hepatitis C virus), and intrinsically disordered (hNlg3cyt: intracellular domain of human NL3) proteins. These comparative investigations suggest that the conformational preferences and local solvation environments of the blocked dipeptides are quite similar to not only those of other short oligopeptides but also those of denatured and natively unfolded proteins.

Phosphorylation alters backbone conformational preferences of serine and threonine peptides

Despite the notion that a control of protein function by phosphorylation works mainly by inducing its conformational changes, the phosphorylation effects on even small peptide conformation have not been fully understood yet. To study its possible effects on serine and threonine peptide conformations, we recently carried out pH- and temperature-dependent circular dichroism (CD) as well as (1)H NMR studies of the phosphorylated serine and threonine peptides and compared them with their unphosphorylated analogs. In the present article, by performing the self-consistent singular value decomposition analysis of the temperature-dependent CD spectra and by analyzing the (3)J(H(N),H(α)) coupling constants extracted from the NMR spectra, the populations of the polyproline II (PPII) and β-strand conformers of the phosphorylated Ser and Thr peptides are determined. As temperature is increased, the β-strand populations of both phosphorylated serine and threonine peptides increase. However, the dependences of PPII/β-strand population ratio on pH are different for these two cases. The phosphorylation of the serine peptide enhances the PPII propensity, whereas that of the threonine peptide has the opposite effect. This suggests that the serine and threonine phosphorylations can alter the backbone conformational propensity via direct but selective intramolecular hydrogen-bonding interactions with the peptide N--H groups. This clearly indicates that the phosphoryl group actively participates in modulating the peptide backbone conformations.