Further developments of protein secondary structure prediction using information theory. New parameters and consideration of residue pairs

We have re-evaluated the information used in the Garnier-Osguthorpe-Robson (GOR) method of secondary structure prediction with the currently available database. The framework of information theory provides a means to formulate the influence of local sequence upon the conformation of a given residue, in a rigorous manner. However, the existing database does not allow the evaluation of parameters required for an exact treatment of the problem. The validity of the approximations drawn from the theory is examined. It is shown that the first-level approximation, involving single-residue parameters, is only marginally improved by an increase in the database. The second-level approximation, involving pairs of residues, provides a better model. However, in this case the database is not big enough and this method might lead to parameters with deficiencies. Attention is therefore given to overcoming this lack of data. We have determined the significant pairs and the number of dummy observations necessary to obtain the best result for the prediction. This new version of the GOR method increases the accuracy of prediction by 7%, bringing the amount of residues correctly predicted to 63% for three states and 68 proteins, each protein to be predicted being removed from the database and the parameters derived from the other proteins. If the protein to be predicted is kept in the database the accuracy goes up to 69.7%.

An algorithm for secondary structure determination in proteins based on sequence similarity

A secondary structure prediction algorithm is proposed on the hypothesis that short homologous sequences of amino acids have the same secondary structure tendencies. Comparisons are made with the secondary structure assignments of Kabsch and Sander from X-ray data [(1983) Biopolymers 22, 2577-2637] and an empirically determined similarity matrix which assigns a sequence similarity score between any two sequences of 7 residues in length. This similarity matrix differs in many respects from that of the Dayhoff substitution matrix [(1978) in: Atlas of Protein Sequence and Structure, (Dayhoff, M.O. ed). vol. 5. suppl. 3, pp. 353-358, National Biochemical Research Foundation, Washington, DC]. This homologue method had a prediction accuracy of 62.2% over 3 states for 61 proteins and 63.6% for a new set of 7 proteins not in the original data base.

Secondary structure prediction: combination of three different methods

A combination of three complementary secondary structure prediction methods is presented. The methods used are the GOR III method, the Homologue method and a new method, the bit pattern method, which is based on hydrophilic/hydrophobic residue patterns. For this purpose a hydropathy scale was developed and is presented here. The combination algorithm (Combine method) was designed to take the best results of each method and use their differences in order to improve the prediction. The combination yields 65.5% correctly predicted residues in three states: alpha-helix (H), beta-strand (E) and aperiodic structure (C) which is an improvement ranging from 2.5 to 6.5% compared with the individual methods when tested with a 67-polypeptide chain database. Seventy-five per cent of the regular secondary structure (H and E) runs are correctly located and beta-sheet runs are much better located by the Combine method in comparison to the other methods.

Computers and viral diseases. Preliminary bioinformatics studies on the design of a synthetic vaccine and a preventative peptidomimetic antagonist against the SARS-CoV-2 (2019-nCoV, COVID-19) coronavirus

This paper concerns study of the genome of the Wuhan Seafood Market isolate believed to represent the causative agent of the disease COVID-19. This is to find a short section or sections of viral protein sequence suitable for preliminary design proposal for a peptide synthetic vaccine and a peptidomimetic therapeutic, and to explore some design possibilities. The project was originally directed towards a use case for the Q-UEL language and its implementation in a knowledge management and automated inference system for medicine called the BioIngine, but focus here remains mostly on the virus itself. However, using Q-UEL systems to access relevant and emerging literature, and to interact with standard publically available bioinformatics tools on the Internet, did help quickly identify sequences of amino acids that are well conserved across many coronaviruses including 2019-nCoV. KRSFIEDLLFNKV was found to be particularly well conserved in this study and corresponds to the region around one of the known cleavage sites of the SARS virus that are believed to be required for virus activation for cell entry. This sequence motif and surrounding variations formed the basis for proposing a specific synthetic vaccine epitope and peptidomimetic agent. The work can, nonetheless, be described in traditional bioinformatics terms, and readily reproduced by others, albeit with the caveat that new data and research into 2019-nCoV is emerging and evolving at an explosive pace. Preliminary studies using molecular modeling and docking, and in that context the potential value of certain known herbal extracts, are also described.

•

Bioinformatics studies are carried out on the COVID-19 virus.

•

A sequence motif KRSFIEDLLFNKV is of particular interest.

