Influence of hydration on the conformational stability and formation of bends in terminally blocked dipeptides

The conformations of 23 terminally blocked dipeptide sequences were examined by conformational energy calculations that included the effects of the aqueous solvent. Starting structures were derived from combinations of minimum-energy conformations of hydrated single residues. Their conformational energies were then minimized using the ECEPP potential (Empirical Conformational Energy Program for Peptides) with hydration included. Short-range interactions dominate in stabilizing the conformations of the hydrated dipeptides. Differences between conformational stabilities of hydrated and unhydrated dipeptides in many cases are due to the competition of solute–water and intramolecular hydrogen bonds. In other cases, perturbation of the hydration shell of the solute by close approach of solute atoms alters conformational preferences. Probabilities of formation of bends were calculated and compared to the corresponding quantities for unhydrated dipeptides and to those calculated from x-ray structures. For bends in dipeptides containing two nonpolar amino acids, computations omitting hydration yield better results. However, better agreement with experimental (x-ray) bend probabilities for dipeptides containing glycine or polar amino acids is obtained only in some sequences when hydration is included. The results are rationalized by the observation that, in proteins, bends containing nonpolar sequences occur on the inside, shielded from the solvent. Bends containing glycine or polar amino acids occur frequently on the surface of the protein, but they are not completely hydrated.     

Solution NMR structure of a D,L-alternating oligonorleucine as a model of beta-helix

beta-Helix structures are of particular interest due to their capacity to form transmembrane channels with different transport properties. However, the relatively large number of beta-helices configurations does not allow a direct conformational analysis of beta-helical oligopeptides. A synthetic alternating D,L-oligopeptide with twelve norleucines (XIIMe) has been used as a model to get insight in the conformational features of beta-helix structures. The spatial configuration of XIIMe in solution has been determined by NMR. An extensive set of distances (nuclear Overhauser effect) and dihedral (J coupling constants) constraints have been included in molecular dynamics calculations. The NMR experimental data and theoretical calculations clearly indicate that the XIIMe adopts a single beta(4.4)-helix-type conformation in nonpolar solvents.     

Stabilization of an alpha-helical conformation in an isolated hexapeptide inhibitor of calmodulin.

The conformational properties of two hexapeptides, Ac-LWRILW-NH(2) and its D-amino acid counterpart Ac-lwrilw-NH(2), identified as calmodulin inhibitors using mixture-based synthetic combinatorial library approaches, have been characterised by NMR and CD spectroscopy. The peptides fold into an alpha-helical conformation in aqueous solution. The observed short- and medium-range nuclear Overhauser effects were consistent with the formation of an alpha-helical structure and a reasonably well-defined set of structures was obtained by using restraints from the NMR data in simulated annealing calculations. Analysis of glycine-substitution analogues demonstrated that all the amino acids that make up the peptide sequence are important for the stabilization of the alpha-helical conformation. The results suggest that a well-defined set of interactions is indispensable to allow alpha-helix formation in this short hexapeptide.     

Conformational behaviour of glycomimetics: NMR and molecular modelling studies of the C-glycoside analogue of the disaccharide methyl beta-D-galactopyranosyl-(1-->3)-beta-D-glucopyranoside

The conformational behaviour of the C-glycoside beta-C-Gal-(1-->3)-beta-Glc-OMe (1) has been studied using a combination of molecular mechanics and NMR spectroscopy (proton-proton coupling constants and nuclear Overhauser effects). It is shown that the C-disaccharide populates two distinctive conformational families in solution, the normal syn-psi conformation, which is the predominating conformation of parent O-glycoside 2, and the anti-psi conformation, which has not been detected for the O-disaccharide.     

Solution conformation of the C-terminal domain of skeletal troponin C. Cation, trifluoperazine and troponin I binding effects

Proton magnetic resonance spectroscopy has been used to study the cation (Mg2+, Ca2+)-dependent conformational states of the C-terminal domain of rabbit skeletal troponin C under a variety of solution conditions. Nuclear Overhauser data and paramagnetic probe observations provide definition of the configuration of this region of troponin C. Comparative study of homologous proteins identify common features of the tertiary structure relevant to the cation binding reaction. Complex formation with troponin I and the drug trifluoperazine is observed to adjust the solution conformation of the C-terminal domain of troponin C. The interactive conformational response to cation coordination and the binding of the drug and troponin I are discussed.     

