Search

The ability of a species to adapt to environmental change is ultimately reflected in its vital rates – i.e. survival and reproductive success of individuals. Together, vital rates determine trends in numbers, commonly monitored using counts of species abundance. Rapid changes in abundance can give rise to concern, leading to calls for research into the biological mechanisms underlying variations in demography. For the northwest European population of Bewick's swan Cygnus columbianus bewickii, there have been major changes in the population trends recorded during nearly five decades of monitoring (1970–2016). The total number of birds increased to a maximum of ca 30 000 in 1995 and subsequently decreased to about 18 000 individuals in 2010. Such large fluctuation in population numbers is rare in long‐lived species and understanding the drivers of this population change is crucial for species management and conservation. Using the integrated population model (IPM) framework, we analysed three demographic datasets in combination: population counts, capture–mark–resightings (CMR) and the proportion of juveniles in winter over a period of ~50 years. We found higher apparent breeding success in the years when the population had a positive growth rate compared to years with a negative growth rate. Moreover, no consistent trend in adult and yearling survival, and an increasing trend in juvenile survival was found. A transient life‐table response experiment showed that apparent breeding success and adult survival contributed most to the variation in population trend. We explored possible explanatory variables for the different demographic rates and found a significant association between juvenile survival both with the water level in lakes during autumn migration, which affects food accessibility for the swans, and with summer temperatures. Such associations are important for understanding the dynamics of species with fluctuating population sizes, and thus for informing management and conservation decisions.

Identifying the environmental drivers of variation in fitness-related traits is a central objective in ecology and evolutionary biology. Temporal fluctuations of these environmental drivers are often synchronized at large spatial scales. Yet, whether synchronous environmental conditions can generate spatial synchrony in fitness-related trait values (i.e., correlated temporal trait fluctuations across populations) is poorly understood. Using data from long-term monitored populations of blue tits (Cyanistes caeruleus, n = 31), great tits (Parus major, n = 35), and pied flycatchers (Ficedula hypoleuca, n = 20) across Europe, we assessed the influence of two local climatic variables (mean temperature and mean precipitation in February–May) on spatial synchrony in three fitness-related traits: laying date, clutch size, and fledgling number. We found a high degree of spatial synchrony in laying date but a lower degree in clutch size and fledgling number for each species. Temperature strongly influenced spatial synchrony in laying date for resident blue tits and great tits but not for migratory pied flycatchers. This is a relevant finding in the context of environmental impacts on populations because spatial synchrony in fitness-related trait values among populations may influence fluctuations in vital rates or population abundances. If environmentally induced spatial synchrony in fitness-related traits increases the spatial synchrony in vital rates or population abundances, this will ultimately increase the risk of extinction for populations and species. Assessing how environmental conditions influence spatiotemporal variation in trait values improves our mechanistic understanding of environmental impacts on populations. birds, climate, clutch size, comparative analysis, fitness-related traits, fledgling number, phenology, spatial synchrony, timing of breeding, weather

Adaptive topography is a central concept in evolutionary biology, describing how the mean fitness of a population changes with gene frequencies or mean phenotypes.We use expected population size as a quantity to be maximized by natural selection to show that selection on pairwise combinations of reproductive traits of collared flycatchers caused by fluctuations in population size generated an adaptive topography with distinct peaks often located at intermediate phenotypes. This occurred because r- and K-selection made phenotypes favored at small densities different from those with higher fitness at population sizes close to the carrying capacity K. Fitness decreased rapidly with a delay in the timing of egg laying, with a densitydependent effect especially occurring among early-laying females. The number of fledglings maximizing fitness was larger at small population sizes than when close to K. Finally, there was directional selection for large fledglings independent of population size. We suggest that these patterns can be explained by increased competition for some limiting resources or access to favorable nest sites at high population densities. Thus, r- and K-selection based on expected population size as an evolutionary maximization criterion may influence life-history evolution and constrain the selective responses to changes in the environment.

The phenology of many species shows strong sensitivity to climate change; however, with few large scale intra-specific studies it is unclear how such sensitivity varies over a species’ range. We document large intra-specific variation in phenological sensitivity to temperature using laying date information from 67 populations of two European songbirds covering a large part of their breeding range. Populations inhabiting deciduous habitats showed stronger phenological sensitivity compared with those in evergreen and mixed habitats. Strikingly, however, the lowest sensitivity was seen in populations that had experienced the greatest change in climate. Therefore, we predict that the strongest phenological advancement will not occur in those populations with the highest sensitivity. Our results show that to effectively assess the impact of climate change on phenology across a species range it will be necessary to account for intra-specific variation in phenological sensitivity, climate change exposure, and the ecological characteristics of a population.

