Computation of the structure of collagen molecule fragments

Step-by-step computations of octapeptide structures was performed for human collagens I and III. It was shown that computational results (coordinates of atoms) practically coincide with X-ray data for the collagen fragment.     

NMR solution structure of a parallel LNA quadruplex

The solution structure of a locked nucleic acid (LNA) quadruplex, formed by the oligomer d(TGGGT), containing only conformationally restricted LNA residues is reported. NMR and CD spectroscopy, as well as molecular dynamics and mechanic calculations, has been used to characterize the complex. The molecule adopts a parallel stranded conformation with a 4-fold rotational symmetry, showing a right-handed helicity and the guanine residues in an almost planar conformation with three well-defined G-tetrads. The thermal stability of Q-LNA has been found to be comparable with that of [r(UGGGU)]₄, while a T_m increment of 20°C with respect to the corresponding DNA quadruplex structure [d(TGGGT)]₄ has been observed. The structural features of the LNA quadruplex reported here may open new perspectives for the biological application of LNAs as novel versatile tools to design aptamer or catalyst oligonucleotides.     

Refinement of herpesvirus B-capsid structure on parallel supercomputers.

Electron cryomicroscopy and icosahedral reconstruction are used to obtain the three-dimensional structure of the 1250-A-diameter herpesvirus B-capsid. The centers and orientations of particles in focal pairs of 400-kV, spot-scan micrographs are determined and iteratively refined by common-lines-based local and global refinement procedures. We describe the rationale behind choosing shared-memory multiprocessor computers for executing the global refinement, which is the most computationally intensive step in the reconstruction procedure. This refinement has been implemented on three different shared-memory supercomputers. The speedup and efficiency are evaluated by using test data sets with different numbers of particles and processors. Using this parallel refinement program, we refine the herpesvirus B-capsid from 355-particle images to 13-A resolution. The map shows new structural features and interactions of the protein subunits in the three distinct morphological units: penton, hexon, and triplex of this T = 16 icosahedral particle.     

Massively parallel interrogation of aptamer sequence, structure and function

Background

Optimization of high affinity reagents is a significant bottleneck in medicine and the life sciences. The ability to synthetically create thousands of permutations of a lead high-affinity reagent and survey the properties of individual permutations in parallel could potentially relieve this bottleneck. Aptamers are single stranded oligonucleotides affinity reagents isolated by in vitro selection processes and as a class have been shown to bind a wide variety of target molecules.

Methodology/Principal Findings

High density DNA microarray technology was used to synthesize, in situ, arrays of approximately 3,900 aptamer sequence permutations in triplicate. These sequences were interrogated on-chip for their ability to bind the fluorescently-labeled cognate target, immunoglobulin E, resulting in the parallel execution of thousands of experiments. Fluorescence intensity at each array feature was well resolved and shown to be a function of the sequence present. The data demonstrated high intra- and inter-chip correlation between the same features as well as among the sequence triplicates within a single array. Consistent with aptamer mediated IgE binding, fluorescence intensity correlated strongly with specific aptamer sequences and the concentration of IgE applied to the array.

Conclusion and Significance

The massively parallel sequence-function analyses provided by this approach confirmed the importance of a consensus sequence found in all 21 of the original IgE aptamer sequences and support a common stem:loop structure as being the secondary structure underlying IgE binding. The microarray application, data and results presented illustrate an efficient, high information content approach to optimizing aptamer function. It also provides a foundation from which to better understand and manipulate this important class of high affinity biomolecules.     

Crystal structure of a naturally occurring parallel right-handed coiled coil tetramer

The crystal structure of a polypeptide chain fragment from the surface layer protein tetrabrachion from Staphylothermus marinus has been determined at 1.8 A resolution. As proposed on the basis of the presence of 11-residue repeats, the polypeptide chain fragment forms a parallel right-handed coiled coil structure. Complementary hydrophobic interactions and complex networks of surface salt bridges result in an extremely thermostable tetrameric structure with remarkable properties. In marked contrast to left-handed coiled coil tetramers, the right-handed coiled coil reveals large hydrophobic cavities that are filled with water molecules. As a consequence, the packing of the hydrophobic core differs markedly from that of a right-handed parallel coiled coil tetramer that was designed on the basis of left-handed coiled coil structures.     

