Validation of an RNA cell cycle progression score for predicting death from prostate cancer in a conservatively managed needle biopsy cohort.

BACKGROUND: The natural history of prostate cancer is highly variable and difficult to predict accurately. Better markers are needed to guide management and avoid unnecessary treatment. In this study, we validate the prognostic value of a cell cycle progression score (CCP score) independently and in a prespecified linear combination with standard clinical variables, that is, a clinical-cell-cycle-risk (CCR) score. METHODS: Paraffin sections from 761 men with clinically localized prostate cancer diagnosed by needle biopsy and managed conservatively in the United Kingdom, mostly between 2000 and 2003. The primary end point was prostate cancer death. Clinical variables consisted of centrally reviewed Gleason score, baseline PSA level, age, clinical stage, and extent of disease; these were combined into a single predefined risk assessment (CAPRA) score. Full data were available for 585 men who formed a fully independent validation cohort. RESULTS: In univariate analysis, the CCP score hazard ratio was 2.08 (95% CI (1.76, 2.46), P<10(-13)) for one unit change of the score. In multivariate analysis including CAPRA, the CCP score hazard ratio was 1.76 (95% CI (1.44, 2.14), P<10(-6)). The predefined CCR score was highly predictive, hazard ratio 2.17 (95% CI (1.83, 2.57), χ(2)=89.0, P<10(-20)) and captured virtually all available prognostic information. CONCLUSIONS: The CCP score provides significant pretreatment prognostic information that cannot be provided by clinical variables and is useful for determining which patients can be safely managed conservatively, avoiding radical treatment.

Prognostic value of a cell cycle progression signature for prostate cancer death in a conservatively managed needle biopsy cohort.

BACKGROUND: The natural history of prostate cancer is highly variable and it is difficult to predict. We showed previously that a cell cycle progression (CCP) score was a robust predictor of outcome in a conservatively managed cohort diagnosed by transurethral resection of the prostate. A greater need is to predict outcome in patients diagnosed by needle biopsy. METHODS: Total RNA was extracted from paraffin specimens. A CCP score was calculated from expression levels of 31 genes. Clinical variables consisted of centrally re-reviewed Gleason score, baseline prostate-specific antigen level, age, clinical stage, and extent of disease. The primary endpoint was death from prostate cancer. RESULTS: In univariate analysis (n=349), the hazard ratio (HR) for death from prostate cancer was 2.02 (95% CI (1.62, 2.53), P<10(-9)) for a one-unit increase in CCP score. The CCP score was only weakly correlated with standard prognostic factors and in a multivariate analysis, CCP score dominated (HR for one-unit increase=1.65, 95% CI (1.31, 2.09), P=3 × 10(-5)), with Gleason score (P=5 × 10(-4)) and prostate-specific antigen (PSA) (P=0.017) providing significant additional contributions. CONCLUSION: For conservatively managed patients, the CCP score is the strongest independent predictor of cancer death outcome yet described and may prove valuable in managing clinically localised prostate cancer.

Patterns of genomic loss of heterozygosity predict homologous recombination repair defects in epithelial ovarian cancer.

BACKGROUND: Defects in BRCA1, BRCA2, and other members of the homologous recombination pathway have potential therapeutic relevance when used to support agents that introduce or exploit double-stranded DNA breaks. This study examines the association between homologous recombination defects and genomic patterns of loss of heterozygosity (LOH). METHODS: Ovarian tumours from two independent data sets were characterised for defects in BRCA1, BRCA2, and RAD51C, and LOH profiles were generated. Publically available data were downloaded for a third independent data set. The same analyses were performed on 57 cancer cell lines. RESULTS: Loss of heterozygosity regions of intermediate size were observed more frequently in tumours with defective BRCA1 or BRCA2 (P=10(-11)). The homologous recombination deficiency (HRD) score was defined as the number of these regions observed in a tumour sample. The association between HRD score and BRCA deficiency was validated in two independent ovarian cancer data sets (P=10(-5) and 10(-29)), and identified breast and pancreatic cell lines with BRCA defects. CONCLUSION: The HRD score appears capable of detecting homologous recombination defects regardless of aetiology or mechanism. This score could facilitate the use of PARP inhibitors and platinum in breast, ovarian, and other cancers.

