Circular proteins in plants - Solution structure of a novel macrocyclic trypsin inhibitor from Momordica cochinchinensis

Much interest has been generated by recent reports on the discovery of circular (i.e. head-to-tail cyclized) proteins in plants. Here we report the three-dimensional structure of one of the newest such circular proteins, MCoTI-II, a novel trypsin inhibitor from Momordica cochinchinensis, a member of the Cucurbitaceae plant family. The structure consists of a small beta -sheet, several turns, and a cystine knot arrangement of the three disulfide bonds. Interestingly, the molecular topology is similar to that of the plant cyclotides (Craik, D. J., Daly, N. L., Bond, T., and Waine, C. (1999) J. Mol. Biol, 294, 1327-1336), which derive from the Rubiaceae and Violaceae plant families, have antimicrobial activities, and exemplify the cyclic cystine knot structural motif as part of their circular backbone. The sequence, biological activity, and plant family of MCoTI-II are all different from known cyclotides. However, given the structural similarity, cyclic backbone, and plant origin of MCoTI-II, we propose that MCoTI-II can be classified as a new member of the cyclotide class of proteins. The expansion of the cyclotides to include trypsin inhibitory activity and a new plant family highlights the importance and functional variability of circular proteins and the fact that they are more common than has previously been believed, Insights into the possible roles of backbone cyclization have been gained by a comparison of the structure of MCoTI-II with the homologous acyclic trypsin inhibitors CMTI-I and EETI-II from the Cucurbitaceae plant family.

Comments on "Exact solution for horizontal redistribution by general similarity"

No abstract

Social structures, genetic structures and dispersal strategies in Australian rabbit (Oryctolagus cuniculus) populations

In this study, the pattern of movement of young male and female rabbits and the genetic structures present in adult male and female populations in four habitats was examined. The level of philopatry in young animals was found to vary between 18-90% for males and 32-95% for females in different populations. It was skewed, with more males dispersing than females in some but not all populations. Analysis of allozyme data using spatial autocorrelation showed that adult females from the same social group, unlike males, were significantly related in four of the five populations studied. Changes in genetic structure and rate of dispersal were measured before and during the recovery of a population that was artificially reduced in size. There were changes in the rate and distance of dispersal with density and sex. Subadults of both sexes moved further in the first year post crash (low density) than in the following years. While the level of dispersal for females was lower than that of the males for the first 3 years, thereafter (high density) both sexes showed similar, low levels of dispersal (20%). The density at which young animals switch behaviour between dispersal and philopatry differed for males and females. The level of genetic structuring in adult females was high in the precrash population, reduced in the first year post crash and undetectable in the second year. Dispersal behaviour of rabbits both affects the genetic structure of the population and changes with conditions. Over a wide range of levels of philopatry, genetic structuring is present in the adult female, but not the male population. Consequently, though genetic structuring is present, it does not lead to inbreeding. More long-distance movements are found in low-density populations, even though vacant warrens are available near birth warrens. The distances moved decreased as density increased. Calculation of the effective population size (N-e) shows that changes in dispersal distance offset changes in density, so that N-e remains constant.

Analytical solution to the zero-inertia problem for surge flow phenomena in nonprismatic channels

