Conformation of the exopolysaccharide of Burkholderia cepacia predicted with molecular mechanics (MM3) using genetic algorithm search

We present a computational conformational analysis of the exopolysaccharide of Burkholderia cepacia, which is believed to play a role in colonization and persistence of B. cepacia in the lungs of cystic fibrosis patients. The repeating unit of the exopolysaccharide is a heptasaccharide with three branches, which cause significant steric restraints. Conformational searches using glygal, an in-house developed software using genetic algorithm search methods, were performed on fragments as well as on the complete repeating unit with wrap-over residues. The force field used for the calculations was MM3(96). The search showed four favored conformations for an isolated repeating unit. However, for a sequence of several repeating units, the calculations indicate a single, well-defined linear conformation.     

A tool for the prediction of structures of complex sugars

In two recent back to back articles(Xia et al., J Chem Theory Comput 3:1620–1628 and 1629–1643, 2007a, b) we have started to address the problem of complex oligosaccharide conformation and folding. The scheme previously presented was based on exhaustive searches in configuration space in conjunction with Nuclear Overhauser Effect (NOE) calculations and the use of a complex rotameric library that takes branching into account. NOEs are extremely useful for structural determination but only provide information about short range interactions and ordering. Instead, the measurement of residual dipolar couplings (RDC), yields information about molecular ordering or folding that is long range in nature. In this article we show the results obtained by incorporation RDC calculations into our prediction scheme. Using this new approach we are able to accurately predict the structure of six human milk sugars: LNF-1, LND-1, LNF-2, LNF-3, LNnT and LNT. Our exhaustive search in dihedral configuration space combined with RDC and NOE calculations allows for highly accurate structural predictions that, because of the non-ergodic nature of these molecules on a time scale compatible with molecular dynamics simulations, are extremely hard to obtain otherwise (Almond et al., Biochemistry 43:5853–5863, 2004). Molecular dynamics simulations in explicit solvent using as initial configurations the structures predicted by our algorithm show that the histo-blood group epitopes in these sugars are relatively rigid and that the whole family of oligosaccharides derives its conformational variability almost exclusively from their common linkage (β-d-GlcNAc-(1→3)-β-d-Gal) which can exist in two distinct conformational states. A population analysis based on the conformational variability of this flexible glycosidic link indicates that the relative population of the two distinct states varies for different human milk oligosaccharides. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Efficient search on energy minima for structure prediction of nucleic acid motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Definition of the phage G13 receptor as structural domains of trisaccharides in Salmonella and Escherichia coli core oligosaccharides