•

Based on the above, synthetic peptides are designed.

•

Preliminary considerations are also given to non-peptide organic molecules.

Markers of risk of asthma death or readmission in the 12 months following a hospital admission for asthma

A case-control study has previously been reported of asthma deaths in people aged 5-45 years who had a hospital admission for asthma (the index admission) in New Zealand during 1981-1987. The study has been re-analysed to examine the association between markers of asthma severity and risk of asthma death or hospital admission; patients prescribed fenoterol were excluded from this re-analysis because of the previously reported interaction between fenoterol, asthma severity, and asthma deaths. The re-analysis included 39 patients who died of asthma during the 12 months after their index admission, 226 patients who had a readmission for asthma during the 12 months after their index admission, and 263 controls chosen from all index admissions. An admission in the previous 12 months was the strongest marker of subsequent risk of death (odds ratio (OR) = 3.5, 95% confidence interval (CI): 1.8-6.9, P less than 0.01), and was also a strong marker of subsequent risk of readmission (OR = 3.0, 95% CI: 2.1-4.2, P less than 0.01); the risk increased with the number of previous admissions. Three or more categories of prescribed asthma drugs was also associated with subsequent death (OR = 1.7, 95% CI: 0.9-3.3, P = 0.13) or readmission (OR = 1.9, 95% CI: 1.3-2.7, P less than 0.01); prescribed oral corticosteroids was only weakly associated with subsequent death (OR = 1.3, 95% CI: 0.6-2.8, P = 0.59), but was more strongly associated with subsequent readmission (OR = 1.9, 95% CI: 1.2-2.8, P less than 0.01).(ABSTRACT TRUNCATED AT 250 WORDS)

The mechanism of folding of globular proteins. Equilibria and kinetics of conformational transitions of penicillinase from Staphylococcus aureus involving a state of intermediate conformation

1. The thermodynamically reversible unfolding and refolding of penicillinase between the native and fully unfolded states were followed by using guanidinium chloride as denaturant. 2. The equilibria, studied by optical rotation, u.v. absorption, viscosity and enzyme activity, show the presence of a state of intermediate conformation, termed state H, which is stable at 20 degrees C in 0.8 M-guanidinium chloride. 3. The physical properties of this state show that it is slightly expanded with an intrinsic viscosity of 8 ml-g-1, that the 13 tyrosine residues, which are distributed through the primary sequence, are maximally exposed to the solvent and that the helix content is the same as that of the native state. 4. The kinetics of the transition between the native state, state H and the fully unfolded state were followed by u.v. absorption and by optical rotation. They are interpreted as showing that state H lies on the folding pathway between the native and fully unfolded states. 5. The transition between the native state and state H exhibits monophasic unfolding kinetics and biphasic refolding kinetics. This indicates that there must be at least two intermediate states in this process, at least one of which lies on the folding pathway which may also involve cul-de-sac paths. 6. The results are discussed in terms of a mechanism involving rapid stabilization of nucleation regions in a moderately compact but internally solvated structure, with 'native format' [Anfinsen (1973) Science 181, 233-230] secondary structure stabilized by tertiary interaction. The final and rate-limiting step in refolding involves shuffling of these structural elements into the native state. 7. This model is discussed in relation to folding in vivo.

COVID-19 Coronavirus spike protein analysis for synthetic vaccines, a peptidomimetic antagonist, and therapeutic drugs, and analysis of a proposed achilles' heel conserved region to minimize probability of escape mutations and drug resistance

This paper continues a recent study of the spike protein sequence of the COVID-19 virus (SARS-CoV-2). It is also in part an introductory review to relevant computational techniques for tackling viral threats, using COVID-19 as an example. Q-UEL tools for facilitating access to knowledge and bioinformatics tools were again used for efficiency, but the focus in this paper is even more on the virus. Subsequence KRSFIEDLLFNKV of the S2′ spike glycoprotein proteolytic cleavage site continues to appear important. Here it is shown to be recognizable in the common cold coronaviruses, avian coronaviruses and possibly as traces in the nidoviruses of reptiles and fish. Its function or functions thus seem important to the coronaviruses. It might represent SARS-CoV-2 Achilles’ heel, less likely to acquire resistance by mutation, as has happened in some early SARS vaccine studies discussed in the previous paper. Preliminary conformational analysis of the receptor (ACE2) binding site of the spike protein is carried out suggesting that while it is somewhat conserved, it appears to be more variable than KRSFIEDLLFNKV. However compounds like emodin that inhibit SARS entry, apparently by binding ACE2, might also have functions at several different human protein binding sites. The enzyme 11β-hydroxysteroid dehydrogenase type 1 is again argued to be a convenient model pharmacophore perhaps representing an ensemble of targets, and it is noted that it occurs both in lung and alimentary tract. Perhaps it benefits the virus to block an inflammatory response by inhibiting the dehydrogenase, but a fairly complex web involves several possible targets.

