Assessing sources of inconsistencies in genotypes and their effects on genome-wide association studies with HapMap samples.

The discordance in results of independent genome-wide association studies (GWAS) indicates the potential for Type I and Type II errors. We assessed the repeatibility of current Affymetrix technologies that support GWAS. Reasonable reproducibility was observed for both raw intensity and the genotypes/copy number variants. We also assessed consistencies between different SNP arrays and between genotype calling algorithms. We observed that the inconsistency in genotypes was generally small at the specimen level. To further examine whether the differences from genotyping and genotype calling are possible sources of variation in GWAS results, an association analysis was applied to compare the associated SNPs. We observed that the inconsistency in genotypes not only propagated to the association analysis, but was amplified in the associated SNPs. Our studies show that inconsistencies between SNP arrays and between genotype calling algorithms are potential sources for the lack of reproducibility in GWAS results.

Analysis of a large structure/biological activity data set using recursive partitioning.

Combinatorial chemistry and high-throughput screening are revolutionizing the process of lead discovery in the pharmaceutical industry. Large numbers of structures and vast quantities of biological assay data are quickly being accumulated, overwhelming traditional structure/activity relationship (SAR) analysis technologies. Recursive partitioning is a method for statistically determining rules that classify objects into similar categories or, in this case, structures into groups of molecules with similar potencies. SCAM is a computer program implemented to make extremely efficient use of this methodology. Depending on the size of the data set, rules explaining biological data can be determined interactively. An example data set of 1650 monoamine oxidase inhibitors exemplifies the method, yielding substructural rules and leading to general classifications of these inhibitors. The method scales linearly with the number of descriptors, so hundreds of thousands of structures can be analyzed utilizing thousands to millions of molecular descriptors. There are currently no methods to deal with statistical analysis problems of this size. An important aspect of this analysis is the ability to deal with mixtures, i.e., identify SAR rules for classes of compounds in the same data set that might be binding in different ways. Most current quantitative structure/activity relationship methods require that the compounds follow a single mechanism. Advantages and limitations of this methodology are presented.

NMR structure of a receptor-bound G-protein peptide.

Heterotrimeric GTP-binding proteins (G proteins) regulate cellular activity by coupling to hormone or sensory receptors. Stimulated receptors catalyse the release of GDP from G protein alpha-subunits and GTP bound to the empty alpha-subunits provides signals that control effectors such as adenylyl cyclases, phosphodiesterases, phospholipases and ion channels. Three cytoplasmic loops of the activated receptor are thought to interact with three sites on the heterotrimeric G protein to provide high-affinity interaction and catalyse G-protein activation. The carboxyl terminus of the alpha-subunit is particularly important for interaction with the receptor. Here we study the structure of part of the active interface between the photon receptor rhodopsin and the G protein transducin, or Gt, using nuclear magnetic resonance. An 11-amino-acid peptide from the C terminus of the alpha-subunit of Gt (alpha t (340-350)) binds to rhodopsin and mimics the G protein in stabilizing its active form, metarhodopsin II. The peptide alpha t (340-350) binds to both excited and unexcited rhodopsin and conformational differences between the two bound forms suggest a mechanism for activation of G proteins by agonist-stimulated receptors. Insight into receptor-catalysed GDP release will have broad application because the GTP/GDP exchange and the intrinsic GTPase activity of GTP-binding proteins constitute a widespread regulatory mechanism.

NMR constrained solution structures for laminin peptide 11. Analogs define structural requirements for inhibition of tumor cell invasion of basement membrane matrix.

Peptide 11, CDPGYIGSR-NH2, is a segment of laminin which blocks tumor cell invasion. A high affinity laminin receptor in tumor cells is thought to be blocked by the carboxyl-terminal YIGSR, and conformational energy calculations suggest that the glycine in YIGSR allows an important conformational bend. We replaced the YIGSR glycine residue in peptide 11 with either D-alanine or L-alanine to allow or disfavor the proposed glycine bend. We found the Gly7-->D-Ala7 analog to be equal to peptide 11 in inhibiting tumor cell invasion of basement membrane matrix. The Gly7-->L-Ala7 analog was much less capable of invasion inhibition. Two-dimensional 1H-1H NMR was used to study the solution conformations of the peptide 11 analogs. NOESY experiments revealed close NH-NH contacts in peptide 11 and the D-Ala7 analog, but not in the L-Ala7 analog. Molecular dynamics generated low energy structures with excellent NOE agreement for peptide 11 and its analogs. Both peptide 11 and the D-Ala7 analog, but not the less active L-Ala7 analog, were predicted to have similar bends around Gly7 or D-Ala7. These results suggest that a bend in the YIGSR region of peptide 11 may be important for the binding of laminin to its metastasis-associated receptor.