Conformational studies on enkephalins using the MOLS technique

Conformational studies of two linear enkephalin molecules, Met-enkephalin and Leu-enkephalin, have been carried using the mutually orthogonal Latin squares (MOLS) technique with the ECEPP/3 force field. This technique was developed recently in our laboratory to perform an unbiased search of the conformational space of peptides and to locate low energy conformations. The present study identified all the folds predicted by other studies, and in addition picked up other energetically favorable structures. The results suggest that the peptide backbone exists as a mixture of folded and unfolded forms (approximately 50% each). The study also provides information on the distribution of the low energy conformations that we have classified on the basis of structural motifs, backbone hydrogen-bonding patterns, and root mean square deviations in atomic positions.     

Enkephalin is a competitive antagonist of cholecystokinin in the gastrointestinal tract, as predicted from prior conformational analysis

Prior calculations based on ECEPP (Empirical Conformational Energies for Peptides Program) of the low energy minima for cholecystokinin (CCK) and Met-enkephalin have demonstrated that significant structural features of these two peptides are identical. This result suggested the possibility that Met-enkephalin, as well as other enkephalin analogues of similar structure, could associate with receptors for CCK. To test this theoretical result, we examined the ability of Met-enkephalin and its analogues to bind to peripheral CCK receptors in the rat gastrointestinal tract; in particular, we measured the ability of the opiate peptide to inhibit the effects of CCK in a physiological assay system which we have previously characterized: CCK-induced contraction of the isolated rat pyloric sphincter. We find that Met-enkephalin is an antagonist of the CCK-8-induced contraction, with a IC₅₀ of 110 nM. Furthermore, antibodies against CCK were found to cross-react with Met-enkephalin and its analogues in a manner which suggests a distinct structure-activity relationship. These experimental results strongly support the theoretical results of conformational analysis showing structural similarity between enkephalin and CCK. They further suggest that enkephalins could modulate the response of CCK systems under physiological conditions.     

Position dependence of amino acid intrinsic helical propensities II: non-charged polar residues: Ser, Thr, Asn, and Gln.

The assumption that the intrinsic alpha-helical propensities of the amino acids are position independent was critical in several helix/coil transition theories. In the first paper of these series, we reported that this is not the case for Gly and nonpolar aliphatic amino acids (Val, Leu, Met, and Ile). Here we have analyzed the helical intrinsic propensities of noncharged polar residues (Ser, Thr, Asn, and Gln) at different positions of a model polyalanine-based peptide. We found that Thr is more favorable (by approximately 0.3 kcal/mol) at positions N1 and N2 than in the helix center, although for Ser, Asn, and Gln the differences are smaller (+/-0.2 kcal/mol), and in many cases within the experimental error. There is a reasonable agreement (+/-0.2 kcal/mol) between the calculated free energies, using the ECEPP/2 force field equipped with a hydration potential, and the experimental data, except at position N1.     

Characterisation of N-terminally extended met-enkephalin Arg6Gly7Leu8 variants in the porcine upper digestive tract

Proenkephalin A-derived peptides are known to occur in the gut, but their precise identity is uncertain. We report here the isolation of N-terminally extended forms of Met-enkephalin Arg6Gly7Leu8 from porcine upper digestive tract monitored by radioimmunoassay. A single major form was identified in pyloric antral muscle and mucosa, but in the duodenum two major forms were detected. Microsequence analysis together with immunological data revealed that the antral mucosal peptide and the most acidic duodenal peptide had identical amino-acid sequences, corresponding to a 5.3 kDa peptide terminating in Met-enkephalin Arg6Gly7Leu8. The data indicate that high-molecular-weight peptides may constitute a major proportion of gut opioid peptide immunoactivity.     

Conformational energy analysis of the substitution of Val for Gly 233 in a functional region of platelet GPIb alpha in platelet-type von Willebrand disease

