Aib-rich peptides containing lactam-bridged side chains as models of the 3(10)-helix.

Aib-rich side chain lactam-bridged oligomers with n =1, 2, 3, were designed and synthesized as putative models of the 3(10)-helix. These peptides were conformationally characterized in aqueous solution containing SDS micelles by CD, NMR, and computer simulations. The lactam bridge between the side chains of L-Glu and L-Lys in (i) and (i+3) positions was introduced in order to enhance the conformational preference toward the right-handed 3(10)-helix. The NMR results clearly indicate that there is an increase of 3(10)-helix formation upon chain elongation. In the dimer and trimer (n = 2 and n = 3, respectively, in the structure reported above) the observed NOE connectivities are compatible with the 3(10)-helical arrangement, confirmed by the temperature coefficients of the amide proton resonances which suggest the presence of a hydrogen-bonded structure. The phi and psi dihedral angles of the structures obtained by molecular dynamics calculations are also compatible with the 3(10)-helix. Identification of the hydrogen-bond pattern indicate that C=O(i)- - -HN(i+3) hydrogen bonds, typical of the 3(10)-helical conformation, are highly probable in all low-energy structures. The CD spectra of these Aib-rich lactam-bridged oligopeptides, obtained in the same solvent system used for NMR experiments, provide important insight into the spectroscopic characteristics of the 3(10)-helix.     

Design of β-Amino Acid with Backbone-Side Chain Interactions: Stabilization of 14/15-Helix in α/β-Peptides

A new C-linked carbo-β-amino acid, (R)-β-Caa((r)), having a carbohydrate side chain with d-ribo configuration, was prepared from d-glucose by inverting the C-3 stereocenter to introduce constraints/interactions. From the NMR studies it was inferred that the new monomer may participate in additional electrostatic interactions, facilitating and enhancing novel folds in oligomeric peptides derived from it. The α/β-peptides, synthesized from alternating l-Ala and (R)-β-Caa((r)), have shown the presence of 14/15-helix by NMR (in CDCl(3), methanol-d(3) and CD(3)CN), CD and MD calculations. The hybrid peptides showed the presence of electrostatic interactions involving the intraresidue amide proton and the C3-OMe, which helped in the stabilization of the NH(i)···CO(i-4) H-bonds and adoption of 14/15-helix. The importance of such additional interactions has been well defined in recent times to stabilize the folding in a variety of peptidic foldamers. These observations suggest and emphasize that the side chain-backbone interactions are crucial in the stabilization of the desired folding propensity. The designed monomer thus enlarges the opportunities for the synthesis of peptides with novel conformations and expands the repertoire of the foldamers.     

Molecular dynamics study of the solvation of an alpha-helical transmembrane peptide by DMSO

A 10-ns molecular dynamics study of the solvation of a hydrophobic transmembrane helical peptide in dimethyl sulfoxide (DMSO) is presented. The objective is to analyze how this aprotic polar solvent is able to solvate three groups of amino acid residues (i.e., polar, apolar, and charged) that are located in a stable helical region of a transmembrane peptide. The 25-residue peptide (sMTM7) used mimics the cytoplasmic proton hemichannel domain of the seventh transmembrane segment (TM7) from subunit a of H(+)-V-ATPase from Saccharomyces cerevisiae. The three-dimensional structure of peptide sMTM7 in DMSO has been previously solved by NMR spectroscopy. The radial and spatial distributions of the DMSO molecules surrounding the peptide as well as the number of hydrogen bonds between DMSO and the side chains of the amino acid residues involved are extracted from the molecular dynamics simulations. Analysis of the molecular dynamics trajectories shows that the amino acid side chains are fully embedded in DMSO. Polar and positively charged amino acid side chains have dipole-dipole interactions with the oxygen atom of DMSO and form hydrogen bonds. Apolar residues become solvated by DMSO through the formation of a hydrophobic pocket in which the methyl groups of DMSO are pointing toward the hydrophobic side chains of the residues involved. The dual solvation properties of DMSO cause it to be a good membrane-mimicking solvent for transmembrane peptides that do not unfold due to the presence of DMSO.     

Structural Investigation of Hybrid Peptide Foldamers Composed of α-Dipeptide Equivalent β-Oxy-δ5-amino Acids