Estimating the twist of beta-strands embedded within a regular parallel beta-barrel structure

The parallel beta-barrel is a recurrent structural motif found in a large variety of different enzymes belonging to the family of alpha/beta-proteins. It has been shown previously that the hyperboloid can be considered as a scaffold describing the parallel beta-barrel structure. To assess restraints on beta-strand twist imposed by a given scaffold geometry, the notion of scaffold twist, Ts, is introduced. From Ts, the beta-strand twist (Tw beta) expected for a given scaffold geometry can be derived and it is verified that this computed twist can be used to identify beta-barrels characterized by good hydrogen bonding. It is noted that Tw beta is only slightly affected for beta-barrels differing in the number (N) of beta-strands, suggesting that restraints on main-chain conformation of beta-strands are not likely to account for the N = 8 invariability observed in natural parallel beta-barrels thereby strengthening previous work rationalizing this constancy.     

A parallel algorithm for estimating the secondary structure in ribonucleic acids

A parallel algorithm for estimating the secondary structure of an RNA molecule is presented in this paper. The mathematical problem to compute an optimal folding based on free-energy minimization is mapped onto a graph planarization problem. In the planarization problem we want to maximize the number of edges in a plane with no two edges crossing each other. To solve a sequence of n bases, n(n — 1)/2 processing elements are used in our algorithm.     

GTfold: Enabling parallel RNA secondary structure prediction on multi-core desktops

ABSTRACT: BACKGROUND: Accurate and efficient RNA secondary structure prediction remains an important open problem in computational molecular biology. Historically, advances in computing technology have enabled faster and more accurate RNA secondary structure predictions. Previous parallelized prediction programs achieved significant improvements in runtime, but their implementations were not portable from niche high-performance computers or easily accessible to most RNA researchers. With the increasing prevalence of multi-core desktop machines, a new parallel prediction program is needed to take full advantage of today's computing technology. FINDINGS: We present here the first implementation of RNA secondary structure prediction by thermodynamic optimization for modern multi-core computers. We show that GTfold predicts secondary structure in less time than UNAfold and RNAfold, without sacrificing accuracy, on machines with four or more cores. CONCLUSIONS: GTfold supports advances in RNA structural biology by reducing the timescales for secondary structure prediction. The difference will be particularly valuable to researchers working with lengthy RNA sequences, such as RNA viral genomes.     

L-alanyl-L-alanyl-L-alanine: parallel pleated sheet arrangement in unhydrated crystal structure, and comparisons with the antiparallel sheet structure

The tripeptide L-alanyl-L-alanyl-L-alanine has been crystallized from a water/dimethylformamide solution in an unhydrated form, with cell dimensions a = 11.849, b = 10.004, c = 9.862 A, beta = 101.30 degrees, monoclinic space group P21 with 4 molecules per cell (2 independent molecules in the asymmetric unit). The structure was determined by direct methods and refined to a discrepancy index R = 0.057. The tri-L-alanine molecules are packed in a parallel pleated sheet arrangement with unusually long amide nitrogen-carbonyl oxygen contacts within sheets. Comparisons are made with the antiparallel pleated sheet structure of tri-L-alanine hemihydrate, previously crystallized from the same solvent system.     

Amyloids of shuffled prion domains that form prions have a parallel in-register beta-sheet structure