Evidence for HTR1A and LHPP as interacting genetic risk factors in major depression.

The HTR1A -1019C>G genotype was associated with major depression in the Utah population. Linkage analysis on Utah pedigrees with strong family histories of major depression including only cases with the HTR1A -1019G allele revealed a linkage peak on chromosome 10 (maximum HLOD=4.4). Sequencing of all known genes in the linkage region revealed disease-segregating single-nucleotide polymorphisms (SNPs) in LHPP. LHPP SNPs were also associated with major depression in both Utah and Ashkenazi populations. Consistent with the linkage evidence, LHPP associations depended on HTR1A genotype. Lhpp or a product of a collinear brain-specific transcript, therefore, may interact with Htr1a in the pathogenesis of major depression.

Impact of local and non-local interactions on thermodynamics and kinetics of protein folding.

To address the question of how the geometry of a protein's native conformation affects its folding and stability, we studied three model 36-mers on a cubic lattice. The native structure of one of these model 36-mers consisted mostly of local contacts, while that of a second consisted mostly of non-local contacts. The third native structure had a typical compact native conformation, and served as our reference. For each protein, the amino acid sequence was designed to have a pronounced energy minimum at its native conformation. We observed dramatic differences in folding, dependent on the presence or absence of non-local contacts. For the proteins with a typical large number of non-local contacts, the folding transition was all-or-none, whereas for the one with mostly local contacts, it was not. Although the maximum rate of folding was similar for all three proteins, we found that under conditions at which each native conformation was stable, the structure with mostly non-local contacts folded two orders of magnitude faster than the one with mostly local contacts. The statistical analysis of protein structure agrees fully with the implications of the theory. We discuss the importance of cooperativity in protein folding for its stability.

Two slow stages in refolding of bovine carbonic anhydrase B are due to proline isomerization.

Kinetics of refolding of bovine carbonic anhydrase B have been studied by the "double-jump" technique (i.e. the dependence of protein refolding on delay time in the unfolded state after fast unfolding). It is shown that two stages (the slow with a relaxation time of t1/2 approximately equal to 120 s and the superslow with t1/2 approximately equal to 600 s) observed during refolding of bovine carbonic anhydrase B are due to trans-cis isomerization of proline residues. The dependences of rate constants of these processes on temperature and on the final denaturant concentration were measured. Activation energies of both processes are the same, Ea = 18(+/- 2) kcal/mol. The rate constants of protein refolding do not depend on the final concentration of urea under native conditions. In addition, the rate of isomerization of essential proline residues in the "molten globule" intermediate state of bovine carbonic anhydrase was measured and found to be equal to that for unstructural polypeptides. The effect of several proline residues on carbonic anhydrase refolding is discussed.

Analysis of missense variation in human BRCA1 in the context of interspecific sequence variation.

INTRODUCTION: Interpretation of results from mutation screening of tumour suppressor genes known to harbour high risk susceptibility mutations, such as APC, BRCA1, BRCA2, MLH1, MSH2, TP53, and PTEN, is becoming an increasingly important part of clinical practice. Interpretation of truncating mutations, gene rearrangements, and obvious splice junction mutations, is generally straightforward. However, classification of missense variants often presents a difficult problem. From a series of 20,000 full sequence tests of BRCA1 carried out at Myriad Genetic Laboratories, a total of 314 different missense changes and eight in-frame deletions were observed. Before this study, only 21 of these missense changes were classified as deleterious or suspected deleterious and 14 as neutral or of little clinical significance. METHODS: We have used a combination of a multiple sequence alignment of orthologous BRCA1 sequences and a measure of the chemical difference between the amino acids present at individual residues in the sequence alignment to classify missense variants and in-frame deletions detected during mutation screening of BRCA1. RESULTS: In the present analysis we were able to classify an additional 50 missense variants and two in-frame deletions as probably deleterious and 92 missense variants as probably neutral. Thus we have tentatively classified about 50% of the unclassified missense variants observed during clinical testing of BRCA1. DISCUSSION: An internal test of the analysis is consistent with our classification of the variants designated probably deleterious; however, we must stress that this classification is tentative and does not have sufficient independent confirmation to serve as a clinically applicable stand alone method.