•

This paper “drills down” into the studies of the author's previous COVID-19 paper.

•

Designing vaccine and drugs must seek to avoid escape mutations.

•

Subsequence KRSFIEDLLFNKV seems recognizable across many coronaviruses.

•

The ACE2 binding domain is a target, but shows variation.

•

A steroid dehydrogenase is argued to remain an interesting model pharmacophore.

PRO-LIGAND: an approach to de novo molecular design. 1. Application to the design of organic molecules

An approach to de novo molecular design, PRO-LIGAND, has been developed that, in the environment of a large, integrated molecular design and simulation system, provides a unified framework for the generation of novel molecules which are either similar or complementary to a specified target. The approach is based on a methodology that has proved to be effective in other studies--placing molecular fragments upon target interaction sites-but incorporates many novel features such as the use of a rapid graph-theoretical algorithm for fragment placing, a generalised driver for structure generation which offers a large variety of fragment assembly strategies to the user and the pre-screening of library fragments. After a detailed description of the relevant modules of the package, PRO-LIGAND's efficacy in aiding rational drug design is demonstrated by its ability to design mimics of methotrexate and potential inhibitors for dihydrofolate reductase and HIV-1 protease.

Analysis of code relating sequences to conformation in globular prtoeins. Theory and application of expected information

1. An information theory analysis of the folding of a globular protein is proposed. 2. The folding is seen as a transfer of information between two messages, the primary sequence and the biologically active conformation. 3. It is shown how the information transferred was estimated by inspection of proteins of known primary sequence and conformation. 4. In this estimation, concerted use of subjective (Bayesian) probabilities leads to a more robust approach which can be employed whether the number of proteins of known sequence and conformation is large or small. 5. Further, it is demonstrated that the problem then becomes a very simple algebraic formulation for information estimates. 6. Finally, it is shown how this process of information theory analysis can be reversed to predict the conformation of a protein by using its primary sequence and the above information estimates obtained from other proteins. 7. The present paper provides the theoretical basis for the derivation and application of a stereochemical alphabet (Robson & Pain, 1974a,c), and for an investigation of the effects of residues on the conformations of their neighbours (Robson & Pain, 1974b).

Community-based asthma care: trial of a "credit card" asthma self-management plan

Although asthma self-management plans are widely recommended as essential in the long-term treatment of adult asthma, there have been few studies examining their use. Our objective was to assess the effect of a "credit card" adult asthma self-management plan in a community experiencing major health problems from asthma, by means of a before and after intervention trial of the efficacy of the "credit card" plan, when introduced through community-based asthma clinics. The participants were 69 Maori people with asthma. The "credit card" plan consisted of written guidelines for the self-management of asthma, based on self-assessment of asthma severity, printed on a plastic card. On one side, management guidelines were based on the interpretation of peak expiratory flow rate (PEFR) recordings, whilst the reverse side was based on symptoms. The outcome measures used were before and after comparison of markers of asthma morbidity and requirement for acute medical treatment; and a structured questionnaire assessing the acceptability and use of the credit card plan. Following the introduction of the plan, the mean PEFR increased from 347 to 389 l.min-1, the percentage of nights woken fell from 30.4 to 16.9%, and the number of days "out of action" fell from 3.8 to 1.7%. The requirements for acute medical treatment also fell during the intervention period. Most participants commented favourably on the content and usefulness of the plan. In the situation of worsening asthma, 28% of subjects found the peak flow side of the card most helpful, 7% the symptoms side, and 48% found both sides equally helpful.(ABSTRACT TRUNCATED AT 250 WORDS)

Fine sediment and particulate organic matter: A review and case study on ridge-to-reef transport, transformations, fates, and impacts on marine ecosystems

