Assessment of the ability to model proteins with leucine-rich repeats in light of the latest structural information

The three-dimensional structures of leucine-rich repeat (LRR) -containing proteins from five different families were previously predicted based on the crystal structure of the ribonuclease inhibitor. using an approach that combined homology-based modeling, structure-based sequence alignment of LRRs, and several rational assumptions. The structural models have been produced based on very limited sequence similarity, which, in general. cannot yield trustworthy predictions. Recently, the protein structures from three of these five families have been determined. In this report we estimate the quality of the modeling approach by comparing the models with the experimentally determined structures. The comparison suggests that the general architecture, curvature, interior/exterior orientations of side chains. and backbone conformation of the LRR structures can be predicted correctly. On the other hand. the analysis revealed that, in some cases. it is difficult to predict correctly the twist of the overall super-helical structure. Taking into consideration the conclusions from these comparisons, we identified a new family of bacterial LRR proteins and present its structural model. The reliability of the LRR protein modeling suggests that it would be informative to apply similar modeling approaches to other classes of solenoid proteins.

Vibration instability in rolling mills: Modeling and experimental results

Models for the occurrence of the vibrational instability during rolling known as third octave chatter are presented and discussed. An analysis of rolling mill chatter was performed for the purpose of identifying characteristics of the vibrations and to determine any dependency on the rolling schedule. In particular, a stability criterion for the critical rolling speed is used to predict the maximum rolling speed without chatter instability on schedules from a 5 stand tandem mill rolling thin steel product. The results correlate well with measurements of critical speed occurring on the mill using a vibration monitor: This research provides significant insights into the chatter phenomena and has been used to investigate control methods for suppression of the instability.

Empirical evidence for ultrametric structure in multi-layer perceptron error surfaces

Combinatorial optimization problems share an interesting property with spin glass systems in that their state spaces can exhibit ultrametric structure. We use sampling methods to analyse the error surfaces of feedforward multi-layer perceptron neural networks learning encoder problems. The third order statistics of these points of attraction are examined and found to be arranged in a highly ultrametric way. This is a unique result for a finite, continuous parameter space. The implications of this result are discussed.

Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda)

ITS2 sequences are used extensively in molecular taxonomy and population genetics of arthropods and other animals yet little is known about the molecular evolution of ITS2. We studied the secondary structure of ITS2 in species from each of the six main lineages of hard ticks (family Ixodidae). The ITS2 of these ticks varied in length from 679 bp in Ixodes scapularis to 1547 bp in Aponomma concolor. Nucleotide content varied also: the ITS2 of ticks from the Prostriata lineage (Ixodes spp.) had 46-49% GC whereas ITS2 sequences of ticks from the Metastriata lineage (all other hard ticks) had 61-62% GC. Despite variation in nucleotide sequence, the secondary structure of the ITS2 of all of these ticks apparently has five domains. Stems 1, 3, 4 and 5 of this secondary structure were obvious in all of the species studied. However, stem 2 was not always obvious despite the fact that it is flanked by highly conserved sequence motifs in the adjacent stems, stems 1 and 3. The ITS2 of hard ticks has apparently evolved mostly by increases and decreases in length of the nucleotide sequences, which caused increases, and decreases in the length of stems of the secondary structure. This is most obvious when stems of the secondary structures of the Prostriata (Ixodes spp.) are compared to those of the Metastriata (all other hard ticks). Increases in the size of the ITS2 may have been caused by replication slippage which generated large repeats, like those seen in Haemaphysalis humerosa and species from the Rhipicepalinae lineage, and the small repeats found in species from the other lineages of ticks.

Molecular analysis of dimethyl sulphide dehydrogenase from Rhodovulum sulfidophilum: its place in the dimethyl sulphoxide reductase family of microbial molybdopterin-containing enzymes

