Effect of Zeitgeber Intensity Reduction on a Simulated Dual-Oscillator Human Circadian System: Classical and Dynamic Analysis

The two-oscillator model of human circadian rhythmicity was analyzed when a zeitgeber relative intensity of 1, 0.5, or 0.1 was introduced into the equations. Fourier analysis was compared with dynamic analysis such as attractor reconstruction or Liapunov exponent calculation. After a 50 or 90% reduction in zeitgeber intensity, the dynamics of the system became equivalent and differed significantly from those of a system with maximal zeitgeber intensity. When 10% aleatory noise was added to the data, the analysis was still applicable, and the results obtained were essentially the same as in the absence of noise. Dynamic analysis could thus provide a distinct classification for periodic data, based on the type of analysis.     

Molecular modeling of the amphipathic helices of the plasma apolipoproteins.

In this paper we propose a classification of the amphipathic helical repeats occurring in the plasma apolipoprotein sequences. It is based upon the calculation of the molecular hydrophobicity potential around the helical segments. The repeats were identified using a new autocorrelation matrix, based upon similarities of hydrophobic and hydrophilic properties of the amino acid residues within the apolipoprotein sequences. The helices were constructed by molecular modeling, the molecular hydrophobicity potential was calculated, and isopotential contour lines drawn around the helices yielded a three-dimensional visualization of the hydrophobicity potential. Two classes of apolipoproteins could be differentiated by comparing the hydrophobic angles obtained by projection of the isopotential contour lines on a plane perpendicular to the long axis of the helix. The isopotential contour lines around apo AI, AIV, and E are more hydrophilic than hydrophobic, whereas they are of similar intensity for apo AII, CI, and CIII. In both cases discoidal lipid-protein complexes are generated, with the amphipathic helices around the edge of the lipid core. The long axis of the helices is oriented parallel to the phospholipid acyl chains and the hydrophilic side of the helix toward the aqueous phase. As a result of the differences in hydrophobicity potential, the contact between the hydrophobic side of the helices and the phospholipid acyl chains is larger for apo AII, CI, and CIII than for the other apolipoproteins. This might account for the greater stability of the discoidal complexes generated between phospholipids and these apoproteins.     

Crystal structure of human immunoglobulin fragment Fab New refined at 2.0 A resolution.

The three-dimensional structure of the human immunoglobulin fragment Fab New (IgG1, lambda) has been refined to a crystallographic R-factor of 16.9% to 2 A resolution. Rms deviations of the final model from ideal geometry are 0.014 A for bond distances and 3.03 degrees for bond angles. Refinement was based on a new X-ray data set including 28,301 reflections with F > 2.5 sigma(F) from 6.0 to 2.0 A resolution. The starting model for the refinement procedure reported here is from the Brookhaven Protein Data Bank entry 3FAB (rev. 1981). Differences between the initial and final models include modified polypeptide-chain folding in the third complementarity-determining region (CDR3) and the third framework region (FR3) of VH and in some exposed loops of CL and CH1. Amino acid sequence changes were determined at a number of positions by inspection of difference electron density maps. The incorporation of amino acid sequence changes results in an improved VH framework model for the "humanization" of monoclonal antibodies.     

Recurrent alpha beta loop structures in TIM barrel motifs show a distinct pattern of conserved structural features.

A systematic survey of seven parallel alpha/beta barrel protein domains, based on exhaustive structural comparisons, reveals that a sizable proportion of the alpha beta loops in these proteins--20 out of a total of 49--belong to either one of two loop types previously described by Thornton and co-workers. Six loops are of the alpha beta 1 type, with one residue between the alpha-helix and beta-strand, and 13 are of the alpha beta 3 type, with three residues between the helix and the strand. Protein fragments embedding the identified loops, and termed alpha beta connections since they contain parts of the flanking helix and strand, have been analyzed in detail revealing that each type of connection has a distinct set of conserved structural features. The orientation of the beta-strand relative to the helix and loop portions is different owing to a very localized difference in backbone conformation. In alpha beta 1 connections, the chain enters the beta-strand via a residue adopting an extended conformation, while in alpha beta 3 it does so via a residue in a near alpha-helical conformation. Other conserved structural features include distinct patterns of side chain orientation relative to the beta-sheet surface and of main chain H-bonds in the loop and the beta-strand moieties. Significant differences also occur in packing interactions of conserved hydrophobic residues situated in the last turn of the helix. Yet the alpha-helix surface of both types of connections adopts similar orientations relative to the barrel sheet surface. Our results suggest furthermore that conserved hydrophobic residues along the sequence of the connections, may be correlated more with specific patterns of interactions made with neighboring helices and sheet strands than with helix/strand packing within the connection itself. A number of intriguing observations are also made on the distribution of the identified alpha beta 1 and alpha beta 3 loops within the alpha/beta-barrel motifs. They often occur adjacent to each other; alpha beta 3 loops invariably involve even numbered beta-strands, while alpha beta 1 loops involve preferentially odd beta-strands; all the analyzed proteins contain at least one alpha beta 3 loop in the first half of the eightfold alpha/beta barrel. Possible origins of all these observations, and their relevance to the stability and folding of parallel alpha/beta barrel motifs are discussed.     

Detection of native-like models for amino acid sequences of unknown three-dimensional structure in a data base of known protein conformations.