The integration and synthesis of the data in different areas of science is drastically slowed and hindered by a lack of standards and networking programmes. Long-term studies of individually marked animals are not an exception. These studies are especially important as instrumental for understanding evolutionary and eco-logical processes in the wild. Furthermore, their number and global distribution provides a unique opportunity to assess the generality of patterns and to address broad-scale global issues (e.g. climate change)., To solve data integration issues and enable a new scale of ecological and evolution-ary research based on long-term studies of birds, we have created the SPI-Birds Network and Database (www.spibirds.org)—a large-scale initiative that connects data from, and researchers working on, studies of wild populations of individually recognizable (usually ringed) birds. Within year and a half since the establishment, SPI-Birds has recruited over 120 members, and currently hosts data on almost 1.5 million individual birds collected in 80 populations over 2,000 cumulative years, and counting., SPI-Birds acts as a data hub and a catalogue of studied populations. It prevents data loss, secures easy data finding, use and integration and thus facilitates collab-oration and synthesis. We provide community-derived data and meta-data stand-ards and improve data integrity guided by the principles of Findable, Accessible, Interoperable and Reusable (FAIR), and aligned with the existing metadata lan-guages (e.g. ecological meta-data language)., The encouraging community involvement stems from SPI-Bird's decentralized ap-proach: research groups retain full control over data use and their way of data management, while SPI-Birds creates tailored pipelines to convert each unique data format into a standard format. We outline the lessons learned, so that other communities (e.g. those working on other taxa) can adapt our successful model. Creating community-specific hubs (such as ours, COMADRE for animal demogra-phy, etc.) will aid much-needed large-scale ecological data integration., The SPI-Birds have been supported by an NWO personal grant (grant number 016.Veni.181.054) to A.C., and a Research Council of Norway grant: 223257 (SFF-III) and 267511 (EVOCLIM).

International audience; Understanding species coexistence has long been a major goal of ecology. Coexistence theory for two competing species posits that intraspecific density dependence should be stronger than interspecific density dependence. Great tits and blue tits are two bird species that compete for food resources and nesting cavities. On the basis of long-term monitoring of these two competing species at sites across Europe, combining observational and manipulative approaches, we show that the strength of density regulation is similar for both species, and that individuals have contrasting abilities to compete depending on their age. For great tits, density regulation is driven mainly by intraspecific competition. In contrast, for blue tits, interspecific competition contributes as much as intraspecific competition, consistent with asymmetric competition between the two species. In addition, including age-specific effects of intra-and interspecific competition in density-dependence models improves predictions of fluctuations in population size by up to three times.

Genomic and proteomic projects are producing a flood of data that all require interpretation which often is best performed based on a three dimensional structure of the molecule(s) involved. These structures can be determined experimentally, or modelled by homology. Because of the complexity of the questions and the heterogeneity of the data, the software used for modelling proteins must become even more versatile. We describe several case studies in which the questions asked, the data, and the requirements on the software all are very different. It is shown how structural knowledge about a protein helps to determine the best place to bind a fluorescent dye. Such dyes are needed to determine protein–protein, protein–DNA interactions or intrinsic fluorescence microscopy. Further, using dyes you can trace molecules in the cell and thus get a handle on subcellular localisation. The first example (OCT-1) involves the search for free amino groups in a protein–DNA complex. The second example (BPTI) is a case, in which the amino acid distribution shows that amino groups are spread all over the structure, so that the natural structure has to be modified to get an answer. The third example (HFE) involves a model built by homology. In this case the amino group distribution can also be predicted. All these studies were performed using the what if software package. This package is available including source code, documentation, etc. See http://www.cmbi.kun.nl/whatif/

Protein stabilization by immobilization has been proposed to be most effective if the protein is attached to the carrier at that region where unfolding is initiated. To probe this hypothesis, we have studied the effects of site-specific immobilization on the thermal stability of mutants of the thermolysin-like protease from Bacillus stearothermophilus (TLP-ste). This enzyme was chosen because previous studies had revealed which parts of the molecule are likely to be involved in the early steps of thermal unfolding. Cysteine residues were introduced by site-directed mutagenesis into various positions of a cysteine-free variant of TLP-ste. The mutant enzymes were immobilized in a site-specific manner onto Activated Thiol-Sepharose. Two mutants (T56C, S65C) having their cysteine in the proposed unfolding region of TLP-ste showed a 9- and 12-fold increase in half-lives at 75 degrees C due to immobilization. The stabilization by immobilization was even larger (33-fold) for the T56C/S65C double mutant enzyme. In contrast, mutants containing cysteines in other parts of the TLP-ste molecule (N181C, S218C, T299C) showed only small increases in half-lives due to immobilization (maximum 2.5-fold). Thus, the stabilization obtained by immobilization was strongly dependent on the site of attachment. It was largest when TLP-ste was fixed to the carrier through its postulated unfolding region. The concept of the unfolding region may be of general use for the design of strategies to stabilize proteins.

The human cDNA encoding prolyl endopeptidase, a cytoplasmic endoprotease which hydrolyses the peptide bond at the C-terminal side of proline, was sequenced. After the isolation of the 3′ terminal fragment of the pep cDNA sequence from a human lymphocyte cDNA library, an approach based on the polymerase chain reaction (PCR) was undertaken to obtain the complete pep cDNA. Overlapping DNA fragments were generated by PCR from cDNA synthesized from human lymphocyte mRNA. The DNA fragments were subcloned and sequenced. The complete cDNA is 2562 nucleotides (nt) in length and contains an open reading frame coding for a protein of 710 amino acids (aa). Comparison of the primary PEP sequences from human lymphocyte and pig brain shows 97% identity. The aa sequence analysis shows homology with bacterial PEPs and with protease II from Escherichia coli. Asp641 probably participates in the active site of PEP.

