Mean geometry of the thiopeptide unit and conformational features of dithiopeptides and polythiopeptides

The mean geometry of the thiopeptide [Ca-N-C(=S)-Ca] unit has been derived from an analysis of X-ray crystal structure data, as well as MM2 and Gaussian 80/82 calculations. The conformational flexibilities of dithiopeptides with glycl- and alanyl-side chains have been investigated by molecular mechanics. Minimum energy conformations were examined using interactive computer graphics molecular modeling techniques. Alanyl-dithiopeptide substitution within an oligopeptide results in considerable restriction of conformational freedom whereas the effect is minimal for glycyl-dithiopeptide substitution. Polyglycyl-thiopeptide adopts a left-handed three or fourfold or right-handed threefold helical structure with favorable interchain C = S...H-N hydrogen bond interactions. A poly-L-alanyl-thiopeptide prefers a left-handed threefold poly-L-proline-like helical structure.     

Modes of binding of 2'-AMP to RNase T1. A computer modeling study.

The modes of binding of adenosine 2'-monophosphate (2'-AMP) to the enzyme ribonuclease (RNase) T1 were determined by computer modelling studies. The phosphate moiety of 2'-AMP binds at the primary phosphate binding site. However, adenine can occupy two distinct sites--(1) The primary base binding site where the guanine of 2'-GMP binds and (2) The subsite close to the N1 subsite for the base on the 3'-side of guanine in a guanyl dinucleotide. The minimum energy conformers corresponding to the two modes of binding of 2'-AMP to RNase T1 were found to be of nearly the same energy implying that in solution 2'-AMP binds to the enzyme in both modes. The conformation of the inhibitor and the predicted hydrogen bonding scheme for the RNase T1-2'-AMP complex in the second binding mode (S) agrees well with the reported x-ray crystallographic study. The existence of the first mode of binding explains the experimental observations that RNase T1 catalyses the hydrolysis of phosphodiester bonds adjacent to adenosine at high enzyme concentrations. A comparison of the interactions of 2'-AMP and 2'-GMP with RNase T1 reveals that Glu58 and Asn98 at the phosphate binding site and Glu46 at the base binding site preferentially stabilise the enzyme-2'-GMP complex.     

Regionalisation of cell fate and morphogenetic movement of the mesoderm during mouse gastrulation

The developmental fate of cells in the epiblast of early-primitive-streak-stage mouse embryos was assessed by studying the pattern of tissue colonisation displayed by lac Z-expressing cells grafted orthotopically to nontransgenic embryos. Results of these fate-mapping experiments revealed that the lateral and posterior epiblast contain cells that will give rise predominantly to mesodermal derivatives. The various mesodermal populations are distributed in overlapping domains in the lateral and posterior epiblast, with the embryonic mesoderm such as heart, lateral, and paraxial mesoderm occupying a more distal position than the extraembryonic mesoderm. Heterotopic grafting of presumptive mesodermal cells results in the grafted cells adopting the fate appropriate to the new site, reflecting a plasticity of cell fate determination before ingression. The first wave of epiblast cells that ingress through the primitive streak are those giving rise to extraembryonic mesoderm. Cells that will form the mesoderm of the yolk sac and the amnion make up a major part of the mesodermal layer of the midprimitive-streak-stage embryo. Cells that are destined for embryonic mesoderm are still found within the epiblast, but some have been recruited to the distal portion of the mesoderm. By the late-primitive-streak-stage, the mesodermal layer contains only the precursors of embryonic mesoderm. This suggests that there has been a progressive displacement of the midstreak mesoderm to extraembryonic sites, which is reminiscent of that occurring in the overlying endodermal tissue. The regionalisation of cell fate in the late-primitive-streak mesoderm bears the same spatial relationship as their ancestors in the epiblast prior to cell ingression. This implies that both the position of the cells in the proximal-distal axis and their proximity to the primitive streak are major determinants for the patterning of the embryonic mesoderm. © 1995 Wiley-Liss, Inc.     

Molecular dynamics simulations of oligosaccharides and their conformation in the crystal structure of lectin-carbohydrate complex: importance of the torsion angle {psi} for the orientation of {alpha}1,6-arm

The conformation of the heptasaccharide Man-1,6-(Man-1,3)(Xyl-ß1,2)-Man-ß,4-GlcNAc₂-ß1,4-(L-Fuc-1,3)-GlcNAc₁,the carbohydrate moiety of Erythrina corallodendron lectin (EcorL),the hexasaccharide Man-1,6-(Man-1,3) (GlcNAc-ß1,4)-Man-ß1,4-GlcNAc-ß1,4-GlcNAcand their disaccharide fragments have been studied by moleculardynamics (MD) simulations for 1000 ps with different initialconformations. In the isolated heptasaccharide, the most frequentlyaccessed conformation during MD has a value of 180° aroundMan-1,6-Man linkage. This conformation is stabilized by theformation of a hydrogen bond between the carbonyl oxygen ofGlcNAc₂ with the O3/O4 hydroxyls of the 1,6-linked mannose residue.The conformation of the heptasaccharide found in the crystalstructure of the EcorL-lactose complex (Shaanan et al., Science,254, 862, 1991), that has a value of 76° around Man-1,6-Manlinkage, is accessed, although less frequently, during MD ofthe isolated oligosaccharide. The ,, = 58°,–134°,–60°conformation around Man-₁,6-Man fragment observed in the crystalstructure of the Lathyrus ochnrs lectin complexed with a biantennaryoctasaccharide (Table I in Homans,S.W., Glycobiology, 3, 551,1993) has also been accessed in the present MD simulations.These values for the 1,6-linkage, which are observed in theprotein-carbohydrate crystal structures and are accessed inthe MD simulations, though occasionally, have not been predictedfrom NMR studies. Furthermore, these different values of leadto significantly different orientations of the 1,6-arm for thesame value of . This contrasts with the earlier predictionsthat only different values of can bring about significant changesin the orientation of the ₁,6-arm. The MD simulations also showthat the effects of bisecting GlcNAc or ß1,2-xyloseare very similar on the ₁,3-arm and slightly different on the₁,6-arm. bisecting GlcNAc carbohydrates glycoprotein lectinsaccharide complex     

