Why twenty amino acid residue types suffice(d) to support all living systems.

It is well known that proteins are built up from an alphabet of 20 different amino acid types. These suffice to enable the protein to fold into its operative form relevant to its required functional roles. For carrying out these allotted functions, there may in some cases be a need for post-translational modifications and it has been established that an additional three types of amino acid have at some point been recruited into this process.

A tensegrity model for hydrogen bond networks in proteins.

Hydrogen-bonding networks in proteins considered as structural tensile elements are in balance separately from any other stabilising interactions that may be in operation. The hydrogen bond arrangement in the network is reminiscent of tensegrity structures in architecture and sculpture. Tensegrity has been discussed before in cells and tissues and in proteins.

The fate of proteins in outer space.

It is well established that any properly conducted biophysical studies of proteins must take appropriate account of solvent. For water-soluble proteins it has been an article of faith that water is largely responsible for stabilizing the fold, a notion that has recently come under increasing scrutiny. Further, there are some instances when proteins are studied experimentally in the absence of solvent, as in matrix-assisted laser desorption/ionization or electrospray mass spectrometry, for example, or in organic solvents for protein engineering purposes.

Melody discrimination and protein fold classification.

One of the greatest challenges in theoretical biophysics and bioinformatics is the identification of protein folds from sequence data. This can be regarded as a pattern recognition problem. In this paper we report the use of a melody generation software where the inputs are derived from calculations of evolutionary information, secondary structure, flexibility, hydropathy and solvent accessibility from multiple sequence alignment data.

Comparison of Algorithms for Prediction of Protein Structural Features from Evolutionary Data.

Proteins have many functions and predicting these is still one of the major challenges in theoretical biophysics and bioinformatics. Foremost amongst these functions is the need to fold correctly thereby allowing the other genetically dictated tasks that the protein has to carry out to proceed efficiently. In this work, some earlier algorithms for predicting protein domain folds are revisited and they are compared with more recently developed methods.

The dipeptide conformations of all twenty amino acid types in the context of biosynthesis.

There have been many studies of dipeptide structure at a high level of accuracy using quantum chemical methods. Such calculations are resource-consuming (in terms of memory, CPU and other computational imperatives) which is the reason why most previous studies were restricted to the two simplest amino-acid residue types, glycine and alanine. We improve on this by extending the scope of residue types to include all 20 naturally occurring residue types.

Prediction of protein structural features from sequence data based on Shannon entropy and Kolmogorov complexity.

While the genome for a given organism stores the information necessary for the organism to function and flourish it is the proteins that are encoded by the genome that perhaps more than anything else characterize the phenotype for that organism. It is therefore not surprising that one of the many approaches to understanding and predicting protein folding and properties has come from genomics and more specifically from multiple sequence alignments. In this work I explore ways in which data derived from sequence alignment data can be used to investigate in a predictive way three different aspects of protein structure: secondary structures, inter-residue contacts and the dynamics of switching between different states of the protein.

The preferred conformation of dipeptides in the context of biosynthesis.

Globular proteins are folded polypeptide structures comprising stretches of secondary structures (helical (α- or 310 helix type), polyproline helix or β-strands) interspersed by regions of less well-ordered structure ("random coil"). Protein fold prediction is a very active field impacting inte alia on protein engineering and misfolding studies. Apart from the many studies of protein refolding from the denatured state, there has been considerable interest in studying the initial formation of peptides during biosynthesis, when there are at the outset only a few residues in the emerging polypeptide.

Current progress in Structure-Based Rational Drug Design marks a new mindset in drug discovery.

The past decade has witnessed a paradigm shift in preclinical drug discovery with structure-based drug design (SBDD) making a comeback while high-throughput screening (HTS) methods have continued to generate disappointing results. There is a deficit of information between identified hits and the many criteria that must be fulfilled in parallel to convert them into preclinical candidates that have a real chance to become a drug. This gap can be bridged by investigating the interactions between the ligands and their receptors.

Protein folding: a problem with multiple solutions.

There is continued interest in predicting the structure of proteins either at the simplest level of identifying their fold class or persevering all the way to an atomic resolution structure. Protein folding methods have become very sophisticated and many successes have been recorded with claims to have solved the native structure of the protein. But for any given protein, there may be more than one solution.

Accelerated simulation of unfolding and refolding of a large single chain globular protein.

We have developed novel strategies for contracting simulation times in protein dynamics that enable us to study a complex protein with molecular weight in excess of 34 kDa. Starting from a crystal structure, we produce unfolded and then refolded states for the protein. We then compare these quantitatively using both established and new metrics for protein structure and quality checking.

Inactivation and reactivation of ribonuclease A studied by computer simulation.

The year 2011 marked the half-centenary of the publication of what came to be known as the Anfinsen postulate, that the tertiary structure of a folded protein is prescribed fully by the sequence of its constituent amino acid residues. This postulate has become established as a credo, and, indeed, no contradictions seem to have been found to date. However, the experiments that led to this postulate were conducted on only a single protein, bovine ribonuclease A (RNAse).

G-protein-coupled receptor dynamics: dimerization and activation models compared with experiment.

Our previously derived models of the active state of the β2-adrenergic receptor are compared with recently published X-ray crystallographic structures of activated GPCRs (G-protein-coupled receptors). These molecular dynamics-based models using experimental data derived from biophysical experiments on activation were used to restrain the receptor to an active state that gave high enrichment for agonists in virtual screening. The β2-adrenergic receptor active model and X-ray structures are in good agreement over both the transmembrane region and the orthosteric binding site, although in some regions the active model is more similar to the active rhodopsin X-ray structures.