Synthesis and conformational studies by 1H- and 13C-NMR spectroscopy of a novel, sterically constrained analogue of thyrotropin-releasing hormone

A novel thyrotropin-releasing hormone (TRH) analogue, [2,4-MePro3]-TRH (2,4-MePro: 2-carboxy-2,4-methanopyrrolidine), has been synthesized using a rapid solid phase peptide synthesis method, and its conformational properties investigated by one- and two-dimensional (2D) nmr spectroscopy and by proton Overhauser measurements. Following a published approach, calibrated interproton Overhauser effects were used together with distance geometry analysis to deduce that the single conformation of the His-2,4-MePro tertiary amide bond is trans in aqueous solution. This conclusion was corroborated by 2D dipolar-correlated (NOESY) spectroscopy. A preferentially extended conformation is indicated by the nmr data, similar to that of TRH. The phi, psi conformational space of 2,4-MePro is, however, significantly different from that of trans proline and the structural consequences of these differences at the C-terminus are discussed. The distribution of histidine side-chain conformations in the TRH analogue was deduced from coupling constants and from the short-range interaction between the imidazole ring and one of the prochiral faces of the 2,4-MePro side chain.     

Reconstruction of NOESY maps. A requirement for a reliable conformational analysis of biomolecules using 2D NMR

The modelling of the conformation of a biomolecule in solution is based mainly on the internuclear distances deduced from measurements of nuclear Overhauser effects (nOe) in NOESY correlation maps. The distances are then used as restraints in the energy minimization procedure, which leads to one or several optimized conformations. A general and safe technique for validating these structures with respect to the experimental data is here proposed: from the internuclear distances, the relaxation matrix can be computed under the assumption of a unique rotational correlation time. By stepwise integration of these relaxation equations, the NOESY maps can be accurately reconstructed for any mixing time. Because multi-spin effects are correctly taken into account, any difference between the experimental and theoretical maps can be easily interpreted in terms of conformation, and possible inconsistencies due to conformational averaging can be pointed out. The technique is illustrated for a bacterial lipopeptide, mycosubtilin, the spectrum of which is completely assigned.     

NMR-derived model of interconverting conformations of an ICAM-1 inhibitory cyclic nonapeptide.

We have produced by phage-display a disulfide-linked cyclic nonapeptide (inhibitory peptide-01, IP01), CLLRMRSIC, that binds to intracellular adhesion molecule-1 (ICAM-1) and blocks binding to its counter-structure, leukocyte functional antigen-1 (LFA-1). As a first step towards improving its pharmacologic properties, we have performed a structural and functional analysis of this peptide inhibitor to determine the features relevant to ICAM-1 binding. We report here the solution model of our initial product, IP01, as derived from two-dimensional nuclear magnetic resonance (NMR) restraints and molecular modeling. Distance and dihedral angle restraints, generated from nuclear Overhauser effect spectroscopy (NOESY) and one-dimensional-NMR experiments respectively, were used to generate an ensemble of structures using distance geometry and simulated annealing. Molecular dynamic simulations produced three interconverting conformational families consistent with the NMR-derived constraints. We describe these conformations and their mechanism of interconversion. Furthermore, we have measured the IC50 s of a series of inhibitors generated from IP01 through alanine substitution of each residue. These results show that the L2-L3-R4-M5-R6 segment is functionally active, conformationally flexible, and contains a beta-turn involving residues R4-S7, while the C1-C9-I8-S7 segment is less functionally-active but adopts a more defined solution conformation, consistent with a scaffolding function. This model will be useful for designing nonpeptide-based organic inhibitors with improved pharmacologic properties.     

High-field NMR and circular dichroism solvent-dependent conformational studies of the bradykinin C-terminal tetrapeptide Ser-Pro-Phe-Arg

The conformational properties of the tetrapeptide Ser1-Pro2-Phe3-Arg4, the C-terminal fragment of the nonapeptide hormone bradykinin, have been studied by circular dichroism and two-dimensional NMR techniques. Measurements of coupling constants, NH temperature dependence rates and nuclear Overhauser effects (performed with rotating frame nuclear Overhauser spectroscopy, ROESY) in H2O and CD3OH/D2O (80/20, v/v) reveal different conformations in the corresponding solvent. In aqueous solution the molecule exists in a random conformation or as an average of several conformations in rapid exchange. In CD3OH/D2O, however, the conformation is well-defined. The backbone of the peptide is extended, and the side-chains of Phe3 and Arg4 exhibit unusual rigidity for a peptide of this size. Evidently, the secondary structure is stabilized by a charge interaction between the guanidino group of Arg4 and the terminal carboxyl group, since experiments at various pH's show clearly that the definition of conformation decreases strongly upon protonation of the carboxyl function. A NH3+(Ser1)-COO-(Arg4) salt bridge, as well as any form of turn stabilized by hydrogen bonds can be ruled out with certainty.     

