Restrained Energy Minimization Research Articles

PANDORA: A Fast, Anchor-Restrained Modelling Protocol for Peptide: MHC Complexes.

Deeper understanding of T-cell-mediated adaptive immune responses is important for the design of cancer immunotherapies and antiviral vaccines against pandemic outbreaks. T-cells are activated when they recognize foreign peptides that are presented on the cell surface by Major Histocompatibility Complexes (MHC), forming peptide:MHC (pMHC) complexes. 3D structures of pMHC complexes provide fundamental insight into T-cell recognition mechanism and aids immunotherapy design. High MHC and peptide diversities necessitate efficient computational modelling to enable whole proteome structural analysis. We developed PANDORA, a generic modelling pipeline for pMHC class I and II (pMHC-I and pMHC-II), and present its performance on pMHC-I here. Given a query, PANDORA searches for structural templates in its extensive database and then applies anchor restraints to the modelling process. This restrained energy minimization ensures one of the fastest pMHC modelling pipelines so far. On a set of 835 pMHC-I complexes over 78 MHC types, PANDORA generated models with a median RMSD of 0.70 Å and achieved a 93% success rate in top 10 models. PANDORA performs competitively with three pMHC-I modelling state-of-the-art approaches and outperforms AlphaFold2 in terms of accuracy while being superior to it in speed. PANDORA is a modularized and user-configurable python package with easy installation. We envision PANDORA to fuel deep learning algorithms with large-scale high-quality 3D models to tackle long-standing immunology challenges.

Membrane-induced structure of novel human tachykinin hemokinin-1 (hHK1).

