PubMed Health⌕ Search

Biomedical subjects

G Vergoten

Publications and source records attributed to G Vergoten.

At least 37 records · Page 2Linked to original sources

A new method to characterize hydrophobic organization of proteins: application to rational protein engineering of barnase.

We present a new algorithm for characterization of protein spatial structure basing on the molecular hydrophobicity potential approach. The method is illustrated by the analysis of three-dimensional structure of barnase and barnase-barstar complex. Current approach enables identification of amino acid residues situated in unfavorable environment (these residues may be "active" for binding), and to map quantitatively hydrophobic, hydrophilic and unfavorable hydrophobic-hydrophilic intra- and inter-molecular contacts involving backbone and side-chain segments of amino acid residues. Calculation of individual contributions of amino acid residues to such contacts permits identification of structurally-important residues. The contact plots obtained with molecular hydrophobicity potential calculations, provide easy rules to choose sites for mutations, which can increase a strength of intra- or inter-molecular hydrophobic interactions. The unfavorable hydrophobic-hydrophilic contact can be mutated to favorable hydrophobic, and already existing weak hydrophobic contact can be strengthen by increasing hydrophobicity of residues in contact. Basing on the analysis of the contact plots, we suggest several mutations of barnase which are supposed to increase intramolecular hydrophobic interactions, and thus might lead to increased stability of the protein. Part of these mutations was studied previously experimentally, and indeed stabilized barnase. The other of predicted mutations were not studied experimentally yet. Several new mutations of barnase and barstar are also proposed to enhance the hydrophobic interactions on their binding interface.

Algorithms↗

Atomic solvation parameters for proteins in a membrane environment. Application to transmembrane alpha-helices.

Several sets of atomic solvation parameters imitating: (i) nonpolar environment of hydrocarbon core of a membrane, (ii) aqueous solution, and (iii) weakly-polar solvents have been developed. The parameters have been incorporated into the ECEPP/2 and CHARMM force fields and employed in non-restrained Monte Carlo and molecular dynamics simulations of membrane-spanning alpha-helical peptides (segment A of bacteriorhodopsin, melittin). Through these simulations, the structure and energetics of the helices have been examined as a function of the solvation term in the potential energy function. For the peptides under study, the set (i) of atomic solvation parameters reveals good retention of the alpha-helical conformation. By contrast, the simulations in vacuum or with the parameters imitating a polar solvent (sets (ii) or (iii)) show fast helix destabilization and tight packing of the structure accompanied by significant decreasing of the surface area accessible to solvent. Increased helical propensity for amino acid residues, population of side-chain rotamers as well as hydrogen-bonding pattern in nonpolar membrane-like environment agree well with available experimental and computational data. The problems related to further applications of the membrane-mimicking sets of atomic solvation parameters to simulations of membrane proteins and peptides are addressed.

Bacteriorhodopsins↗

Hydrophobic organization of alpha-helix membrane bundle in bacteriorhodopsin.

The hydrophobic organization of the intramembrane alpha-helical bundle in bacteriorhodopsin (BRh) was assessed based on a new approach to characterization of spatial hydrophobic properties of transmembrane (TM) alpha-helical peptides. The method employs two independent techniques: Monte Carlo simulations of nonpolar solvent around TM peptides and analysis of molecular hydrophobicity potential on their surfaces. The results obtained by the two methods agree with each other and permit precise hydrophobicity mapping of TM peptides. Superimposition of such data on the experimentally derived spatial model of the membrane moiety together with 2D maps of hydrophobic hydrophilic contacts provide considerable insight into the hydrophobic organization of BRh. The helix bundle is stabilized to a large extent by hydrophobic interactions between helices--neighbors in the sequence of BRh, by long-range interactions in helix pairs C-E, C-F, and C-G, and by nonpolar contracts between retinal and helices C, D, E, F. Unlike globular proteins, no polar contacts between residues distantly separated in the sequence of BRh were found in the bundle. One of the most striking results of this study is the finding that the hydrophobic organization of BRh is significantly different from those in bacterial photoreaction centers. Thus, TM alpha-helices in BRh expose their most nonpolar sides to the bilayer as well as to the neighboring helices and to the interior of the bundle. Some of them contact lipids with their relatively hydrophilic surfaces. No correlation was found between disposition of the most hydrophobic and the most variable sides of the TM helices.

Bacteriorhodopsins↗

Recognition of transmembrane alpha-helical segments with environmental profiles.