Determination of three-dimensional solution structure of waglerin I,a toxin from Trimeresurus wagleri,using 2D-NMR and molecular dynamics simulation

The solution conformation of a synthetic snake venom toxin waglerin I, has been determined by using proton nuclear magnetic resonance spectroscopy. By y a combination of various two-dimensional NMR techniques, the ¹H-NMR spectrum of waglerin I was completely assigned. A set of 247 interproton distance restraints was derived from nuclear Overhauser enhancement (NOE) measurements. These NOE constraints, in addition to the 2 dihedral angle restraints (from coupling constant measurements) and 7 ω torsion angle restraints for prolines, formed the basis of three-dimensional structure determined by molecular dynamics techniques. The 19 structures that were obtained satisfy the experimental restraints, and display small deviation from idealized covalent geometry. Analysis of converged structures indicates that the toxin has no special secondary structure. In the solution structure of waglerin I, the central ring region is well defined but the N- and C-termini possesses more disorder.     

PEPT1 as a paradigm for membrane carriers that mediate electrogenic bidirectional transport of anionic,cationic, and neutral substrates

The capability for electrogenic inward transport of substrates that carry different net charge is a phenomenon observed in a variety of membrane-solute transporters but is not yet understood. We employed the two-electrode voltage clamp technique combined with intracellular pH recordings and the giant patch technique to assess the selectivity for bidirectional transport and the underlying stoichiometries in proton to substrate flux coupling for electrogenic transfer of selected anionic, cationic, and neutral dipeptides by the intestinal peptide transporter PEPT1. Anionic dipeptides such as Gly-Asp and Asp-Gly are transported in their neutral and negatively charged forms with high and low affinities, respectively. The positive transport current obtained with monoanionic substrates results from the cotransport of two protons. Cationic dipeptides can be transported in neutral and positively charged form, resulting in an excess transport current as compared with neutral substrates. However, binding and transport of cationic dipeptides shows a pronounced selectivity for the position of charged side chains demonstrating that the binding domain of PEPT1 is asymmetric, both in its inward and outward facing conformation. The simultaneous presence of identically charged substrates on both membrane surfaces generates outward and, unexpectedly, enhanced inward transport currents probably by increasing the turnover rate.     

Invertible micellar polymer assemblies for delivery of poorly water-soluble drugs

Strategically designed amphiphilic invertible polymers (AIPs) are capable of (i) self-assembling into invertible micellar assemblies (IMAs) in response to changes in polarity of environment, polymer concentration, and structure, (ii) accommodating (solubilizing) substances that are otherwise insoluble in water, and (iii) inverting their molecular conformation in response to changes in the polarity of the local environment. The unique ability of AIPs to invert the molecular conformation depending on the polarity of the environment can be a decisive factor in establishing the novel stimuli-responsive mechanism of solubilized drug release that is induced just in response to a change in the polarity of the environment. The IMA capability to solubilize lipophilic drugs and deliver and release the cargo molecules by conformational inversion of polymer macromolecules in response to a change of the polarity of the environment was demonstrated by loading IMA with a phytochemical drug, curcumin. It was demonstrated that four sets of micellar vehicles based on different AIPs were capable of delivering the curcumin from water to an organic medium (1-octanol) by means of unique mechanism: AIP conformational inversion in response to changing polarity from polar to nonpolar. The IMAs are shown to be nontoxic against human cells up to a concentration of 10 mg/L. On the other hand, the curcumin-loaded IMAs are cytotoxic to breast carcinoma cells at this concentration, which confirms the potential of IMA-based vehicles in controlled delivery of poorly water-soluble drug candidates and release by means of this novel stimuli-responsive mechanism.     

An ellipticine derivative (oxazolopyridocarbazolium) 3' linked to tetrathymidylate stacks intramolecularly with the nearest thymine at low concentration and head-to-tail intermolecularly at high concentration

