Trials of the beta model for complex inheritance.

Theoretical advantages of nonparametric logarithm of odds to map polygenic diseases are supported by tests of the beta model that depends on a single logistic parameter and is the only model under which paternal and maternal transmissions to sibs of specified phenotypes are independent. Although it does not precisely describe recurrence risks in monozygous twins, the beta model has greater power to detect family resemblance or linkage than the more general delta model which describes the probability of 0, 1, or 2 alleles identical by descent (ibd) with two parameters. Available data on ibd in sibs are consistent with the beta model, but not with the equally parsimonious but less powerful gamma model that assumes a fixed probability of 1/2 for 1 allele ibd. Additivity of loci on the liability scale is not disproven. A simple equivalence extends the beta model to multipoint analysis.

A general additive distance with time-reversibility and rate variation among nucleotide sites.

As additivity is a very useful property for a distance measure, a general additive distance is proposed under the stationary time-reversible (SR) model of nucleotide substitution or, more generally, under the stationary, time-reversible, and rate variable (SRV) model, which allows rate variation among nucleotide sites. A method for estimating the mean distance and the sampling variance is developed. In addition, a method is developed for estimating the variance-covariance matrix of distances, which is useful for the statistical test of phylogenies and molecular clocks. Computer simulation shows (i) if the sequences are longer than, say, 1000 bp, the SR method is preferable to simpler methods; (ii) the SR method is robust against deviations from time-reversibility; (iii) when the rate varies among sites, the SRV method is much better than the SR method because the distance is seriously underestimated by the SR method; and (iv) our method for estimating the sampling variance is accurate for sequences longer than 500 bp. Finally, a test is constructed for testing whether DNA evolution follows a general Markovian model.

From residue matching patterns to protein folding topographies: General model and bovine pancreatic trypsin inhibitor

A coarse-grained model for protein-folding dynamics is introduced based on a discretized representation of torsional modes. The model, based on the Ramachandran map of the local torsional potential surface and the class (hydrophobic/polar/neutral) of each residue, recognizes patterns of both torsional conformations and hydrophobic-polar contacts, with tolerance for imperfect patterns. It incorporates empirical rates for formation of secondary and tertiary structure. The method yields a topological representation of the evolving local torsional configuration of the folding protein, modulo the basins of the Ramachandran map. The folding process is modeled as a sequence of transitions from one contact pattern to another, as the torsional patterns evolve. We test the model by applying it to the folding process of bovine pancreatic trypsin inhibitor, obtaining a kinetic description of the transitions between the contact patterns visited by the protein along the dominant folding pathway. The kinetics and detailed balance make it possible to invert the result to obtain a coarse topographic description of the potential energy surface along the dominant folding pathway, in effect to go backward or forward between a topological representation of the chain conformation and a topographical description of the potential energy surface governing the folding process. As a result, the strong structure-seeking character of bovine pancreatic trypsin inhibitor and the principal features of its folding pathway are reproduced in a reasonably quantitative way.

A general two-process model describes the hydrogen exchange behavior of RNase A in unfolding conditions.

When NMR hydrogen exchange was used previously to monitor the kinetics of RNase A unfolding, some peptide NH protons were found to show EX2 exchange (detected by base catalysis) in addition to the expected EX1 exchange, whose rate is limited by the kinetic unfolding process. In earlier work, two groups showed independently that a restricted two-process model successfully fits published hydrogen exchange rates of native RNase A in the range 0-0.7 M guanidinium chloride. We find that this model predicts properties that are very different from the observed properties of the EX2 exchange reactions of RNase A in conditions where guanidine-induced unfolding takes place. The model predicts that EX2 exchange should be too fast to measure by the technique used, whereas it is readily measurable. Possible explanations for the contradiction are considered here, and we show that removing the restriction from the earlier two-process model is sufficient to resolve the contradiction; instead of specifying that exchange caused by global unfolding occurs by the EX2 mechanism, we allow it to occur by the general mechanism, which includes both the EX1 and EX2 cases. It is logical to remove this restriction because global unfolding of RNase A is known to give rise to EX1 exchange in these unfolding conditions. Resolving the contradiction makes it possible to determine whether populated unfolding intermediates contribute to the EX2 exchange, and this question is considered elsewhere. The results and simulations indicate that moderate or high denaturant concentrations readily give rise to EX1 exchange in native proteins. Earlier studies showed that hydrogen exchange in native proteins typically occurs by the EX2 mechanism but that high temperatures or pH values above 7 may give rise to EX1 exchange. High denaturant concentrations should be added to the list of variables likely to cause EX1 exchange.

A general model of invariant chain association with class II major histocompatibility complex proteins.

The binding of invariant chain to major histocompatibility complex (MHC) proteins is an important step in processing of MHC class II proteins and in antigen presentation. The question of how invariant chain can bind to all MHC class II proteins is central to understanding these processes. We have employed molecular modeling to predict the structure of class II-associated invariant chain peptide (CLIP)-MHC protein complexes and to ask whether the predicted mode of association could be general across all MHC class II proteins. CLIP fits identically into the MHC class II alleles HLA-DR3, I-Ak, I-Au, and I-Ad, with a consistent pattern of hydrogen bonds, contacts, and hydrophobic burial and without bad contacts. Our model predicts the burial of CLIP residues Met-91 and Met-99 in the deep P1 and P9 anchor pockets and other detailed interactions, which we have compared with available data. The predicted pattern of I-A allele-specific effects on CLIP binding is very similar to that observed experimentally by alanine-scanning mutations of CLIP. Together, these results indicate that CLIP may bind in a single, general way across products of MHC class II alleles.