Domains in folding of model proteins.

By means of Monte Carlo simulation, we investigated the equilibrium between folded and unfolded states of lattice model proteins. The amino acid sequences were designed to have pronounced energy minimum target conformations of different length and shape. For short fully compact (36-mer) proteins, the all-or-none transition from the unfolded state to the native state was observed. This was not always the case for longer proteins. Among 12 designed sequences with the native structure of a fully compact 48-mer, a simple all-or-none transition was observed in only three cases. For the other nine sequences, three states of behavior-the native, denatured, and intermediate states-were found. The contiguous part of the native structure (domain) was conserved in the intermediate state, whereas the remaining part was completely unfolded and structureless. These parts melted separately from each other.

Formation of unique structure in polypeptide chains. Theoretical investigation with the aid of a replica approach.

A replica approach analogous to that used in spin glass systems is implemented to study the configurational space of a heteropolymeric model of protein with a quenched, disordered sequence of links in the limit of a large number of link types. It is shown that there exists a threshold value of chain heterogeneity which separates two qualitatively different types of behavior. For a low degree of heterogeneity the protein globule is like a homopolymer in a collapsed state without definite chain folds: an exponentially large number of folds make a significant contribution to the partition function in this regime. After the threshold heterogeneity has been overcome, the chain freezes drastically but without latent heat; few (approx. 1) frozen states with definite chain folds are thermodynamically dominant in this state. The relation of these results to thermodynamic aspects of protein folding is discussed.

[Helix-coil transition in the simplest model of large native RNA. I. Consideration of only native helices]. / Perekhod spiral'-klubok v prosteishei modeli bol'shikh prirodnykh RNK. I. Uchet tol'ko nativnykh spiralei.

A simple model of secondary structure of large native RNAs is proposed. The recurrence approach is described allowing to investigate helix-coil transition for this model by computer simulations. "Broken like" behavior of molecules of geometrical size has been observed upon changing the stability of helices at smooth changing of the helix fraction. Such behavior may be explained by the fact that helices formed by interactions between remote parts of the chain melt cooperatively by the mechanism of the second order phase transition. At the same time, helices formed by interactions between proximate segments in the chain melt practically independent of each other.

[Helix-coil transition in the simplest models of large native RNA. II. Consideration of nonspecific interactions]. / Perekhod spiral'-klubok v prosteishei modeli bol'shikh prirodnykh RNK. II. Uchet nespetsificheskikh vzaimodeistvii.

A simple model of secondary structure of large native RNAs is proposed which makes allowance for the formation of week non-specific helices in addition to the formation of strong helices. It is shown that in this model, the same as in the one investigated earlier [1], cooperative melting of helices formed by interactions between remote parts of chain by the mechanism of the second order phase transition is observed upon changing the stability of helices.

Genomic DNA pooling for whole-genome association scans in complex disease: empirical demonstration of efficacy in rheumatoid arthritis.

A pragmatic approach that balances the benefit of a whole-genome association (WGA) experiment against the cost of individual genotyping is to use pooled genomic DNA samples. We aimed to determine the feasibility of this approach in a WGA scan in rheumatoid arthritis (RA) using the validated human leucocyte antigen (HLA) and PTPN22 associations as test loci. A total of 203 269 single-nucleotide polymorphisms (SNPs) on the Affymetrix 100K GeneChip and Illumina Infinium microarrays were examined. A new approach to the estimation of allele frequencies from Affymetrix hybridization intensities was developed involving weighting for quality signals from the probe quartets. SNPs were ranked by z-scores, combined from United Kingdom and New Zealand case-control cohorts. Within a 1.7 Mb HLA region, 33 of the 257 SNPs and at PTPN22, 21 of the 45 SNPs, were ranked within the top 100 associated SNPs genome wide. Within PTPN22, individual genotyping of SNP rs1343125 within MAGI3 confirmed association and provided some evidence for association independent of the PTPN22 620W variant (P=0.03). Our results emphasize the feasibility of using genomic DNA pooling for the detection of association with complex disease susceptibility alleles. The results also underscore the importance of the HLA and PTPN22 loci in RA aetiology.

