Glycoproteins: glycan presentation and protein-fold stability.

Glycosylation of proteins has been shown to play a role in a variety of cellular events. Thanks to recent advances in obtaining conformational constraints across glycosidic linkages, structural characterisation of glycoproteins has improved considerably. It is now becoming apparent that N-glycosylation of a folded protein can have a significant stabilising effect on large regions of the backbone structure.     

Conformation of 3'CMP bound to RNase A using TrNOESY

The conditions for accurately determining distance constraints from TrNOESY data on a small ligand (3'CMP) bound to a small protein (RNase A, <14 kDa) are described. For small proteins, normal TrNOESY conditions of 10:1 ligand:protein or greater can lead to inaccurate structures for the ligand-bound conformation due to the contribution of the free ligand to the TrNOESY signals. By using two ligand:protein ratios (2:1 and 5:1), which give the same distance constraints, a conformation of 3'CMP bound to RNase A was determined (glycosidic torsion angle, chi=-166 degrees ; pseudorotational phase angle, 0 degrees < or = P < or =36 degrees ). Ligand-protein NOESY cross peaks were also observed and used to dock 3'CMP into the binding pocket of the apo-protein (7rsa). After energy minimization, the conformation of the 3'CMP:RNase A complex was similar to the X-ray structure (1 rpf) except that a C3'-endo conformation for the ribose ring (rather than C2'-exo conformation) was found in the TrNOESY structure.     

An evaluation of the combined use of nuclear magnetic resonance and distance geometry for the determination of protein conformations in solution

An evaluation of the potential of nuclear magnetic resonance (n.m.r.) as a means of determining polypeptide conformation in solution is performed with the aid of a new distance geometry program which is capable of computing complete spatial structures for small proteins from n.m.r. data. Ten sets of geometric constraints which simulate the results available from n.m.r. experiments of varying precision and completeness were extracted from the crystal structure of the basic pancreatic trypsin inhibitor, and conformers consistent with these constraints were computed. Comparison of these computed structures with each other and with the original crystal structure shows that it is possible to determine the global conformation of a polypeptide chain from the distance constraints which are available from n.m.r. experiments. The results obtained with the different data sets also provide a standard by which the quality of protein structures computed from n.m.r. data can be evaluated when no crystal structure is available, and indicate directions in which n.m.r. experiments for protein structure determination could be further improved.     

Application of restrained minimization, simulated annealing and molecular dynamics simulations for the conformational analysis of oligosaccharides.

The purpose of the present study was to determine the confidence with which the small number of 1H NMR nuclear Overhauser effect (NOE) distance constraints measurable across glycosidic linkages in oligosaccharides could be used for solution conformational analysis. This was assessed by use of these constraints in restrained molecular mechanical minimization of the tetrasaccharide Gal beta 1----4(Fuc alpha 1----3)Glc-NAc beta 1----3Gal, a model compound of the Lewis-X antigenic determinant. This presents a particularly severe test case in view of extreme resonance overlap and a dearth of inter-residue distance constraints. It is concluded that these constraints, when used in conventional restrained minimization, result in the generation of 'virtual conformations' and local minima about glycosidic linkages. However, these restraints are nevertheless found to be useful in the initial stages of a conformational analysis strategy involving restrained minimization combined with dynamical simulated annealing to define more accurately the global minimum energy configuration, together with molecular dynamics simulation to explore conformational mobility about this minimum. Theoretical ROE values calculated over the time course of the MD simulation, using a formalism appropriate for the time scale of the internal motion, are compared with those obtained experimentally in the oligosaccharide.     

Conformational analysis of a dinucleotide photodimer with the aid of the genetic algorithm.

The solution structure of the photodimer cis,syn-dUp[]dT is derived with the aid of the genetic algorithm. The conformational space available for the molecule is sampled efficiently using the computer program DENISE and tested against a set of constraints available from nmr experiments. The dominant conformation in solution found with this approach can be described by the following combinations of sugar-phosphate backbone torsion angles: epsilon(t), zeta(t), alpha(+), beta(-ac), and gamma(t). The conformation of the sugars and glycosidic torsion angles are S type and syn, respectively. The cyclobutane ring and pyrimidines are puckered. In addition, other conformations that exist in equilibrium with the first are found. It is concluded that the cyclobutane-pyrimidine system is rigid, whereas the sugar-phosphate backbone is flexible. The solution structures are compared with the crystal structure of the strongly related cyano-ethyl ester of cis,syn-dTp[]dT.     