Due to their equivalent lengths, δ-amino acids can serve as surrogates of α-dipeptides. However, δ-amino acids with proteinogenic side chains have not been well studied because of synthetic difficulties and because of their insolubility in organic solvents. Recently we reported the spontaneous supramolecular gelation of δ-peptides composed of β(O)-δ⁵-amino acids. Here, we report the incorporation of β(O)-δ⁵-amino acids as guests into the host α-helix, α,γ-hybrid peptide 12-helix and their single-crystal conformations. In addition, we studied the solution conformations of hybrid peptides composed of 1:1 alternating α and β(O)-δ⁵-amino acids. In contrast to the control α-helix structures, the crystal structure of peptides with β(O)-δ⁵-amino acids exhibit α-helical conformations consisting of both 13- and 10-membered H-bonds. The α,δ-hybrid peptide adopted mixed 13/11-helix conformation in solution with alternating H-bond directionality. Crystal-structure analysis revealed that the α,γ⁴-hybrid peptide accommodated the guest β(O)-δ⁵-amino acid without significant deviation to the overall helix folding. The results reported here emphasize that β(O)-δ⁵-amino acids with proteinogenic side chains can be accommodated into regular α-helix or 12-helix as guests without much deviation of the overall helix folding of the peptides.     

Dipole-induced self-assembly of helical beta-peptides

In this work, the interactions between beta-peptides are investigated for helix-forming peptides using molecular simulation. The role of electrostatic interactions in the self-assembly of these peptides is studied by calculating the dipole moment of various 14-helical beta-peptides using molecular dynamics simulations. The stability of a beta-peptide that is known to form a liquid crystalline phase is determined by calculating the potential of mean force using the expanded ensemble density of states method. This peptide is found to form a mechanically stable 14-helix in an implicit solvent model. The interaction between two of these peptides is examined by calculating the potential of mean force between the two peptides in implicit solvent. The peptides are shown to favorably associate in an end-to-end manner, driven largely by dipolar interactions. In order to understand the possible structures that form when many peptides interact in solution, a coarse-grained model is developed. Brownian dynamics simulations of the coarse-grained model at intermediate concentrations (1-50 mM) are performed, and the aggregation behavior is understood by calculating the diffusivity and the radial distribution function. An analysis of the resultant structures reveals that the coarse-grained model of the peptide leads to the formation of spherical clusters.     

The balance between side-chain and backbone-driven association in folding of the α-helical influenza A transmembrane peptide

The correct balance between attractive, repulsive and peptide hydrogen bonding interactions must be attained for proteins to fold correctly. To investigate these important contributors, we sought a comparison of the folding between two 25-residues peptides, the influenza A M2 protein transmembrane domain (M2TM) and the 25-Ala (Ala₂₅). M2TM forms a stable α-helix as is shown by circular dichroism (CD) experiments. Molecular dynamics (MD) simulations with adaptive tempering show that M2TM monomer is more dynamic in nature and quickly interconverts between an ensemble of various α-helical structures, and less frequently turns and coils, compared to one α-helix for Ala₂₅. DFT calculations suggest that folding from the extended structure to the α-helical structure is favored for M2TM compared with Ala₂₅. This is due to CH⋯O attractive interactions which favor folding to the M2TM α-helix, and cannot be described accurately with a force field. Using natural bond orbital (NBO) analysis and quantum theory atoms in molecules (QTAIM) calculations, 26 CH⋯O interactions and 22 NH⋯O hydrogen bonds are calculated for M2TM. The calculations show that CH⋯O hydrogen bonds, although individually weaker, have a cumulative effect that cannot be ignored and may contribute as much as half of the total hydrogen bonding energy, when compared to NH⋯O, to the stabilization of the α-helix in M2TM. Further, a strengthening of NH⋯O hydrogen bonding interactions is calculated for M2TM compared to Ala₂₅. Additionally, these weak CH⋯O interactions can dissociate and associate easily leading to the ensemble of folded structures for M2TM observed in folding MD simulations.     

The role of cation-pi interactions in biomolecular association. Design of peptides favoring interactions between cationic and aromatic amino acid side chains.