The solution conformation of a tetrathymidylate linked through an ester bond to an ellipticine derivative oxazolopyridocarbazolium (OPC) at the 3' position was investigated using one- and two-dimensional nmr experiments. Since the total electric charge of the OPC ring may influence self-association, we first determined the pKa of the oxazole cyclic acidic function. Nuclear Overhauser effect spectroscopy experiments showed that, at low concentration, the OPC stacks intramolecularly with the nearest thymine at the 3' end. At highest concentration, however, the OPC rings are self-associated. The stacking constant was calculated using 1H chemical shift dilution experiment. The conformational model suggested by P-nmr was tested by molecular mechanics computations.     

Metropolis Monte Carlo calculations of DNA structure using internal coordinates and NMR distance restraints: An alternative method for generating a high-resolution solution structure

Summary A new method, a restrained Monte Carlo (rMC) calculation, is demonstrated for generating high-resolution structures of DNA oligonucleotides in solution from interproton distance restraints and bounds derived from complete relaxation matrix analysis of two-dimensional nuclear Overhauser effect (NOE) spectral peak intensities. As in the case of restrained molecular dynamics (rMD) refinement of structures, the experimental distance restraints and bounds are incorporated as a pseudo-energy term (or penalty function) into the mathematical expression for the molecular energy. However, the use of generalized helical parameters, rather than Cartesian coordinates, to define DNA conformation increases efficiency by decreasing by an order of magnitude the number of parameters needed to describe a conformation and by simplifying the potential energy profile. The Metropolis Monte Carlo method is employed to simulate an annealing process. The rMC method was applied to experimental 2D NOE data from the octamer duplex d(GTA-TAATG)·d(CATTATAC). Using starting structures from different locations in conformational space (e.g. A-DNA and B-DNA), the rMC calculations readily converged, with a root-mean-square deviation (RMSD) of <0.3 Å between structures generated using different protocols and starting structures. Theoretical 2D NOE peak intensities were calculated for the rMC-generated structures using the complete relaxation matrix program CORMA, enabling a comparison with experimental intensities via residual indices. Simulation of the vicinal proton coupling constants was carried out for the structures generated, enabling a comparison with the experimental deoxyribose ring coupling constants, which were not utilized in the structure determination in the case of the rMC simulations. Agreement with experimental 2D NOE and scalar coupling data was good in all cases. The rMC structures are quite similar to that refined by a traditional restrained MD approach (RMSD<0.5 Å) despite the different force fields used and despite the fact that MD refinement was conducted with additional restraints imposed on the endocyclic torsion angles of deoxyriboses. The computational time required for the rMC and rMD calculations is about the same. A comparison of structural parameters is made and some limitations of both methods are discussed with regard to the average nature of the experimental restraints used in the refinement.Abbreviations MC Monte Carlo - rMC restrained Monte Carlo - MD molecular dynamics - rMD restrained molecular dynamics - DG distance geometry - EM energy minimization - 2D NOE two-dimensional nuclear Overhauser effect - DQF-COSY double-quantum-filtered correlation spectroscopy - RMSD root-mean-square deviation To whom correspondence should be addressed.     

Nmr and molecular dynamics study of four carbocyclic muramyl dipeptide analogues

We have performed a conformational analysis of the carbocyclic muramyl dipeptide analogues (1′R, 2′R)- and (1′S, 2′S)-N-[2-(2′-acetamidocyclohexyIoxy)acetyl]-L Ala-D -iGln (-D -Glu) utilizing ¹H-nmr spectroscopy and nuclear Overhauser effect restrained molecular dynamics. Intramolecular H bonding for all four diastereoisomers is suggested by the Ala-NH temperature coefficients. Distance restraints were obtained by NOE spectroscopy and rotating frame NOE spectroscopy experiments. Structures with low potential energy and high agreement with NOE data were sought by restrained molecular dynamics. The ring configuration was found to induce conformational preferences. The β-like turn characterized by the intramolecular C₁₀ H-bond Ala-NH-acetamido-CO is preferred with the (1′S, 2′S), but appears to be less stable with (1′R, 2′R), diastereoisomers. Calculations show the double β-turn proposed for muramyl dipeptide [S. Fermandjian, B. Perly, M. Level, and P. Lefrancier (1987) Carbohydrate Research, Vol. 162, pp. 23–32] to have higher potential energy. © 1993 John Wiley & Sons, Inc.     

Conformational analysis of a IgG1 hinge peptide derivative in solution determined by NMR spectroscopy and refined by restrained molecular dynamics simulations.

