The solution structure of C1-T1, a two-domain proteinase inhibitor derived from a circular precursor protein from Nicotiana alata

A two-domain portion of the proteinase inhibitor precursor from Nicotiana alata (NaProPI) has been expressed and its structure determined by NMR spectroscopy. NaProPI contains six almost identical 53 amino acid repeats that fold into six highly similar domains; however, the sequence repeats do nut coincide with the structural domains. Five of the structural domains comprise the C-terminal portion of one repeat and the N-terminal portion of the next. The sixth domain contains the C-terminal portion of the sixth repeat and the N-terminal portion of the first repeat. Disulphide bonds link these C and N-terminal fragments to generate the clasped-bracelet fold of NaProPI. The three-dimensional structure of NaProPI is not known, but it is conceivable that adjacent domains in NaProPI interact to generate the circular bracelet with the N and C termini in close enough proximity to facilitate formation of the disulphide bonds that form the clasp The expressed protein, examined in the current study, comprises residues 25-135 of NaProPI and encompasses the first two contiguous structural domains, namely the chymotrypsin inhibitor C1 and the trypsin inhibitor T1, joined by a five-residue linker, and is referred to as C1-T1. The tertiary structure of each domain in C1-T1 is identical to that found in the isolated inhibitors. However, no nuclear Overhauser effect contacts are observed between the two domains and the five-residue linker adopts an extended conformation. The absence of interactions between the domains indicates that adjacent domains do not specifically interact to drive the circularisation of NaProPI. These results are in agreement with recent data which describe similar PI precursors from other members of the Solanaceae having two, three, or four repeats. The lack of strong interdomain association is likely to be important for the function of individual inhibitors by ensuring that there is no masking of reactive sites upon release from the precursor. (C) 2001 Academic Press.

Three-dimensional structure of RTD-1, a cyclic antimicrobial defensin from rhesus macaque leukocytes

Most mammalian defensins are cationic peptides of 29-42 amino acids long, stabilized by three disulfide bonds. However, recently Tang et al. (1999, Science 286, 498-502) reported the isolation of a new defensin type found in the leukocytes of rhesus macaques. In contrast to all the other defensins found so far, rhesus theta defensin-1 (RTD-1) is composed of just 18 amino acids with the backbone cyclized through peptide bonds. Antibacterial activities of both the native cyclic peptide and a linear form were examined, showing that the cyclic form was 3-fold more active than the open chain analogue [Tang et al. (1999) Science 286, 498-502]. To elucidate the three-dimensional structure of RTD-1 and its open chain analogue, both peptides were synthesized using solid-phase peptide synthesis and tert-butyloxycarbonyl chemistry. The structures of both peptides in aqueous solution were determined from two-dimensional H-1 NMR data recorded at 500 and 750 MHz. Structural constraints consisting of interproton distances and dihedral angles were used as input for simulated-annealing calculations and water refinement with the program CNS. RTD-1 and its open chain analogue oRTD-1 adopt very similar structures in water. Both comprise an extended beta -hairpin structure with turns at one or both ends. The turns are well defined within themselves and seem to be flexible with respect to the extended regions of the molecules. Although the two strands of the beta -sheet are connected by three disulfide bonds, this region displays a degree of flexibility. The structural similarity of RTD-1 and its open chain analogue oRTD-1, as well as their comparable degree of flexibility, support the theory that the additional charges at the termini of the open chain analogue rather than overall differences in structure or flexibility are the cause for oRTD-1's lower antimicrobial activity. In contrast to numerous other antimicrobial peptides, RTD-1 does not display any amphiphilic character, even though surface models of RTD-1 exhibit a certain clustering of positive charges. Some amide protons of RTD-1 that should be solvent-exposed in monomeric beta -sheet structures show low-temperature coefficients, suggesting the possible presence of weak intermolecular hydrogen bonds.

Getting the adrenaline going: Crystal structure of the adrenaline-synthesizing enzyme PNMT