Dimethyl sulphide dehydrogenase catalyses the oxidation of dimethyl sulphide to dimethyl sulphoxide (DMSO) during photoautotrophic growth of Rhodovulum sulfidophilum . Dimethyl sulphide dehydrogenase was shown to contain bis (molybdopterin guanine dinucleotide)Mo, the form of the pterin molybdenum cofactor unique to enzymes of the DMSO reductase family. Sequence analysis of the ddh gene cluster showed that the ddhA gene encodes a polypeptide with highest sequence similarity to the molybdop-terin-containing subunits of selenate reductase, ethylbenzene dehydrogenase. These polypeptides form a distinct clade within the DMSO reductase family. Further sequence analysis of the ddh gene cluster identified three genes, ddhB , ddhD and ddhC . DdhB showed sequence homology to NarH, suggesting that it contains multiple iron-sulphur clusters. Analysis of the N-terminal signal sequence of DdhA suggests that it is secreted via the Tat secretory system in complex with DdhB, whereas DdhC is probably secreted via a Sec-dependent mechanism. Analysis of a ddhA mutant showed that dimethyl sulphide dehydrogenase was essential for photolithotrophic growth of Rv. sulfidophilum on dimethyl sulphide but not for chemo-trophic growth on the same substrate. Mutational analysis showed that cytochrome c (2) mediated photosynthetic electron transfer from dimethyl sulphide dehydrogenase to the photochemical reaction centre, although this cytochrome was not essential for photoheterotrophic growth of the bacterium.

Analysis of genetic diversity within Australian lucerne cultivars and implications for future genetic improvement

Lucerne (Medicago sativa L.) is autotetraploid, and predominantly allogamous. This complex breeding structure maximises the genetic diversity within lucerne populations making it difficult to genetically discriminate between populations. The objective of this study was to evaluate the level of random genetic diversity within and between a selection of Australian-grown lucerne cultivars, with tetraploid M. falcata included as a possible divergent control source. This diversity was evaluated using random amplified polymorphic DNA (RAPDs). Nineteen plants from each of 10 cultivars were analysed. Using 11 RAPD primers, 96 polymorphic bands were scored as present or absent across the 190 individuals. Genetic similarity estimates (GSEs) of all pair-wise comparisons were calculated from these data. Mean GSEs within cultivars ranged from 0.43 to 0.51. Cultivar Venus (0.43) had the highest level of intra-population genetic diversity and cultivar Sequel HR (0.51) had the lowest level of intra-population genetic diversity. Mean GSEs between cultivars ranged from 0.31 to 0.49, which overlapped with values obtained for within-cultivar GSE, thus not allowing separation of the cultivars. The high level of intra- and inter-population diversity that was detected is most likely due to the breeding of synthetic cultivars using parents derived from a number of diverse sources. Cultivar-specific polymorphisms were only identified in the M. falcata source, which like M. sativa, is outcrossing and autotetraploid. From a cluster analysis and a principal components analysis, it was clear that M. falcata was distinct from the other cultivars. The results indicate that the M. falcata accession tested has not been widely used in Australian lucerne breeding programs, and offers a means of introducing new genetic diversity into the lucerne gene pool. This provides a means of maximising heterozygosity, which is essential to maximising productivity in lucerne.

Analysis of the wind field and heat budget in an alpine lake basin during summertime fair weather conditions

Observational data collected in the Lake Tekapo hydro catchment of the Southern Alps in New Zealand are used to analyse the wind and temperature fields in the alpine lake basin during summertime fair weather conditions. Measurements from surface stations, pilot balloon and tethersonde soundings, Doppler sodar and an instrumented light aircraft provide evidence of multi-scale interacting wind systems, ranging from microscale slope winds to mesoscale coast-to-basin flows. Thermal forcing of the winds occurred due to differential heating as a consequence of orography and heterogeneous surface features, which is quantified by heat budget and pressure field analysis. The daytime vertical temperature structure was characterised by distinct layering. Features of particular interest are the formation of thermal internal boundary layers due to the lake-land discontinuity and the development of elevated mixed layers. The latter were generated by advective heating from the basin and valley sidewalls by slope winds and by a superimposed valley wind blowing from the basin over Lake Tekapo and up the tributary Godley Valley. Daytime heating in the basin and its tributary valleys caused the development of a strong horizontal temperature gradient between the basin atmosphere and that over the surrounding landscape, and hence the development of a mesoscale heat low over the basin. After noon, air from outside the basin started flowing over mountain saddles into the basin causing cooling in the lowest layers, whereas at ridge top height the horizontal air temperature gradient between inside and outside the basin continued to increase. In the early evening, a more massive intrusion of cold air caused rapid cooling and a transition to a rather uniform slightly stable stratification up to about 2000 m agl. The onset time of this rapid cooling varied about 1-2 h between observation sites and was probably triggered by the decay of up-slope winds inside the basin, which previously countered the intrusion of air over the surrounding ridges. The intrusion of air from outside the basin continued until about mid-night, when a northerly mountain wind from the Godley Valley became dominant. The results illustrate the extreme complexity that can be caused by the operation of thermal forcing processes at a wide range of spatial scales.