We present an approach which can be used to identify native-like folds in a data base of protein conformations in the absence of any sequence homology to proteins in the data base. The method is based on a knowledge-based force field derived from a set of known protein conformations. A given sequence is mounted on all conformations in the data base and the associated energies are calculated. Using several conformations and sequences from the globin family we show that the native conformation is identified correctly. In fact the resolution of the force field is high enough to discriminate between a native fold and several closely related conformations. We then apply the procedure to several globins of known sequence but unknown three dimensional structure. The homology of these sequences to globins of known structures in the data base ranges from 49 to 17%. With one exception we find that for all globin sequences one of the known globin folds is identified as the most favorable conformation. These results are obtained using a force field derived from a data base devoid of globins of known structure. We briefly discuss useful applications in protein structural research and future development of our approach.     

Three-dimensional structure of the quinoprotein methylamine dehydrogenase from Paracoccus denitrificans determined by molecular replacement at 2.8 A resolution.

The three-dimensional structure of the quinoprotein methylamine dehydrogenase from Paracoccus dentrificans (PD-MADH) has been determined at 2.8 A resolution by the molecular replacement method combined with map averaging procedures, using data collected from an area detector. The structure of methylamine dehydrogenase from Thio-bacillus versutus, which contains an "X-ray" sequence, was used as the starting search model. MADH consists of 2 heavy (H) and 2 light (L) subunits related by a molecular 2-fold axis. The H subunit is folded into seven four-stranded beta segments, forming a disk-shaped structure, arranged with pseudo-7-fold symmetry. A 31-residue elongated tail exists at the N-terminus of the H subunit in MADH from T. versutus but is partially digested in this crystal form of MADH from P. denitrificans, leaving the H subunit about 18 residues shorter. Each L subunit contains 127 residues arranged into 10 beta-strands connected by turns. The active site of the enzyme is located in the L subunit and is accessible via a hydrophobic channel between the H and L subunits. The redox cofactor of MADH, tryptophan tryptophylquinone is highly unusual. It is formed from two covalently linked tryptophan side chains at positions 57 and 107 of the L subunit, one of which contains an orthoquinone.     

Molecular motors in axonal transport

Neurons require a large amount of intracellular transport. Cytoplasmic polypeptides and membrane-bounded organelles move from the perikaryon, down the length of the axon, and to the synaptic terminals. This movement occurs at distinct rates and is termed axonal transport. Axonal transport is divided into the slow transport of cytoplasmic proteins including glycolytic enzymes and cytoskeletal structures and the fast transport of membrane-bounded organelles along linear arrays of microtubules. The polypeptide compositions of the rate classes of axonal transport have been well characterized, but the underlying molecular mechanisms of this movement are less clear. Progress has been particularly slow toward understanding force-generation in slow transport, but recent developments have provided insight into the molecular motors involved in fast axonal transport. Recent advances in the cellular and molecular biology of one fast axonal transport motor, kinesin, have provided a clearer understanding of organelle movement along microtubules. The availability of cellular and molecular probes for kinesin and other putative axonal transport motors have led to a reevaluation of our understanding of intracellular motility.     

Stochastic simulation of clonal growth in the tall goldenrod,Solidago altissima

Summary As clonal plants grow they move through space. The movement patterns that result can be complex and difficult to interpret without the aid of models. We developed a stochastic simulation model of clonal growth in the tall goldenrod, Solidago altissima. Our model was calibrated with field data on the clonal expansion of both seedlings and established clones, and model assumptions were verified by statistical analyses.When simulations were based on empirical distributions with long rhizome lengths, there was greater dispersal, less leaf overlap, and less spatial aggregation than when simulations were based on distributions with comparatively short rhizome lengths. For the field data that we utilized, variation in rhizome lengths had a greater effect than variation for either branching angles or rhizome initiation points (see text). We also found that observed patterns of clonal growth in S. altissima did not cause the formation of fairy rings. However, simulations with an artificial distribution of branching angles demonstrate that fairy rings can result solely from a plant's clonal morphology.Stochastic simulation models that incorporated variation in rhizome lengths, branching angles, and rhizome initiation points produced greater dispersal and less leaf overlap than deterministic models. Thus, variation for clonal growth parameters may increase the efficiency of substrate exploration by increasing the area covered and by decreasing the potential for intraclonal competition. We also demonstrated that ramet displacements were slightly, but consistently lower in stochastic simulation models than in random-walk models. This difference was due to the incorporation of details on rhizome bud initiation into stochastic simulation models, but not random-walk models. We discuss the advantages and disadvantages of deterministic, stochastic simulation, and random-walk models of clonal growth.     

Kinetic and molecular orbital studies on the rate of oxidation of monosubstituted phenols and anilines by lactoperoxidase compound II in comparison with the case of horseradish peroxidase

The second-order rate constant (k4) for the oxidation of monosubstituted phenols and anilines by lactoperoxidase compound II was examined by Chance's method [B. Chance, Arch. Biochem. Biophys. 71 (1957), 130–136]. When the electronic states of these substrates were calculated by an ab initio molecular orbital method, it was found that the log k4 value correlates well with the highest occupied molecular orbital (HOMO) energy level but not with the net charge or frontier electron density. These results are essentially similar to those reported previously in the case of horseradish peroxidase [J. Sakurada, R. Sekiguchi, K. Sato, and T. Hosoya, Biochemistry 29 (1990), 4093–4098], showing some dissimilar features which are considered to reflect the structural difference between the two enzymes.Abbreviations HOMO highest occupied molecular orbital - HRP horseradish peroxidase - LPO lactoperoxidase (EC 1.11.1.7) - LUMO lowest unoccupied molecular orbital     

Weak field influence on biomolecular changes

Model equations for the kinetics of the synthesis and decay of molecular aggregates are used to show the high sensitivity of equilibrium concentrations of high-molecular aggregates to external radiation. This phenomenon is used to explain the effects of low-intensity microwave fields on the functioning of biological systems. The experimental results on the influence of SHF-radiation on ferricyanide reduction by erythrocytes are interpreted in detail.