The [URE3] and [PSI (+)] prions of Saccharomyces cerevisiae are self-propagating amyloid forms of Ure2p and Sup35p, respectively. The Q/N-rich N-terminal domains of each protein are necessary and sufficient for the prion properties of these proteins, forming in each case their amyloid cores. Surprisingly, shuffling either prion domain, leaving amino acid content unchanged, does not abrogate the ability of the proteins to become prions. The discovery that the amino acid composition of a polypeptide, not the specific sequence order, determines prion capability seems contrary to the standard folding paradigm that amino acid sequence determines protein fold. The shuffleability of a prion domain further suggests that the beta-sheet structure is of the parallel in-register type, and indeed, the normal Ure2 and Sup35 prion domains have such a structure. We demonstrate that two shuffled Ure2 prion domains capable of being prions form parallel in-register beta-sheet structures, and our data indicate the same conclusion for a single shuffled Sup35 prion domain. This result confirms our inference that shuffleability indicates parallel in-register structure.     

Branched oligonucleotide-intercalator conjugate forming a parallel stranded structure inhibits HIV-1 integrase.

Integration of a DNA copy of the HIV-1 genome into chromosomal DNA of infected cells is a key step of viral replication. Integration is carried out by integrase, a viral protein which binds to both ends of viral DNA and catalyses reactions of the 3'-end processing and strand transfer. A 3'-3' branched oligonucleotide functionalised by the intercalator oxazolopyridocarbazole at each 5'-end was found to inhibit integration in vitro. We show that both a specific (G,A) sequence and the OPC intercalating agent contribute to the capability of the branched oligonucleotide to form a parallel stranded structure responsible for the inhibition.     

The Yersinia adhesin YadA collagen-binding domain structure is a novel left-handed parallel beta-roll

The crystal structure of the recombinant collagen-binding domain of Yersinia adhesin YadA from Yersinia enterocolitica serotype O:3 was solved at 1.55 A resolution. The trimeric structure is composed of head and neck regions, and the collagen binding head region is a novel nine-coiled left-handed parallel beta-roll. Before the beta-roll, the polypeptide loops from one monomer to the rest, and after the beta-roll the neck region does the same, making the transition from the globular head region to the narrower stalk domain. This creates an intrinsically stable 'lock nut' structure. The trimeric form of YadA is required for collagen binding, and mutagenesis of its surface residues allowed identification of a putative collagen-binding surface. Furthermore, a new structure-sequence motif for YadA beta-roll was used to identify putative YadA-head-like domains in a variety of human and plant pathogens. Such domains may therefore be a common bacterial strategy for avoiding host response.     

Crystal structure of pullulanase: evidence for parallel binding of oligosaccharides in the active site

The crystal structures of Klebsiella pneumoniae pullulanase and its complex with glucose (G1), maltose (G2), isomaltose (isoG2), maltotriose (G3), or maltotetraose (G4), have been refined at around 1.7-1.9A resolution by using a synchrotron radiation source at SPring-8. The refined models contained 920-1052 amino acid residues, 942-1212 water molecules, four or five calcium ions, and the bound sugar moieties. The enzyme is composed of five domains (N1, N2, N3, A, and C). The N1 domain was clearly visible only in the structure of the complex with G3 or G4. The N1 and N2 domains are characteristic of pullulanase, while the N3, A, and C domains have weak similarity with those of Pseudomonas isoamylase. The N1 domain was found to be a new type of carbohydrate-binding domain with one calcium site (CBM41). One G1 bound at subsite -2, while two G2 bound at -1 approximately -2 and +2 approximately +1, two G3, -1 approximately -3 and +2 approximately 0', and two G4, -1 approximately -4 and +2 approximately -1'. The two bound G3 and G4 molecules in the active cleft are almost parallel and interact with each other. The subsites -1 approximately -4 and +1 approximately +2, including catalytic residues Glu706 and Asp677, are conserved between pullulanase and alpha-amylase, indicating that pullulanase strongly recognizes branched point and branched sugar residues, while subsites 0' and -1', which recognize the non-reducing end of main-chain alpha-1,4 glucan, are specific to pullulanase and isoamylase. The comparison suggested that the conformational difference around the active cleft, together with the domain organization, determines the different substrate specificities between pullulanase and isoamylase.     