Platelet-type von Willebrand disease (PT-vWD) is an autosomal dominant bleeding disorder in which patient platelets exhibit an abnormally increased binding of circulating von Willebrand factor (vWF). We have recently shown that this abnormality is associated with a point mutation resulting in substitution of Val for Gly 233 in platelet membrane glycoprotein Ib alpha (GPIb alpha), a major component of the platelet GPIb/IX receptor for vWF. To investigate the effect of this substitution on the three-dimensional structure of this region of the protein, we have generated the allowed (low energy) conformations of the region of the GPI alpha protein containing residues 228-238 (with 5 residues on either side of the critical residue 233) with Gly 233 (wild type) and Val 233 (PT-vWD) using the computer program ECEPP (Empirical Conformational Energies of Peptides Program). The wild-type sequence is Tyr-Val-Trp-Lys-Gln-Gly-Val-Asp-Val-Lys-Ala. We find that the Gly 233-containing peptide can exist in two low energy conformers. The lowest energy conformer is a structure containing a beta-turn at Gln 232-Gly 233 while the alternative conformation is an amphipathic helical structure. Only the amphipathic helical structure is allowed for the Val 233-containing peptide which contains a hydrophobic 'face' consisting of Val 229, Val 233 and Val 236 and another hydrophilic surface composed of such residues as Lys 231 and Asp 235. No such surfaces exist for the lowest energy bend conformer for the Gly 233-containing peptide, but do exist in the higher energy helical structure. The amphipathic surfaces in the 228-238 region of the Val 233-containing GPIb alpha protein may associate strongly with complementary surfaces during vWF binding to the GPIb/IX receptor complex and may help explain heightened association of vWF with this receptor in PT-vWD.     

Enkephalin degradation in mouse brain studied by a new H.P.L.C. method: further evidence for the involvement of carboxydipeptidase

The degradation of the Met-enkephalin by peptidases present in striatum extracts was followed by a new H.P.L.C. method, based on ligand exchange chromatography. An aminopeptidase activity is present in both the soluble and the particulate fractions of mouse striatum. In addition we confirmed the presence, mainly in the particulate fraction, of a carboxydipeptidase releasing TyrGlyGly and PheMet from Met-enkephalin. The tetrapeptide GlyGlyPheMet also appears to be cleaved by the same enzyme activity.     

Activated conformations of the ras-gene-encoded p21 protein. 1. An energy-refined structure for the normal p21 protein complexed with GDP.

A complete three-dimensional structure for the ras-gene-encoded p21 protein with Gly 12 and Gln 61, bound to GDP, has been constructed in four stages using the available alpha-carbon coordinates as deposited in the Brookhaven National Laboratories Protein Data Bank. No all-atom structure has been made available despite the fact that the first crystallographic structure for the p21 protein was reported almost four years ago. In the p21 protein, if amino acid substitutions are made at any one of a number of different positions in the amino acid sequence, the protein becomes permanently activated and causes malignant transformation of normal cells or, in some cell lines, differentiation and maturation. For example, all amino acids except Gly and Pro at position 12 result in an oncogenic protein; all amino acids except Gln, Glu and Pro at position 61 likewise cause malignant transformation of cells. We have constructed our all-atom structure of the non-oncogenic protein from the x-ray structure in order to determine how oncogenic amino acid substitutions affect the three-dimensional structure of this protein. In Stage 1 we generated a poly-alanine backbone (except at Gly and Pro residues) through the alpha-carbon structure, requiring the individual Ala, Pro or Gly residues to conform to standard amino acid geometry and to form trans-planar peptide bonds. Since no alpha-carbon coordinates for residues 60-65 have been determined, these residues were modeled by generating them in the extended conformation and then subjecting them to molecular dynamics using the computer application DISCOVER and energy minimization using DISCOVER and the ECEPP (Empirical Conformational Energies for Peptides Program). In Stage 2, the positions of residues that are homologous to corresponding residues of bacterial elongation factor Tu (EF-Tu) to which p21 bears an overall 40% sequence homology, were determined from their corresponding positions in a high-resolution structure of EF-Tu. Non-homologous loops were taken from the structure generated in Stage 1 and were placed between the appropriate homologous segments so as to connect them. In Stage 3, all bad contacts that occurred in this resulting structure were removed, and the coordinates of the alpha-carbon atoms were forced to superimpose as closely as possible on the corresponding atoms of the reference (x-ray) structure. Then the side chain positions of residues of the non-homologous loop regions were modeled using a combination of molecular dynamics and energy minimization using DISCOVER and ECEPP respectively.(ABSTRACT TRUNCATED AT 400 WORDS)     