Background: Adrenaline is localized to specific regions of the central nervous system (CNS), but its role therein is unclear because of a lack of suitable pharmacologic agents. Ideally, a chemical is required that crosses the blood-brain barrier, potently inhibits the adrenaline-synthesizing enzyme PNMT, and does not affect other catecholamine processes. Currently available PNMT inhibitors do not meet these criteria. We aim to produce potent, selective, and CNS-active PNMT inhibitors by structure-based design methods. The first step is the structure determination of PNMT. Results: We have solved the crystal structure of human PNMT complexed with a cofactor product and a submicromolar inhibitor at a resolution of 2.4 Angstrom. The structure reveals a highly decorated methyltransferase fold, with an active site protected from solvent by an extensive cover formed from several discrete structural motifs. The structure of PNMT shows that the inhibitor interacts with the enzyme in a different mode from the (modeled) substrate noradrenaline. Specifically, the position and orientation of the amines is not equivalent. Conclusions: An unexpected finding is that the structure of PNMT provides independent evidence of both backward evolution and fold recruitment in the evolution of a complex enzyme from a simple fold. The proposed evolutionary pathway implies that adrenaline, the product of PNMT catalysis, is a relative newcomer in the catecholamine family. The PNMT structure reported here enables the design of potent and selective inhibitors with which to characterize the role of adrenaline in the CNS. Such chemical probes could potentially be useful as novel therapeutics.

Detection and description of non-linear interdependence in normal multichannel human EEG data

Objectives: This study examines human scalp electroencephalographic (EEG) data for evidence of non-linear interdependence between posterior channels. The spectral and phase properties of those epochs of EEG exhibiting non-linear interdependence are studied. Methods: Scalp EEG data was collected from 40 healthy subjects. A technique for the detection of non-linear interdependence was applied to 2.048 s segments of posterior bipolar electrode data. Amplitude-adjusted phase-randomized surrogate data was used to statistically determine which EEG epochs exhibited non-linear interdependence. Results: Statistically significant evidence of non-linear interactions were evident in 2.9% (eyes open) to 4.8% (eyes closed) of the epochs. In the eyes-open recordings, these epochs exhibited a peak in the spectral and cross-spectral density functions at about 10 Hz. Two types of EEG epochs are evident in the eyes-closed recordings; one type exhibits a peak in the spectral density and cross-spectrum at 8 Hz. The other type has increased spectral and cross-spectral power across faster frequencies. Epochs identified as exhibiting non-linear interdependence display a tendency towards phase interdependencies across and between a broad range of frequencies. Conclusions: Non-linear interdependence is detectable in a small number of multichannel EEG epochs, and makes a contribution to the alpha rhythm. Non-linear interdependence produces spatially distributed activity that exhibits phase synchronization between oscillations present at different frequencies. The possible physiological significance of these findings are discussed with reference to the dynamical properties of neural systems and the role of synchronous activity in the neocortex. (C) 2002 Elsevier Science Ireland Ltd. All rights reserved.

Air-water flows down stepped chutes: turbulence and flow structure observations

Interactions between turbulent waters and atmosphere may lead to strong air-water mixing. This experimental study is focused on the flow down a staircase channel characterised by very strong flow aeration and turbulence. Interfacial aeration is characterised by strong air-water mixing extending down to the invert. The size of entrained bubbles and droplets extends over several orders of magnitude, and a significant number of bubble/droplet clusters was observed. Velocity and turbulence intensity measurements suggest high levels of turbulence across the entire air-water flow. The increase in turbulence levels, compared to single-phase flow situations, is proportional to the number of entrained particles. (C) 2002 Elsevier Science Ltd. All rights reserved.

Social structure in migrating humpback whales (Megaptera novaeangliae)

Although largely solitary, humpback whales exhibit a number of behaviours where individuals co-operate with one another, for example during bubble net feeding. Such cases could be due to reciprocal altruism brought on by exceptional circumstances, for example the presence of abundant shoaling fish. An alternative explanation is that these behaviours have evolved through kin selection. With little restriction to either communication or movement, diffuse groups of relatives could maintain some form of social organization without the need to travel in tight-nit units. To try to distinguish between these hypotheses, we took advantage of the fact that migrating humpback whales often swim together in small groups. If kin selection is important in humpback whale biology, these groups should be enriched for relatives. Consequently, we analysed biopsy samples from 57 groups of humpback whales migrating off Eastern Australia in 1992. A total of 142 whales were screened for eight microsatellite markers. Mitochondrial DNA sequences (371 bp) were also used to verify and assist kinship identification. Our data add support to the notion that mothers travel with their offspring for the first year of the calf's life. However, beyond the presence of mother-calf/yearling pairs, no obvious relatedness pattern was found among whales sampled either in the same pod or on the same day. Levels of relatedness did not vary between migratory phases (towards or away from the breeding ground), nor between the two sexes considered either overall or in the north or south migrations separately. These findings suggest that, if any social organization does exist, it is formed transiently when needed rather than being a constant feature of the population, and hence is more likely based on reciprocal altruism than kin selection.