Cation-pi interactions between amino acid side chains are increasingly being recognized as important structural and functional features of proteins and other biomolecules. Although these interactions have been found in static protein structures, they have not yet been detected in dynamic biomolecular systems. We determined, by (1)H NMR spectroscopic titrations, the energies of cation-pi interactions of the amino acid derivative AcLysOMe (1) with AcPheOEt (2) and with AcTyrOEt (3) in aqueous and three organic solvents. The interaction energy is substantial; it ranges from -2.1 to -3.4 kcal/mol and depends only slightly on the dielectric constant of the solvent. To assess the effects of auxiliary interactions and structural preorganization on formation of cation-pi interactions, we studied these interactions in the association of pentapeptides. Upon binding of the positively-charged peptide AcLysLysLysLysLysNH(2) (5) to the negatively-charged partner AcAspAspXAspAspNH(2) (6), in which X is Leu (6a), Tyr (6b), and Phe (6c), multiple interactions occur. Association of the two pentapeptides is dynamic. Free peptides and their complex are in fast exchange on the NMR time-scale, and 2D (1)H ROESY spectra of the complex of the two pentapeptides do not show intermolecular ROESY peaks. Perturbations of the chemical shifts indicated that the aromatic groups in peptides 6b and 6c were affected by the association with 5. The association constants K(A) for 5 with 6a and with 6b are nearly equal, (4.0 +/- 0.7) x 10(3) and (5.0 +/- 1.0) x 10(3) M(-)(1), respectively, while K(A) for 5 with 6c is larger, (8.3 +/- 1.3) x 10(3) M(-)(1). Molecular-dynamics (MD) simulations of the pentapeptide pairs confirmed that their association is dynamic and showed that cation-pi contacts between the two peptides are stereochemically possible. A transient complex between 5 and 6 with a prominent cation-pi interaction, obtained from MD simulations, was used as a template to design cyclic peptides C(X) featuring persistent cation-pi interactions. The cyclic peptide C(X) had a sequence in which X is Tyr, Phe, and Leu. The first two peptides do, but the third does not, contain the aromatic residue capable of interacting with a cationic Lys residue. This covalent construct offered conformational stability over the noncovalent complexes and allowed thorough studies by 2D NMR spectroscopy. Multiple conformations of the cyclic peptides C(Tyr) and C(Phe) are in slow exchange on the NMR time-scale. In one of these conformations, cation-pi interaction between Lys3 and Tyr9/Phe9 is clearly evident. Multiple NOEs between the side chains of residues 3 and 9 are observed; chemical-shift changes are consistent with the placement of the side chain of Lys3 over the aromatic ring. In contrast, the cyclic peptide C(Leu) showed no evidence for close approach of the side chains of Lys3 and Leu9. The cation-pi interaction persists in both DMSO and aqueous solvents. When the disulfide bond in the cyclic peptide C(Phe) was removed, the cation-pi interaction in the acyclic peptide AC(Phe) remained. To test the reliability of the pK(a) criterion for the existence of cation-pi interactions, we determined residue-specific pK(a) values of all four Lys side chains in all three cyclic peptides C(X). While NOE cross-peaks and perturbations of the chemical shifts clearly show the existence of the cation-pi interaction, pK(a) values of Lys3 in C(Tyr) and in C(Phe) differ only marginally from those values of other lysines in these dynamic peptides. Our experimental results with dynamic peptide systems highlight the role of cation-pi interactions in both intermolecular recognition at the protein-protein interface and intramolecular processes such as protein folding.     

Numerical Simulation of the Effect of Solvent Viscosity on the Motions of a β‐Peptide Heptamer

This report examines the effect of a decrease in solvent viscosity on the simulated folding behaviour of a β‐peptide heptamer in methanol. Simulations of the molecular dynamics of the heptamer H‐β³‐HVal‐β³‐HAla‐β³‐HLeu‐(S,S)‐β³‐HAla(αMe)‐β³‐HVal‐β³‐HAla‐β³‐HLeu‐OH in methanol, with an explicit representation of the methanol molecules, were performed for 80 ns at various solvent viscosities. The simulations indicate that at a solvent viscosity of one third of that of methanol, only the dynamic aspects of the folding process are altered, and that the rate of folding is increased. At a viscosity of one tenth of that of methanol, insufficient statistics are obtained within the 80 ns period. We suggest that 80 ns is an insufficient time to reach conformational equilibrium at very low viscosity because the dependence of the folding rate of a β‐peptide on solvent viscosity has two regimes; a result that was observed in another computational study for α‐peptides.     

Do Valine Side Chains Have an Influence on the Folding Behavior of β‐Substituted β‐Peptides?