A method for assessing the environmental properties of membrane-spanning alpha-helical peptides in proteins has been proposed. The algorithm employs a set of environmental preference parameters derived for amino acid residues based on the analysis of the 3-D structures of membrane domains in bacteriorhodopsin and photoreaction centers Rhodopseudomonas viridis and Rhodobacter sphaeroides. The resulting 3-D-1-D scores for transmembrane segments are significantly different from those derived for alpha-helices in globular proteins. The parameters obtained have been used to construct environmental profiles for membrane alpha-helices in bacteriorhodopsin and photoreaction centers. The profiles successfully recognize their own sequences in several specially designed large databases. The method has been applied to several membrane proteins with unknown spatial structures. Most of their membrane-spanning peptides were efficiently recognized by the profiles. The predicted environment of the residues in the membrane segments fits the experimental data well. The approach is independent of any homology data and can be employed to delineate the membrane segments of a protein with environmental characteristics close to those of bacteriorhodopsin and photoreaction centers. The alignment of these segments with the reference profiles provides a considerable amount of data about their lipid and protein exposure.

Algorithms↗

Amino acid residue: is it structural or functional?

A new approach is suggested for delineating the structural and functional amino acid residues in proteins with known three-dimensional structure, basing on the involvement of residues in intramolecular hydrophobic and hydrophilic interactions and additional information about the conservativity of the residues. The approach is applied to the families of homologous neurotoxins and cardiotoxins. The results obtained concerning the role of amino acid residues in both families of toxins accord well with the similarity of their fold, but different mechanisms of action. Current approach can be used for detailed characterization of protein spatial structures, as well as for rational protein engineering.

Amino Acid Sequence↗

Characterization of a new class of inhibitors of the recombinant human liver UDP-glucuronosyltransferase, UGT1*6.

The inhibitory effect of a series of novel structurally related compounds on the human UDP-glucuronosyltransferase UGT1*6 stably expressed in a V79 cell line was investigated. The inhibitors contain a lipophilic N-acyl phenylaminoalcohol residue and a uridine moiety connected by a spacer varying for each compound. The effects of these compounds on the glucuronidation reaction measured with 4-methylumbelliferone as substrate were determined. The best inhibitor of the series, D-DPMSU, had an IC50 of 39 microM in the assay conditions. Low Ki values were found toward both UDP-glucuronic acid and 4-methylumbelliferone (17 and 21 microM, respectively). The inhibition was competitive toward both substrates. A similar strong and competitive inhibitory effect was observed with two other inhibitors, DHPASU and DHPASiU. Another compound, D-DPASiU, showed a pure competitive inhibition towards UDP-glucuronic acid, but a non-competitive inhibition towards the acceptor substrate. These data and the optimization of the structures of the inhibitors by molecular modeling suggest that D-DPMSU and DHPASiU compounds may be transition state analog inhibitors of the recombinant UGT1*6 enzyme.

Binding Sites↗

Detailed assessment of spatial hydrophobic and electrostatic properties of 2D NMR-derived models of neurotoxin II.

2D NMR-derived spatial structures of neurotoxin II (NtII) and several homologous toxins in solution were assessed by comparison with their own amino acid sequences using a three-dimensional (3D) profile method. 3D profiles of all the toxin models match the sequences well and, therefore, the method of 3D profile was demonstrated to work correctly for these well-resolved NMR structures in aqueous solution. At the same time, the profile window plots reveal low scores in the bottom tip of loop II (residues 22-34 in NtII) and in beta-strand of loop III (residues 49-52). Some residues in the first poor-scoring region are of functional importance being involved in binding with nicotinic acetylcholine receptor (AChR). Furthermore, the second segment participates in intermolecular hydrogen bonding upon dimerization of postsynaptic neurotoxins in solution resulting in increasing of the 3D-1D score for residues at the interface between monomers. Therefore, the 3D profile method can be useful for detection functionally-important regions in well-resolved protein structures.

Algorithms↗

Molecular interactions between human lactotransferrin and the phytohemagglutinin-activated human lymphocyte lactotransferrin receptor lie in two loop-containing regions of the N-terminal domain I of human lactotransferrin.