The Three-dimensional Solution Structure of NaD1, a New Floral Defensin from Nicotiana alata and its Application to a Homology Model of the Crop Defense Protein alfAFP

NMR spectroscopy and simulated annealing calculations have been used to determine the three-dimensional structure of NaD1, a novel antifungal and insecticidal protein isolated from the flowers of Nicotiana alata. NaD1 is a basic, cysteine-rich protein of 47 residues and is the first example of a plant defensin from flowers to be characterized structurally. Its three-dimensional structure consists of an a-helix and a triple-stranded anti-parallel beta-sheet that are stabilized by four intramolecular disulfide bonds. NaD1 features all the characteristics of the cysteine-stabilized up motif that has been described for a variety of proteins of differing functions ranging from antibacterial insect defensins and ion channel-perturbing scorpion toxins to an elicitor of the sweet taste response. The protein is biologically active against insect pests, which makes it a potential candidate for use in crop protection. NaD1 shares 31% sequence identity with alfAFP, an antifungal protein from alfalfa that confers resistance to a fungal pathogen in transgenic potatoes. The structure of NaD1 was used to obtain a homology model of alfAFP, since NaD1 has the highest level of sequence identity with alfAFP of any structurally characterized antifungal defensin. The structures of NaD1 and alfAFP were used in conjunction with structure - activity data for the radish defensin Rs-AFP2 to provide an insight into structure-function relationships. In particular, a putative effector site was identified in the structure of NaD1 and in the corresponding homology model of alfAFP. (C) 2002 Elsevier Science Ltd. All rights reserved.

Secondary structure of the third extracellular loop responsible for ligand selectivity of a mammalian gonadotropin-releasing hormone receptor

The extracellular loop 3 (ECL3) of the mammalian gonadotropin-releasing hormone receptor (GnRH-R) contains an acidic amino acid (Glu(301) in the mouse GnRH-R,) that confers agonist selectivity for Are in mammalian GnRH. It is proposed that a specific conformation of ECL3 is necessary to orientate the carboxyl side chain of the acidic residue for interaction with Arg(8) of GnRH, which is supported by decreased affinity for Arg(8) GnRH but not Gln(8) GnRH when an adjacent Pro is mutated to Ala. To probe the structural contribution of the loop domain to the proposed presentation of the carboxyl side chain, we synthesized a model peptide (CGPEMLNRVSEPGC) representing residues 293-302 of mouse ECL3, where Cys and Gly residues are added symmetrically at the N and C termini, respectively, allowing the introduction of a disulfide bridge to simulate the distances at which the ECL3 is tethered to the transmembrane domains 6 and 7 of the receptor. The ability of the ECL3 peptide to bind GnRH with low affinity was demonstrated by its inhibition of GnRH stimulation of inositol phosphate production in cells expressing the GnRH-R. The CD bands of the ECL3 peptides exhibited a superposition of predominantly unordered structure and partial contributions from beta-sheet structure. Likewise, the analysis of the amide I and amide III bands from micro-Raman and FT Raman experiments revealed mainly unordered conformations of the cyclic and of the linear peptide. NMR data demonstrated the presence of a beta-hairpin among an ensemble of largely disordered structures in the cyclic peptide. The location of the turn linking the two strands of the hairpin was assigned to the three central residues L-296, N-297, and R-298. A small population of structured species among an ensemble of predominantly random coil conformation suggests that the unliganded receptor represents a variety of structural conformers, some of which have the potential to make contacts with the ligand. We propose a mechanism of receptor activation whereby binding of the agonist to the inactive receptor state induces and stabilizes a particular structural state of the loop domain, leading to further conformational rearrangements across the transmembrane domain and signal propagating interaction with G proteins. Interaction of the Glu(301) of the receptor with Arg(8) of GnRH induces a folded configuration of the ligand. Our proposal thus suggests that conformational changes of both ligand and receptor result from this interaction.