The influence of valine side chains on the folding/unfolding equilibrium and, in particular, on the 3₁₄‐helical propensity of β³‐peptides were investigated by means of molecular‐dynamics (MD) simulation. To that end, the valine side chains in two different β³‐peptides were substituted by leucine side chains. The resulting four peptides, of which three have never been synthesized, were simulated for 150 to 200 ns at 298 and 340 K, starting from a fully extended conformation. The simulation trajectories obtained were compared with respect to structural preferences and folding behavior. All four peptides showed a similar folding behavior and were found to predominantly adopt 3₁₄‐helical conformations, irrespective of the presence of valine side chains. No other well‐defined conformation was observed at significant population in any of the simulations. Our results imply that β³‐peptides show a structural preference for 3₁₄‐helices independent of the branching nature of the side chains, in contrast to what has been previously proposed on the basis of circular‐dichroism (CD) measurements.     

Evaluating force field accuracy with long-time simulations of a β-hairpin tryptophan zipper peptide

We have combined graphics processing unit-accelerated all-atom molecular dynamics with parallel tempering to explore the folding properties of small peptides in implicit solvent on the time scale of microseconds. We applied this methodology to the synthetic β-hairpin, trpzip2, and one of its sequence variants, W2W9. Each simulation consisted of over 8 μs of aggregated virtual time. Several measures of folding behavior showed good convergence, allowing comparison with experimental equilibrium properties. Our simulations suggest that the intramolecular interactions of tryptophan side chains are responsible for much of the stability of the native fold. We conclude that the ff99 force field combined with ff96 φ and ψ dihedral energies and an implicit solvent can reproduce plausible folding behavior in both trpzip2 and W2W9.     

Solvent polarity controls the helical conformation of short peptides rich in Calpha-tetrasubstituted amino acids

The two peptides, rich in C(alpha)-tetrasubstituted amino acids, Ac-[Aib-L-(alphaMe)Val-Aib](2)-L-His-NH(2) (1) and Ac-[Aib-L-(alphaMe)Val-Aib](2)-O-tBu (2 a) are prevalently helical. They present the unique property of changing their conformation from the alpha- to the 3(10)-helix as a function of the polarity of the solvent: alpha in more polar solvents, 3(10) in less polar ones. Conclusive evidence of this reversible change of conformation is reported on the basis of the circular dichroism (CD) spectra and a detailed two-dimensional NMR analysis in two solvents (trifluoroethanol and methanol) refined with molecular dynamics calculations. The X-ray diffractometric analysis of the crystals of both peptides reveals that they assume a prevalent 3(10)-helix conformation in the solid state. This conformation is practically superimposable on that obtained from the NMR analysis of 1 in methanol. The NMR results further validate the reported CD signature of the 3(10)-helix and the use of the CD technique for its assessment.     

Chain-length-dependent helical motifs and self-association of beta-peptides with constrained side chains

Homo-oligomers constructed by using trans-2-aminocyclohexanecarboxylic acid monomers without protecting groups were studied. Both ab initio theory and NMR measurements showed that the tetramer tends to adopt a 10-helix motif, while the pentamer and hexamer form the known 14-helix. It was concluded that the conformationally constrained backbone is flexible enough to afford both 10-helical and 14-helical motifs, this observation in turn providing evidence of the true folding process. Self-association of the helical units was also detected, and the results of variable-temperature diffusion NMR measurements strongly suggested the presence of helical bundles in methanol solution.     

Relationship between side chain structure and 14-helix stability of beta3-peptides in water

Folded polymers are used in Nature for virtually every vital process. Nonnatural folded polymers, or foldamers, have the potential for similar versatility, and the design and refinement of such molecules is of considerable current interest. Here we report a complete and systematic analysis of the relationship between side chain structure and the 14-helicity of a well-studied class of foldamers, beta(3)-peptides, in water. Our experimental results (1) verify the importance of macrodipole stabilization for maintaining 14-helix structure, (2) provide comprehensive evidence that beta(3)-amino acids branched at the first side chain carbon are 14-helix-stabilizing, (3) suggest a novel role for side chain hydrogen bonding as an additional stabilizing force in beta(3)-peptides containing beta(3)-homoserine or beta(3)-homothreonine, and (4) demonstrate that diverse functionality can be incorporated into a stable 14-helix. Gas- and solution-phase calculations and Monte Carlo simulations recapitulate the experimental trends only in the context of oligomers, yielding insight into the mechanisms behind 14-helix folding. The 14-helix propensities of beta(3)-amino acids differ starkly from the alpha-helix propensities of analogous alpha-amino acids. This contrast informs current models for alpha-helix folding, and suggests that 14-helix folding is governed by different biophysical forces than is alpha-helix folding. The ability to modulate 14-helix structure through side chain choice will assist rational design of 14-helical beta-peptide ligands for macromolecular targets.     