Conformational energy analysis of the substitution of Val for Gly 233 in a functional region of platelet GPIbα in platelet-type von Willebrand disease

Platelet-type von Willebrand disease (PT-vWD) is an autosomal dominant bleeding disorder in which patient platelets exhibit an abnormally increased binding of circulating von Willebrand factor (vWF). We have recently shown that this abnormality is associated with a point mutation resulting in substitution of Val for Gly 233 in platelet membrane glycoprotein Ibα (GPIbga), a major component of the platelet (GPIb/IX receptor for vWF. To investigate the effect of this substitution on the three-dimensional structure of this region of the protein, we have generated the allowed (low energy) conformations of the region of the GPIα protein containing residues 228–238 (with 5 residues on either side of the critical residue 233) with Gly 233 (wild type) and Val 233 (PT-vWD) using the computer program ECEPP (Empirical Conformational Energies of Peptides Program). The wild-type sequence is Tyr-Val-Trp-Lys-Gln-Gly-Val-Asp-Val-Lys-Ala. We find that the Gly 233-containing peptide can exist in two low energy conformers. The lowest energy conformer is a structure containing a β-turn at Gln 232-Gly 233 while the alternative conformation is an amphipathic helical structure. Only the amphipathic helical structure is allowed for the Val 233-containing peptide which contains a hydrophobic ‘face’ consisting of Val 229, Val 233 and Val 236 and another hydrophilic surface composed of such residues as Lys 231 and Asp 235. No such surfaces exist for the lowest energy bend conformer for the Gly 233-containing peptide, but do exist in the higher energy helical structure. The amphiphatic surfaces in the 228–238 region of the Val 233-containing GPIbα protein may associate strongly with complementary surfaces during vWF binding to the GPIb/IX receptor complex and may help explain heightened association of vWF with this receptor in PT-vWD.     

Preferred conformations of cyclic Ac-Cys-Pro-Xaa-Cys-NHMe peptides: a model for chain reversal and active site of disulfide oxidoreductase

The conformational study on cyclic Ac-Cys-Pro-Xaa-Cys-NHMe (Ac-CPXC-NHMe; X=Ala, Val, Leu, Aib, Gly, His, Phe, Tyr, Asn and Ser) peptides has been carried out using the Empirical Conformational Energy Program for Peptides, version 3 (ECEPP/3) force field and the hydration shell model in the unhydrated and hydrated states. This work has been undertaken to investigate structural implications of the CPXC sequence as the chain reversal for the initiation of protein folding and as the motif for active site of disulfide oxidoreductases. The backbone conformation DAAA is commonly the most feasible for cyclic CPXC peptides in the hydrated state, which has a type I beta-turn at the Pro-Xaa sequence. The proline residue and the hydrogen bond between backbones of two cystines as well as the formation of disulfide bond appear to play a role in stabilizing this preferred conformation of cyclic CPXC peptides. However, the distributions of backbone conformations and beta-turns may indicate that the cyclic CPXC peptide seems to exist as an ensemble of beta-turns and coiled conformations in aqueous solution. The intrinsic stability of the cyclic CPXC motif itself for the active conformation seems to play a role in determining electrochemical properties of disulfide oxidoreductases.     

Conformational analysis of the 20-residue membrane-bound portion of melittin by conformational space annealing

The conformational space of the 20-residue membrane-bound portion of melittin has been investigated extensively with the conformational space annealing (CSA) method and the ECEPP/3 (Empirical Conformational Energy Program for Peptides) algorithm. Starting from random conformations, the CSA method finds that there are at least five different classes of conformations, within 4 kcal/mol, which have distinct backbone structures. We find that the lowest energy conformation of this peptide from previous investigations is not the global minimum-energy conformation (GMEC); but it belongs to the second lowest energy class of the five classes found here. In four independent runs, one conformation is found repeatedly as the lowest energy conformation of the peptide (two of the four lowest energy conformations are identical; the other two have essentially identical backbone conformations but slightly different side-chain conformations). We propose this conformation, whose energy is lower than that found previously by 1.9 kcal/mol, as the GMEC of the ECEPP/3 force field. The structure of the proposed GMEC is less helical and more compact than the previous one. It appears that the CSA method can find several classes of conformations of a 20-residue peptide starting from random conformations utilizing only its amino acid sequence information. The proposed GMEC has also been found with a modified electrostatically driven Monte Carlo method [D. R. Ripoll, A. Liwo, and H.A. Scheraga (1998) “New Developments of the Electrostatically Driven Monte Carlo Method: Test on the Membrane-Bound Portion of Melittin,” Biopolymers, Vol. 46, pp. 117–126]. © 1998 John Wiley & Sons, Inc. Biopoly 46: 103–115, 1998     

A statistical analysis of side-chain conformations in proteins: Comparison with ECEPP predictions

A comparison of the statistical distributions of side-chain conformations of 17 amino acids (Gly, Ala, and Pro excluded), observed in 63 nonhomologous globular proteins (covering 10,832 residues), is made with similar distributions calculated from the low-energy conformational states for the same amino acids (blocked with acetyl and N-methylamide groups at the N- and C-termini, respectively) obtained by Vásquezet al. [(1983),Macromolecules 16, 1043–1049 using the ECEPP/2 force field. Those residues (i) with linear side chains (Arg, Lys, Met, Cys, Ser), or those that are unbranched through the-carbon atom (Glu, Gln) show good agreement, whereas (ii) those with side chains that are branched at C or C show poor agreement with ECEPP calculations. A possible explanation for this is shown to be the greater tendency for side-chain atoms in class (ii) to interact with the backbone and/or adjacent side chains. Accordingly, ECEPP/3 calculations, carried out after elongating the backbone chain of the model peptide unit (by adding three Ala residues on each side of the central residue, and then blocking the termini as before), result in distributions that are often closer to the observed side-chain distributions. The implications of these results for the relative importance of short-range versus long-range interactions in determining protein structure are discussed.     

Conformational solution studies of neuropeptide gamma using CD and NMR spectroscopy.

Neuropeptide gamma is one of the largest members of the tachykinin family of peptides, exhibiting strong agonistic activity towards the NK-2 tachykinin receptor. This peptide was synthesized by the solid-phase method using the Fmoc chemistry. Circular-dichroism spectroscopy (CD) investigations of this peptide were performed in phosphate buffer, in the presence of sodium dodecylsulphate (SDS) micelles and trifluoroethanol (TFE) solutions and in DMSO-d6 using the 2D NMR technique in conjunction with two different theoretical approaches. The first assumes multiconformational equilibrium of the peptide studied characterized by the values of statistical weights of low-energy conformations. These calculations were performed using three different force fields ECEPP/3, AMBER4.1 and CHARMM (implemented in the X-PLOR program). The second method incorporates interproton distance and dihedral angle constraints into the starting conformation using the Simulated Annealing algorithm (X-PLOR program). The CD experiments revealed that although the peptide studied is flexible in polar solvents, a tendency to adopt a helical structure was observed in the hydrophobic environment. The NMR data (NOE effects) indicate a helical or reverse structure in the Ile7-His12 fragment of the peptide studied in DMSO-d6 solution. The results obtained cannot be interpreted in terms of a single conformation. Most of the conformations obtained with the ECEPP/3 force field possess a high content of a helical structure. None of the conformers, obtained with the AMBER4.1 and CHARMM force fields, can be considered as the dominant one. In all conformations several beta-turns were detected and in some cases gamma-turns were also found. But in fact, it is rather difficult to select the position of the secondary element(s) present in the structure of NPgamma in solution. All conformers calculated with the X-PLOR program (with using NMR derived distance and torsion angle constraints) are stabilized by several beta-turns. Common structural motives are a type IV beta-turn in the Gln6-His12 fragment. All conformations obtained using two approaches adopt very similar turn shapes in the middle region of molecule and a random structure on the N- and C-terminal fragments.     

Folding of the villin headpiece subdomain from random structures. Analysis of the charge distribution as a function of pH

The structure of the 36 residue villin headpiece subdomain is investigated with the electrostatically driven Monte Carlo method. The ECEPP/3 (Empirical Conformational Energy Program for Peptides) force field, plus two different continuum solvation models, were used to describe the conformational energy of the chain with both blocked and unblocked N and C termini. A statistical analysis of an ensemble of ab initio generated conformations was carried out, based on a comparison with a set of ten native-like structures derived from published experimental data, by using rigid geometry and NMR-derived constraints obtained at pH 3.7. The ten native-like structures satisfy the NMR-derived constraints. The whole ensemble of conformations of the terminally unblocked villin headpiece sub-domain, generated by using ECEPP/3 with a continuum solvation model, were subsequently evaluated at pH 3.7 with a potential function that includes ECEPP/3 combined with a fast multigrid boundary element method. At pH 3.7, the lowest-energy conformation found during the conformational search satisfies approximately 70% of both the distance and the dihedral-angle constraints, and possesses the characteristic packing of three phenylalanine residues that constitute the main part of the hydrophobic core of the molecule. On the other hand, computations at pH 3.7 and pH 7.0 for the ten native-like structures satisfying the NMR-derived constraints indicate a substantial change in the charge distribution for each type of amino acid residue with the change in pH. The results of this study provide a basis to understand the effect of the interactions, such as hydrophobicity, charge-charge interaction and solvent polarization, on the stability of this small alpha-helical protein.     

Computational study of the Trp-cage miniprotein based on the ECEPP/3 force field

Using a newly developed Monte Carlo global optimization method called basin paving, we have performed an ab initio computation for the structure of Trp-cage based on the ECEPP/3 force field in vacuo. The lowest energy minimum has been located. Its corresponding configuration is comparable to the native structure of Trp-cage (PDB code 1L2Y) with a backbone root mean square deviation of 2.24 A.     

The multiple-minima problem in the conformational analysis of polypeptides. III. An Electrostatically Driven Monte Carlo Method: Tests on enkephalin

The three-dimensional conformation of Met-enkephalin, corresponding to the lowest minimum of the empirical potential energy function ECEPP/2 (empirical conformational energy program for peptides), has been determined using a new algorithm, viz. the Electrostatically Driven Monte Carlo Method. This methodology assumes that a polypeptide or protein molecule is driven toward the native structure by the combined action of electrostatic interactions and stochastic conformational changes associated with thermal movements. These features are included in the algorithm that produces a Monte Carlo search in the conformational hyperspace of the polypeptide, using electrostatic predictions and a random sampling technique to locate low-energy conformations. In addition, we have incorporated an alternative mechanism that allows the structure to escape from some conformational regions representing metastable local energy minima and even from regions of the conformational space with great stability. In 33 test calculations on Met-enkephalin, starting from arbitrary or completely random conformations, the structure corresponding to the global energy minimum was found inall the cases analyzed, with a relatively small search of the conformational space. Some of these starting conformations wereright orleft-handed -helices, characterized by good electrostatic interactions involving their backbone peptide dipoles; nevertheless, the procedure was able to convert such locally stable structures to the global-minimum conformation.On leave from the National University of San Luis, Faculty of Sciences and Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Instituto de Matemática Aplicada, San Luis, Ejército de los Andes 950, 5700 San Luis, Argentina.     

New general approach for determining the solution structure of a ligand bound weakly to a receptor: structure of a fibrinogen Aalpha-like peptide bound to thrombin (S195A) obtained using NOE distance constraints and an ECEPP/3 flexible docking program.

A new approach incorporating flexible docking simulations and NMR data is presented for calculating the bound conformation of a ligand that interacts weakly with an enzyme. This approach consists of sampling directly the conformation of a flexible ligand inside a receptor active site containing surrounding flexible loops. To make this sampling efficient, a ligand-growing procedure has been adopted. Optimization of the ECEPP/3-plus-NOE constraint function is carried out by using a collective variable Monte Carlo minimization technique. Numerous energy minimizations are made possible for such a large system by using a Bezier splines energy grid technique. This new flexible docking approach was applied to determine the structure of a fibrinogen Aalpha-like peptide (7DFLAEGGGVRGPRV20) bound to an active site mutant of thrombin [thrombin(S195A)]. Structure calculations of the bound ligand, using 2D-transferred NOESY distance constraints in the DIANA program, showed that the N-terminal portion of the peptide (D7-R16) involves a chain reversal, whereas the C-terminal portion (G17-V20) adopts a fold that exists in several different orientations. In addition, the ECEPP/3 flexible docking package was used to assess the conformational variability of the ligand and surrounding 60D-insertion loop of thrombin. Amino acid residues (17-20) of the peptide interact with a region of the enzyme that exhibits broad specificity, with a preferred direction between the 60D-insertion loop and Pro37 of thrombin.     

Differential accumulation of catecholamines, proenkephalin- and chromogranin A-derived peptides in the medium after chronic nicotine stimulation of cultured bovine adrenal chromaffin cells

We have used an antiserum to a synthetic peptide fragment of bovine chromogranin A (ChrgA)[Tyr0] bovine ChrgA (306-313): YLSKEWEDA, together with antibodies to proenkephalin-derived peptides, to measure the release of immunoreactive peptides from nicotine-stimulated cultured bovine adrenal chromaffin cells. Over a period of 6 hr the accumulation of YLSKEWEDA immunoreactivity and Met-enkephalin Arg6Gly7Leu8 (MERGL) immunoreactivity in the medium of 10 microM nicotine-stimulated cells was shown to be biphasic; the initial phase occurred in the first 15-30 min and the second phase reached a peak after 4 hr. In contrast, catecholamine release occurred monophasically over the initial 15-30 min. Investigation of the second phase of peptide accumulation revealed that it was due in part to nicotine-evoked exocytosis and in part to extracellular processing of high molecular weight precursor proteins.     

A dipeptidyl carboxypeptidase in brain synaptic membranes active in the metabolism of enkephalin containing peptides

A dipeptidyl carboxypeptidase activity has been localized in synaptic plasma membranes which have been prepared from isolated rat brain cortical synaptosomes. The specificity of this proteolytic activity towards various synthetic and biological active peptides is compared to the peptidase activities of intact synaptosomes. In contrast to the synaptosomal peptidases which are capable of cleaving all peptide bonds of Met-enkephalin-Arg6-Phe7 the peptidase activity associated with the synaptic plasma membrane exclusively hydrolyses a dipeptide from the carboxyl terminus of all hepta- and hexapeptides tested. The fact that this dipeptidyl carboxypeptidase does not cleave the Gly3-Phe4 peptide bond of Met-enkephalin suggests that this enzyme is different from "enkephalinase". The synaptic membrane dipeptidyl carboxypeptidase is inhibited by metal chelating agents and thiols but is not affected by compounds known to inhibit serine proteases, thermolysin and "enkephalinase".     

Comparison of the low energy conformations of an oncogenic and a non-oncogenic p21 protein, neither of which binds GTP or GDP

Oncogenic p21 protein, encoded by theras-oncogene, that causes malignant transformation of normal cells and many human tumors, is almost identical in sequence to its normal protooncogene-encoded counterpart protein, except for the substitution of arbitrary amino acids for the normally occurring amino acids at critical positions such as Gly 12 and Gin 61. Since p21 is normally activated by the binding of GTP in place of GDP, it has been postulated that oncogenic forms must retain bound GTP for prolonged time periods. However, two multiply substituted p21 proteins have been cloned, neither of which binds GDP or GTP. One of these mutant proteins with Val for Gly 10, Arg for Gly 12, and Thr for Ala 59 causes cell transformation, while the other, similar protein with Gly 10, Arg 12, Val for Gly 13 and Thr 59 does not transform cells. To define the critical conformational changes that occur in the p21 protein that cause it to become oncogenic, we have calculated the low energy conformations of the two multiply substituted mutant p21 proteins using a new adaptation of the electrostatically driven Monte Carlo (EDMC) technique, based on the program ECEPP. We have used this method to explore the conformational space available to both proteins and to compute the average structures for both using statistical mechanical averaging. Comparison of the average structures allows us to detect the major differences in conformation between the two proteins. Starting structures for each protein were calculated using the recently deposited x-ray crystal coordinates for the p21 protein, that was energy-refined using ECEPP, and then perturbed using the EDMC method to compute its average structure. The specific amino acid substitutions for both proteins were then generated into the lowest energy structure generated by this procedure, subjected to energy minimization and then to full EDMC perturbations. We find that both mutant proteins exhibit major differences in conformation in specific regions, viz., residues 35–47, 55–78, 81–93, 96–110, 115–126, and 123–134, compared with the EDMC-refined x-ray structure of the wild-type protein. These regions have been found to be the most flexible in the p21 protein bound to GDP from prior molecular dynamics calculations (Dykeset al., 1993). Comparison of the EDMC-average structure of the transforming mutant with that of the nontransforming mutant reveals major structural differences at residues 10–16, 32–40, and 60–68. These structural differences appear to be the ones that are critical in activation of the p21 protein. Analysis of the correlated motions of the different regions of the two mutant proteins reveals that changes in the conformation of regions in the carboxyl half of the protein are caused by changes in conformation around residues 10–16 and are transmitted by means of residues around Gln 61. The latter region therefore constitutes a molecular switch unit, in agreement with conclusions from prior work.On leave from the Department of Chemistry, University of Gdask, ul. Sobieskiego 18, 80-952 Gdask, Poland.     

A new force-field program for the calculation of glycopeptides and its application to a heptacosapeptide-decasaccharide of immunoglobulin G1. Importance of 1-6-glycosidic linkages in carbohydrate.peptide interactions

Energetically favored conformations of glycopeptide 1 were calculated using the newly developed force-field program, GEGOP (geometry of glycopeptides). The three-dimensional structure of glycopeptide 1, which is part of the Fc fragment of IgG1, has been calculated. 1 contains 27 amino acid residues from Pro291 to Lys317 and a biantennary decasaccharide N-linked to Asn297. The conformations of the peptide and the carbohydrate parts are shown to be mutually dependent. Single glycosyl residues of 1 exhibit interaction energies of up to -31.8 kJ/mol with the peptide portion. Generally, only a few of the glycosyl residues of the oligosaccharide moiety express significant interaction energies with the peptide part. No easy prediction is possible of glycosyl residues which exhibit favorable interaction energies. However, in all of the calculated structures, the glycosyl residues of the 1-6-linked branches show strong attractive forces for the peptide part. 1-6-glycosidically linked branches can adopt a larger number of conformations than other linkages due to their high flexibility which allows more favorable interactions with proteins. We developed the GEGOP program in order to be able to study the preferred conformations of large glycopeptides. The program is based on the GESA (geometry of saccharides) program and utilizes the HSEA (hard sphere exo anomeric) force field for the carbohydrate part and the ECEPP/2 (empirical conformation energy program for peptides) force field [Némethy, G., Pottle, M. S. & Scheraga, H. A. (1983) J. Phys. Chem. 87, 1883-1887] for the peptide part. The GEGOP program allows the simultaneous relaxation of all rotational degrees of freedom of these glycoconjugates during the energy optimization process. Thus, mutual interactions between glycosyl and amino acid residues can be studied in detail.