Studies documenting the effects of land-derived suspended particulate matter (SPM, i.e., particulate organic matter and mineral sediment) on marine ecosystems are typically disconnected from terrestrial studies that determine their origin, transport and fate. This study reviews sources, transport, transformations, fate and effects of SPM along the 'ridge-to-reef' continuum. We show that some of the SPM can be transported over long distances and transformed into large and easily resuspendible organic-rich sediment flocs. These flocs may lead to prolonged reductions in water clarity, impacting upon coral reef, seagrass and fish communities. Using the Great Barrier Reef (NE Australia) as a case study, we identify the latest research tools to determine thresholds of SPM exposure, allowing for an improved appreciation of marine risk. These tools are used to determine ecologically-relevant end-of-basin load targets and reliable marine water quality guidelines, thereby enabling enhanced prioritisation and management of SPM export from ridge-to-reef.

Bioinformatics studies on a function of the SARS-CoV-2 spike glycoprotein as the binding of host sialic acid glycans

SARS-CoV and SARS-CoV-2 do not appear to have functions of a hemagglutinin and neuraminidase. This is a mystery, because sugar binding activities appear essential to many other viruses including influenza and even most other coronaviruses in order to bind to and escape from the glycans (sugars, oligosaccharides or polysaccharides) characteristic of cell surfaces and saliva and mucin. The S1 N terminal Domains (S1-NTD) of the spike protein, largely responsible for the bulk of the characteristic knobs at the end of the spikes of SARS-CoV and SARS-CoV-2, are here predicted to be “hiding” sites for recognizing and binding glycans containing sialic acid. This may be important for infection and the ability of the virus to locate ACE2 as its known main host cell surface receptor, and if so it becomes a pharmaceutical target. It might even open up the possibility of an alternative receptor to ACE2. The prediction method developed, which uses amino acid residue sequence alone to predict domains or proteins that bind to sialic acids, is naïve, and will be advanced in future work. Nonetheless, it was surprising that such a very simple approach was so useful, and it can easily be reproduced in a very few lines of computer program to help make quick comparisons between SARS-CoV-2 sequences and to consider the effects of viral mutations.

•

This paper extends the studies of the author's previous SARS-CoV-2 papers.

•

Designing vaccine and drugs must seek to avoid escape mutations.

•

Strangely, SARS-CoV and SARS-CoV-2 appear to lack sialic acid binding functions.

•

Sequence motifs are found, but they require a simple prediction method.

Monte Carlo simulation of water behavior around the dipeptide N-acetylalanyl-N-methylamide

Applications of Monte Carlo and molecular dynamics computer simulation techniques indicate that they are potentially powerful tools for understanding biological systems at the molecular level. The Monte Carlo technique can be used to study the solvent structure around a small peptide and the effect of the aqueous environment on the conformational equilibria of the peptide.

The mechanism of folding of globular proteins. Suitability of a penicillinase from Staphylococcus Aureus as a model for refolding studies

1. A homogeneous preparation of penicillinase (penicillin amido-beta-lactamhydrolase, EC 3.5.2.6) was isolated and purified from cultures of Staphylococcus aureus by a simple two-stage procedure. 2. The native protein contains 20-30% helix as determined by optical-rotatory-dispersion and circular-dichroism measurements. Some 54(+/-5)% of the 13 tyrosine residues are exposed to solvent molecules of diameter 0.44 and 0.94 nm. 3. Conditions that allow full recovery of enzymic activity and native conformation from the fully unfolded state in 4M-guanidinium chloride were defined. 4. Refolding of the protein was shown to be inhibited by intermolecular interaction, by small changes in ionization and by low concentrations (0.025 M) of phenol.

Refined models for computer calculations in protein engineering. Calibration and testing of atomic potential functions compatible with more efficient calculations

A reappraisal has been made of interatomic potential functions for protein structure calculations using the all-atom approximation (except CH, CH2 and CH3, which are treated as "united atoms"). Some key problems are identified and treated. The potential functions are somewhat novel in form and consistent with more efficient and robust folding algorithms. In addition, the potentials are calibrated for for the rigid geometry approximation, since use of fixed standard bond lengths and valence angles (and fixed trans planar peptide groups) reduces the number of conformational variables and saves a great deal of computer time. Though these algorithms demand the use of potential functions of this special type, these functions can be readily implemented in more classical programs for the conformational analysis of proteins. They are calibrated or tested against a large body of experimental data, including extended basis set ab initio, quantum mechanical calculations, nuclear magnetic resonance spectroscopic data and dipole moment data for di- and oligopeptides, characteristic ratio data for random coil homopolypeptides, extensive data from peptide solubility studies, and experimental structures of polyalanine fibres and globular proteins. This paper will form the basis of a further report, which will include investigations of how water might be more realistically represented subject to the computing power available.