A major incentive in inhibitor research is that control of limited proteolysis constitutes a valuable pharmacological tool. Protease inhibitors have proved to be successful in influencing pathogenesis in many experimental models but a breakthrough to use in human therapy has mainly been restricted to aprotinin and angiotensin converting enzyme (ACE) inhibitors. However, the success of ACE inhibitors as pharmacological tools in hypertension has proved to be a strong stimulant for new protease inhibitor approaches to drug therapy. While emphasis in the search for next generations of ACE inhibitors may move from the circulation renin-angiotensin system to the local tissue systems, including heart, brain and genital tract, persistent and insightful design of renin inhibitors has already yielded highly specific molecules with potent activities in several in vivo models. The development of orally effective long-acting inhibitors will finally allow an evaluation to be made of their therapeutic profile with regard to the family of ACE inhibitors. The close relationship between renin and HIV-1 protease presents an exceptional opportunity for transfer of the knowledge acquired in renin inhibitor development during the past decade, to an accelerated generation of specific HIV-1 protease inhibitors as effective agents in treatment of AIDS. The self-assembly of 2 identical monomers into a symmetrical structure in HIV-1 protease is not only an elegant way to create an active enzyme while encoding a minimal amount of genetic information, but is also in concordance with the bilobular active-site found in mammalian aspartic proteases.(ABSTRACT TRUNCATED AT 250 WORDS)

Many GPCR models have been built over the years for many different purposes, of which drug-design undoubtedly has been the most frequent one. The release of the structure of bovine rhodopsin in August 2000 enabled us to analyze models built before that period to learn things for the models we build today. We conclude that the GPCR modeling field is riddled with "common knowledge". Several characteristics of the bovine rhodopsin structure came as a big surprise, and had obviously not been predicted, which led to large errors in the models. Some of these surprises, however, could have been predicted if the modelers had more rigidly stuck to the rule that holds for all models, namely that a model should explain all experimental facts, and not just those facts that agree with the modeler's preconceptions.

The hydrophobic S1' subsite is one of the major determinants of the substrate specificity of thermolysin and related M4 family proteases. In the thermolysin-like protease (TLP) produced by Bacillus stearothermophilus (TLP-ste), the hydrophobic S1' subsite is mainly formed by Phe130, Phe133, Val139 and Leu202. In the present study, we have examined the effects of replacing Leu202 by smaller (Gly, Ala, Val) and larger (Phe, Tyr) hydrophobic residues. The mutational effects showed that the wild-type S1' pocket is optimal for binding leucine side chains. Reduction of the size of residue 202 resulted in a higher efficiency towards substrates with Phe in the P1' position. Rather unexpectedly, the Leu202--Phe and Leu202--Tyr mutations, which were expected to decrease the size of the S1' subsite, resulted in a large increase in activity towards dipeptide substrates with Phe in the P1' position. This is probably due to the fact that 202Phe and 202Tyr adopt a second possible rotamer that opens up the subsite compared to Leu202, and also favours interactions with the substrate. To validate these results, we constructed variants of thermolysin with changes in the S1' subsite. Thermolysin and TLP-ste variants with identical S1' subsites were highly similar in terms of their preference for Phe vs. Leu in the P1' position.

A comparison is made between the distribution of residue preferences, three dimensional nearest neighbour contacts, preferred rotamers, helix-helix crossover angles and peptide bond angles in three sets of proteins: a non-redundant set of accurately determined globular protein structures, a set of four-helix bundle structures and a set of membrane protein structures. Residue preferences for the latter two sets may reflect overall helix stabilising propensities but may also highlight differences arising out of the contrasting nature of the solvent environments in these two cases. The results bear out the expectation that there may be differences between residue type preferences in membrane proteins and in water soluble globular proteins. For example, the beta-branched residue types valine and isoleucine are considerably more frequently encountered in membrane helices. Likewise, glycine and proline. residue types normally associated with 'helix-breaking' propensity are found to be relatively more common in membrane helices. Three dimensional nearest neighbour contacts along the helix, preferred rotamers, and peptide bond angles are very similar in the three sets of proteins as far as can be ascertained within the limits of the relatively low resolution of the membrane proteins dataset. Crossing angles for helices in the membrane protein set resemble the four helix bundle set more than the general non-redundant set, but in contrast to both sets they have smaller crossing angles consistent with the dual requirements for the helices to form a compact structure while having to span the membrane. In addition to the pairwise packing of helices we investigate their global packing and consider the question of helix supercoiling in helix bundle proteins.

A statistical analysis is reported of 1,200 of the 1,404 nuclear magnetic resonance (NMR)-derived protein and nucleic acid structures deposited in the Protein Data Bank (PDB) before 1999. Excluded from this analysis were the entries not yet fully validated by the PDB and the more than 100 entries that contained95% of the expected hydrogens. The aim was to assess the geometry of the hydrogens in the remaining structures and to provide a check on their nomenclature. Deviations in bond lengths, bond angles, improper dihedral angles, and planarity with respect to estimated values were checked. More than 100 entries showed anomalous protonation states for some of their amino acids. Approximately 250,000 (1.7%) atom names differed from the consensus PDB nomenclature. Most of the inconsistencies are due to swapped prochiral labeling. Large deviations from the expected geometry exist for a considerable number of entries, many of which are average structures. The most common causes for these deviations seem to be poor minimization of average structures and an improper balance between force-field constraints for experimental and holonomic data. Some specific geometric outliers are related to the refinement programs used. A number of recommendations for biomolecular databases, modeling programs, and authors submitting biomolecular structures are given.