Modeling of phase equilibrium and vapor adsorption on carbon black based on a combination of a lattice theory and equation of state

A thermodynamic approach is developed in this paper to describe the behavior of a subcritical fluid in the neighborhood of vapor-liquid interface and close to a graphite surface. The fluid is modeled as a system of parallel molecular layers. The Helmholtz free energy of the fluid is expressed as the sum of the intrinsic Helmholtz free energies of separate layers and the potential energy of their mutual interactions calculated by the 10-4 potential. This Helmholtz free energy is described by an equation of state (such as the Bender or Peng-Robinson equation), which allows us a convenient means to obtain the intrinsic Helmholtz free energy of each molecular layer as a function of its two-dimensional density. All molecular layers of the bulk fluid are in mechanical equilibrium corresponding to the minimum of the total potential energy. In the case of adsorption the external potential exerted by the graphite layers is added to the free energy. The state of the interface zone between the liquid and the vapor phases or the state of the adsorbed phase is determined by the minimum of the grand potential. In the case of phase equilibrium the approach leads to the distribution of density and pressure over the transition zone. The interrelation between the collision diameter and the potential well depth was determined by the surface tension. It was shown that the distance between neighboring molecular layers substantially changes in the vapor-liquid transition zone and in the adsorbed phase with loading. The approach is considered in this paper for the case of adsorption of argon and nitrogen on carbon black. In both cases an excellent agreement with the experimental data was achieved without additional assumptions and fitting parameters, except for the fluid-solid potential well depth. The approach has far-reaching consequences and can be readily extended to the model of adsorption in slit pores of carbonaceous materials and to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Modeling of gas adsorption equilibrium over a wide range of pressure: A thermodynamic approach based on equation of state

A thermodynamic approach based on the Bender equation of state is suggested for the analysis of supercritical gas adsorption on activated carbons at high pressure. The approach accounts for the equality of the chemical potential in the adsorbed phase and that in the corresponding bulk phase and the distribution of elements of the adsorption volume (EAV) over the potential energy for gas-solid interaction. This scheme is extended to subcritical fluid adsorption and takes into account the phase transition in EAV The method is adapted to gravimetric measurements of mass excess adsorption and has been applied to the adsorption of argon, nitrogen, methane, ethane, carbon dioxide, and helium on activated carbon Norit R I in the temperature range from 25 to 70 C. The distribution function of adsorption volume elements over potentials exhibits overlapping peaks and is consistently reproduced for different gases. It was found that the distribution function changes weakly with temperature, which was confirmed by its comparison with the distribution function obtained by the same method using nitrogen adsorption isotherm at 77 K. It was shown that parameters such as pore volume and skeleton density can be determined directly from adsorption measurements, while the conventional approach of helium expansion at room temperature can lead to erroneous results due to the adsorption of helium in small pores of activated carbon. The approach is a convenient tool for analysis and correlation of excess adsorption isotherms over a wide range of pressure and temperature. This approach can be readily extended to the analysis of multicomponent adsorption systems. (C) 2002 Elsevier Science (USA).

Synthesis and structural analysis of the N-terminal domain of the thyroid hormone-binding protein transthyretin

Transthyretin (TTR) is a 55 kDa protein responsible for the transport of thyroid hormones and retinol in human serum. Misfolded forms of the protein are implicated in the amyloid diseases familial amyloidotic polyneuropathy and senile systemic amyloidosis. Its folding properties and stabilization by ligands are of current interest due to their importance in understanding and combating these diseases. To assist in such studies we developed a method for the solid phase synthesis of the monomeric unit of a TTR analogue and its folding to form a functional 55 kDa tetramer. The monomeric unit of the protein was chemically synthesized in three parts, comprising amino acid residues 151, 5499 and 102127, and ligated using chemoselective thioether ligation chemistry. The synthetic protein was folded and assembled to a tetrameric structure in the presence of the TTRs native ligand, thyroxine, as shown by gel filtration chromatography, native gel electrophoresis, TTR antibody recognition and thyroid hormone binding. In the current study the solution structure of the first of these fragment peptides, TTR(151) is examined to determine its intrinsic propensity to form beta-sheet structure, potentially involved in amyloid fibril formation by TTR. Despite the presence of extensive beta-structure in the native form of the protein, the Nterminal fragment adopts an essentially random coil conformation in solution.

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.