Approximate Bayesian Computation

Approximate Bayesian computation (ABC) constitutes a class of computational methods rooted in Bayesian statistics. In all model-based statistical inference, the likelihood function is of central importance, since it expresses the probability of the observed data under a particular statistical model, and thus quantifies the support data lend to particular values of parameters and to choices among different models. For simple models, an analytical formula for the likelihood function can typically be derived. However, for more complex models, an analytical formula might be elusive or the likelihood function might be computationally very costly to evaluate. ABC methods bypass the evaluation of the likelihood function. In this way, ABC methods widen the realm of models for which statistical inference can be considered. ABC methods are mathematically well-founded, but they inevitably make assumptions and approximations whose impact needs to be carefully assessed. Furthermore, the wider application domain of ABC exacerbates the challenges of parameter estimation and model selection. ABC has rapidly gained popularity over the last years and in particular for the analysis of complex problems arising in biological sciences (e.g., in population genetics, ecology, epidemiology, and systems biology).

This is a “Topic Page” article for PLOS Computational Biology.

     

Sequence and structure analysis of parallel beta helices: implication for constructing amyloid structural models

Increasing evidence suggests that amyloids and parallel beta helices may share similar motifs. A systemic analysis of beta helices is performed to examine their sequence and structural characteristics. Ile prefers to occur in beta strands. In contrast, Pro is disfavored, compatible with the underlying assumption in Pro-scanning mutagenesis. Cys, Asn, and Phe form significant homostacking (identical amino acid interactions). Asn is highly conserved in the high-energy, left-handed alpha-helical conformation, where it frequently forms amide stacking. Based on the observed prominent stacking of chemically similar residues in parallel beta helices, we propose that within the "cross-beta" framework, amyloids with longer peptide chains may have common structural features of in-register, parallel alignment, with the side chains forming identical amino acid ladders. The requirement of ladder formation limits the combinations of side chain interactions. Such a limit combined with environmental conditions (e.g., pH, concentration) could be a major reason for the ability of most polypeptides to form amyloids.     

Computation of 3D structure and distribution of helical parameters for D-period segments of human fibrillar collagen III

The structures of D-period segments of collagen (234 amino residues or ～1/4 of whole length) are established by methods of molecular mechanics and geometry analysis. Each D-period segment proves to have a unique spatial structure. The distributions of local helical parameters along the molecule are calculated. It is found that a second hydrogen bond is formed in every case when the second residue in the tripeptide G-X-Y is an amino acid. With such a combined H-bond network, all the peptide CO groups of glycines and of the third residues in tripeptides have quasi-equivalent positions on the surface of the collagen molecule. The local deformations of the polyproline II helix in the triple complex give rise to the observed modulation of structure at the macromolecular level, which may be important for the mutual orientation of collagen molecules during fibrillogenesis.     

Helix packing of leucine-rich peptides: a parallel leucine ladder in the structure of Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe.

The packing of peptide helices in crystals of the leucine-rich decapeptide Boc-Aib-Leu-Aib-Aib-Leu-Leu-Leu-Aib-Leu-Aib-OMe provides an example of ladder-like leucylleucyl interactions between neighboring molecules. The peptide molecule forms a helix with five 5----1 hydrogen bonds and two 4----1 hydrogen bonds near the C terminus. Three head-to-tail NH ... O = C hydrogen bonds between helices form continuous columns of helices in the crystal. The helicial columns associate in an antiparallel fashion, except for the association of Leu ... Leu side chains, which occurs along the diagonal of the cell where the peptide helices are parallel. The peptide, with formula C56H102N10O13, crystallizes in space group P2(1)2(1)2(1) with Z = 4 and cell parameters a = 16.774(3) A, b = 20.032(3) A and c = 20.117(3) A; overall agreement factor R = 10.7% for 2014 data with magnitude of F(obs) greater than 3 sigma (F); resolution 1.0 A.