De Novo Left‐Handed Synthetic Peptidomimetic Foldamers

The development of peptidomimetic helical foldamers with a wide repertoire of functions is of significant interest. Herein, we report the X‐ray crystal structures of a series of homogeneous l ‐sulfono‐γ‐AA foldamers and elucidate their folding conformation at the atomic level. Single‐crystal X‐ray crystallography revealed that this class of oligomers fold into unprecedented dragon‐boat‐shaped and unexpected left‐handed helices, which are stabilized by the 14‐hydrogen‐bonding pattern present in all sequences. These l ‐sulfono‐γ‐AApeptides have a helical pitch of 5.1 Å and exactly four side chains per turn, and the side chains lie perfectly on top of each other along the helical axis. 2D NMR spectroscopy, computational simulations, and CD studies support the folding conformation in solution. Our results provide a structural basis at the atomic level for the design of novel biomimetics with a precise arrangement of functional groups in three dimensions.     

Establishing effective simulation protocols for beta- and alpha/beta-peptides. II. Molecular mechanical (MM) model for a cyclic beta-residue

All-atom molecular mechanical (MM) force field parameters are developed for a cyclic beta-amino acid, amino-cyclo-pentane-carboxylic acid (ACPC), using a multi-objective evolutionary algorithm. The MM model is benchmarked using several short, ACPC-containing alpha/beta-peptides in water and methanol with SCC-DFTB (self consistent charge-density functional tight binding)/MM simulations as the reference. Satisfactory agreements are found between the MM and SCC-DFTB/MM results regarding the distribution of key dihedral angles for the tetra-alpha/beta-peptide in water. For the octa-alpha/beta-peptide in methanol, the MM and SCC-DFTB/MM simulations predict the 11- and 14/15-helical form as the more stable conformation, respectively; however, the two helical forms are very close in energy (2-4 kcal/mol) at both theoretical levels, which is also the conclusion from recent NMR experiments. As the first application, the MM model is applied to an alpha/beta-pentadeca-peptide in water with both explicit and implicit solvent models. The stability of the peptide is sensitive to the starting configuration in the explicit solvent simulations due to their limited length ( approximately 10-40 ns). Multiple ( approximately 20 x 20 ns) implicit solvent simulations consistently show that the 14/15-helix is the predominant conformation of this peptide, although substantially different conformations are also accessible. The calculated nuclear Overhauser effect (NOE) values averaged over different trajectories are consistent with experimental data, which emphasizes the importance of considering conformational heterogeneity in such comparisons for highly dynamical peptides.     

Helical Folding of Conjugated Oligo(phenyleneethynylene): Chain‐Length Dependence,Solvent Effects,and Intermolecular Assembly

As a representative folding system that features a conjugated backbone, a series of monodispersed (o‐phenyleneethynylene)‐alt‐(p‐phenyleneethynylene) (PE) oligomers of varied chain length and different side chains were studied. Molecules with the same backbone but different side‐chain structures were shown to exhibit similar helical conformations in respectively suitable solvents. Specifically, oligomers with dodecyloxy side chains folded into the helical structure in apolar aliphatic solvents, whereas an analogous oligomer with tri(ethylene glycol) (Tg) side chains adopted the same conformation in polar solvents. The fact that the oligomers with the same backbone manifested a similar folded conformation independent of side chains and the nature of the solvent confirmed the concept that the driving force for folding was the intramolecular aromatic stacking and solvophobic interactions. Although all were capable of inducing folding, different solvents were shown to bestow slightly varied folding stability. The chain‐length dependence study revealed a nonlinear correlation between the folding stability with backbone chain length. A critical size of approximately 10 PE units was identified for the system, beyond which folding occurred. This observation corroborated the helical nature of the folded structure. Remarkably, based on the absorption and emission spectra, the effective conjugation length of the system extended more effectively under the folded state than under random conformations. Moreover, as evidenced by the optical spectra and dynamic light‐scattering studies, intermolecular association took place among the helical oligomers with Tg side chains in aqueous solution. The demonstrated ability of such a conjugated foldamer in self‐assembling into hierarchical supramolecular structures promises application potential for the system.     

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.     