Conformational space of alpha 1-4 glycosidic linkage: a molecular dynamics study

Molecular Dynamics simulations have been carried out for 100 ps on crystal structure of beta-cyclodextrin in vacuo and with explicit inclusion of solvent at constant pressure and constant temperature using the GROMOS MD algorithm, with a time step of 0.005 ps. The conformational space of the glycosidic linkage was studied by calculating two virtual dihedrals connecting the successive glucose units for the 2000 structures saved during the two simulations. Three preferred regions for alpha 1-4 glycosidic linkage were found in both the simulations. The use of these virtual dihedral angles in representing the glycosidic linkage is also brought out from these studies.     

The solution structure of the intramolecular photoproduct of d(TpA) derived with the use of NMR and a combination of distance geometry and molecular dynamics.

One and two dimensional NMR techniques have been used together with molecular modelling to obtain the solution structure for the photoproduct d(TpA)*. The NMR data confirm that the cyclobutane linkage is formed between the bonds thymine C6-C5 and adenine C5-C6. The 2D NOE data are used as constraints in a distance geometry calculation. The structures obtained show a trans-syn cyclobutane linkage and the glycosidic angles are SYN and ANTI for thymidine and deoxyadenosine, respectively. The coupling constant data are used to check the backbone torsion angles of the obtained structures. Typical torsion angles are a gamma+ and beta t for the deoxyadenosine residue. A free molecular dynamics simulation of a trans-syn d(TpA) photoproduct confirmed all these structural characteristics.     

A statistical analysis of N- and O-glycan linkage conformations from crystallographic data

We have generated a database of 639 glycosidic linkage structures by an exhaustive survey of the available crystallographic data for isolated oligosaccharides, glycoproteins, and glycan-binding proteins. For isolated oligosaccharides there is relatively little crystallographic data available. A much larger number of glycoprotein and glycan-binding protein structures have now been solved in which two or more linked monosaccharides can be resolved. In the majority of these cases, only a few residues can be seen. Using the 639 glycosidic linkage structures, we have identified one or more distinct conformers for all the linkages. The O5-C1-O-C(x)' torsion angles for all these distinct conformers appear to be determined chiefly by the exo-anomeric effect. The Manalpha1-6Man linkage appears to be less restrained than the others, showing a wide degree of dispersion outside the ranges of the defined conformers. The identification of distinct conformers for glyco-sidic linkages allows "average" glycan structures to be modeled and also allows the easy identification of distorted glycosidic linkages. Such an analysis shows that the interactions between IgG Fc and its own N-linked glycan result in severe distortion of the terminal Galbeta1-4GlcNAc linkage only, indicating the strong interactions that must be present between the Gal residue and the protein surface. The applicability of this crystallographic based analysis to glycan structures in solution is discussed. This database of linkagestructures should be a very useful reference tool in three-dimensional structure determinations.     

Relationship between protein structure and geometrical constraints.

We evaluate to what extent the structure of proteins can be deduced from incomplete knowledge of disulfide bridges, surface assignments, secondary structure assignments, and additional distance constraints. A cost function taking such constraints into account was used to obtain protein structures using a simple minimization algorithm. For small proteins, the approximate structure could be obtained using one additional distance constraint for each amino acid in the protein. We also studied the effect of using predicted secondary structure and surface assignments. The constraints used in this approach typically may be obtained from low-resolution experimental data. When using a cost function based on distances, half of the resulting structures will be mirrored, because the resulting structure and its mirror image will have the same cost. The secondary structure assignments were therefore divided into chirality constraints and distance constraints. Here we report that the problem of mirrored structures, in some cases, can be solved by using a chirality term in the cost function.     

Solution conformation of a pectin fragment disaccharide using molecular modelling and nuclear magnetic resonance