Fluorescein isothiocyanate derivatization of the human lactotransferrin on Lys-264 inhibits the binding of the protein of human PHA-activated lymphocytes [Legrand, D., Mazurier, J., Maes, P., Rochard, E., Montreuil, J., & Spik, G. (1991) Biochem. J. 276, 733-738], indicating that part of the receptor-binding site is located in the N-terminal domain I of lactotransferrin. In the present study, a 6-kDa peptide (residues 4-52) was isolated from the N-terminal lobe of human lactotransferrin which inhibited the binding of the protein to its cell receptor. In addition, lactotransferrin was derivatized using sulfosuccinimidyl 2-(p-azidosalicylamido)ethyl-1,3'-dithiopropionate (SASD) and sulfosuccinimidyl 6-((4'-azido-2'-nitrophenyl)amino)hexanoate (sulfo-SANPAH), two heterobifunctional reagents generally used for receptor-ligand cross-linking. The azide group of these two reagents was inactivated by photolysis, and only the succinimidyl ester group was allowed to react with lysine residues of the protein. The binding of the derivatized lactotransferrins to the human lymphocyte receptor was assayed. SASD, which binds to Lys-74, was able to inhibit the binding of lactotransferrin to the cell receptor, in contrast to Lys-281-binding sulfo-SANPAH. Molecular modeling showed the position of SASD, sulfo-SANPAH, and fluorescein molecules at the surface of the protein and suggested that SASD and fluorescein could mask residues 4-6 and two loop-containing regions of human lactotransferrin (residues 28-34 and 38-45). The comparison of the primary and tertiary structures of human lactotransferrin and serotransferrin, which bind to specific cell receptors, shows that the above-mentioned regions, which are likely involved in protein-receptor interactions, possess specific structural features.

Amino Acid Sequence↗

Application of three-dimensional molecular hydrophobicity potential to the analysis of spatial organization of membrane domains in proteins: I. Hydrophobic properties of transmembrane segments of Na+, K(+)-ATPase.

A new computer-aided molecular modeling approach based on the concept of three-dimensional (3D) molecular hydrophobicity potential has been developed to calculate the spatial organization of intramembrane domains in proteins. The method has been tested by calculating the arrangement of membrane-spanning segments in the photoreaction center of Rhodopseudomonas viridis and comparing the results obtained with those derived from the X-ray data. We have applied this computational procedure to the analysis of interhelical packing in membrane moiety of Na+, K(+)-ATPase. The work consists of three parts. In Part I, 3D distributions of electrostatic and molecular hydrophobicity potentials on the surfaces of transmembrane helical peptides were computed and visualized. The hydrophobic and electrostatic properties of helices are discussed from the point of view of their possible arrangement within the protein molecule. Interlocation of helical segments connected with short extramembrane loops found by means of optimization of their hydrophobic/hydrophilic contacts is considered in Part II. The most probable 3D model of packing of helical peptides in the membrane domain of Na+, K(+)-ATPase is discussed in the final part of the work.

Amino Acid Sequence↗

Specific histamine release capacity of peptides selected from the modelized Der p I protein, a major allergen of Dermatophagoides pteronyssinus.

Dust mite allergens are considered as a major cause of allergic disease and as a risk factor for asthma. Der p I, a 222 amino-acid residue globular glycoprotein, is one of the major allergens from Dermatophagoides pteronyssinus (Dpt) mites. In this study, we have used predictive conventional algorithms (i.e. hydrophilicity, mobility, accessibility) and a three-dimensional model of Der p I derived from comparison to actinidin and papain to select continuous amino acid sequences as potential B cell epitopes. Four peptides, 52-71, 117-133, 176-187, 188-199 were synthesized. Their antigenic reactivity was investigated, mainly by measuring their capacity to induce in vitro histamine release. Results indicated that only Dpt-sensitive patients react specifically to Der p I-derived peptides and more frequently to 52-71 and 117-133. For each peptide, the intensity of response was dependent on the patient tested and on the peptide concn. The capacity of peptides to induce histamine release was demonstrated to be correlated with the serum level of anti-Der p I IgE (r = 0.86; p less than 10(-2)). Taken together these data emphasize, in Dpt-sensitive patients, the heterogeneity of the specific response to synthetic Der p I-derived peptides and underline the possible variety of epitopes belonging to the allergen Der p I.

Allergens↗

Molecular dynamics simulations of a monofucosylated biantennary glycan of the N-acetyllactosamine type: the human lactotransferrin glycan.