PPT-C encoded hemokinin-1(hHK-1) of Homo sapiens (TGKASQFFGLM) is a structurally distinct neuropeptide among the tachykinin family that participate in the NK-1 receptor downstream signaling processes. Subsequently, signal transduction leads to execution of various effector functions which includes aging, immunological, and central nervous system (CNS) regulatory actions. However the conformational pattern of ligand receptor binding is unclear. The three-dimensional structure of the hemokinin-1 in aqueous and micellar environment has been studied by one and two-dimensional proton nuclear magnetic resonance (2D 1H-NMR spectroscopy) and distance geometry calculations. Data shows that hemokinin-1 was unstructured in aqueous environment; anionic detergent SDS induces α-helix formation. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (CYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. Helical conformation is observed from residue K3-M11. The conformational range of the peptide revealed by NMR studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared to that of Substance P potent NK1 agonist. Thus the report provides a structural insight to study hHK-1-NK1 interaction that is essential for hHK1 based signaling events.

Structural Insights into Human Hemokinin1 -Nk1Receptor Interactions

PPT-C encoded hemokinin-1 (TGKASQFFGLM) of Homo sapiens is a structurally distinct neuropeptide among the tachykinin family that participate in the NK-1 receptor downstream signaling processes. Subsequently, signal transduction leads to execution of various effecter functions which includes aging, immunological and central nervous system (CNS) regulatory actions. However the conformational pattern of ligand receptor binding is unclear. The three-dimensional structure of the hemokinin-1 in aqueous and micellar environment has been studied by CD Spectroscopy and two-dimensional proton nuclear magnetic resonance (2D 1H-NMR spectroscopy) and distance geometry calculations. CD data shows that hemokinin-1 was unstructured in aqueous environment; anionic detergent SDS induces α-helix formation. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (CYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. The results show that in aqueous environment hemokinin-1 lacks a definite secondary structure. The structure is well defined in presence of dodecyl phosphocholine micelles. The conformational range of the peptide revealed by NMR and CD studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared to that of Substance P, Neurokinin A and Neurokinin B, potent NK1, NK2 and NK3 agonists respectively. Thus the report provides a structural insight to study hHK-1-NK1 interaction that is essential for hHK1 based signaling events.

Conformational Characterization of Phosphorylated Analogs of a Peptide Mimetic of the Fourth Cytoplasmic Loop of CB1Cannabinoid Receptor

Phosphorylation is an important regulatory mechanism in receptor signaling. It has been shown that phosphorylation of the CB1 receptor disrupts modulation of ion channels by the receptor and CB1 cannabinoid receptor intracellular C-terminal tail domain (amino acids 401-417) is critical for G(i/o) protein coupling. Further, the proximal portion of the C-terminus of the cannabinoid CB1 receptor is a primary determinant for G-protein activation. The synthetic peptide fragment of the C-terminal juxtamembrane region (CB1 401-417) referred here as IL4 peptide mimicked the receptor's response of inhibiting adenylate cyclase. In the present study, we have used phosphorylated analogues of IL4 peptide to analyze the effect of phosphorylation on the conformation of the peptide and G-protein activation using NMR Spectroscopy. Unambiguous proton NMR assignments have been carried out with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (DYANA) to generate a family of structures which have been refined using restrained energy minimization and dynamics. The conformational range of the phosphorylated IL4 peptide revealed by NMR studies has been analyzed in terms of characteristic secondary structural features. The results obtained provide insight into the mechanism by which the peptide activates G-proteins, as a first step in signal transduction.

Bridging the gap between single-template and fragment based protein structure modeling using Spanner

Background: As the coverage of experimentally determined protein structures increases, fragment-based structural modeling approaches are expected to play an ever more important role in structural modeling. Here we introduce a structural modeling method by which an initial structural template can be extended by the addition of structural fragments to more closely match an aligned query sequence. A database of pro-tein fragments indexed by their internal coordinates was created and a novel methodology for their retrieval was implemented. After fragment selection and assembly, sidechains are replaced and the all-atom model is refined by restrained energy minimization. We implemented the proposed method in the program Span-ner and benchmarked it using a previously published set of 367 immunoglobulin (Ig) loops, 206 historical query-template pairs and alignments from the Critical Assessment of protein Structure Prediction (CASP) experiment, and 217 structural alignments between remotely homologous query-template pairs. The con-straint-based modeling software MODELLER and previously reported results for RosettaAntibody, were used as references. Results: The error in the modeled structures was assessed by root-mean square deviation (RMSD) from the native structure, as a function of the query-template sequence identity. For the Ig benchmark set, for which a single fragment was used to model each loop, the average RMSD for Spanner (3 +/- 1.5 A) was found to lie midway between that of MODELLER (4 +/- 2 A) and RosettaAntibody (2 +/- 1 A). For the CASP and structural alignment benchmarks, for which gaps represent a small fraction of the modeled residues, the difference between Spanner and MODELLER were much smaller then the standard deviations of either program. The Spanner web server and source code are available at http://sysimm.ifrec.osaka-u.ac.jp/Spanner/. Conclusions: For typical homology modeling, Spanner is at least as good, on average as the template-free constraint-driven approach used by MODELLER. The Ig model results suggest that when gap regions represent a significant fraction of the alignment, Spanner’s efficient use of fragment libraries, along with local sequence and secondary structural information, significantly improve model accuracy without a dra-matic increase in computational cost.

Structure and Dynamics of Double Helical DNA in Torsion Angle Hyperspace: A Molecular Mechanics Approach

Analysis of the conformational space populated by the torsion angles and the correlation between the conformational energy and the sequence of DNA are important for fully understanding DNA structure and function. Presence of seven variable torsion angles about single covalent bonds in DNA main chain puts a big challenge for such analysis. We have carried out restrained energy minimization studies for four representative dinucleosides, namely d(ApA):d(TpT), d(CpG):d(CpG), d(GpC):d(GpC) and d(CpA):d(TpG) to determine the energy hyperspace of DNA in context to the values of the torsion angles and the structural properties of the DNA conformations populating the favorable regions of this energy hyperspace. The torsion angles were manipulated by constraining their values at the reference points and then performing energy minimization. The energy minima obtained on the potential energy contour plots mostly correspond to the conformations populated in crystal structures of DNA. Some novel favorable conformations that are not present in crystal structure data are also found. The plots also suggest few low energy routes for conformational transitions or the associated energy barrier heights. Analyses of base pairing and stacking possibility reveal structural changes accompanying these transitions as well as the flexibility of different base steps towards variations in different torsion angles.

The NMR Structures of the Major Intermediates of the Two-domain Tick Carboxypeptidase Inhibitor Reveal Symmetry in Its Folding and Unfolding Pathways

There is a lack of experimental structural information about folding intermediates of multidomain proteins. Tick carboxypeptidase inhibitor (TCI) is a small, disulfide-rich protein consisting of two domains that fold and unfold autonomously through the formation of two major intermediates, IIIa and IIIb. Each intermediate contains three native disulfide bonds in one domain and six free cysteines in the other domain. Here we have determined the NMR structures of these two intermediates trapped and isolated at acidic pH in which they are stable and compared their structures with that of the native protein analyzed under the same conditions. Both IIIa and IIIb were found to contain a folded region that corresponds to the N- and C-terminal domains of TCI, respectively, with structures very similar to the corresponding regions of the native protein. The remainder of the polypeptide chains of the intermediates was shown to be unfolded in a random coil conformation. Solvent exchange measurements further indicated that the two protein domains are not completely independent, but affect each other in terms of dynamics and stability, in agreement with reported inhibitory activity data. The derived results provide structural evidence for symmetric TCI folding and unfolding mechanisms that converge in IIIa and IIIb and reveal the structural basis that accounts for the strong and simultaneous accumulation of both intermediates. Altogether, this work has important implications for a better understanding of the folding mechanisms of multidomain, disulfide-rich proteins.

Structural Characterization of Neurokinin-3 Receptor Selective Peptide Agonist Scyliorhinin II Bound to DPC Micelles

Scyliorhinin II, a cyclic Tachykinin peptide, is a potent NK3 receptor agonist. The pharmacology of NK3 receptor is least characterized out of the three tachykinin receptor subtypes cloned and characterized for Tachykinins. To understand the structural basis of peptide-receptor interaction, the three-dimensional structure of the Scyliorhinin II in aqueous and micellar environments has been studied by two-dimensional proton nuclear magnetic resonance (2D 1H-NMR spectroscopy) and distance geometry calculations. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (DYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. The results show that in an aqueous environment, Scyliorhinin II lacks a definite secondary structure. The structure is well-defined in presence of dodecyl phosphocholine micelles. The global fold of Scyliorhinin II bound to DPC micelles consists of a well-defined helix in the C-terminal region from residue 12–18 and a series of turns towards N-terminus. The structure is further stabilized by disulfide bond between Cys7 and Cys13. The conformational range of the peptide revealed by NMR and CD studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared with those of Substance P, Neurokinin A and Neurokinin B, potent NK1, NK2, and NK3 agonists, respectively.

Unique helical conformation of the fourth cytoplasmic loop of the CB 1 cannabinoid receptor in a negatively charged environment

The proximal portion of the C-terminus of the CB 1 cannabinoid receptor is a primary determinant for G-protein activation. A 17 residue proximal C-terminal peptide (rodent CB1 401–417), the intracellular loop 4 (IL4) peptide, mimicked the receptor’s G-protein activation domain. Because of the importance of the cationic amino acids to G-protein activation, the three-dimensional structure of the IL4 peptide in a negatively charged sodium dodecyl sulfate (SDS) micellar environment has been studied by two-dimensional proton nuclear magnetic resonance (2D 1H NMR) spectroscopy and distance geometry calculations. Unambiguous proton NMR assignments were carried out with the aid of correlation spectroscopy (DQF-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints were used in torsion angle dynamics algorithm for NMR applications (DYANA) to generate a family of structures which were refined using restrained energy minimization and dynamics. In water, the IL4 peptide prefers an extended conformation, whereas in SDS micelles, 3 10-helical conformation is induced. The predominance of 3 10-helical domain structure in SDS represents a unique difference compared with structure in alternative environments, which can significantly impact global electrostatic surface potential on the cytoplasmic surface of the CB 1 receptor and might influence the signal to the G-proteins.

Close Packing of Helices 3 and 12 of 16 S rRNA Is Required for the Normal Ribosome Function

The along-groove packing motif is a quasi-reciprocal arrangement of two RNA double helices in which a backbone of each helix is closely packed within the minor groove of the other helix. At the center of the inter-helix contact, a GU base pair in one helix packs against a Watson-Crick base pair in the other helix. Here, based on in vivo selection from a combinatorial gene library of 16 S rRNA and on functional characterization of the selected clones, we demonstrate that the normal ribosome performance requires that helices 3 and 12 be closely packed. In some clones the Watson-Crick and GU base pairs exchange in their positions between the two helices, which affects neither the quality of the helix packing, nor the ribosome function. On the other hand, perturbations in the close packing usually lead to a substantial drop in the ribosome activity. The functionality of the clones containing such perturbations may depend on the presence of particular elements in the vicinity of the area of contact between helices 3 and 12. Such cases do not exist in natural 16 S rRNA, and their selection enriches our knowledge of the constraints imposed on the structure of ribosomal RNA in functional ribosomes.

Three-Dimensional Structure of Neuropeptide K Bound to Dodecylphosphocholine Micelles

Neuropeptide K (NPK), an N-terminally extended form of neurokinin A (NKA), represents the most potent and longest lasting vasodepressor and cardiomodulatory tachykinin reported thus far. NPK has been shown to have high selectivity for the NK2 receptor. Because the micelle-associated structure may be relevant to the NPK-receptor interaction, the three-dimensional structure of the NPK in aqueous and micellar environments has been studied by two-dimensional proton nuclear magnetic resonance (2D (1)H NMR spectroscopy) and distance geometry calculations. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (DQF-COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The interproton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (DYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. The results show that in an aqueous environment NPK lacks a definite secondary structure, although some turn-like elements are present in the N terminus. The structure is well-defined in the presence of dodecylphosphocholine micelles. The global fold of NPK bound to DPC micelles consists of two well-defined helices from residues 9 to 18 and residues 27 to 33 connected by a noncanonical beta turn. The N terminus of the peptide is characterized by a 3(10) helix or a series of dynamic beta turns. The conformational range of the peptide revealed by NMR and circular dichroism (CD) studies has been analyzed in terms of characteristic secondary features. The observed conformational features have been further compared to a NKA and neuropeptide gamma (NPgamma) potent endogenous agonist for the NK2 receptor.

Experimental constraints on quaternary structure in Alzheimer's beta-amyloid fibrils.

We describe solid-state nuclear magnetic resonance (NMR) measurements on fibrils formed by the 40-residue beta-amyloid peptide associated with Alzheimer's disease (Abeta(1-40)) that place constraints on the identity and symmetry of contacts between in-register, parallel beta-sheets in the fibrils. We refer to these contacts as internal and external quaternary contacts, depending on whether they are within a single molecular layer or between molecular layers. The data include (1) two-dimensional 13C-13C NMR spectra that indicate internal quaternary contacts between side chains of L17 and F19 and side chains of I32, L34, and V36, as well as external quaternary contacts between side chains of I31 and G37; (2) two-dimensional 15N-13C NMR spectra that indicate external quaternary contacts between the side chain of M35 and the peptide backbone at G33; (3) measurements of magnetic dipole-dipole couplings between the side chain carboxylate group of D23 and the side chain amine group of K28 that indicate salt bridge interactions. Isotopic dilution experiments allow us to make distinctions between intramolecular and intermolecular contacts. On the basis of these data and previously determined structural constraints from solid-state NMR and electron microscopy, we construct full molecular models using restrained molecular dynamics simulations and restrained energy minimization. These models apply to Abeta(1-40) fibrils grown with gentle agitation. We also present evidence for different internal quaternary contacts in Abeta(1-40) fibrils grown without agitation, which are morphologically distinct.

SOLUTION STRUCTURE OF IMAZALIL/β-CYCLODEXTRIN INCLUSION COMPLEX

The effectiveness of the Imazalil fungicide is hampered by its water solubility, intrinsically very low: preparation of Imazalil/β-Cyclodextrins inclusion complexes in equimolar ratio resulted in an efficient and convenient water soluble formulation, thus resulting in a simple and non expensive preparation of the fungicide. The structural characterization in solution of this inclusion complex has been elucidated by means of NMR spectroscopy and computational methods (Restrained Molecular Dynamics and Restrained Energy Minimization). In the resulting structures the aromatic ring of Imazalil resulted deeply inserted into the CD cavity, whereas the allyl group and the imidazole ring were excluded from it, but still in close contact with the larger rim. The inclusion mechanism is also demonstrated.

Solution conformations of amphidinolide H

Solution conformations of amphidinolide H ( 1), a 26-membered macrolide exhibiting potent cytotoxic and antitumor activity, in CDCl 3 and DMSO- d 6 were investigated on the basis of NMR data, distance geometry calculation, and restrained energy minimization. Three-dimensional conformations in CDCl 3 were suggested to be close to the X-ray structure of 1, while those in DMSO- d 6 were indicated to be different from both those in CDCl 3 and the X-ray structure.

A Solid-Phase Approach to the Phallotoxins: Total Synthesis of [Ala7]-Phalloidin

Herein we report a solid-phase synthetic approach to [Ala7]-phalloidin (1). Prior syntheses of phallotoxins were carried out using solution-phase routes that required large scale and preclude library production. The route presented here consists of solution-phase preparation of key orthogonally protected amino acid building blocks, followed by a solid-phase peptide synthesis sequence, featuring two resin-bound macro-cyclization reactions. The final product mixture was composed of two atropisomeric compounds, one designated "natural" (1) and the other designated "non-natural" (1'). The structures of these species were modeled using restrained energy minimization with NMR-derived restraints.

Membrane-induced structure of the mammalian tachykinin Neuropeptide gamma

Neuropeptide gamma (NPγ) is a neurokinin-2 (NK-2) receptor selective agonist, which plays an important role in mediation of asthma and elicits a wide range of biological responses like bronchoconstriction, vasodepression and regulation of endocrine functions. The structure determination of this peptide agonist is important in understanding the molecular basis of peptide ligand recognition by the receptor and for rational drug design. In the present study we report the solution structure of NPγ characterized by circular dichroism (CD) spectropolarimetry and 2D 1H NMR spectroscopy in both aqueous and membrane mimetic solvents. Effect of calcium ions on the conformation of NPγ was also studied using CD spectropolarimetry. Sequence-specific resonance assignments of protons have been made with the aid of correlation spectroscopy experiments and nuclear Overhauser effect spectroscopy experiments. The distance constraints obtained from the NMR data have been utilized to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that in water NPγ prefers to be in an extended chain conformation whereas a helical conformation is induced in the central core and the C-terminal region of the peptide (K13–M21) in the presence of perdeuterated dodecylphosphocholine micelles, a membrane model system. A type II′ β turn from H9 to R11 precedes the helical core in the C-terminus of NPγ. N-terminus of NPγ also displays some degree of order and a possible turn structure. Conformation adopted by NPγ in presence of lipid micelles represents a structural motif typical of NK-2 selective agonists and is similar to that observed for Neurokinin A in hydrophobic environment. The observed conformational features have been correlated to the binding ability and biological activity of NPγ.

Three Dimensional Structure of Mammalian Tachykinin Peptide Neurokinin B Bound to Lipid Micelles

Neurokinin B (NKB), a decapeptide of mammalian origin exhibits a variety of biological activities such as regulatory functions in reproduction, pre-eclampsia and neuroprotection in Alzheimer's disease. In order to gain insight into structure-function relationship, three- dimensional structure of NKB has been investigated using CD spectropolarimetry and two-dimensional proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy in aqueous and membrane mimetic solvents. Unambiguous NMR assignments of resonances have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and Nuclear Overhauser Effect Spectroscopy (NOESY) experiments. Distance constraints obtained from the NMR data have been used to generate a family of structures, which have been refined using restrained energy minimization and dynamics. Our data show that a helical structure is induced in NKB, in presence of perdeuterated dodecyl phosphocholine (DPC) micelles, a membrane model system. Further, the conformation adopted by NKB in presence of DPC micelles represents a structural motif typical of neurokinin-3 selective agonists.

Three-Dimensional Structure of the Mammalian Tachykinin Peptide Neurokinin A Bound to Lipid Micelles

The solution structure of NKA, a decapeptide of mammalian origin, has been characterized by CD spectropolarimetry and 2D proton nuclear magnetic resonance (2D 1H-NMR) spectroscopy in both aqueous and membrane mimetic solvents. Unambiguous NMR assignments of protons have been made with the aid of correlation spectroscopy (DQF-COSY and TOCSY) experiments and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The distance constraints obtained from the NMR data have been utilized to generate a family of structures, which have been refined using restrained energy minimization and dynamics. These data show that in water NKA prefers to be in an extended chain conformation whereas a helical conformation is induced in the central core and the C-terminal region (D4-M10) of the peptide in the presence of perdeuterated dodecylphosphocholine (DPC) micelles, a membrane model system. Though less defined the N-terminus also displays some degree of order and a possible turn structure. The conformation adopted by NKA in the presence of DPC micelles represents a structural motif typical of neurokinin-2 selective agonists and is similar to that reported for eledoisin in hydrophobic environment.

The solution structure of the oxidized bovine microsomal cytochrome b5 mutant V61H

Using 1488 NOE constraints, 19 stereo-specific assignments, 13 pairs of H-bond constraints, and 140 pseudo-contact shift constraints, a family of 35 structures of bovine microsomal cytochrome b 5 mutant V61H has been obtained through the program PSEUDYANA. The family has been further refined by restrained energy minimization to give a family of final structures. The RMSD values of final structures with respect to the average structure are 0.45 ± 0.11 and 0.96 ± 0.10 Å for backbone and heavy atoms, respectively. The final Δ χ ax and Δ χ rh values are 2.34 × 10 −32 and −0.67 × 10 −32 m 3, respectively. The comparisons between the solution structures of mutant V61H and WT cytochrome b 5, and X-ray structure of the mutant V61H show that the global folding of the molecule in solution is unchanged and the side-chain of His61 deviates from the heme pocket and extends into the solvent like in its crystal structure. However, the helices around the heme pocket undergo outward global displacement while their local conformations are well maintained. Meanwhile, the heme ring shows a little off the heme pocket, which accounts for the lower stability of the mutant. Additionally, the axial ligand rings counterclockwise rotate around His39 N–Fe axis due to the mutation, which is confirmed by variation of the hyperfine shifts of the heme protons of V61H compared to those of WT cytochrome b 5.

Comparison of the solution structures of angiotensin I & II. Implication for structure-function relationship.

Conformational analysis of angiotensin I (AI) and II (AII) peptides has been performed through 2D 1H-NMR spectroscopy in dimethylsulfoxide and 2,2,2-trifluoroethanol/H2O. The solution structural models of AI and AII have been determined in dimethylsulfoxide using NOE distance and 3JHNHalpha coupling constants. Finally, the AI family of models resulting from restrained energy minimization (REM) refinement, exhibits pairwise rmsd values for the family ensemble 0.26 +/- 0.13 A, 1.05 +/- 0.23 A, for backbone and heavy atoms, respectively, and the distance penalty function is calculated at 0.075 +/- 0.006 A2. Comparable results have been afforded for AII ensemble (rmsd values 0.30 +/- 0.22 A, 1.38 +/- 0.48 A for backbone and heavy atoms, respectively; distance penalty function is 0.029 +/- 0.003 A2). The two peptides demonstrate similar N-terminal and different C-terminal conformation as a consequence of the presence/absence of the His9-Leu10 dipeptide, which plays an important role in the different biological function of the two peptides. Other conformational variations focused on the side-chain orientation of aromatic residues, which constitute a biologically relevant hydrophobic core and whose inter-residue contacts are strong in dimethylsulfoxide and are retained even in mixed organic-aqueous media. Detailed analysis of the peptide structural features attempts to elucidate the conformational role of the C-terminal dipeptide to the different binding affinity of AI and AII towards the AT1 receptor and sets the basis for understanding the factors that might govern free- or bound-depended AII structural differentiation.