The hinge region links the antigen binding Fab part to the constant Fc domain in immunoglobulins. For the hinge peptide derivative [AcThr(OtBu)-Cys-Pro-Pro-Cys-Pro-Ala-ProNH2]2 the assignment of the 1H and 13C resonances was achieved by two-dimensional nmr techniques: total correlation spectroscopy (TOCSY), nuclear Overhauser enhancement spectroscopy (NOESY), rotating frame nuclear Overhauser enhancement spectroscopy (ROESY), heteronuclear multiple quantum coherence (HMQC) transfer, and a HSQC (modified Overbodenhausen experiment) with high resolution in F1, which was several times folded in F1 but still phase correctable. Conformational relevant parameters (78 nuclear Overhauser effect distance restraints, 3JHH for prochiral assignments, temperature gradients) were determined by nmr and served as input data for molecular dynamics (MD) structure refinement. A simulated model compound corresponding to the [Cys-Pro-Pro-Cys]2 core elongated by the peptide chains in the Fab and Fc direction served as a starting structure for the final MD run. The conformation calculated in in vacuo does not agree with the C2 symmetry required from nmr data, but the structure obtained by a water simulation fulfills the requirement. Here the core of the hinge peptide derivative adopts a polyproline II double helix as in the x-ray structure of IgG1. Hence, segments responsible for the internal flexibility are located outside the core as confirmed by the flexibility of the solvent exposed C termini.     

Synthetic and conformational studies on dehydroalanine-containing model peptides

Three model dipeptides containing a dehydroalanine residue (delta Ala) at the C-terminal, Boc-X-delta Ala-NHCH3 [X = Ala, Val, and Phe,] have been synthesized and their solution conformations investigated by 1H-NMR, IR, and CD spectroscopy. NMR studies on these peptides in CDCl3 clearly indicate that the NH group of dehydroalanine is involved in an intramolecular hydrogen bond. This conclusion is supported by IR studies also. Nuclear Overhauser effect (NOE) studies are also accommodative of an inverse gamma-turn-type of conformation that is characterized by conformational angles of phi approximately -70 degrees and psi approximately +70 degrees around the X residue, and a C alpha i + 1 H-Ni + 2H interproton distance of 2.5 A. It appears that unlike dehydrophenylalanine or dehydroleucine, which tend to stabilize beta-turn type of structures occupying the i + 2 position of the turn, dehydroalanine favors the formation of an inverse gamma-turn, centered at the preceding L-residue in such solvents as CDCl3 and (CD3)2SO. A comparison of solution conformation of Boc Val-delta Ala-NHCH3 with the corresponding saturated analogue, Boc-Val-Ala-NHCH3, is also presented and shows that dehydroalanine is responsible for inducing the turn structure. It may be possible to design peptides with different preferred conformations using the suitable dehydroamino acid.     

Solution conformation of a peptide fragment representing a proposed RNA-binding site of a viral coat protein studied by two-dimensional NMR

The first 25 amino acids of the coat protein of cowpea chlorotic mottle virus are essential for binding the encapsidated RNA. Although an alpha-helical conformation has been predicted for this highly positively charged N-terminal region [Argos, P. (1981) Virology 110, 55-62; Vriend, G., Verduin, B. J. M., & Hemminga, M. A. (1986) J. Mol. Biol. 191, 453-460], no experimental evidence for this conformation has been presented so far. In this study, two-dimensional proton NMR experiments were performed on a chemically synthesized pentacosapeptide containing the first 25 amino acids of this coat protein [Ten Kortenaar, P. B. W., Krüse, J., Hemminga, M. A., & Tesser, G. I. (1986) Int. J. Pept. Protein Res. 27, 401-413]. All resonances could be assigned by a combined use of two-dimensional correlated spectroscopy and nuclear Overhauser enhancement spectroscopy carried out at four different temperatures. Various NMR parameters indicate the presence of a conformational ensemble consisting of helical structures rapidly converting into more extended states. Differences in chemical shifts and nuclear Overhauser effects indicate that lowering the temperature induces a shift of the dynamic equilibrium toward more helical structures. At 10 degrees C, a perceptible fraction of the conformational ensemble consists of structures with an alpha-helical conformation between residues 9 and 17, likely starting with a turnlike structure around Thr9 and Arg10. Both the conformation and the position of this helical region agree well with the secondary structure predictions mentioned above.