Conformational changes are known to play a crucial role in the function of the bacterial GroE chaperonin system. Here, results are presented from an essential dynamics analysis of known experimental structures and from computer simulations of GroEL using the CONCOORD method. The results indicate a possible direct form of inter-ring communication associated with internal fluctuations in the nucleotide-binding domains upon nucleotide and GroES binding that are involved in the allosteric mechanism of GroEL. At the level of conformational transitions in entire GroEL rings, nucleotide-induced structural changes were found to be distinct and in principle uncoupled from changes occurring upon GroES binding. However, a coupling is found between nucleotide-induced conformational changes and GroES-mediated transitions, but only in simulations of GroEL double rings, and not in simulations of single rings. This provides another explanation for the fact that GroEL functions a double ring system.

Starting from a cysteine-free mutant of the thermolysin-like neutral protease from Bacillus stearothermophilus, cysteines were introduced into different positions on the enzyme surface by site-directed mutagenesis. The mutant enzymes were immobilized via the SH-groups to Activated Thiol-Sepharose 4B and their thermostabilities were compared to those of the soluble enzymes. The results showed that the effects of immobilization on stability strongly depend on the site of attachment. Binding of the enzyme via engineered cysteines in the critical unfolding region between residues 56 and 69 led to a considerable increase of thermal stability, whereas the immobilization via a cysteine introduced remote from the unfolding region yielded less stabilization. An extremely strong stabilization was obtained upon binding via T56C where the half-life at 75 °C was increased by the factor of 24.

A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances measured in experimental structures and the strength of the interatomic interaction. Structural differences between generated structures are similar to those obtained from experiment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large number of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints.

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.

A family of G-protein-coupled receptors is believed to mediate the recognition of odor molecules. In order to identify potential ligand-binding residues, we have applied correlated mutation analysis to receptor sequences from the rat. This method identifies pairs of sequence positions where residues remain conserved or mutate in tandem, thereby suggesting structural or functional importance. The analysis supported molecular modeling studies in suggesting several residues in positions that were consistent with ligand-binding function. Two of these positions, dominated by histidine residues, may play important roles in ligand binding and could confer broad specificity to mammalian odor receptors. The presence of positive (overdominant) selection at some of the identified positions provides additional evidence for roles in ligand binding. Higher-order groups of correlated residues were also observed. Each group may interact with an individual ligand determinant, and combinations of these groups may provide a multi-dimensional mechanism for receptor diversity.

An extracellular lipase (triacylglycerol acylhydrolase EC 3.1.1.3), produced by the fungus Rhizopus javanicus was purified to homogeneity using an expeditious two-step isolation method. The enzyme, with a molecular mass of 36 kDa and a specific activity of 9260 microequivalent of fatty acid released per minute and mg under standard conditions, consists of three isoforms with isoelectric points of 7.8, 7.7, and 7.1, respectively. The purified lipase was digested using chemical and enzymatical procedures: CNBr cleavage, partial acid hydrolysis, and proteolytic cleavage by means of trypsin. Amino-acid sequencing of the resulting peptides indicates that the three lipases from Rhizopus javanicus, Rhizopus niveus and Rhizopus delemar are produced as identical proenzymes but processed differently. These Rhizopus lipases show 54% identity with the lipase from Rhizomucor miehei. Using the structure of the Rhizomucor miehei lipase, the molecular model of Rhizopus javanicus lipase was constructed. Both enzymes are alpha/beta type proteins with a central 8-stranded mixed beta-pleated sheet and have a remarkably similar distribution of hydrophobic amino acids at their surface. The tryptophan in the center of the helical lid covering the active site of Rhizomucor miehei lipase is mutated into an alanine, indicating that it is not essential for the proper movement of the helical lid.

The progesterone-binding protein uteroglobin has been expressed in Escherichia coli in an unfused, soluble form. Like mature uteroglobin from rabbit endometrium (UG), the E.coli produced uteroglobin (UG1) dimerizes in vitro, forms an antiparallel dimer with Cys3-Cys69' and Cys69-Cys3' disulfide bonds and binds progesterone under reducing conditions. In order to analyze the dimerization and the reduction dependence of progesterone binding in more detail, we separately replaced cysteine 3 and cysteine 69 by serines. Under reducing conditions, both uteroglobin variants (UG1-3Ser and UG1-69Ser) bind progesterone with the same affinity as the wild-type suggesting that both cysteine residues are not directly involved in progesterone binding. In contrast to the wild-type protein, both cysteine variants also bind progesterone with high affinity in the absence of reducing agents. In addition, UG1-3Ser and UG1-69Ser both form covalently linked homodimers. Thus, unnatural Cys69-69' and Cys3-3' disulfide bonds exist in UG1-3Ser and UG1-69Ser, respectively. These data together with computer models based on X-ray diffraction data strongly support the idea that progesterone reaches its binding site located in an internal hydrophobic cavity via a hydrophobic tunnel along helices 1 and 4. Under non-reducing conditions the tunnel is closed by two disulfide bridges (Cys3-Cys69' and Cys69-Cys3') that lie in the most flexible region of the dimer. Reduction or replacement of a cysteine residue enables conformational changes that open the channel allowing progesterone to enter.

