A method for α-helical integral membrane protein fold prediction

Integral membrane proteins (of the α-helical class) are of central importance in a wide variety of vital cellular functions. Despite considerable effort on methods to predict the location of the helices, little attention has been directed toward developing an automatic method to pack the helices together. In principle, the prediction of membrane proteins should be easier than the prediction of globular proteins: there is only one type of secondary structure and all helices pack with a common alignment across the membrane. This allows all possible structures to be represented on a simple lattice and exhaustively enumerated. Prediction success lies not in generating many possible folds but in recognizing which corresponds to the native. Our evaluation of each fold is based on how well the exposed surface predicted from a multiple sequence alignment fits its allocated position. Just as exposure to solvent in globular proteins can be predicted from sequence variation, so exposure to lipid can be recognized by variable-hydrophobic (variphobic) positions. Application to both bacteriorhodopsin and the eukaryotic rhodopsin/opsin families revealed that the angular size of the lipid-exposed faces must be predicted accurately to allow selection of the correct fold. With the inherent uncertainties in helix prediction and parameter choice, this accuracy could not be guaranteed but the correct fold was typically found in the top six candidates. Our method provides the first completely automatic method that can proceed from a scan of the protein sequence databanks to a predicted three-dimensional structure with no intervention required from the investigator. Within the limited domain of the seven helix bundle proteins, a good chance can be given of selecting the correct structure. However, the limited number of sequences available with a corresponding known structure makes further characterization of the method difficult. © 1994 John Wiley & Sons, Inc.     

Monte carlo simulations of protein folding. II. Application to protein A,ROP, and crambin

The hierarchy of lattice Monte Carlo models described in the accompanying paper (Kolinski, A., Skolnick, J. Monte Carlo simulations of protein folding. I. Lattice model and interaction scheme. Proteins 18:338–352, 1994) is applied to the simulation of protein folding and the prediction of 3-dimensional structure. Using sequence information alone, three proteins have been successfully folded: the B domain of staphylococcal protein A, a 120 residue, monomeric version of ROP dimer, and crambin. Starting from a random expanded conformation, the model proteins fold along relatively well-defined folding pathways. These involve a collection of early intermediates, which are followed by the final (and rate-determining) transition from compact intermediates closely resembling the molten globule state to the native-like state. The predicted structures are rather unique, with native-like packing of the side chains. The accuracy of the predicted native conformations is better than those obtained in previous folding simulations. The best (but by no means atypical) folds of protein A have a coordinate rms of 2.25 Å from the native Cα trace, and the best coordinate rms from crambin is 3.18 Å. For ROP monomer, the lowest coordinate rms from equivalent Cαs of ROP dimer is 3.65 Å. Thus, for two simple helical proteins and a small α/β protein, the ability to predict protein structure from sequence has been demonstrated. © 1994 John Wiley & Sons, Inc.     

Energy minimization method using automata network for sequence and side-chain conformation prediction from given backbone geometry

Globular proteins have high packing densities as a result of residue side chains in the core achieving a tight, complementary packing. The internal packing is considered the main determinant of native protein structure. From that point of view, we present here a method of energy minimization using an automata network to predict a set of amino acid sequences and their side-chain conformations from a desired backbone geometry for de novo design of proteins. Using discrete side-chain conformations, that is, rotamers, the sequence generation problem from a given backbone geometry becomes one of combinatorial problems. We focused on the residues composing the interior core region and predicted a set of amino acid Sequences and their side-chain conformations only from a given backbone geometry. The kinds of residues were restricted to six hydrophobic amino acids (Ala, Ile, Met, Leu, Phe, and Val) because the core regions are almost always composed of hydrophobic residues. The obtained sequences were well packed as was the native sequence. The method can be used for automated sequence generation in the de novo design of proteins. © 1994 Wiley-Liss, Inc.     

RFLP analysis for species separation in the generaBipolaris andCurvularia

Restriction fragment length polymorphism (RFLP) of the total DNA ofBipolaris andCurvularia species was analysed using arbitrarily chosen genomic clones of DNA fromCurvularia lunata andBipolaris maydis as probes. Clear differences among species in both genera, resulting in different banding positions, were obtained with some probe-enzyme combinations. Intraspecific polymorphism in banding positions with these probe-enzyme combinations was slight. These analyses allow discrimination between the species. DNA fingerprinting with intrageneric probes is a potentially useful tool for species separation and identification inBipolaris andCurvularia when coupled with another characteristic such as conidial morphology.Curvularia aeria comb. nov. was proposed forCurvularia lunata var.aeria on the basis of differences in RFLP banding patterns and differences in conidial morphology.     

Discrimination of folding types of globular proteins based on average distance maps constructed from their sequences