Molecular dynamics simulations of high-mannose oligosaccharides

Conformations of several high-mannose-type oligosaccharidesthat are generated during the biosynthetic degradation of Man₉GlcNAc₂to Man₅GlcNAc₂ have been studied by molecular dynamics (MD).Simulations were performed on NCI-FCRDC's Cray Y-MP 8D/8128supercomputer using Biosym's CVFF force field for 1000 Ps withdifferent initial conformations. The conformations of the two1,3- and the two 1,6-linkages in each oligomannose were different,suggesting that deriving oligosaccharide conformations basedon the conformational preferences of the constituent disaccharidefragments will not always yield correct results. Unlike otheroligomannoses, Man₉GlcNAc₂ appears to take more than one distinctconformation around the core 1,6-linkage. These various conformationsmay play an important role in determining the processing pathways.Using the data on the preferred conformations of these oligomannosesand the available experimental results, possible pathways forprocessing Man₉GlcNAc₂ to Man₅GlcNAc₂ by 1,2-linkage-specificmannosidases have been proposed. Conformational analysis ofMan₅GlcNAc₂ indicates that the addition of ß1,2-GlcNActo the 1,3-linked core mannose, besides serving as a prerequisitefor mannosidase II action as suggested earlier, may also preventthe removal of 1,3-mannose. The MD simulations also suggestthat the processing of the precursor oligosaccharide duringAsn-linked complex and hybrid glycan biosynthesis proceeds ina well-defined pathway involving more than one 1,2-linkage-specificmannosidase. Knowledge of the conformation of the processingintermediates obtained from the present study can be used todesign highly specific substrate analogues to inhibit a particularmannosidase, thereby blocking one processing pathway withoutinterfering with the others. carbohydrates conformation glycosidase inhibitors mannosidase oligosaccharide processing     

Functional domains of bovine β-1,4 galactosyltransferase

A number of N- and C-terminal deletion and point mutants of bovine -1,4 galactosyltransferase (-1,4GT) were expressed inE. coli to determine the binding regions of the enzyme that interact withN-acetylglucosamine (NAG) and UDP-galactose. The N-terminal truncated forms of the enzyme between residues 1–129, do not show any significant difference in the apparentK _ms toward NAG or linear oligosaccharide acceptors e.g. for chitobiose and chitotriose, or for the nucleotide donor UDP-galactose. Deletion or mutation of Cys 134 results in the loss of enzymatic activity, but does not affect the binding properties of the protein either to NAG- or UDP-agarose. From these columns the protein can be eluted with 15mm NAG and 50mm EDTA, like the enzymatically active protein, TL-GT129, that contains residues 130–402 of bovine -1,4GT. Also the N-terminus fragment, TL-GT129NAG, that contains residues 130–257 of the -1,4GT, binds to, and elutes with 15mm NAG and 50mm EDTA from the NAG-agarose column as efficiently as the enzymatically active TL-GT129. Unlike TL-GT129, the TL-GT129NAG binds to UDP-columns less efficiently and can be eluted from the column with only 15mm NAG. The C-terminus fragment GT-257UDP, containing residues 258–402 of -1,4GT, binds tightly to both NAG- and UDP-agarose columns. A small fraction, 5–10% of the bound protein, can be eluted from the UDP-agarose column with 50mm EDTA alone. The results show that the binding behaviour of N- and C-terminal fragments of -1,4GT towards the NAG- and UDP-agarose columns differ, the former binds preferentially to NAG-columns, while the latter binds to UDP-agarose columns via Mn²⁺.     

Molecular dynamics simulations of alpha2 --> 8-linked disialoside: conformational analysis and implications for binding to proteins.

Computational methods have played a key role in elucidating the various three-dimensional structures of oligosaccharides. Such structural information, together with other experimental data, leads to a better understanding of the role of oligosaccharide in various biological processes. The disialoside Neu5Ac-alpha2-->8-Neu5Ac appears as the terminal glycan in glycoproteins and glycolipids, and is known to play an important role in various events of cellular communication. Neurotoxins such as botulinum and tetanus require Neu5Ac-alpha2 --> 8-Neu5Ac for infecting the host. Glycoconjugates containing this disialoside and the enzymes catalyzing their biosynthesis are also regulated during cell growth, development, and differentiation. Unlike other biologically relevant disaccharides that have only two linkage bonds, the alpha2-->8-linked disialoside has four: C2-O, O-C8', C8'-C7', and C7'-C6'. The present report describes the results from nine 1 ns MD simulations of alpha2-->8-linked disialoside (Neu5Ac-alpha2-->8-Neu5Ac); simulations were run using GROMOS96 by explicitly considering the solvent molecules. Conformations around the O-C8' bond are restricted to the +sc/+ap regions due to stereochemical reasons. In contrast, conformations around the C2-O and C8'-C7' bonds were found to be largely unrestricted and all the three staggered regions are accessible. The conformations around the C7'-C6' bond were found to be in either the -sc or the anti region. These results are in excellent agreement with the available NMR and potential energy calculation studies. Overall, the disaccharide is flexible and adopts mainly two ensembles of conformations differing in the conformation around the C7'-C6' bond. The flexibility associated with this disaccharide allows for better optimization of intermolecular contacts while binding to proteins and this may partially compensate for the loss of conformational entropy that may be incurred due to disaccharide's flexibility.