In the present study, the conformational behaviour of methylated pectic disaccharide 4-O-alpha-D-galactopyranurosyl 1-O-methyl-alpha-D-galactopyranuronic 6,6'-dimethyl diester 1 has been completely characterized through combined n.m.r. and molecular modelling studies. The 1H-1H n.O.e. across the glycosidic bond was measured by both steady-state and transient 1D and 2D experiments. In parallel, the complete conformational analysis of the disaccharide has been achieved with the MM3 molecular mechanics method. The conformation of the pyranose ring is confirmed by the excellent agreement between the experimental and calculated intracyclic scalar coupling constants. The iso-energy contours displayed on the 'relaxed' map indicate an important flexibility about the glycosidic linkage. There is no significant influence of the methoxyl group on the conformational behaviour of the disaccharide. The theoretical n.m.r. data were calculated taking into account all the accessible conformations and using the averaging methods appropriate for slow internal motions. 3JC-H coupling constants were calculated using an equation suitable for C-O-C-H segments. The agreement between experimental and theoretical data is excellent. Within the potential energy surface calculated for the disaccharide, several conformers can be identified. When these conformations are extrapolated to a regular polymer structure, they generate pectins with right- and left-handed chirality along with a two-fold helix. These different types of helical structure are the result of small changes in conformation, without any drastic variation of the fibre repeat.     

Molecular cloning, purification, and characterization of a novel polyMG-specific alginate lyase responsible for alginate MG block degradation in Stenotrophomas maltophilia KJ-2

A gene for a polyMG-specific alginate lyase possessing a novel structure was identified and cloned from Stenotrophomas maltophilia KJ-2 by using PCR with homologous nucleotide sequences-based primers. The recombinant alginate lyase consisting of 475 amino acids was purified on Ni-Sepharose column and exhibited the highest activity at pH 8 and 40?°C. Interestingly, the recombinant alginate lyase was expected to have a similar catalytic active site of chondroitin B lyase but did not show chondroitin lyase activity. In the test of substrate specificity, the recombinant alginate lyase preferentially degraded the glycosidic bond of polyMG-block than polyM-block and polyG-block. The chemical structures of the degraded alginate oligosaccharides were elucidated to have mannuronate (M) at the reducing end on the basis of NMR analysis, supporting that KJ-2 polyMG-specific alginate lyase preferably degraded the glycosidic bond in M-G linkage than that in G-M linkage. The KJ-2 polyMG-specific alginate lyase can be used in combination with other alginate lyases for a synergistic saccharification of alginate.     

Reconstruction of Natural RNA Sequences from RNA Shape,Thermodynamic Stability,Mutational Robustness,and Linguistic Complexity by Evolutionary Computation

Abstract

The process of designing novel RNA sequences by inverse RNA folding, available in tools such as RNAinverse and InfoRNA, can be thought of as a reconstruction of RNAs from secondary structure. In this reconstruction problem, no physical measures are considered as additional constraints that are independent of structure, aside of the goal to reach the same secondary structure as the input using energy minimization methods. An extension of the reconstruction problem can be formulated since in many cases of natural RNAs, it is desired to analyze the sequence and structure of RNA molecules using various physical quantifiable measures. In prior works that used secondary structure predictions, it has been shown that natural RNAs differ significantly from random RNAs in some of these measures. Thus, we relax the problem of reconstructing RNAs from secondary structure into reconstructing RNAs from shapes, and in turn incorporate physical quantities as constraints. This allows for the design of novel RNA sequences by inverse folding while considering various physical quantities of interest such as thermodynamic stability, mutational robustness, and linguistic complexity. At the expense of altering the number of nucleotides in stems and loops, for example, physical measures can be taken into account. We use evolutionary computation for the new reconstruction problem and illustrate the procedure on various natural RNAs.     

Sequencing of Parts in Robotic Cells

This paper considers scheduling problems in robotic cells that produce a set of part types on several machines served by one robot. We study the problem of sequencing parts of different types in a cell to minimize the production cycle time when the sequence of the robot moves is given. This problem is NP-hard for most of the one-unit robot move cycles in a robotic cell with more than two machines and producing more than two part types. We first give a mathematical formulation to the problem, and then propose a branch-and-bound algorithm to solve it. The bounding scheme of this algorithm is based on relaxing, for all of the machines except two, the constraints that a machine should be occupied by a part for a period at least as long as the processing time of the part. The lower bound obtained in this way is tight. This relaxation allows us to overcome the complexity of the problem. The lower bound can be computed using the algorithm of Gilmore and Gomory. Computational experiments on part sequencing problems in three-machine robotic cells are given.     

Conformation of the glycotripeptide repeating unit of antifreeze glycoprotein of polar fish as determined from the fully assigned proton n.m.r. spectrum