It has been shown that probable portions which form contacts in a protein can be predicted by means of an average distance map (ADM) as well as regular structures (-helices and -turns) defined as short-range compact regions (Kikuchiet al., 1988a,c). In this paper, we analyze the occurrence of those portions and short-range compact regions on ADMs for various proteins regarding their folding types. We have found out that each folding type of proteins shows characteristic distribution of such parts on ADMS. We also discuss the possibility of the prediction of folding types of proteins by ADMs.     

Prediction of near-term climate change impacts on UK wheat quality and the potential for adaptation through plant breeding

Wheat is a major crop worldwide, mainly cultivated for human consumption and animal feed. Grain quality is paramount in determining its value and downstream use. While we know that climate change threatens global crop yields, a better understanding of impacts on wheat end-use quality is also critical. Combining quantitative genetics with climate model outputs, we investigated UK-wide trends in genotypic adaptation for wheat quality traits. In our approach, we augmented genomic prediction models with environmental characterisation of field trials to predict trait values and climate effects in historical field trial data between 2001 and 2020. Addition of environmental covariates, such as temperature and rainfall, successfully enabled prediction of genotype by environment interactions (G × E), and increased prediction accuracy of most traits for new genotypes in new year cross validation. We then extended predictions from these models to much larger numbers of simulated environments using climate scenarios projected under Representative Concentration Pathways 8.5 for 2050–2069. We found geographically varying climate change impacts on wheat quality due to contrasting associations between specific weather covariables and quality traits across the UK. Notably, negative impacts on quality traits were predicted in the East of the UK due to increased summer temperatures while the climate in the North and South-west may become more favourable with increased summer temperatures. Furthermore, by projecting 167,040 simulated future genotype–environment combinations, we found only limited potential for breeding to exploit predictable G × E to mitigate year-to-year environmental variability for most traits except Hagberg falling number. This suggests low adaptability of current UK wheat germplasm across future UK climates. More generally, approaches demonstrated here will be critical to enable adaptation of global crops to near-term climate change.     

Genomic Analysis of Borrelia japonica Sp. Nov. Isolated from Ixodes ovatus in Japan

Genetic characteristics of 12 Borrelia isolates from the tick, Ixodes ovatus, I. persulcatus, and the rodent, Apodemus speciosus ainu, in Japan were compared to members of the three genospecies of Borrelia burgdorferi sensu lato; B. burgdorferi sensu stricto, B. garinii and group VS461. The methods used in this study were the quantitative microplate DNA hybridization assay and restriction fragment length polymorphism (RFLP) analyses of the flagellin structural genes and the 16S rRNA genes. The six isolates from I. persulcatus and A. speciosus ainu were identified as genospecies B. garinii using RFLP analysis of the flagellin and 16S rRNA genes. In contrast, RFLP analysis of the six isolates from I. ovatus indicated that they were different from the three reported genospecies. DNA homology studies confirmed the RFLP results. The six isolates from I. ovatus had DNA homologies ranging from 85 to 99%, whereas DNA relatedness of the I. ovatus isolate with strains belonging to the three genospecies was 50 to 64%. These results suggest that the strains isolated from I. ovatus in Japan differ from the three genospecies and should be classified as a new genospecies of B. burgdorferi sensu lato. We propose that strains isolated from I. ovatus should be classified as B. japonica sp. nov.     

Structures of randomly generated mutants of T4 lysozyme show that protein stability can be enhanced by relaxation of strain and by improved hydrogen bonding via bound solvent.

The structures of three mutants of bacteriophage T4 lysozyme selected using a screen designed to identify thermostable variants are described. Each of the mutants has a substitution involving threonine. Two of the variants, Thr 26-->Ser (T26S) and Thr 151-->Ser (T151S), have increased reversible melting temperatures with respect to the wild-type protein. The third, Ala 93-->Thr (A93T), has essentially the same stability as wild type. Thr 26 is in the wall of the active-site cleft. Its replacement with serine results in the rearrangement of nearby residues, most notably Tyr 18, suggesting that the increase in stability may result from the removal of strain. Thr 151 in the wild-type structure is far from the active site and appears to sterically prevent the access of solvent to a preformed binding site. In the mutant, the removal of the methyl group allows access to the solvent binding site and, in addition, the Ser 151 hydroxyl rotates to a new position so that it also contributes to solvent binding. Residue 93 is in a highly exposed site on the surface of the molecule, and presumably is equally solvent exposed in the unfolded protein. It is, therefore, not surprising that the substitution Ala 93-->Thr does not change stability. The mutant structures show how chemically similar mutations can have different effects on both the structure and stability of the protein, depending on the structural context. The results also illustrate the power of random mutagenesis in obtaining variants with a desired phenotype.(ABSTRACT TRUNCATED AT 250 WORDS)