Influence of point mutations on protein structure: probability of a neutral mutation.

We investigate the probability X that a random mutation (i.e. the substitution in a random site of one amino acid residue by randomly chosen residue) will be a neutral one, i.e. it will not lead to a change in structure. Using a random energy model for the description of protein energy "levels" we show that this probability depends only on stiffness of a chain which is characterized by the number of conformations per peptide bond gamma. The result is X approximately gamma -8. The application of this result for protein engineering experiments and for possible scenario of neutral evolution is discussed.

Engineering of stable and fast-folding sequences of model proteins.

The statistical mechanics of protein folding implies that the best-folding proteins are those that have the native conformation as a pronounced energy minimum. We show that this can be obtained by proper selection of protein sequences and suggest a simple practical way to find these sequences. The statistical mechanics of these proteins with optimized native structure is discussed. These concepts are tested with a simple lattice model of a protein with full enumeration of compact conformations. Selected sequences are shown to have a native state that is very stable and kinetically accessible.

A new approach to the design of stable proteins.

We propose a simple algorithm to design a sequence which fits a given protein structure with a given energy. The algorithm is a modification of the Metropolis Monte Carlo scheme in sequence space with an evolutionary temperature which sets the energy scale. There is a one to one correspondence between this optimization scheme and the Ising model of ferromagnetism. This analogy implies that the design algorithm does not encounter multiple-minima problems and is very fast. The algorithm is tested by 'predicting' the primary structures of four proteins. In each case the calculated primary structures had statistically significant homology with the natural structures.

Why do protein architectures have Boltzmann-like statistics?

A theoretical study has shown that the occurrence of various structural elements in stable folds of random copolymers is exponentially dependent on the own energy of the element. A similar occurrence-on-energy dependence is observed in globular proteins from the level of amino acid conformations to the level of overall architectures. Thus, the structural features stabilized by many random sequences are typical of globular proteins while the features rarely observed in proteins are those which are stabilized by only a minor part of the random sequences.

Implications of thermodynamics of protein folding for evolution of primary sequences.

Natural proteins exhibit essentially two-state thermodynamics, with one stable fold that dominates thermodynamically over a vast number of possible folds, a number that increases exponentially with the size of the protein. Here we address the question of whether this feature of proteins is a rare property selected by evolution or whether it is in fact true of a significant proportion of all possible protein sequences. Using statistical procedures developed to study spin glasses, we show that, given certain assumptions, the probability that a randomly synthesized protein chain will have a dominant fold (which is the global minimum of free energy) is a function of temperature, and that below a critical temperature the probability rapidly increases as the temperature decreases. Our results suggest that a significant proportion of all possible protein sequences could have a thermodynamically dominant fold.

Evolution-like selection of fast-folding model proteins.

We propose an algorithm providing sequences of model proteins with rapid folding into a given target (native) conformation. This algorithm is applied to a chain of 27 residues on a cubic lattice. It generates sequences with folding 2 orders of magnitude faster than that of the practically random starting sequence. Thermodynamic analysis shows that the increase in speed is matched by an increase in stability: the evolved sequences are much more stable in their native conformation than the initial random sequence. The unfolding temperature for evolved sequences is slightly higher than the simulation temperature, bearing direct correspondence to the relatively low stability of real proteins.

How the first biopolymers could have evolved.

In this work, we discuss a possible origin of the first biopolymers with stable unique structures. We suggest that at the prebiotic stage of evolution, long organic polymers had to be compact to avoid hydrolysis and had to be soluble and thus must not be exceedingly hydrophobic. We present an algorithm that generates such sequences for model proteins. The evolved sequences turn out to have a stable unique structure, into which they quickly fold. This result illustrates the idea that the unique three-dimensional native structures of first biopolymers could have evolved as a side effect of nonspecific physicochemical factors acting at the prebiotic stage of evolution.