Multi-Overlap Simulations for Transitions between Reference Configurations

We introduce a new procedure to construct weight factors, which flatten the probability density of the overlap with respect to some pre-defined reference configuration. This allows one to overcome free energy barriers in the overlap variable. Subsequently, we generalize the approach to deal with the overlaps with respect to two reference configurations so that transitions between them are induced. We illustrate our approach by simulations of the brainpeptide Met-enkephalin with the ECEPP/2 energy function using the global-energy-minimum and the second lowest-energy states as reference configurations. The free energy is obtained as functions of the dihedral and the root-mean-square distances from these two configurations. The latter allows one to identify the transition state and to estimate its associated free energy barrier.Comment: 12 pages, (RevTeX), 14 figures, Phys. Rev. E, submitte

Parallelization of Markov chain generation and its application to the multicanonical method

We develop a simple algorithm to parallelize generation processes of Markov chains. In this algorithm, multiple Markov chains are generated in parallel and jointed together to make a longer Markov chain. The joints between the constituent Markov chains are processed using the detailed balance. We apply the parallelization algorithm to multicanonical calculations of the two-dimensional Ising model and demonstrate accurate estimation of multicanonical weights.Comment: 15 pages, 5 figures, uses elsart.cl

Secondary-Structure Design of Proteins by a Backbone Torsion Energy

We propose a new backbone-torsion-energy term in the force field for protein systems. This torsion-energy term is represented by a double Fourier series in two variables, the backbone dihedral angles phi and psi. It gives a natural representation of the torsion energy in the Ramachandran space in the sense that any two-dimensional energy surface periodic in both phi and psi can be expanded by the double Fourier series. We can then easily control secondary-structure-forming tendencies by modifying the torsion-energy surface. For instance, we can increase/decrease the alpha-helix-forming-tendencies by lowering/raising the torsion-energy surface in the alpha-helix region and likewise increase/decrease the beta-sheet-forming tendencies by lowering/raising the surface in the beta-sheet region in the Ramachandran space. We applied our approach to AMBER parm94 and AMBER parm96 force fields and demonstrated that our modifications of the torsion-energy terms resulted in the expected changes of secondary-structure-forming-tendencies by performing folding simulations of alpha-helical and beta-hairpin peptides.Comment: 13 pages, (Revtex4), 5 figure

Temperature dependence of ESR intensity for the nanoscale molecular magnet V15

The electron spin resonance (ESR) of nanoscale molecular magnet ${\rm V}_{15}$ is studied. Since the Hamiltonian of ${\rm V}_{15}$ has a large Hilbert space and numerical calculations of the ESR signal evaluating the Kubo formula with exact diagonalization method is difficult, we implement the formula with the help of the random vector technique and the Chebyshev polynominal expansion, which we name the double Chebyshev expansion method. We calculate the temperature dependence of the ESR intensity of ${\rm V}_{15}$ and compare it with the data obtained in experiment. As another complementary approach, we also implement the Kubo formula with the subspace iteration method taking only important low-lying states into account. We study the ESR absorption curve below $100{\rm K}$ by means of both methods. We find that side peaks appear due to the Dzyaloshinsky-Moriya interaction and these peaks grows as temperature decreases.Comment: 9 pages, 4 figures. To appear in J. Phys. Soc. Jpn. Supp

Identification of the first ATRIP-deficient patient and novel mutations in ATR define a clinical spectrum for ATR-ATRIP Seckel Syndrome

A homozygous mutational change in the Ataxia-Telangiectasia and RAD3 related (ATR) gene was previously reported in two related families displaying Seckel Syndrome (SS). Here, we provide the first identification of a Seckel Syndrome patient with mutations in ATRIP, the gene encoding ATR-Interacting Protein (ATRIP), the partner protein of ATR required for ATR stability and recruitment to the site of DNA damage. The patient has compound heterozygous mutations in ATRIP resulting in reduced ATRIP and ATR expression. A nonsense mutational change in one ATRIP allele results in a C-terminal truncated protein, which impairs ATR-ATRIP interaction; the other allele is abnormally spliced. We additionally describe two further unrelated patients native to the UK with the same novel, heterozygous mutations in ATR, which cause dramatically reduced ATR expression. All patient-derived cells showed defective DNA damage responses that can be attributed to impaired ATR-ATRIP function. Seckel Syndrome is characterised by microcephaly and growth delay, features also displayed by several related disorders including Majewski (microcephalic) osteodysplastic primordial dwarfism (MOPD) type II and Meier-Gorlin Syndrome (MGS). The identification of an ATRIP-deficient patient provides a novel genetic defect for Seckel Syndrome. Coupled with the identification of further ATR-deficient patients, our findings allow a spectrum of clinical features that can be ascribed to the ATR-ATRIP deficient sub-class of Seckel Syndrome. ATR-ATRIP patients are characterised by extremely severe microcephaly and growth delay, microtia (small ears), micrognathia (small and receding chin), and dental crowding. While aberrant bone development was mild in the original ATR-SS patient, some of the patients described here display skeletal abnormalities including, in one patient, small patellae, a feature characteristically observed in Meier-Gorlin Syndrome. Collectively, our analysis exposes an overlapping clinical manifestation between the disorders but allows an expanded spectrum of clinical features for ATR-ATRIP Seckel Syndrome to be define

