Confidence Intervals on Ratios of Variance Components for the Unbalanced Two Factor Nested Model

The model considered in this article is the two-factor nested unbalanced variance component model: for p = 1, 2, …, P; q = 1, 2, …, Q_p; and r = 1, 2, …, R_pq. The random variables Y_pqr are observable. The constant μ is an unknown parameter, and A_p, B_pq and C_pqr are (unobservable) normal and independently distributed random variables with zero means and finite variances σ²_A, σ²_B, and σ²_C, respectively. Approximate confidence intervals on ?_A and ?_B using unweighted means are derived, where The performance of these approximate confidence intervals are evaluated using computer simulation. The simulated results indicate that these proposed confidence intervals perform satisfactorily and can be used in applied problems.     

Preferential inhibition of acetylcholinesterase molecular forms in rat brain

The effect of eight different acetylcholinesterase inhibitors (AChEIs) on the activity of acetylcholinesterase (AChE) molecular forms was investigated. Aqueous-soluble and detergent-soluble AChE molecular forms were separated from rat brain homogenate by sucrose density sedimentation. The bulk of soluble AChE corresponds to globular tetrameric (G₄), and monomeric (G₁) forms. Heptylphysostigmine (HEP) and diisopropylfluorophosphate were more selective for the G₁ than for the G₄ form in aqueous-soluble extract. Neostigmine showed slightly more selectivity for the G₁ form both in aqueous- and detergent-soluble extracts. Other drugs such as physostigmine, echothiophate, BW284C51, tetrahydroaminoacridine, and metrifonate inhibited both aqueous- and detergent-soluble AChE molecular forms with similar potency. Inhibition of aqueous-soluble AChE by HEP was highly competitive with Triton X-100 in a gradient, indicating that HEP may bind to a detergent-sensitive non-catalytic site of AChE. These results suggest a differential sensitivity among AChE molecular forms to inhibition by drugs through an allosteric mechanism. The application of these properties in developing AChEIs for treatment of Alzheimer disease is considered.Special issue dedicated to Dr. Morris H. Aprison.     

Molecular mimicry and the autoimmune response to the peripheral nerve myelin PO glycoprotein

In the Lewis rat immunisation with the myelin PO glycoprotein can induce an inflammatory demyelinating disease of the peripheral nervous system, experimental allergic neuritis (EAN), which has many clinical and histopathological parallels with the human disease the Guillain-Barre syndrome. In view of the reported association of GBS with a number of infectious agents we have investigated whether molecular mimicry may occur between microbial antigens and the PO protein that could possibly trigger a similar pathogenic autoimmune response in man. A computer search of the available protein sequence data bases identified several absolute sequence homologies between PO and viral proteins that involve five or more consecutive amino acid residues. Four of these sequence homologies involved viral pathogens previously associated with the Guillain-Barre syndrome, namely Epstein-Barr virus (EBV), cytomegalovirus (CMV), Varicella zoster virus (VZV) and human immunodeficiency virus I (HIV I). Although, sequence homologies were also found between viral peptides and the neuritogenic determinants of PO, residues 56–71 and 180–199, these homologies proved incapable of eliciting EAN in the Lewis rat. These observations are discussed with reference to the role that molecular mimicry between T cell epitopes on pathogen derived antigens and the PO protein may play in the pathogenesis of the Guillain-Barre syndrome.Abbreviations EAN Experimental allergic neuritis - EAE experimental allergic encephalomyelitis - PNS peripheral nervous system - CNS central nervous system - MBP myelin basic protein - GBS Guillain Barre syndrome - CFA complete Freund's Adjuvant - LPC lysophosphatidyl choline - VZV Varicella zoster virus - CMV cytomegalovirus - EBV Epstein Barr virus - HIV I human immunodeficiency virus I Special issue dedicated to Dr. Alan N. Davison     

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.     

Multiple protein sequence alignment from tertiary structure comparison: assignment of global and residue confidence levels.

An algorithm is presented for the accurate and rapid generation of multiple protein sequence alignments from tertiary structure comparisons. A preliminary multiple sequence alignment is performed using sequence information, which then determines an initial superposition of the structures. A structure comparison algorithm is applied to all pairs of proteins in the superimposed set and a similarity tree calculated. Multiple sequence alignments are then generated by following the tree from the branches to the root. At each branchpoint of the tree, a structure-based sequence alignment and coordinate transformations are output, with the multiple alignment of all structures output at the root. The algorithm encoded in STAMP (STructural Alignment of Multiple Proteins) is shown to give alignments in good agreement with published structural accounts within the dehydrogenase fold domains, globins, and serine proteinases. In order to reduce the need for visual verification, two similarity indices are introduced to determine the quality of each generated structural alignment. Sc quantifies the global structural similarity between pairs or groups of proteins, whereas Pij' provides a normalized measure of the confidence in the alignment of each residue. STAMP alignments have the quality of each alignment characterized by Sc and Pij' values and thus provide a reproducible resource for studies of residue conservation within structural motifs.     

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.