Linkage and Association Analysis of Radiation Damage Repair Genes XRCC3 and XRCC5 with Nevus Density in Adolescent Twins

Previous studies have shown that a deficiency in DNA damage repair is associated with increased cancer risk, and exposure to UV radiation is a major risk factor for the development of malignant melanoma. High density of common nevi (moles) is a major risk factor for cutaneous melanoma. A nevus may result from a mutation in a single UV-exposed melanocyte which failed to repair DNA damage in one or more critical genes. XRCC3 and XRCC5 may have an effect on nevus count through their function as components of DNA repair processes that may be involved directly or indirectly in the repair of DNA damage due to UV radiation. This study aims to test the hypothesis that the frequency of flat or raised moles is associated with polymorphism at or near these DNA repair genes, and that certain alleles are associated with less efficient DNA repair, and greater nevus density. Twins were recruited from schools in south eastern Queensland and were examined close to their 12th birthday. Nurses examined each individual and counted all moles on the entire body surface. A 10cM genome scan of 274 families (642 individuals) was performed and microsatellite polymorphisms in XRCC3 and adjacent to XRCC5 were also typed. Linkage and association of nevus count to these loci were tested simultaneously using a structural-equation modeling approach implemented in MX. There is weak evidence for linkage of XRCC5 to a QTL influencing raised mole count, and also weak association. There is also weak evidence for association between flat mole count and XRCC3. No tests were significant after correction for testing multiple alleles, nor were any of the tests for total association significant. If variation in XRCC3 or XRCC5 influences UV sensitivity, and indirectly affects nevus density, then the effects are small.

Direction of Causation Modeling Between Cross-Sectional Measures of Parenting and Psychological Distress in Female Twins

Under certain conditions, cross-sectional analysis of cross-twin intertrait correlations can provide important information about the direction of causation (DOC) between two variables. A community-based sample of Australian female twins aged 18 to 45 years was mailed an extensive Health and Lifestyle Questionnaire (HLQ) that covered a wide range of personality and behavioral measures. Included were self-report measures of recent psychological distress and perceived childhood environment (PBI). Factor analysis of the PBI yielded three interpretable dimensions: Coldness, Overprotection, and Autonomy. Univariate analysis revealed that parental Overprotection and Autonomy were best explained by additive genetic, shared, and nonshared environmental effects (ACE), whereas the best-fitting model for PBI Coldness and the three measures of psychological distress (Depression, Phobic Anxiety, and Somatic Distress) included only additive genetic and nonshared environmental effects (AE). A common pathway model best explained the covariation between (1) the three PBI dimensions and (2) the three measures of psychological distress. DOC modeling between latent constructs of parenting and psychological distress revealed that a model which specified recollected parental behavior as the cause of psychological distress provided a better fit than a model which specified psychological distress as the cause of recollected parental behavior. Power analyses and limitations of the findings are discussed.

NMR analysis of covalent intermediates in thiamin diphosphate enzymes

Enzymic catalysis proceeds via intermediates formed in the course of substrate conversion. Here, we directly detect key intermediates in thiamin diphosphate (ThDP)-dependent enzymes during catalysis using H-1 NMR spectroscopy. The quantitative analysis of the relative intermediate concentrations allows the determination of the microscopic rate constants of individual catalytic steps. As demonstrated for pyruvate decarboxylase (PDC), this method, in combination with site-directed mutagenesis, enables the assignment of individual side chains to single steps in catalysis. In PDC, two independent proton relay systems and the stereochemical control of the enzymic environment account for proficient catalysis proceeding via intermediates at carbon 2 of the enzyme-bound cofactor. The application of this method to other ThDP-dependent enzymes provides insight into their specific chemical pathways.