Molecular dynamics simulations were carried out to explore the conformational flexibility of the antennae of N-linked glycans. They were performed over 200 ps in vacuo on the complete disialylated monofucosylated biantennary glycan of the N-acetyllactosaminic type. Starting from a bird-conformation, the 3-D-structure evolved through 9 successive transitional states to a new, compact and energetically favorable conformation which had never been previously described. In this conformation, both antennae are organized in two coplanar loops rolled in a contrary direction and oriented perpendicularly to the plane of the di-N-acetyl chitobiose residue leading to a 'lobster conformation'. All the glycosidic linkages of the disialylated monofucosylated biantennary glycan, except the Fuc(alpha 1-6)GlcNAc(beta 1-), were modified by a phase transition. Particularly, the Man(beta 1-4) GlcNAc(beta 1-) linkage, which was previously described by NMR and X-ray diffraction as a rigid one, was involved in numerous conformational changes during 83 ps, even before the first transition phase. The freedom of mobility of the torsional angles of the Man(alpha 1-6)Man(beta 1-) linkage was limited, under these simulation conditions, to the angle psi which took three values: 30 degrees, 90 degrees and 180 degrees. Moreover, from 150 ps up to the end of the simulation, the value of the torsional angle omega of the NeuAc(alpha 2-6)Gal(beta 1-) linkage of the alpha-1,6-antenna continuously swung between 60 degrees and -60 degrees. Finally, we observed that the values of the torsional angles of the three linkages: NeuAc(alpha 2-6)Gal(beta 1-), Gal(beta 1-4)GlcNAc(beta 1-) and GlcNAc(beta 1-2)Man(beta 1-) of each of the two antennae were different, demonstrating their asymmetric conformation.

Amino Sugars↗

Effect of Urey-Bradley-Shimanouchi force field on the harmonic dynamics of proteins.

A normal mode analysis of bovine pancreatic trypsin inhibitor is carried out by using a Urey-Bradley-Shimanouchi potential energy function. The density of vibrational states, the magnitudes, and time scales of the atomic fluctuations are compared with experimental and theoretical results obtained by the more commonly used potential energy functions. The atomic fluctuations of Lys-15 are subject to extensive considerations as this residue is buried in the trypsin specificity pocket. It is found that Arg-17 is likely to be of importance in order to understand the way BPTI binds on trypsin.

Acetamides↗

Molecular modeling of a disialylated monofucosylated biantennary glycan of the N-acetyllactosamine type.

The conformations of a disialylated monofucosylated biantennary glycan of the N-acetyllactosamine type were analysed using the Tripos 5.3 force field from the Sybyl software currently used for molecular modelling. The conformation of each glycosidic linkage was calculated when included in oligosaccharide structures of up to 5 units and the influence of the glycosidic environment on the overall structure was measured. The study clearly shows that the conformation of a branched glycan cannot result from the simple addition of the different low energy conformers of each of the glycosidic linkages constituting the glycan structure. The asymmetrical conformation of the two antennae was demonstrated. The lowest energy conformations of the overall glycan structure were built and classified into 5 main models: the Y, T, bird and broken wing conformations already described and a new one called the 'back folded wing conformation.'

Amino Sugars↗

[Molecular modelling of glycans: three-dimensional structure and protein fraction interaction. The example of rabbit sero-transferrin].

On the basis of experimental data and of computer calculations using the Tripos 5.3 force field in order to examine the three-dimensional structures which are sterically feasible and the conformations which are energetically the most favourable, we have designed a program of molecular modelling of biantennary glycans of the N-acetyllactosaminic type (complex type). We demonstrate that, in absence of any interaction with the protein, a high number of glycan conformations exists which can be classified into five basic conformations, four of which have already been described. In fact, in addition to the Y-, T-, "bird" and "broken-wing" conformations, a "back-folded wing" conformation is energetically feasible. In contrast, the glycan linked to the protein is immobilized into only one conformation: the "broken-wing" conformation. Forming a bridge between the two lobes of the peptide chain, it probably contributes to the maintenance of the protein in a biologically active conformation.

Animals↗

Raman spectroscopy and normal vibrations of peptides. Characteristic normal modes of a type-II beta turn.

Raman and infrared spectra have been recorded on the peptide Gly-LPro-Gly-Gly in its zwitterionic form and on its protected homolog Z-Gly-LPro-Gly-Gly-OMe (Z = N-benzyloxycarbonyl and OMe = methyl ester). Both peptides adopt a type-II beta-turn conformation into appropriate solvents, that is in dimethylsulfoxide and chloroform respectively. A normal coordinate treatment has been performed on Gly-LPro-Gly-Gly using a modified Urey-Bradley force field and the characteristic normal modes of vibrations. Their corresponding frequencies have been extracted from experimental investigations.

Amides↗

Normal modes of vibration of a model peptide adopting the C7 C5 structure. Application toi angiotensin II.

A normal coordinate treatment was performed on the C7C5 conformation using a modified Urey-Bradley force field refined from previous studies on beta-turns. The predicted frequencies were compared with the experimental ones obtained on the peptide hormone human angiotensin II. The existence of both a beta-turn and a C7C5 conformation in solid and in aqueous solution is discussed.

Angiotensin II↗