Effect of hydrophobic amino acids on the conformational change of decapeptides in micellar environments

A series of peptides containing various hydrophobic amino acids [methionine (Met), leucine (Leu), norleucine (Nle), phenylalanine (Phe), 2-aminooctanoic acid (Aoc), and 2-aminodecanoic acid (Ade)] were synthesized and their conformations were studied using circular dichroism (CD) spectroscopy in different solvents such as water, methanol, and aqueous solution of ammonium tetradecanesulfonate. Peptides containing hydrophobic amino acids with linear side chains formed -sheets in water and methanol. Electrostatic interaction between the charged side chain (lysine) and a micelle consisting of an anionic surfactant, ammonium tetradecanesulfonate, is necessary for the formation of -helices in micellar environments. The conformational transition from -helix to -sheet structure required moderate hydrophobicity and linear side chains. This conformational transition depended on the surfactant concentration.     

Interpreting NMR data for beta-peptides using molecular dynamics simulations

NMR is one of the most used techniques to resolve structure of proteins and peptides in solution. However, inconsistencies may occur due to the fact that a polypeptide may adopt more than one conformation. Since the NOE distance bounds and (3)J-values used in such structure determination represent a nonlinear average over the total ensemble of conformers, imposition of NOE or (3)J-value restraints to obtain one unique conformation is not an appropriate procedure in such cases. Here, we show that unrestrained MD simulation of a solute in solution using a high-quality force field yields a conformational ensemble that is largely compatible with the experimental NMR data on the solute. Four 100 ns MD simulations of two forms of a nine-residue beta-peptide in methanol at two temperatures produced conformational ensembles that were used to interpret the NMR data on this molecule and resolve inconsistencies between the experimental NOEs. The protected and unprotected forms of the beta-peptide adopt predominantly a 12/10-helix in agreement with the qualitative interpretation of the NMR data. However, a particular NOE was not compatible with this helix indicating the presence of other conformations. The simulations showed that 3(14)()-helical structures were present in the ensemble of the unprotected form and that their presence correlates with the fulfillment of the particular NOE. Additionally, all inter-hydrogen distances were calculated to compare NOEs predicted by the simulations to the ones observed experimentally. The MD conformational ensembles allowed for a detailed and consistent interpretation of the experimental data and showed the small but specific conformational differences between the protected and unprotected forms of the peptide.     

Interactions of peptide mimics of hyaluronic acid with the receptor for hyaluronan mediated motility (RHAMM)

Summary Using the hyaluronic acid (HA) binding region of the receptor for hyaluronan-mediated motility (RHAMM) as a model, a molecular perspective for peptide mimicry of the natural ligand was established by comparing the interaction sites of HA and unnatural peptide–ligands to RHAMM. This was accomplished by obtaining a series of octapeptide–ligands through screening experiments that bound to the HA binding domains of RHAMM (amino acids 517–576) and could be displaced by HA. These molecules were computationally docked onto a three-dimensional NMR based model of RHAMM. The NMR model showed that RHAMM(517–576) was a set of three helices, two of which contained the HA binding domains (HABDs) flanking a central groove. The structure was stabilized by hydrophobic interactions from four pairs of Val and Ile side chains extending into the groove. The presence of solvent exposed, positively charged side chains spaced 11Å apart matched the spacing of negative charges on HA. Docking experiments using flexible natural and artificial ligands demonstrated that HA and peptide–mimetics preferentially bound to the second helix that contains HABD-2. Three salt bridges between HA carboxylates and Lys548, Lys553 and Lys560 and two hydrophobic interactions involving Val538 and Val559 were predicted to stabilize the RHAMM-HA complex. The high affinity peptides and HA utilized the same charged residues, with additional contacts to other basic residues. However, hydrophobic contacts do not contribute to affinity for peptide ligand-RHAMM complexes. These results offer insight into how selectivity is achieved in the binding of HA to RHAMM, and how peptide competitors may compete for binding with HA on a single hyaladherin.