Modelling High-Dimensional Data by Mixtures of Factor Analyzers

We focus on mixtures of factor analyzers from the perspective of a method for model-based density estimation from high-dimensional data, and hence for the clustering of such data. This approach enables a normal mixture model to be fitted to a sample of n data points of dimension p, where p is large relative to n. The number of free parameters is controlled through the dimension of the latent factor space. By working in this reduced space, it allows a model for each component-covariance matrix with complexity lying between that of the isotropic and full covariance structure models. We shall illustrate the use of mixtures of factor analyzers in a practical example that considers the clustering of cell lines on the basis of gene expressions from microarray experiments. (C) 2002 Elsevier Science B.V. All rights reserved.

Systematic characterization of mutations in yeast acetohydroxyacid synthase: Interpretation of herbicide-resistance data

Acetohydroxyacid synthase (AHAS, EC 4.1.3.18) catalyses the first step in branched-chain amino acid biosynthesis and is the target for sulfonylurea and imidazolinone herbicides, which act as potent and specific inhibitors. Mutants of the enzyme have been identified that are resistant to particular herbicides. However, the selectivity of these mutants towards various sulfonylureas and imidazolinones has not been determined systematically. Now that the structure of the yeast enzyme is known, both in the absence and presence of a bound herbicide, a detailed understanding of the molecular interactions between the enzyme and its inhibitors becomes possible. Here we construct 10 active mutants of yeast AHAS, purify the enzymes and determine their sensitivity to six sulfonylureas and three imidazolinones. An additional three active mutants were constructed with a view to increasing imidazolinone sensitivity. These three variants were purified and tested for their sensitivity to the imidazolinones only. Substantial differences are observed in the sensitivity of the 13 mutants to the various inhibitors and these differences are interpreted in terms of the structure of the herbicide-binding site on the enzyme.

Numerical simulations of wind and temperature structure within an Alpine lake basin, Lake Tekapo, New Zealand

High-resolution numerical model simulations have been used to study the local and mesoscale thermal circulations in an Alpine lake basin. The lake (87 km(2)) is situated in the Southern Alps, New Zealand and is located in a glacially excavated rock basin surrounded by mountain ranges that reach 3000 m in height. The mesoscale model used (RAMS) is a three-dimensional non-hydrostatic model with a level 2.5 turbulence closure scheme. The model demonstrates that thermal forcing at local (within the basin) and regional (coast-to-basin inflow) scales drive the observed boundary-layer airflow in the lake basin during clear anticyclonic summertime conditions. The results show that the lake can modify (perturb) both the local and regional wind systems. Following sunrise, local thermal circulations dominate, including a lake breeze component that becomes embedded within the background valley wind system. This results in a more divergent flow in the basin extending across the lake shoreline. However, a closed lake breeze circulation is neither observed nor modelled. Modelling results indicate that in the latter part of the day when the mesoscale (coast-to-basin) inflow occurs, the relatively cold pool of lake air in the basin can cause the intrusion to decouple from the surface. Measured data provide qualitative and quantitative support for the model results.

The pyramidal neuron in occipital, temporal and prefrontal cortex of the owl monkey (Aotus trivirgatus): regional specialization in cell structure

Recent studies have revealed marked regional variation in pyramidal cell morphology in primate cortex. In particular, pyramidal cells in human and macaque prefrontal cortex (PFC) are considerably more spinous than those in other cortical regions. PFC pyramidal cells in the New World marmoset monkey, however, are less spinous than those in man and macaques. Taken together, these data suggest that the pyramidal cell has become more branched and more spinous during the evolution of PFC in only some primate lineages. This specialization may be of fundamental importance in determining the cognitive styles of the different species. However, these data are preliminary, with only one New World and two Old World species having been studied. Moreover, the marmoset data were obtained from different cases. In the present study we investigated PFC pyramidal cells in another New World monkey, the owl monkey, to extend the basis for comparison. As in the New World marmoset monkey, prefrontal pyramidal cells in owl monkeys have relatively few spines. These species differences appear to reflect variation in the extent to which PFC circuitry has become specialized during evolution. Highly complex pyramidal cells in PFC appear not to have been a feature of a common prosimian ancestor, but have evolved with the dramatic expansion of PFC in some anthropoid lineages.