Surge flow phenomena. e.g.. as a consequence of a dam failure or a flash flood, represent free boundary problems. ne extending computational domain together with the discontinuities involved renders their numerical solution a cumbersome procedure. This contribution proposes an analytical solution to the problem, It is based on the slightly modified zero-inertia (ZI) differential equations for nonprismatic channels and uses exclusively physical parameters. Employing the concept of a momentum-representative cross section of the moving water body together with a specific relationship for describing the cross sectional geometry leads, after considerable mathematical calculus. to the analytical solution. The hydrodynamic analytical model is free of numerical troubles, easy to run, computationally efficient. and fully satisfies the law of volume conservation. In a first test series, the hydrodynamic analytical ZI model compares very favorably with a full hydrodynamic numerical model in respect to published results of surge flow simulations in different types of prismatic channels. In order to extend these considerations to natural rivers, the accuracy of the analytical model in describing an irregular cross section is investigated and tested successfully. A sensitivity and error analysis reveals the important impact of the hydraulic radius on the velocity of the surge, and this underlines the importance of an adequate description of the topography, The new approach is finally applied to simulate a surge propagating down the irregularly shaped Isar Valley in the Bavarian Alps after a hypothetical dam failure. The straightforward and fully stable computation of the flood hydrograph along the Isar Valley clearly reflects the impact of the strongly varying topographic characteristics on the How phenomenon. Apart from treating surge flow phenomena as a whole, the analytical solution also offers a rigorous alternative to both (a) the approximate Whitham solution, for generating initial values, and (b) the rough volume balance techniques used to model the wave tip in numerical surge flow computations.

Analytical solution for tidal propagation in a coupled semi-confined/phreatic coastal aquifer

An analytical solution is derived for tidal fluctuations in a coupled coastal aquifer system consisting of a semi-confined aquifer, a thin semi-permeable layer and a phreatic aquifer. Based on the solution, we study the interactions (via leakage) between the confined and unconfined aquifers in response to tides. The results show that, under certain conditions, leakage from the confined aquifer can affect considerably the tidal water table fluctuation in the phreatic aquifer and vice versa. Ignoring these effects could lead to errors in estimating aquifer properties based on tidal signals. (C) 2002 Elsevier Science Ltd. All rights reserved.

A director theory for visco-elastic folding instabilities in multilayered rock

A model for finely layered visco-elastic rock proposed by us in previous papers is revisited and generalized to include couple stresses. We begin with an outline of the governing equations for the standard continuum case and apply a computational simulation scheme suitable for problems involving very large deformations. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered beam under compression. We analyse folding up to 40% shortening. The standard continuum solution becomes unstable for extreme values of the shear/normal viscosity ratio. The instability is a consequence of the neglect of the bending stiffness/viscosity in the standard continuum model. We suggest considering these effects within the framework of a couple stress theory. Couple stress theories involve second order spatial derivatives of the velocities/displacements in the virtual work principle. To avoid C-1 continuity in the finite element formulation we introduce the spin of the cross sections of the individual layers as an independent variable and enforce equality to the spin of the unit normal vector to the layers (-the director of the layer system-) by means of a penalty method. We illustrate the convergence of the penalty method by means of numerical solutions of simple shears of an infinite layer for increasing values of the penalty parameter. For the shear problem we present solutions assuming that the internal layering is oriented orthogonal to the surfaces of the shear layer initially. For high values of the ratio of the normal-to the shear viscosity the deformation concentrates in thin bands around to the layer surfaces. The effect of couple stresses on the evolution of folds in layered structures is also investigated. (C) 2002 Elsevier Science Ltd. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Estimation of crystalline phase present in the glucose crystal-solution mixture by water activity measurement

The amount of crystalline fraction present in monohydrate glucose crystal-solution mixture up to 110% crystal in relation to solution (crystal:solution=110:100) was determined by water activity measurement. It was found that the water activity had a strong linear correlation (R-2=0.994) with the amount of glucose present above saturation. Difference in the water activities of the crystal-solution mixture (a(w1)) and the supersaturated solution (a(w2)) by re-dissolving the crystalline fraction allowed calculation of the amount of crystalline phase present (DeltaG) in the mixture by an equation DeltaG=846.97(a(w1)-a(w2)). Other methods such as Raoult's, Norrish and Money-Born equations were also tested for the prediction of water activity of supersaturated glucose solution. (C) 2003 Swiss Society of Food Science and Technology. Published by Elsevier Science Ltd. All rights reserved.

A comparison of the self-association behavior of the plant cyclotides kalata B1 and kalata B2 via analytical ultracentrifugation

