Heavier-than-air flying machines are impossible.

Many G protein-coupled receptor (GPCR) models have been built over the years. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn more about the models we build today. We conclude that the GPCR modelling field is riddled with 'common knowledge' similar to Lord Kelvin's remark in 1895 that "heavier-than-air flying machines are impossible", and we summarize what we think are the (im)possibilities of modelling GPCRs using the coordinates of bovine rhodopsin as a template. Associated WWW pages: www.gpcr.org/articles/2003_mod

Correlated mutation analyses on very large sequence families.

The 'omics era' (the era of genomics, proteomics, and so forth) is marked by a flood of data that need to be interpreted to become useful information. Thanks to genome sequencing projects, large numbers of sequence families with more than a thousand members each are now available. Novel analytical techniques are needed to deal with this avalanche of sequence data. Sequence entropy is a measure of the information present in an alignment, whereas sequence variability represents the mutational flexibility at a particular position. Entropy versus variability plots can reveal the roles of groups of residues in the overall function of a protein. Such roles can be as part of the main active site, part of a modulator binding site, or transduction of a signal between those sites. Residues that are involved in a common function tend to stay conserved as a group, but when they mutate, they tend to mutate together. Correlated mutation analysis can detect groups of residue positions that show this behaviour. The combination of entropy, variability and correlation is a powerful tool to convert sequence data into useful information. This analysis can, for example, detect the key residues involved in cooperativity in globins, the switch regions in ras-like proteins and the calcium binding and signalling residues in serine proteases. We have extrapolated from these three classes of structurally and functionally well-described proteins to G-protein-coupled receptors (GPCRs). We can detect the residues in the main functional site in GPCRs that are responsible for G-protein coupling, the residues in the endogenous agonist binding site, and the residues in between that transduce the signal to and fro between these sites. The results are discussed in the light of a simple two-step evolutionary model for the development of functional proteins.

Full window stereo.

Visualisation is the bioinformaticist's most important tool for the study of macromolecules, and being able to see molecules in stereo is a crucial aspect. Stereo vision is based on the principle that each eye is presented with the best possible image of what it would have seen if the object was really there in 3D. The simplest approach to stereo vision is to display the right eye picture on the right half of the screen and the left eye picture on the left half while using a mirror system to ensure that each eye sees what it is supposed to see. More expensive workstations use hardware to alternately display the left and right eye pictures while synchronously blocking the transparency in the right or left lens of the special glasses worn by the user. We present here some simple software that uses inexpensive hardware, originally designed for the computer game industry, to make full screen stereo available on Linux-based PCs. The quality of the stereo vision is similar to the top-of-the-line graphics workstations that are capable of quad-buffering. This stereo option has been incorporated in the XII based version of WHAT IF (Vriend, G. J. Mol. Graphics 1990, 8, 52-56), but the stereo source code is freely available and can easily be incorporated in other visualization packages.

ACGT and vicilin core sequences in a promoter domain required for seed-specific expression of a 2S storage protein gene are recognized by the opaque-2 regulatory protein.

The expression of Brazil nut storage albumin genes is highly regulated during seed development. Several sequences in the promoter of one of these genes show homologies with the target sites of the maize O2 bZIP regulatory protein. We therefore asked whether the O2 protein would recognize these promoter sequences. We show that the O2 protein binds to three different sequences (F1, F2 and F3). F1 and F3 are hybrid C/G and A/G boxes, respectively, that are homologous to the O2-binding site of a maize alpha-zein gene. F2 is a new O2-binding sequence related to the O2 target sites of the Coix alpha-coxin, the maize b-32 genes and the AP-1 pseudopalindrome. Molecular modelling showed that an Asn and a Ser in the 02 DNA binding domain make different base-specific contacts with each operator. 5' Promoter deletions of the be2S1 gene showed that the domain containing the O2 target sites F1 and F2 is required for detectable reporter gene expression in transgenic tobacco seeds. Moreover, the homologous coix O2 protein was shown to in situ transactivate the promoter region encompassing the three O2-binding sites F1, F2 and F3. Thus, these sites may be in vivo regulatory sequences mediating activation by bZIP regulatory proteins.

Analysis of the black-eyed pea trypsin and chymotrypsin inhibitor-alpha-chymotrypsin complex.

The black-eyed pea trypsin and chymotrypsin inhibitor (BTCI) is a member of the Bowman-Birk protease inhibitor (BBI) family. The three-dimensional model of the BTCI-chymotrypsin complex was built based on the homology to Bowman-Birk inhibitors with known structures. An extensive theoretical and experimental study of these known structures has been performed. The model confirms the ideas about Bowman-Birk inhibitor structure-function relations and agrees well with our experimental data (circular dichroism, IR and light scattering). The electrostatic potentials at the enzyme-inhibitor contact surface reveal a pattern of complementary electrostatic potentials from which mutations can be inferred that could give these inhibitors an altered specificity.

The use of position-specific rotamers in model building by homology.

In this study we concentrate on replacing side chains as a subtask of model building by homology. Two problems arise. How to determine potential low energy rotamers? And how to avoid the combinatorial explosion that results from the combination of many residues for which multiple good rotamers are predicted? We attempt to solve these problems by choosing position-specific rather than generalized rotamers and by sorting the residues that have to be modelled as a function of their freedom in rotamer space. The practical advantages of our method are the quality of the models for cases of high backbone similarity, the small amount of human intervention needed, and the fact that the method automatically estimates the reliability with which each residue has been modeled. Other methods described in this issue are probably more suitable if large backbone rearrangements or loop insertions and deletions need to be modeled.

Structure-function relationship for the highly toxic crotoxin from Crotalus durissus terrificus.

The three-dimensional structure of the highly toxic crotoxin from Crotalus durissus terrificus was modelled based on sequence analysis and the refined structure of calcium-free phospholipase of Crotalus atrox venom. Small-angle x-ray scattering experiments were performed on aqueous solutions of crotoxin. The radial distribution function derived from these scattering experiments and the one calculated from the model structure are in good agreement. Crotoxin consists of a basic and an acidic subunit. The model strongly suggests that the overall folding motif of phospholipases has been preserved in both subunits. The basic domain has an intact active site. The residues that are expected to contact the lipid tails of the phospholipid are different from other phospholipases, but they are all hydrophobic. The acidic domain consists of three independent chains interconnected by disulfide bonds. Compared to other phospholipases the active site for the greater part has been preserved in this domain, but it is not very well shielded from solvent. Most residues normally in contact with the lipid tails of the phospholipid are missing, which might explain the acidic subunit's lack of phospholipase activity. A homology between the third chain of the acidic domain and neurophysins suggests that the acidic domain may act as a chaperone for the basic domain.