Mechanical Behavior of the Lamellar Structure in Semi-Crystalline Polymers

We have employed molecular dynamics simulations to study the behavior of virtual polymeric materials under an applied uniaxial tensile load. Through computer simulations, one can obtain experimentally inaccessible information about phenomena taking place at the molecular and microscopic levels. Not only can the global material response be monitored and characterized along time, but the response of macromolecular chains can be followed independently if desired. The computer-generated materials were created by emulating the step-wise polymerization, resulting in self-avoiding chains in 3D with controlled degree of orientation along a certain axis. These materials represent a simplified model of the lamellar structure of semi-crystalline polymers,being comprised of an amorphous region surrounded by two crystalline lamellar regions. For the simulations, a series of materials were created, varying i) the lamella thickness, ii) the amorphous region thickness, iii) the preferential chain orientation, and iv) the degree of packing of the amorphous region. Simulation results indicate that the lamella thickness has the strongest influence on the mechanical properties of the lamella-amorphous structure, which is in agreement with experimental data. The other morphological parameters also affect the mechanical response, but to a smaller degree. This research follows previous simulation work on the crack formation and propagation phenomena, deformation mechanisms at the nanoscale, and the influence of the loading conditions on the material response. Computer simulations can improve the fundamental understanding about the phenomena responsible for the behavior of polymeric materials, and will eventually lead to the design of knowledge-based materials with improved properties.

Semi-automatic 3D Segmentation Of Costal Cartilage In CT Data From Pectus Excavatum Patients

One of the current frontiers in the clinical management of Pectus Excavatum (PE) patients is the prediction of the surgical outcome prior to the intervention. This can be done through computerized simulation of the Nuss procedure, which requires an anatomically correct representation of the costal cartilage. To this end, we take advantage of the costal cartilage tubular structure to detect it through multi-scale vesselness filtering. This information is then used in an interactive 2D initialization procedure which uses anatomical maximum intensity projections of 3D vesselness feature images to efficiently initialize the 3D segmentation process. We identify the cartilage tissue centerlines in these projected 2D images using a livewire approach. We finally refine the 3D cartilage surface through region-based sparse field level-sets. We have tested the proposed algorithm in 6 noncontrast CT datasets from PE patients. A good segmentation performance was found against reference manual contouring, with an average Dice coefficient of 0.75±0.04 and an average mean surface distance of 1.69±0.30mm. The proposed method requires roughly 1 minute for the interactive initialization step, which can positively contribute to an extended use of this tool in clinical practice, since current manual delineation of the costal cartilage can take up to an hour.

Grid structure impact in sparse point representation of derivatives

In the Sparse Point Representation (SPR) method the principle is to retain the function data indicated by significant interpolatory wavelet coefficients, which are defined as interpolation errors by means of an interpolating subdivision scheme. Typically, a SPR grid is coarse in smooth regions, and refined close to irregularities. Furthermore, the computation of partial derivatives of a function from the information of its SPR content is performed in two steps. The first one is a refinement procedure to extend the SPR by the inclusion of new interpolated point values in a security zone. Then, for points in the refined grid, such derivatives are approximated by uniform finite differences, using a step size proportional to each point local scale. If required neighboring stencils are not present in the grid, the corresponding missing point values are approximated from coarser scales using the interpolating subdivision scheme. Using the cubic interpolation subdivision scheme, we demonstrate that such adaptive finite differences can be formulated in terms of a collocation scheme based on the wavelet expansion associated to the SPR. For this purpose, we prove some results concerning the local behavior of such wavelet reconstruction operators, which stand for SPR grids having appropriate structures. This statement implies that the adaptive finite difference scheme and the one using the step size of the finest level produce the same result at SPR grid points. Consequently, in addition to the refinement strategy, our analysis indicates that some care must be taken concerning the grid structure, in order to keep the truncation error under a certain accuracy limit. Illustrating results are presented for 2D Maxwell's equation numerical solutions.