Copy Number Alteration and Uniparental Disomy Analysis Categorizes Japanese Papillary Thyroid Carcinomas into Distinct Groups

The aim of the present study was to investigate chromosomal aberrations in sporadic Japanese papillary thyroid carcinomas (PTCs), concomitant with the analysis of oncogene mutational status. Twenty-five PTCs (11 with BRAFV600E, 4 with RET/PTC1, and 10 without mutation in HRAS, KRAS, NRAS, BRAF, RET/PTC1, or RET/PTC3) were analyzed using Genome-Wide Human SNP Array 6.0 which allows us to detect copy number alteration (CNA) and uniparental disomy (UPD), also referred to as copy neutral loss of heterozygosity, in a single experiment. The Japanese PTCs showed relatively stable karyotypes. Seven cases (28%) showed CNA(s), and 6 (24%) showed UPD(s). Interestingly, CNA and UPD were rarely overlapped in the same tumor; the only one advanced case showed both CNA and UPD with a highly complex karyotype. Thirteen (52%) showed neither CNA nor UPD. Regarding CNA, deletions tended to be more frequent than amplifications. The most frequent and recurrent region was the deletion in chromosome 22; however, it was found in only 4 cases (16%). The degree of genomic instability did not depend on the oncogene status. However, in oncogene-positive cases (BRAFV600E and RET/PTC1), tumors with CNA/UPD were less frequent (5/15, 33%), whereas tumors with CNA/UPD were more frequent in oncogene-negative cases (7/10, 70%), suggesting that chromosomal aberrations may play a role in the development of PTC, especially in oncogene-negative tumors. These data suggest that Japanese PTCs may be classified into three distinct groups: CNA+, UPD+, and no chromosomal aberrations. BRAFV600E mutational status did not correlate with any parameters of chromosomal defects

Dedifferentiation of Human Primary Thyrocytes into Multilineage Progenitor Cells without Gene Introduction

While identification and isolation of adult stem cells have potentially important implications, recent reports regarding dedifferentiation/reprogramming from differentiated cells have provided another clue to gain insight into source of tissue stem/progenitor cells. In this study, we developed a novel culture system to obtain dedifferentiated progenitor cells from normal human thyroid tissues. After enzymatic digestion, primary thyrocytes, expressing thyroglobulin, vimentin and cytokeratin-18, were cultured in a serum-free medium called SAGM. Although the vast majority of cells died, a small proportion (∼0.5%) survived and proliferated. During initial cell expansion, thyroglobulin/cytokeratin-18 expression was gradually declined in the proliferating cells. Moreover, sorted cells expressing thyroid peroxidase gave rise to proliferating clones in SAGM. These data suggest that those cells are derived from thyroid follicular cells or at least thyroid-committed cells. The SAGM-grown cells did not express any thyroid-specific genes. However, after four-week incubation with FBS and TSH, cytokeratin-18, thyroglobulin, TSH receptor, PAX8 and TTF1 expressions re-emerged. Moreover, surprisingly, the cells were capable of differentiating into neuronal or adipogenic lineage depending on differentiating conditions. In summary, we have developed a novel system to generate multilineage progenitor cells from normal human thyroid tissues. This seems to be achieved by dedifferentiation of thyroid follicular cells. The presently described culture system may be useful for regenerative medicine, but the primary importance will be as a tool to elucidate the mechanisms of thyroid diseases

Enhanced and effective conformational sampling of protein molecular systems for their free energy landscapes

Protein folding and protein–ligand docking have long persisted as important subjects in biophysics. Using multicanonical molecular dynamics (McMD) simulations with realistic expressions, i.e., all-atom protein models and an explicit solvent, free-energy landscapes have been computed for several systems, such as the folding of peptides/proteins composed of a few amino acids up to nearly 60 amino-acid residues, protein–ligand interactions, and coupled folding and binding of intrinsically disordered proteins. Recent progress in conformational sampling and its applications to biophysical systems are reviewed in this report, including descriptions of several outstanding studies. In addition, an algorithm and detailed procedures used for multicanonical sampling are presented along with the methodology of adaptive umbrella sampling. Both methods control the simulation so that low-probability regions along a reaction coordinate are sampled frequently. The reaction coordinate is the potential energy for multicanonical sampling and is a structural identifier for adaptive umbrella sampling. One might imagine that this probability control invariably enhances conformational transitions among distinct stable states, but this study examines the enhanced conformational sampling of a simple system and shows that reasonably well-controlled sampling slows the transitions. This slowing is induced by a rapid change of entropy along the reaction coordinate. We then provide a recipe to speed up the sampling by loosening the rapid change of entropy. Finally, we report all-atom McMD simulation results of various biophysical systems in an explicit solvent