The interaction between phage G13 and different bacterial and synthetic oligosaccharides has been studied using equilibrium dialysis inhibition. The results, and conformational analysis of the oligosaccharides, make us conclude that the phage G13 carbohydrate receptor is a conformational domain involving three sugar residues. The following trisaccharide elements contain the domain: alpha-D-Galp-(1----3)-[alpha-D-Galp-(1----6)]-alpha-D-Glcp, alpha-D-Manp-(1----3)-[alpha-D-Manp-(1----6)-alpha-D-Manp , and alpha-D-Glcp-(1----3)-[L-gly-alpha-D-man-Hepp-(1----7)]-L-gly-alph a-D- man-Hepp. Thus two structures, either a hexose substituted with alpha-D-glycopyranosyl groups in the 3- and 6-positions, or a heptose substituted with such groups in the 3- and 7-positions are functional G13 binding sites. Such domains are present in several cores of lipopolysaccharides from Salmonella and Escherichia coli species. Some cores, e.g. those from S. typhimurium chemotypes Ra, Rb1 and Rb2, contain two such domains. The identification of two G13 receptor domains within different core saccharides could explain the broad host range of this phage.     

Efficient Search on Energy Minima for Structure Prediction of Nucleic Acid Motifs

Abstract

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during con-formational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

The conformations of the O-specific polysaccharides of Shigella dysenteriae type 4 and Escherichia coli O159 studied with molecular mechanics (MM3) filtered systematic search

The branched O-antigens of Escherichia coli O159 and Shigella dysenteriae type 4 are structurally related and are known to show cross-reactivity with antibodies. In the present study, conformational analyses were performed on these two O-antigens using molecular mechanics MM3(96) with filtered systematic search. The results show very strong steric restrictions for the trisaccharide at the branch point of the E. coli O159 antigen, especially for the β-d-GlcNAc-(1 → 3)-β-d-GlcNAc linkage of the main chain. For the type 4 O-antigen the calculations show essentially a single conformation with respect to the α-d-GlcNAc-(1 → 3)-α-d-GlcNAc linkage of the main chain and three different favoured conformations for the fucose branch. Consecutive repeating units of the S. dysenteriae type 4 and E. coli O159 O-antigens form linear extended chains with significant flexibility between the branches. Comparative calculations carried out with the SWEET server indicate that our method of filtered systematic search is a superior method in the case of branched, constrained oligosaccharides. Based on the results of the MM3 calculations, we propose that the common epitope explaining the cross-reactivity comprises the fucose branch, the downstream GlcNAc and part of the uronic acid.     

Pharmacophoric pattern matching in files of three-dimensional chemical structures: use of bounded distance matrices for the representation and searching of conformationally flexible molecules

This paper discusses the use of bounded distance matrices for the representation of conformationally flexible threedimensional (3D) molecules. It is shown that pharmacophoric pattern searches of databases of flexible 3D molecules represented in this way can be carried out using screen and geometric searching algorithms that are analogous to those used for searching databases of rigid 3D structures. Molecules matching a query pattern after the geometric search must then undergo a final conformational search to determine whether they can, in fact, adopt a conformation that matches the query. An analysis of this three-stage searching procedure shows that searching databases of flexible 3D molecules is extremely demanding of computational resources.     

A computational comparison of the atomic models of the actomyosin interface

Several atomic models of the actomyosin interface have been proposed based on the docking together of their component structures using electron microscopy and resonance energy-transfer measurements. Although these models are in approximate agreement in the location of the binding interfaces when myosin is tightly bound to actin, their relationships to molecular docking simulations based on computational free-energy calculations are investigated here. Both rigid-docking and flexible-docking conformational search strategies were used to identify free-energy minima at the interfaces between atomic models of myosin and actin. These results suggest that the docking model produced by resonance energy-transfer data is closer to a free-energy minimum at the interface than are the available atomic models based on electron microscopy. The conformational searches were performed using both scallop and chicken skeletal muscle myosins and identified similarly oriented actin-binding interfaces that serve to validate that these models are at the global minimum. These results indicate that the existing docking models are close to but not precisely at the lowest-energy initial contact site for strong binding between myosin and actin that should represent an initial contact between the two proteins; therefore, conformational changes are likely to be important during the transition to a strongly bound complex.     

Conformational sampling using high-temperature molecular dynamics

High-temperature molecular dynamics as a method for conformational search was explored on the antigen combining site of McPC 603, a phosphorylcholine binding immunoglobulin. Simulations at temperatures of 500, 800, and 1500 K were run for 111.5, 101.7, and 76.3 ps, respectively. The effectiveness of the search was assessed using a variety of methods. For the shorter hypervariable loops, molecular dynamics explored an appreciable fraction of the conformational space as evidenced by a comparison to a simple theoretical model of the size of the conformational space. However, for the longer loops and the antigen combining site as a whole, the simulation times were too short for a complete search. The simulations at 500 and 800 K both generated conformations that minimized to energies 200 kcal/mole lower than the crystal structure. However, the 1500 K simulation produced higher energy structures, even after minimization; in addition, this highest temperature run had many cis-trans peptide isomerizations. This suggests that 1500 K is too high a temperature for unconstrained conformational sampling. Comparison of the results of high temperature molecular dynamics with a direct conformational search method, [R. E. Bruccoleri & M. Karplus (1987) Biopolymers 26, 137-168]. showed that the two methods did not overlap much in conformational space. Simple geometric measures of the conformational space indicated that the direct method covered more space than molecular dynamics at the lower temperature, but not at 1500 K. The results suggest that high-temperature molecular dynamics can aid in conformational searches.     

Empirical solvation models in the context of conformational energy searches: application to bovine pancreatic trypsin inhibitor.

Continuum solvation models that estimate free energies of solvation as a function of solvent accessible surface area are computationally simple enough to be useful for predicting protein conformation. The behavior of three such solvation models has been examined by applying them to the minimization of the conformational energy of bovine pancreatic trypsin inhibitor. The models differ only with regard to how the constants of proportionality between free energy and surface area were derived. Each model was derived by fitting to experimentally measured equilibrium solution properties. For two models, the solution property was free energy of hydration. For the third, the property was NMR coupling constants. The purpose of this study is to determine the effect of applying these solvation models to the nonequilibrium conformations of a protein arising in the course of global searches for conformational energy minima. Two approaches were used: (1) local energy minimization of an ensemble of conformations similar to the equilibrium conformation and (2) global search trajectories using Monte Carlo plus minimization starting from a single conformation similar to the equilibrium conformation. For the two models derived from free energy measurements, it was found that both the global searches and local minimizations yielded conformations more similar to the X-ray crystallographic structures than did searches or local minimizations carried out in the absence of a solvation component of the conformational energy. The model derived from NMR coupling constants behaved similarly to the other models in the context of a global search trajectory. For one of the models derived from measured free energies of hydration, it was found that minimization of an ensemble of near-equilibrium conformations yielded a new ensemble in which the conformation most similar to the X-ray determined structure PTI4 had the lowest total free energy. Despite the simplicity of the continuum solvation models, the final conformation generated in the trajectories for each of the models exhibited some of the characteristics that have been reported for conformations obtained from molecular dynamics simulations in the presence of a bath of explicit water molecules. They have smaller root mean square (rms) deviations from the experimentally determined conformation, fewer incorrect hydrogen bonds, and slightly larger radii of gyration than do conformations derived from search trajectories carried out in the absence of solvent.     

A theoretical study of pentacyclo-undecane cage peptides of the type [Ac-X-Y-NHMe].

The conformational preferences of peptides of the type, Ac-X-Y-NHMe, where X and Y = Ala, cage and Pro, were studied by means of computational techniques within the framework of a molecular mechanics approach. For each of the eight peptide analogues, extensive conformational searches were carried out using molecular dynamics (MD) and simulated annealing (SA) protocols in an iterative fashion. Both results are in good agreement and complement each other. The conformational search indicates that the cage residue restricts the conformational freedom of the dipeptide considerably in comparison with the other model residues used. This study revealed that proline exhibits a greater tendency in promoting reverse-turn characteristics in comparison to the cage peptides, which show promising beta-turn characteristics. It was also found that 300-500 K is not sufficient to overcome rotational barriers for cage peptides. In all cases, the low-energy conformers have a tendency to form bent structures.     

SimShiftDB; local conformational restraints derived from chemical shift similarity searches on a large synthetic database

We present SimShiftDB, a new program to extract conformational data from protein chemical shifts using structural alignments. The alignments are obtained in searches of a large database containing 13,000 structures and corresponding back-calculated chemical shifts. SimShiftDB makes use of chemical shift data to provide accurate results even in the case of low sequence similarity, and with even coverage of the conformational search space. We compare SimShiftDB to HHSearch, a state-of-the-art sequence-based search tool, and to TALOS, the current standard tool for the task. We show that for a significant fraction of the predicted similarities, SimShiftDB outperforms the other two methods. Particularly, the high coverage afforded by the larger database often allows predictions to be made for residues not involved in canonical secondary structure, where TALOS predictions are both less frequent and more error prone. Thus SimShiftDB can be seen as a complement to currently available methods.     

Profiling the morphological distribution of O-linked oligosaccharides

The morphological distribution of oligosaccharides is determined in the egg jelly surrounding Xenopus laevis eggs. This biological system is used to illustrate a method for readily identifying and quantifying oligosaccharides in specific tissues. The extracellular matrix surrounding X. laevis eggs consists of a vitelline envelope and a jelly coat. The jelly coat contains three morphologically distinct layers designated J1, J2, and J3 from the innermost to the outermost and is composed of 9-11 distinct glycoproteins. Each jelly layer is known to have specific functions in the fertilization of the egg. We developed a rapid method to separate and identify the oligosaccharides from X. laevis egg jelly layers. Identification was based on the retention times in high-performance liquid chromatography (porous graphitized carbon column), exact masses, and tandem mass spectrometry. Over 40 neutral and 30 sulfated oligosaccharides were observed in the three jelly layers. Neutral oligosaccharide structures from different jelly layers were both unique and overlapping, while sulfated oligosaccharides were detected only in layers J1 and J2. Neutral oligosaccharides unique to jelly layer J3 and the combined layers J1+J2 had similar core structures and similar residues. However, differences between these two sets of unique oligosaccharides were also observed and were primarily due to the branching carbohydrate moieties rather than the core structures.     

Studies on the conformational behaviour of GlcNAc-Man3-GlcNAc2 oligosaccharides using molecular dynamics simulations

Three-dimensional structures of the natural substrate unit for the enzyme N-acetylglucosamine-transferase II, GlcNAc-Man3-GlcNAc2, were investigated by molecular modelling methods. Molecular dynamics (MD) and molecular mechanics calculations on two hexasaccharides, namely GlcNAc-Man3-GlcNAc2-Asn and GlcNAc-Man3-GlcNAc2-OMe were performed by the Biosym/MSI software using the CVFF and CFF95 force fields in vacuum. The MD simulations were calculated for 3 ns at different simulation temperatures and for two values of dielectric constant, = 1 and = 4. From each 3 ns trajectory, 3050 structures have been optimized. The local minima obtained have been clustered into families exhibiting similar values of glycosidic torsional angles phi, psi, and omega. The influence of the simulation conditions and force fields used on the conformational behaviour and structure of the title oligosaccharides is discussed.     

Tackling the conformational sampling of larger flexible compounds and macrocycles in pharmacology and drug discovery

Computational conformational sampling underpins much of molecular modeling and design in pharmaceutical work. The sampling of smaller drug-like compounds has been an active area of research. However, few studies have tested in details the sampling of larger more flexible compounds, which are also relevant to drug discovery, including therapeutic peptides, macrocycles, and inhibitors of protein–protein interactions. Here, we investigate extensively mainstream conformational sampling methods on three carefully curated compound sets, namely the ‘Drug-like’, larger ‘Flexible’, and ‘Macrocycle’ compounds. These test molecules are chemically diverse with reliable X-ray protein-bound bioactive structures. The compared sampling methods include Stochastic Search and the recent LowModeMD from MOE, all the low-mode based approaches from MacroModel, and MD/LLMOD recently developed for macrocycles. In addition to default settings, key parameters of the sampling protocols were explored. The performance of the computational protocols was assessed via (i) the reproduction of the X-ray bioactive structures, (ii) the size, coverage and diversity of the output conformational ensembles, (iii) the compactness/extendedness of the conformers, and (iv) the ability to locate the global energy minimum. The influence of the stochastic nature of the searches on the results was also examined. Much better results were obtained by adopting search parameters enhanced over the default settings, while maintaining computational tractability. In MOE, the recent LowModeMD emerged as the method of choice. Mixed torsional/low-mode from MacroModel performed as well as LowModeMD, and MD/LLMOD performed well for macrocycles. The low-mode based approaches yielded very encouraging results with the flexible and macrocycle sets. Thus, one can productively tackle the computational conformational search of larger flexible compounds for drug discovery, including macrocycles.     

Application of Enhanced Sampling Monte Carlo Methods for High-Resolution Protein-Protein Docking in Rosetta

The high-resolution refinement of docked protein-protein complexes can provide valuable structural and mechanistic insight into protein complex formation complementing experiment. Monte Carlo (MC) based approaches are frequently applied to sample putative interaction geometries of proteins including also possible conformational changes of the binding partners. In order to explore efficiency improvements of the MC sampling, several enhanced sampling techniques, including temperature or Hamiltonian replica exchange and well-tempered ensemble approaches, have been combined with the MC method and were evaluated on 20 protein complexes using unbound partner structures. The well-tempered ensemble method combined with a 2-dimensional temperature and Hamiltonian replica exchange scheme (WTE-H-REMC) was identified as the most efficient search strategy. Comparison with prolonged MC searches indicates that the WTE-H-REMC approach requires approximately 5 times fewer MC steps to identify near native docking geometries compared to conventional MC searches.     

Characterization of N-linked oligosaccharides in chorion peroxidase of Aedes aegypti mosquito

A peroxidase is present in the chorion of Aedes aegypti eggs and catalyzes chorion protein cross-linking during chorion hardening, which is critical for egg survival in the environment. The unique chorion peroxidase (CPO) is a glycoprotein. This study deals with the N-glycosylation site, structures, and profile of CPO-associated oligosaccharides using mass spectrometric techniques and enzymatic digestion. CPO was isolated from chorion by solubilization and several chromatographic methods. Mono-saccharide composition was analyzed by HPLC with fluorescent detection. Our data revealed that carbohydrate (D-mannose, N-acetyl D-glucosamine, D-arabinose, N-acetyl D-galactosamine, and L-fucose) accounted for 2.24% of the CPO molecular weight. A single N-glycosylation site (Asn328-Cys- Thr) was identified by tryptic peptide mapping and de novo sequencing of native and PNGase A-deglycosylated CPO using matrix-assisted laser/desorption/ionization time-of-flight mass spectrometry (MALDI/TOF/MS) and liquid chromatography/tandem mass spectrometry (LC/MS/MS). The Asn328 was proven to be a major fully glycosylated site. Potential tryptic glycopeptides and profile were first assessed by MALDI/TOF/MS and then by precursor ion scanning during LC/MS/MS. The structures of N-linked oligosaccharides were elucidated from the MS/MS spectra of glycopeptides and exoglycosidase sequencing of PNGase A-released oligosaccharides. These CPO-associated oligosaccharides had dominant Man3GlcNAc2 and Man3 (Fuc) GlcNAc2 and high mannose-type structures (Man(4-8)GlcNAc2). The truncated structures, Man2GlcNAc2 and Man2 (Fuc) GlcNAc2, were also identified. Comparison of CPO activity and Stokes radius between native and deglycosylated CPO suggests that the N-linked oligosaccharides influence the enzyme activity by stabilizing its folded state.     

Factors influencing accuracy of computer-built models: a study based on leucine zipper GCN4 structure.

A three-dimensional model of the leucine zipper GCN4 built from its amino acid sequence had been reported previously by us. When the two alternative x-ray structures of the GCN4 dimer became available, the root mean square (r.m.s.) shifts between our model and the structures were determined as approximately 2.7 A on all atoms. These values are similar to the r.m.s. shift of 2.8 A between the two GCN4 structures in the different crystal forms (C2 and P2(1)2(1)2(1)). CONGEN conformational searches were run to better understand the conditions that may determine the preference of different conformers in different environments and to test the sensitivity of our current modeling techniques. With a judicious choice of CONGEN search parameters, the backbone r.m.s. deviation improved to 0.8 A and 2.5 A on all atoms. The side-chain conformations of Val and Leu at the helical interface were well reproduced (1.2 A r.m.s.), and the large side-chain misplacements occurred with only a small number of charged amino acids and a tyrosine. Inclusion of the crystal environment (C2 symmetry), as a passive background, into the side-chain conformational search further improved the accuracy of the model to an r.m.s. deviation of 2.1 A. Conformational searches carried out in the two different crystal environments and employing the AMBER protein/DNA forcefield, as implemented in CONGEN, gave the r.m.s. values of 2.2 A (for the C2 symmetry) and 2.5 A (for the P2(1)2(1)2(1) symmetry). In the C2 symmetry crystal, as much as 40% of the surface of each dimer was involved in crystal contacts with other dimers, and the charged residues on the surface often interacted with immobilized water molecules. Thus, occasional large r.m.s. deviations between the model and the x-ray side chains were due to specific conditions that did not occur in solution.     

Determination of the conformation of Lewis blood group oligosaccharides by simulation of two-dimensional nuclear Overhauser data

Through control of both the nmr probe temperature and of the solvent viscosity, phase-sensitive two-dimensional 1H nuclear Overhauser data (NOESY) at 300 and 500 MHz are obtained with excellent signal-to-noise ratios for Lewis blood group penta- and hexasaccharides isolated from human milk. Relatively long mixing times are required to produce measurable NOE intensities in these oligosaccharides, which makes a full relaxation matrix analysis necessary. By measurements of selective T1 for a few isolated 1H resonances, it was possible to generate a simulation of the complete NOESY spectrum at arbitrary mixing time for comparison with the experimental data. From an exhaustive search of the conformational space, it was found that only a small range of glycosidic dihedral angles of the nonreducing terminal Lewis blood group determinant fragments of the milk oligosaccharides LNF-2 and LND-1 produce simulated spectra agreeing within experimental error to the data. Conformational energy calculations reveal that each of these conformations is also one of minimum energy. It is concluded that the Lewis(a) and Lewis(b) oligosaccharides adopt relatively compact rigid structures in solution, as shown by the observation of cross peaks between protons in nonadjacent residues. Like the blood group A and H oligosaccharides, there exists only a small dependence of the conformation for Lewis(a) and Lewis(b) oligosaccharides on solvent. The apparent lack of dependence of conformation of these oligosaccharides on DMSO in D2O suggests that modification of solvent viscosity with mixtures of DMSO:D2O may provide a useful general strategy of NOESY studies of oligosaccharides.