The recently discovered cyclotides kalata B1 and kalata B2 are miniproteins containing a head-to-tail cyclized backbone and a cystine knot motif, in which disulfide bonds and the connecting backbone segments form a ring that is penetrated by the third disulfide bond. This arrangement renders the cyclotides extremely stable against thermal and enzymatic decay, making them a possible template onto which functionalities can be grafted.We have compared the hydrodynamic properties of two prototypic cyclotides, kalata B1 and kalata B2, using analytical ultracentrifugation techniques. Direct evidence for oligomerization of kalata B2 was shown by sedimentation velocity experiments in which a method for determining size distribution of polydisperse molecules in solution was employed. The shape of the oligomers appears to be spherical. Both sedimentation velocity and equilibrium experiments indicate that in phosphate buffer kalata B1 exists mainly as a monomer, even at millimolar concentrations. In contrast, at 1.6 mM, kalata B2 exists as an equilibrium mixture of monomer (30%), tetramer (42%), octamer (25%), and possibly a small proportion of higher oligomers. The results from the sedimentation equilibrium experiments show that this self-association is concentration dependent and reversible. We link our findings to the three-dimensional structures of both cyclotides, and propose two putative interaction interfaces on opposite sides of the kalata B2 molecule, one involving a hydrophobic interaction with the Phe(6), and the second involving a charge-charge interaction with the Asp(25) residue. An understanding of the factors affecting solution aggregation is of vital importance for future pharmaceutical application of these molecules.

Seedling responses of four Australian tree species to toxic concentrations of manganese in solution culture

Little is known about the responses of Australian plants to excess metal, including Mn. It is important to remedy this lack of information so that knowledgeable decisions can be made about managing Mn contaminated sites where inhabited by Australian vegetation. Acacia holosericea, Melaleuca leucadendra, Eucalyptus crebra and Eucalyptus camaldulensis were grown in dilute solution culture for 10 weeks. The seedlings ( 42 days old) were exposed to six Mn treatments viz., 1, 8, 32, 128, 512 and 2048 muM. The order of tolerance to toxic concentrations of Mn was A. holosericea congruent to = E. crebra < M. leucadendra < E. camaldulensis, the critical external concentrations being approximately 5.1, 5.0, 21 and 330 muM, respectively. The critical tissue Mn concentrations for the youngest fully expanded leaf and total shoots were, respectively, 265 and 215 mug g(-1) DM for A. holosericea, 445 and 495 mug g(-1) DM for M. leucadendra, 495 and 710 mug g(-1) DM for E. crebra and 7230 and 6510 mug g(-1) DM for E. camaldulensis. The high tolerance of E. camaldulensis ( as opposed to the sensitivity of E. crebra) to excess Mn raises concern about fauna feeding on the plant and is consistent with hypotheses suggesting the Eucalyptus subgenus Symphomyrtus is particularly tolerant of stress, including excess Mn. The results from this paper provide the first comprehensive combination of growth responses, critical external concentrations, critical tissue concentrations and plant toxicity symptoms for three important Australian genera, viz., Eucalyptus, Acacia and Melaleuca, for use in the management of Mn toxic sites.

Stability property of matchings is a natural solution concept in coalitional market games

Stability of matchings was proved to be a new cooperative equilibrium concept in Sotomayor (Dynamics and equilibrium: essays in honor to D. Gale, 1992). That paper introduces the innovation of treating as multi-dimensional the payoff of a player with a quota greater than one. This is done for the many-to-many matching model with additively separable utilities, for which the stability concept is defined. It is then proved, via linear programming, that the set of stable outcomes is nonempty and it may be strictly bigger than the set of dual solutions and strictly smaller than the core. The present paper defines a general concept of stability and shows that this concept is a natural solution concept, stronger than the core concept, for a much more general coalitional game than a matching game. Instead of mutual agreements inside partnerships, the players are allowed to make collective agreements inside coalitions of any size and to distribute his labor among them. A collective agreement determines the level of labor at which the coalition operates and the division, among its members, of the income generated by the coalition. An allocation specifies a set of collective agreements for each player.