The structure of the enzymatically active subunit of human plasma carboxypeptidase N was determined by computer aided model building by homology using the structural coordinates from carboxypeptidase A. The active site of carboxypeptidase N has been well conserved in comparison with carboxypeptidase A. Differences in substrate specificity can be explained by the comparison of energetically favorable binding sites for different atomic probe groups.

The structure of the enzymatically active subunit of human plasma carboxypeptidase N was determined by computer aided model building by homology using the structural coordinates from carboxypeptidase A. The active site of carboxypeptidase N has been well conserved in comparison with carboxypeptidase A. Differences in substrate specificity can be explained by the comparison of energetically favorable binding sites for different atomic probe groups.

Triosephosphate isomerase (TIM) is a dimeric glycolytic enzyme. TIM from Trypanosoma brucei brucei has been crystallized at pH 7.0 in 2.4 M-ammonium sulphate. The well-diffracting crystals have one dimer per asymmetric unit. The structure has been refined at 1.83 A resolution with an R-factor of 18.3% for all data between 6 A and 1.83 A (37,568 reflections). The model consists of 3778 protein atoms and 297 solvent atoms. Subunit 1 is involved in considerably more crystal contacts than subunit 2. Correlated with these differences in crystal packing is the observation that only in the active site of subunit 2 is a sulphate ion bound. Furthermore, significant differences with respect to structure and flexibility are observed in three loops near the active site. In particular, there is a 7 A positional difference of the tip of the flexible loop (loop 6) when comparing subunit 1 and subunit 2. Also, the neighbouring loops (loop 5 and loop 7) have significantly different conformations and flexibility. In subunit 1, loop 6 is in an "open" conformation, in subunit 2, loop 6 is in an "almost closed" conformation. Only in the presence of a phosphate-containing ligand, such as glycerol-3-phosphate, does loop 6 take up the "closed" conformation. Loop 6 and loop 7 (and also to some extent loop 5) are rather flexible in the almost closed conformation, but well defined in the open and closed conformations. The closing of loop 6 (167 to 180), as observed in the almost closed conformation, slightly changes the main-chain conformation of the catalytic glutamate, Glu167, leading to a change of the chi 1 angle of this residue from approximately -60 degrees to approximately 60 degrees and the weakening of the hydrogen bonds between its polar side-chain atoms and Ser96. In the closed conformation, in the presence of glycerol-3-phosphate, the main-chain atoms of Glu167 remain in the same position as in the almost closed conformation, but the side-chain has rotated around the CA-CB bond changing chi 1 from approximately 60 degrees to approximately -60 degrees. In this new position the hydrogen bonding to Ser96 is completely lost and also a water-mediated salt bridge between OE2(Glu167) and NE(Arg99) is lost. Comparison of the two independently refined subunits, showed that the root-mean-square deviation for all 249 CA atoms is 0.9 A; for the CA atoms of the beta-strands this is only 0.2 A. The average B-factor for all subunit 1 and subunit 2 atoms is 20 A2 and 25 A2, respectively.(ABSTRACT TRUNCATED AT 400 WORDS)

The thermostability of the B. subtilis neutral protease was studied under various conditions. At elevated temperatures the enzyme was inactivated as a result of autolysis. The rate of inactivation did not depend on the enzyme concentration and the enzyme was most stable near its pH optimum. The rate of inactivation was unaffected by the presence of a second protease during the incubation at high temperatures. The results indicate that the rate of thermal inactivation of the neutral protease is determined by the kinetics of local unfolding processes that precede autolysis rather than by the catalytic rate of the autodigestion reaction or an irreversible unfolding step.

The oxidation of alkanes to alkanols by Pseudomonas oleovorans involves a three-component enzyme system: alkane hydroxylase, rubredoxin and rubredoxin reductase. Alkane hydroxylase and rubredoxin are encoded by the alkBFGHJKL operon, while previous studies indicated that rubredoxin reductase is most likely encoded on the second alk cluster: the alkST operon. In this study we show that alkT encodes the 41 x 10(3) Mr rubredoxin reductase, on the basis of a comparison of the expected amino acid composition of AlkT and the previously established amino acid composition of the purified rubredoxin reductase. The alkT sequence revealed significant similarities between AlkT and several NAD(P)H and FAD-containing reductases and dehydrogenases. All of these enzymes contain two ADP binding sites, which can be recognized by a common beta alpha beta-fold or fingerprint, derived from known structures of cofactor binding enzymes. By means of this amino acid fingerprint we were able to determine that one ADP binding site in rubredoxin reductase (AlkT) is located at the N terminus and is involved in FAD binding, while the second site is located in the middle of the sequence and is involved in the binding of NAD or NADP. In addition, we derived from the sequences of FAD binding reductases a second amino acid fingerprint for FAD binding, and we used this fingerprint to identify a third amino acid sequence in AlkT near the carboxy terminus for binding of the flavin moiety of FAD. On the basis of the known architecture and relative spatial orientations of the NAD and FAD binding sites in related dehydrogenases, a model for part of the tertiary structure of AlkT was developed.