With its simple glycotripeptide repeating structure the antifreeze glycoprotein of polar fish may be an especially simple conformational mode for mucin glycoproteins with similar but more complex structures. The fully assigned proton n.m.r. spectrum confirms the anomeric configurations of the hexapyranosidic sugars of the side chains and the coupling constants of the alpha GalNAc and the beta Gal residues show both to be in the expected 4C1 chair conformation. The assignment of a single resonance for each proton of the (Ala-Thr-Ala)n repeat unit coupled with the observation of long range nuclear Overhauser effects (n.O.e.) implies a three-fold repeating conformation. The resonances of the two alanines are distinct and can be assigned to their correct positions in the peptide sequence by n.O.e. observed at the amide proton resonances on saturation of the alpha proton signals. The amide proton coupling constants of all three peptide residues are similar and imply a limited range of peptide backbone torsion angles, phi CN. The large n.O.e. which has been observed between the amide proton and the alpha proton of the residue preceding it in the sequence implies large positive values for the peptide dihedral angle, psi CC. Limits are placed on possible values of side chain dihedral angles by the observation of the coupling constant between the alpha and beta protons of the threonyl residue. The observation of n.O.e. between the anomeric proton of GalNAc and the threonyl side chain protons gives information on the conformation of the alpha glycosidic linkage between the disaccharide and the peptide. n.O.e. observed between the protons of the beta glycosidic linkage indicates the conformation of the disaccharide and the large amide proton coupling constant of the GalNAc residue shows a trans proton relationship. The spectroscopically derived data have been combined with conformational energy calculations to give a conformational model for antifreeze glycoprotein in which the hydrophobic surfaces of the disaccharide side chains are wrapped closely against a three-fold left handed helical peptide backbone. The hydrophilic sides of the disaccharides are aligned so that they may bind to the ice crystal face, which is perpendicular to the fast growth axis inhibiting normal crystal growth.     

Reconstruction of natural RNA sequences from RNA shape, thermodynamic stability, mutational robustness, and linguistic complexity by evolutionary computation

The process of designing novel RNA sequences by inverse RNA folding, available in tools such as RNAinverse and InfoRNA, can be thought of as a reconstruction of RNAs from secondary structure. In this reconstruction problem, no physical measures are considered as additional constraints that are independent of structure, aside of the goal to reach the same secondary structure as the input using energy minimization methods. An extension of the reconstruction problem can be formulated since in many cases of natural RNAs, it is desired to analyze the sequence and structure of RNA molecules using various physical quantifiable measures. In prior works that used secondary structure predictions, it has been shown that natural RNAs differ significantly from random RNAs in some of these measures. Thus, we relax the problem of reconstructing RNAs from secondary structure into reconstructing RNAs from shapes, and in turn incorporate physical quantities as constraints. This allows for the design of novel RNA sequences by inverse folding while considering various physical quantities of interest such as thermodynamic stability, mutational robustness, and linguistic complexity. At the expense of altering the number of nucleotides in stems and loops, for example, physical measures can be taken into account. We use evolutionary computation for the new reconstruction problem and illustrate the procedure on various natural RNAs.     

Kinetics of hydrolysis of chitin/chitosan oligomers in concentrated hydrochloric acid

The kinetics of hydrolysis in concentrated hydrochloric acid (12.07 M) of the fully N-acetylated chitin tetramer (GlcNAc(4)) and the fully N-deacetylated chitosan tetramer (GlcN(4)) were followed by determining the amounts of the lower DP oligomers as a function of time. A theoretical model was developed to simulate the kinetics of hydrolysis of the three different glycosidic linkages in the tetramers. The model uses two different rate constants for the hydrolysis of the glycosidic bonds in the oligomers, assuming that the glycosidic bond next to one of the end residues are hydrolysed faster than the two other glycosidic linkages. The two rate constants were estimated by fitting model data to experimental results. The results show that the hydrolysis of the tetramers is a nonrandom process as the glycosidic bonds next to one of the end residues are hydrolysed 2.5 and 2.0 times faster as compared to the other glycosidic linkages in the fully N-acetylated and fully N-deacetylated tetramer, respectively. From previous results on other oligomers and the reaction mechanism, it is likely that the glycosidic bond that is hydrolysed fastest is the one next to the nonreducing end. The absolute values for the rate constants for the hydrolysis of the glycosidic linkages in GlcNAc(4) were found to be 50 times higher as compared to the glycosidic linkages in GlcN(4), due to the catalytic role of the N-acetyl group and the presence of the positively charged amino-group on the N-deacetylated sugar residue.     

Multiple RNA structure alignment