On dating stages in prebiotic chemical evolution.

The notion that RNA must have had a unique and decisive role in the development of life needs hardly be questioned. However, the chemical complexity and other properties of RNA, such as high solubility in water and vulnerability to degradation, make it improbable that RNA could have had an early presence in the development of life on Earth or on any comparable telluric planet. Rather, the task of origin of life research must surely be to identify those chemical processes which could have taken place on Earth that could accumulate the complexity and rich molecular information content needed to sustain primitive life, and ultimately give rise to RNA.

Membrane-spanning peptides and the origin of life.

An explanation is given as to why membrane-spanning peptides must have been the first "information-rich" molecules in the development of life. These peptides are stabilised in a lipid bilayer membrane environment and they are preferentially made from the simplest, and likewise oldest, of the amino acids that survive today. Transmembrane peptides can exercise functions that are essential for biological systems such as signal transduction and material transport across membranes.

Residues from transmembrane helices 3 and 5 participate in leukotriene B4 binding to BLT1.

Leukotrienes are inflammatory mediators that bind to seven transmembrane, G-protein-coupled receptors (GPCRs). Here we examine residues from transmembrane helices 3 and 5 of the leukotriene B4 (LTB4) receptor BLT1 to elucidate how these residues are involved in ligand binding. We have selected these residues on the basis of (1) amino acid sequence analysis, (2) receptor binding and activation studies with a variety of leukotriene-like ligands and recombinant BLT1 receptors, (3) previously published recombinant BLT1 mutants, and (4) a computed model of the active structure of the BLT1 receptor.

Did life begin on the beach?

Water is one of the prerequisites of life. Further requirements are the existence of a system of interacting organic molecules capable of capturing and converting the supply of external energy and elaborating the replicating function that is needed for propagation. None of this would be possible without the existence of some means of concentrating, selecting, and then containing these mutually interacting substances in proximity to one another, i.

Location and nature of the residues important for ligand recognition in G-protein coupled receptors.

The overall structure of the biogenic amine subclass of the G-protein-coupled receptors, and of their ligand binding sites, is discussed with the aim of highlighting the major structural features of these receptors that are responsible for ligand recognition. A comparison is made between biogenic amine receptors, peptide receptors of the rhodopsin class, and the secretin receptors which all have peptide ligands. The question of where the peptide ligands bind, whether at extracellular sites or within the transmembrane helix bundle, is discussed.

Toward the active conformations of rhodopsin and the beta2-adrenergic receptor.

Using sets of experimental distance restraints, which characterize active or inactive receptor conformations, and the X-ray crystal structure of the inactive form of bovine rhodopsin as a starting point, we have constructed models of both the active and inactive forms of rhodopsin and the beta2-adrenergic G-protein coupled receptors (GPCRs). The distance restraints were obtained from published data for site-directed crosslinking, engineered zinc binding, site-directed spin-labeling, IR spectroscopy, and cysteine accessibility studies conducted on class A GPCRs. Molecular dynamics simulations in the presence of either "active" or "inactive" restraints were used to generate two distinguishable receptor models.

De novo generation of molecular structures using optimization to select graphs on a given lattice.

A recurrent problem in organic chemistry is the generation of new molecular structures that conform to some predetermined set of structural constraints that are imposed in an endeavor to build certain required properties into the newly generated structure. An example of this is the pharmacophore model, used in medicinal chemistry to guide de novo design or selection of suitable structures from compound databases. We propose here a method that efficiently links up a selected number of required atom positions while at the same time directing the emergent molecular skeleton to avoid forbidden positions.

Recognition of privileged structures by G-protein coupled receptors.

Privileged structures are ligand substructures that are widely used to generate high-affinity ligands for more than one type of receptor. To explain this, we surmised that there must be some common feature in the target proteins. For a set of class A GPCRs, we found a good correlation between conservation patterns of residues in the ligand binding pocket and the privileged structure fragments in class A GPCR ligands.

Ab initio modelling of the N-terminal domain of the secretin receptors.

G protein coupled receptors of the secretin family are activated by peptide hormones of about 30 residues in length. There is considerable sequence homology within both the hormone and receptor families. The receptors possess in addition to the integral membrane domain a characteristic extracellular domain of about 120 residues in length, having conserved cysteine residues, which are involved in disulphide bridge formation, and tryptophanes, which have been shown to be critical for hormone binding.

Essential motions in a fungal lipase with bound substrate, covalently attached inhibitor and product.

As an aid to understanding the influence of dynamic fluctuations during esterolytic catalysis, we follow protein flexibility at three different steps along the catalytic pathway from substrate binding to product clearance via a covalently attached inhibitor, which represents a transition-state mimic. We have applied a classical approach, using molecular dynamics simulations to monitor protein dynamics in the nanosecond regime. We filter out small amplitude fluctuations and focus on the anharmonic contributions to the overall dynamics.

Prediction of protein residue contacts with a PDB-derived likelihood matrix.

Proteins with similar folds often display common patterns of residue variability. A widely discussed question is how these patterns can be identified and deconvoluted to predict protein structure. In this respect, correlated mutation analysis (CMA) has shown considerable promise.