The interaction of the oligonucleotides (Ap)8A and (A-T)5 with empty capsids of the coat protein of cowpea chlorotic mottle virus (CCMV) has been studied with 500 MHz 1H nuclear magnetic resonance. It is found that these oligonucleotides specifically bind to the arginine and lysine residues of the N-terminal arm of the protein. Upon this binding, immobilization of part of the N-terminal arm occurs. In addition, secondary structure predictions and energy calculations have been performed on the N-terminal arm. These calculations were carried out as a function of the charges on the arginine and lysine side-chains. For free coat protein, where the arginine and lysine side-chains are charged, the arm is found in a random-coil conformation. In the neutralized state, as for the coat protein in the virus, the arm adopts an alpha-helical conformation. The results support a previously published model for the assembly of CCMV, in which a random-coil to alpha-helix conformational transition, induced by neutralizing the arginine and lysine side-chains, plays an essential role.

In a previous paper [I] we have shown that the protein region in cowpea chlorotic mottle virus (CCMV) involved in protein-RNA interaction is very mobile in the absence of RNA. Upon binding RNA, immobilization of this protein region occurred. From a comparison with other systems, e.g. tobacco mosaic virus [2-41, tomato bushy stunt virus [5], southern bean mosaic virus [6], histones [7] and protamines [8] and using protein secondary structure predictions [9] we have suggested a general model for the protein-RNA interaction in simple plant viruses. This model is based on a random coil to a-helix transition for the protein region, interacting with RNA upon binding. In this paper we extend our work to three other spherical plant viruses. We used brome mosaic virus (BMV), belladonna mottle virus (BdMV) and cowpea mosaic virus (CPMV) since the stability of these three viruses depends to a different extent upon protein-RNA and protein-protein interactions [lo]. The results obtained confirm our previous ideas about protein-RNA interaction in simple plant viruses and permit a further retinement of our model.

Secondary structural elements of the phosphatidylcholine-transfer protein from bovine liver have been predicted from its primary structure with the aid of two computerized methods. The predicted alpha-helix and beta-strand content have been compared with the values derived from circular dichroism spectra. The hydrophobicity profile (Rose plot) of the protein indicated that the supposed lipid-binding site occurs in the most hydrophobic region. The predicted secondary structural elements have been folded in a tentative model of the protein molecule according to its hydrophobicity profile.

One of the two SH groups of the coat protein of Cowpea chlorotic mottle virus (CCMV) was spin-labeled with 3-maleimido-2,2,5,5-tetramethyl-l-pyrrolodinyloxyl. Both electron spin resonance and saturation-transfer electron spin resonance spectra were recorded of intact CCMV and CCMV in which the protein subunits are cross-linked with glutaraldehyde. From an analysis of these spectra it follows that the spin label performs a local anisotropic motion with respect to the protein to which it is attached. It is also concluded that no internal protein mobility is present on the ST-ESR timescale.

Using synchrotron radiation, the X-ray diffraction intensities of crystals of p-hydroxy-benzoate hydroxylase, complexed with the substrate p-hydroxybenzoate, were measured to a resolution of 1.9 A. Restrained least-squares refinement alternated with rebuilding in electron density maps yielded an atom model of the enzyme-substrate complex with a crystallographic R-factor of 15.6% for 31,148 reflections between 6.0 and 1.9 A. A total of 330 solvent molecules was located. In the final model, only three residues have deviating phi-psi angle combinations. One of them, the active site residue Arg44, has a well-defined electron density and may be strained to adopt this conformation for efficient catalysis. The mode of binding of FAD is distinctly different for the different components of the coenzyme. The adenine ring is engaged in three water-mediated hydrogen bonds with the protein, while making only one direct hydrogen bond with the enzyme. The pyrophosphate moiety makes five water-mediated versus three direct hydrogen bonds. The ribityl and ribose moieties make only direct hydrogen bonds, in all cases, except one, with side-chain atoms. The isoalloxazine ring also makes only direct hydrogen bonds, but virtually only with main-chain atoms. The conformation of FAD in p-hydroxybenzoate hydroxylase is strikingly similar to that in glutathione reductase, while the riboflavin-binding parts of these two enzymes have no structural similarity at all. The refined 1.9 A structure of the p-hydroxybenzoate hydroxylase-substrate complex was the basis of further refinement of the 2.3 A structure of the enzyme-product complex. The result was a final R-factor of 16.7% for 14,339 reflections between 6.0 and 2.3 A and an improved geometry. Comparison between the complexes indicated only small differences in the active site region, where the product molecule is rotated by 14 degrees compared with the substrate in the enzyme-substrate complex. During the refinements of the enzyme-substrate and enzyme-product complexes, the flavin ring was allowed to bend or twist by imposing planarity restraints on the benzene and pyrimidine ring, but not on the flavin ring as a whole. The observed angle between the benzene ring and the pyrimidine ring was 10 degrees for the enzyme-substrate complex and 19 degrees for the enzyme-product complex. Because of the high temperature factors of the flavin ring in the enzyme-product complex, the latter value should be treated with caution. Six out of eight peptide residues near the flavin ring are oriented with their nitrogen atom pointing towards the ring.(ABSTRACT TRUNCATED AT 400 WORDS)

Cowpea chlorotic mottle virus (CCMV) is a spherical plantvirus of the group of bromoviruses. An interesting feature of bromoviruses is the increase in hydrodynamic volume when the pH is raised from 5.0-7.5 at low ionic strength (

Proteins are key to many different processes in and around living cells, they are the molecular mechanisms that facilitate life. Therefore, protein function and structure have been of great interest to researchers over the past decades. A better understanding of how proteins work is important, as it will help us understand the various diseases caused by their malfunction and help us apply proteins in various applications, such as industrial processes. The research on proteins has resulted in an ever-increasing amount of data related to proteins.Protein information systems were developed to make sense out of all the data available for specific protein classes. The work described in this thesis makes use of the 3DM platform to store and order the enormous amount of protein data that is nowadays available. The 3DM platform integrates heterogeneous datasets for complete protein superfamilies to generate molecular class specific information systems (3DM systems).The main objective of this thesis is to facilitate a better understanding of protein function through data integration- and analysis-tools to support both DNA diagnostics and protein engineering efforts. To that end, chapter 2 and chapter 3 focus on the collection and integration of mutational data from scientific literature, followed by protein superfamily analysis to reveal predictive patterns. In chapter 4, mutational data mined for whole families are combined with a superfamily alignment analysis to assign function to networks of residues that appear to have co-evolved. And in chapter 5, the data of a 3DM system is used to train a machine learning model for the prediction of pathogenicity of novel genetic variants.In chapter 2 a method for the automated collection and integration of mutation data from scientific literature is described. This text-mining tool, Mutator, was designed to automatically identify and retrieve relevant literature for a complete protein family.Mutator was tested on the alpha-amylase protein family which contains a protein that, when mutated, can lead to Fabry’s disease. The Fabry disease mutation database was constructed which contained 30% more mutation data when compared to the manually curated HGMD. To differentiate between natural variation and pathogenic mutations, alignment patterns in the alpha-amylase superfamily were investigated and Validator was developed, a tool to evaluate the potential pathogenicity of individual mutations.In a follow up study, chapter 3, nine disease-related proteins across two different protein families are investigated. An improved version of Mutator was developed and used to collect mutational data for these protein families. A comparison to the manually curated HGMD showed that Mutator extracts far more mutation data, both in terms of unique mutations extracted, as well as number of references per mutation.These mutation datasets were analyzed in the context of the protein families to identify patterns to assess pathogenicity for novel variants. This analysis showed that pathogenic mutations would more often be observed in the structurally conserved core of the family than in structurally variable regions. A significant overlap was observed of protein positions with pathogenic mutations between very distant homologs, indicating that data can be transferred between equivalent residue positions even across large phylogenetic distances. Protein superfamily information systems that use structural alignments, such as 3DM, can facilitate the transfer of data across equivalent residues of structurally conserved positions.Chapter 4 describes a web-based tool for the analysis of co-evolving residues in protein super-family alignments. This tool, CorNet, combines networks of co-evolving positions with literature data to shed light on the function of these residues. When keywords, such as “specificity”, are mentioned in the sentences that contain mutational data for an alignment position, these were used to annotate the alignment position of a network.For six different protein families, a list of literature keywords was used to annotate network positions and to calculate an enrichment score of any keyword to be specifically associated with a network of co-evolving positions. Enriched keywords indicate the function behind evolutionary pressures that caused networks of residues to co-evolve. It was shown that an assigned function can be transferred from annotated to non-annotated positions in a network with a mutagenesis study that targeted only the non-annotated positions of a network enriched for the keyword “enantioselectivity”.In chapter 5, the integrated data of 3DM systems is utilized to train a model to predict variant pathogenicity with machine learning techniques. For three genes related to Long QT syndrome and Brugada syndrome, both cardiac arrhythmia syndromes, datasets of both benign and pathogenic variants were collected. 3DM was used to generate much heterogeneous data to describe each variant. Using this data, a pathogenicity predictor was generated to separate disease-causing from non-disease causing variants. This predictor outperformed frequently used genome-wide prediction tools and even gene-family specific predictors.In the general discussion an outlook is given on the scale-up of protein information systems to cover the complete structural space and even the whole human exome. This thesis provides a variety of avenues to apply integrated data to gain insight in proteins and support protein related research and development.

As the active component of many biological systems, proteins are of great interest to life scientists. Proteins are used in a large number of different applications such as the production of precursors and compounds, for bioremediation, as drug targets, to diagnose patients suffering from genetic disorders, etc. Many research projects have therefore focused on the characterization of proteins and on improving the understanding of the functional and mechanistic properties of proteins. Studies have examined folding mechanisms, reaction mechanisms, stability under stress, effects of mutations, etc. All these research projects have resulted in an enormous amount of available data in lots of different formats that are difficult to retrieve, combine, and use efficiently. The main topic of this thesis is the 3DM platform that was developed to generate Molecular Class Specific Information Systems (3DM systems) for protein superfamilies. These superfamily systems can be used to collect and interlink heterogeneous data sets based on structure based multiple sequence alignments. 3DM systems can be used to integrate protein, structure, mutation, reaction, conservation, correlation, contact, and many other types of data. Data is visualized using websites, directly in protein structures using YASARA, and in literature using Utopia Documents. 3DM systems contain a number of modules that can be used to analyze superfamily characteristics namely Comulator for correlated mutation analyses, Mutator for mutation retrieval, and Validator for mutant pathogenicity prediction. To be able to determine the characteristics of subsets of proteins and to be able to compare the characteristics of different subsets a powerful filtering mechanism is available. 3DM systems can be used as a central knowledge base for projects in protein engineering, DNA diagnostics, and drug design. The scientific and technical background of the 3DM platform is described in the first two chapters. Chapter 1 describes the scientific background, starting with an overview of the foundations of the 3DM platform. Alignment methods and tools for both structure and sequence alignments, and the techniques used in the 3DM modules are described in detail. Alternative methods are also described with the advantages and disadvantages of the various strategies. Chapter 2 contains a technical description of the implementation of the 3DM platform and the 3DM modules. A schematic overview of the database used to store the data is provided together with a description of the various tables and the steps required to create new 3DM systems. The techniques used in the Comulator, Mutator and Validator modules of the 3DM platforms are discussed in more detail. Chapter 3 contains a concise overview of the 3DM platform, its capabilities, and the results of protein engineering projects using 3DM systems. Thirteen 3DM systems were generated for superfamilies such as the PEPM/ICL and Nuclear Receptors. These systems are available online for further examination. Protein engineering studies aimed at optimizing substrate specificity, enzyme activity, or thermostability were designed targeting proteins from these superfamilies. Preliminary results of drug design and DNA diagnostics projects are also included to highlight the diversity of projects 3DM systems can be applied to. Project HOPE: a biomedical tool to predict the effect of a mutation on the structure of a protein is described in chapter 4. Project HOPE is developed at the Radboud University Nijmegen Medical Center under supervision of H. Venselaar. Project HOPE employs webservices to optimally reuse existing databases and computing facilities. After selection of a mutant in a protein, data is collected from various sources such as UniProt and PISA. A homology model is created to determine features such as contacts and side-chain accessibility directly in the structure. Using a decision tree, the available data is evaluated to predict the effects of the mutation on the protein. Chapter 5 describes Comulator: the 3DM module for correlated mutation analyses. Two positions in an alignment correlate when they co-evolve, that is they mutate simultaneously or not at all. Comulator uses a statistical coupling algorithm to calculate correlated mutation analyses. Correlated mutations are visualized using heatmaps, or directly in protein structures using YASARA. Analyses of correlated mutations in various superfamilies showed that positions that correlate are often found in networks and that the positions in these networks often share a common function. Using these networks, mutants were predicted to increase the specificity or activity of proteins. Mutational studies confirmed that correlated mutation analyses are a valuable tool for rational design of proteins. Mutator, the text mining tool used to incorporate mutations into 3DM systems is described in chapter 6. Mutator was designed to automatically retrieve mutations from literature and store these mutations in a 3DM system. A PubMed search using keywords from the 3DM system is used to preselect articles of interest. These articles are retrieved from the internet, converted to text, and parsed for mutations. Mutations are then grounded to proteins and stored in a 3DM database. Mutation retrieval was tested on the alpha-amylase superfamily as this superfamily contains the enzyme involved in Fabry’s disease: an x linked lysosomal storage disease. Compared to existing mutant databases, such as the HGMD and SwissProt, Mutator retrieved 30% more mutations from literature. A major problem in DNA diagnostics is the differentiation between natural variants and pathogenic mutations. To distinguish between pathogenic mutations and natural variation in proteins the Validator modules was added to 3DM. Validator uses the data available in a 3DM system to predict the pathogenicity of a mutant using, for example, the residue conservation of the mutants alignment position, side-chain accessibility of the mutant in the structure, and the number of mutations found in literature for the alignment position. Mutator and Validator can be used to study mutants found in disorder related genes. Although these tools are not the definitive solution for DNA diagnostics they can hopefully be used to increase our understanding of the molecular basis of disorders. Chapter 7 and 8 describe applied research projects using 3DM systems containg proteins of potential commercial interest. A 3DM system for the a/b-beta hydrolases superfamily is described in chapter 7. This superfamily consists of almost 20,000 proteins with a diverse range of functions. Superfamily alignments were generated for the common beta-barrel fold shared by all superfamily members, and for five distinct subtypes within the superfamily. Due to the size and functional diversity of the superfamily, there is a lot of potential for industrial application of superfamily members. Chapter 8 describes a study focusing on a sucrose phosphorylase enzyme from the a-amylase superfamily. This enzyme can be potentially used in an industrial setting for the transfer of glucose to a wide variety of molecules. The aim of the study was to increase the stability of the protein at higher temperatures. A combination of rational design using a 3DM system, and in-depth study of the protein structure, led to a series of mutations that resulted in more than doubling the half-life of the protein at 60°C. 3DM systems have been successfully applied in a wide range of protein engineering and DNA diagnostics studies. Currently, 3DM systems are applied most successfully in project studying a single protein family or monogenetic disorder. In the future, we hope to be able to apply 3DM to more complex scenarios such as enzyme factories and polygenetic disorders by combining multiple 3DM systems for interacting proteins.