Ribonucleic Acid (RNA) structures can be viewed as a special kind of strings where characters in a string can bond with each other. The question of aligning two RNA structures has been studied for a while, and there are several successful algorithms that are based upon different models. In this paper, by adopting the model introduced in Wang and Zhang,(19) we propose two algorithms to attack the question of aligning multiple RNA structures. Our methods are to reduce the multiple RNA structure alignment problem to the problem of aligning two RNA structure alignments. Meanwhile, we will show that the framework of sequence center star alignment algorithm can be applied to the problem of multiple RNA structure alignment, and if the triangle inequality is met in the scoring matrix, the approximation ratio of the algorithm remains to be 2-2(over)n, where n is the total number of structures.     

Conformational analysis and molecular modelling of the branching point of amylopectin

The conformational properties of 6(2) alpha-D-glucosylmaltotriose have been studied using energy calculations that include van der Waals interactions, hydrogen bond stabilization, exo-anometric effect and torsional potential contributions. The calculations focused mainly on the conformational properties displayed at the alpha(1----6) linkage within the tetrasaccharide for which the conformational space is reported. The tetrasaccharide molecule was then considered as a model compound of the branching point in amylopectin. From molecular modelling, some basic structural features associated with branching were clearly established. It was found that, among the low energy arrangements, the side chain would fold back onto the main backbone, thereby producing dense three-dimensional structures in which a 'parallel' arrangement is achieved. The branching between two strands of the double helix, as found in the crystalline moiety of A and B starches, was further investigated. It was found that one particular set of conformations about the glycosidic linkages in the two different strands, could result in an arrangement such that strands could be connected through an alpha(1----6) glycosidic linkage, with a minimum of distortion. The three-dimensional features derived from the molecular modelling agree with the physical properties and mode of biogenesis within the starch granule; they are in accord with a 'cluster' type of structure.     

Systematic study of nuclear Overhauser effects vis-à-vis local helical parameters, sugar puckers, and glycosidic torsions in B DNA: insensitivity of NOE to local transitions in B DNA oligonucleotides due to internal structural compensations.

A method has been developed to solve structures of DNA oligomers in solution from the experimental NOE data. The method is a combination of two approaches: (1) full matrix NOESY simulations and (2) conformational calculations of DNA double helix based on generalized helical parameters. The process of the refinement of a solution structure does not involve NMR-derived interproton distance constraints; rather it consists of a direct fitting of a structure to the experimental NOE data, a weighted sum of energy, and R factor being under minimization. A helical parameters-based generation of DNA forms makes it possible to organize the search for the optimal structure more effectively, systematically varying starting conformations. The method has been used to calculate a structure for the self-complementary DNA hexamer GGATCC, which is consistent with the available experimental data. The structure belongs to the B family of forms, although the local structural heterogeneity is very strong. Sugar puckers vary from O4'-exo to C3'-exo; helical steps are open with different magnitudes toward the minor groove. Next, we have addressed the question of how uniquely the structure is defined by the existing NMR data. Different structural parameters have been systematically varied, and their effect on individual NOE's and the R factor has been studied. Two energetically conjugated parameters, sugar puckers and glycosidic angles, can be determined very reliably, because of the strong dependences of the intraresidue H6/H8 to H2'/H2'/H3' NOE's. In contrast, the local helical conformation of DNA and the geometry of base pairs proved to be underdetermined by the existing NOE information, because the effect of any helical parameter on interproton distances can be compensated by the concerted changes in other parameters.     

A study of nurse scheduling in Japan

Scheduling nurses to staff shifts is a major problem in hospitals. The necessity of maintaining a certain level of service and skill in the makeup of every shift, while balancing the workload among the nurses involved, is incredibly difficult. It is often impossible to develop a schedule which satisfies all the requirements despite the time and resources spent in the effort. This paper summarizes all our published research on nurse scheduling to date. The difficulties realized by our two investigations in Japan are shown first, together with a resulting scheduling problem. The nurse scheduling model based on the results is then described. In this model, all constraints are divided into two essentially different types; that which maintains a certain level of skill for each shift ('shift constraints') and that which concerns the workload for each nurse ('nurse constraints'). By classifying the constraints in this manner, we can determine what is affected by a specific constraint when the constraint is not satisfied. We developed efficient algorithms while taking advantage of the structure of this model. Finally, it is shown that our algorithm can solve this problem for a 2-shift system efficiently.