Cross-Regulation between Oncogenic BRAFV600E Kinase and the MST1 Pathway in Papillary Thyroid Carcinoma

BACKGROUND:The BRAF(V600E) mutation leading to constitutive signaling of MEK-ERK pathways causes papillary thyroid cancer (PTC). Ras association domain family 1A (RASSF1A), which is an important regulator of MST1 tumor suppressor pathways, is inactivated by hypermethylation of its promoter region in 20 to 32% of PTC. However, in PTC without RASSF1A methylation, the regulatory mechanisms of RASSF1A-MST1 pathways remain to be elucidated, and the functional cooperation or cross regulation between BRAF(V600E) and MST1,which activates Foxo3,has not been investigated. METHODOLOGY/PRINCIPAL FINDINGS:The negative regulators of the cell cycle, p21 and p27, are strongly induced by transcriptional activation of FoxO3 in BRAF(V600E) positive thyroid cancer cells. The FoxO3 transactivation is augmented by RASSF1A and the MST1 signaling pathway. Interestingly, introduction of BRAF(V600E)markedly abolished FoxO3 transactivation and resulted in the suppression of p21 and p27 expression. The suppression of FoxO3 transactivation by BRAF(V600E)is strongly increased by coexpression of MST1 but it is not observed in the cells in which MST1, but not MST2,is silenced. Mechanistically, BRAF(V600E)was able to bind to the C-terminal region of MST1 and resulted in the suppression of MST1 kinase activities. The induction of the G1-checkpoint CDK inhibitors, p21 and p27,by the RASSF1A-MST1-FoxO3 pathway facilitates cellular apoptosis, whereas addition of BRAF(V600E) inhibits the apoptotic processes through the inactivation of MST1. Transgenic induction of BRAF(V600E)in the thyroid gland results in cancers resembling human papillary thyroid cancers. The development of BRAF(V600E)transgenic mice with the MST1 knockout background showed that these mice had abundant foci of poorly differentiated carcinomas and large areas without follicular architecture or colloid formation. CONCLUSIONS/SIGNIFICANCE:The results of this study revealed that the oncogenic effect of BRAF(V600E) is associated with the inhibition of MST1 tumor suppressor pathways, and that the activity of RASSF1A-MST1-FoxO3 pathways determines the phenotypes of BRAF(V600E) tumors

Self-renewal and chemotherapy resistance of p75NTR positive cells in esophageal squamous cell carcinomas

<p>Abstract</p> <p>Background</p> <p>p75^NTRhas been used to isolate esophageal and corneal epithelial stem cells. In the present study, we investigated the expression of p75^NTRin esophageal squamous cell carcinoma (ESCC) and explored the biological properties of p75^NTR+cells.</p> <p>Methods</p> <p>p75^NTRexpression in ESCC was assessed by immunohistochemistry. p75^NTR+and p75^NTR-cells of 4 ESCC cell lines were separated by fluorescence-activated cell sorting. Differentially expressed genes between p75^NTR+and p75^NTR-cells were determined by real-time quantitative reverse transcription-PCR. Sphere formation assay, DDP sensitivity assay, ⁶⁴copper accumulation assay and tumorigenicity analysis were performed to determine the capacity of self-renewal, chemotherapy resistance and tumorigenicity of p75^NTR+cells.</p> <p>Results</p> <p>In ESCC specimens, p75^NTRwas found mainly confined to immature cells and absent in cells undergoing terminal differentiation. The percentage of p75^NTR+cells was 1.6%–3.7% in Eca109 and 3 newly established ESCC cell lines. The expression of Bmi-1, which is associated with self-renewal of stem cells, was significantly higher in p75^NTR+cells. p63, a marker identified in keratinocyte stem cells, was confined mainly to p75^NTR+cells. The expression of CTR1, which is associated with cisplatin (DDP)-resistance, was significantly decreased in p75^NTR+cells. Expression levels of differentiation markers, such as involucrin, cytokeratin 13, β1-integrin and β4-integrin, were lower in p75^NTR+cells. In addition, p75^NTR+cells generated both p75^NTR+and p75^NTR-cells, and formed nonadherent spherical clusters in serum-free medium supplemented with growth factors. Furthermore, p75^NTR+cells were found to be more resistant to DDP and exhibited lower ⁶⁴copper accumulation than p75^NTR-cells.</p> <p>Conclusion</p> <p>Our results demonstrated that p75^NTR+cells possess some characteristics of CSCs, namely, self-renewal and chemotherapy resistance. Chemotherapy resistance of p75^NTR+cells may probably be attributable to decreased expression of CTR1.</p