Population genetic structure of a colonising, triploid weed, Hieracium lepidulum

Understanding the breeding system and population genetic structure of invasive weed species is important for biocontrol, and contributes to our understanding of the evolutionary processes associated with invasions. Hieracium lepidulum is an invasive weed in New Zealand, colonising a diverse range of habitats including native Nothofagus forest, pine plantations, scrubland and tussock grassland. It is competing with native subalpine and alpine grassland and herbfield vegetation. H. lepidulum is a triploid, diplosporous apomict, so theoretically all seed is clonal, and there is limited potential for the creation of variation through recombination. We used intersimple sequence repeats (ISSRs) to determine the population genetic structure of New Zealand populations of H. lepidulum. ISSR analysis of five populations from two regions in the South Island demonstrated high intrapopulation genotypic diversity, and high interpopulation genetic structuring; PhiST = 0.54 over all five populations. No private alleles were found in any of the five populations, and allelic differentiation was correlated to geographic distance. Cladistic compatibility analysis indicated that both recombination and mutation were important in the creation of genotypic diversity. Our data will contribute to any biocontrol program developed for H. lepidulum. It will also be a baseline data set for future comparisons of genetic structure during the course of H. lepidulum invasions.

PRO-LIGAND: an approach to de novo molecular design. 3. A genetic algorithm for structure refinement

Recently, the development of computer programs which permit the de novo design of molecular structures satisfying a set of steric and chemical constraints has become a burgeoning area of research and many operational systems have been reported in the literature. Experience with PRO-LIGAND-the de novo design methodology embodied in our in-house molecular design and simulation system PRO-METHEUS-has suggested that the addition of a genetic algorithm (GA) structure refinement procedure can 'add value' to an already useful tool. Starting with the set of designed molecules as an initial population, the GA can combine features from both high- and low-scoring structures and, over a number of generations, produce individuals of better score than any of the starting structures. This paper describes how we have implemented such a procedure and demonstrates its efficacy in improving two sets of molecules generated by different de novo design projects.

Biochemical effects of the hypoglycaemic compound pent--4-enoic acid and related non-hypoglycaemic fatty acids

1. The effects of the hypoglycaemic compound, pent-4-enoic acid, and of four structurally related non-hypoglycaemic compounds (pentanoic acid, pent-2-enoic acid, cyclopropanecarboxylic acid and cyclobutanecarboxylic acid), on the oxidation of saturated fatty acids by rat liver mitochondria were determined. 2. The formation of (14)CO(2) from [1-(14)C]palmitate was strongly inhibited by 0.01mm-pent-4-enoic acid. 3. The inhibition of oxygen uptake was less than that of (14)CO(2) formation, presumably because fumarate was used as a sparker. 4. The oxidation of [1-(14)C]-butyrate, -octanoate or -laurate was not strongly inhibited by 0.01mm-pent-4-enoic acid. 5. The other four non-hypoglycaemic compounds did not inhibit the oxidation of any saturated fatty acid when tested at 0.01mm concentration, though they all inhibited strongly at 10mm. 6. The oxidation of [1-(14)C]-myristate and -stearate, but not of [1-(14)C]decanoate, was strongly inhibited by 0.01mm-pent-4-enoic acid. 7. The oxidation of [1-(14)C]palmitate was about 50% carnitine-dependent under the experimental conditions used. 8. The percentage inhibition of [1-(14)C]palmitate oxidation by pent-4-enoic acid was the same whether carnitine was present or not. 9. Acetoacetate formation from saturated fatty acids was inhibited by 0.1mm-cyclopropanecarboxylic acid to a greater extent than their oxidation. 10. The other compounds tested inhibited acetoacetate formation from saturated fatty acids proportionately to the inhibition of oxidation. 11. Possible mechanisms for the inhibition of long-chain fatty acid oxidation by pent-4-enoic acid are discussed. 12. There was a correlation between the ability to inhibit long-chain fatty acid oxidation and hypoglycaemic activity in this series of compounds.