Brownian dynamics simulations of protein folding: Access to milliseconds time scale and beyond

Protein folding occurs on a time scale ranging from milliseconds to minutes for a majority of proteins. Computer simulation of protein folding, from a random configuration to the native structure, is nontrivial owing to the large disparity between the simulation and folding time scales. As an effort to overcome this limitation, simple models with idealized protein subdomains, e.g., the diffusion–collision model of Karplus and Weaver, have gained some popularity. We present here new results for the folding of a four-helix bundle within the framework of the diffusion–collision model. Even with such simplifying assumptions, a direct application of standard Brownian dynamics methods would consume 10,000 processor-years on current supercomputers. We circumvent this difficulty by invoking a special Brownian dynamics simulation. The method features the calculation of the mean passage time of an event from the flux overpopulation method and the sampling of events that lead to productive collisions even if their probability is extremely small (because of large free-energy barriers that separate them from the higher probability events). Using these developments, we demonstrate that a coarse-grained model of the four-helix bundle can be simulated in several days on current supercomputers. Furthermore, such simulations yield folding times that are in the range of time scales observed in experiments.

A maternal diet high in n − 6 polyunsaturated fats alters mammary gland development, puberty onset, and breast cancer risk among female rat offspring

We hypothesized that feeding pregnant rats with a high-fat diet would increase both circulating 17β-estradiol (E2) levels in the dams and the risk of developing carcinogen-induced mammary tumors among their female offspring. Pregnant rats were fed isocaloric diets containing 12% or 16% (low fat) or 43% or 46% (high fat) of calories from corn oil, which primarily contains the n − 6 polyunsaturated fatty acid (PUFA) linoleic acid, throughout pregnancy. The plasma concentrations of E2 were significantly higher in pregnant females fed a high n − 6 PUFA diet. The female offspring of these rats were fed with a laboratory chow from birth onward, and when exposed to 7,12-dimethylbenz(a)anthracene had a significantly higher mammary tumor incidence (60% vs. 30%) and shorter latency for tumor appearance (11.4 ± 0.5 weeks vs. 14.2 ± 0.6 weeks) than the offspring of the low-fat mothers. The high-fat offspring also had puberty onset at a younger age, and their mammary glands contained significantly higher numbers of the epithelial structures that are the targets for malignant transformation. Comparable changes in puberty onset, mammary gland morphology, and tumor incidence were observed in the offspring of rats treated daily with 20 ng of E2 during pregnancy. These data, if extrapolated to humans, may explain the link among diet, early puberty onset, mammary parenchymal patterns, and breast cancer risk, and indicate that an in utero exposure to a diet high in n − 6 PUFA and/or estrogenic stimuli may be critical for affecting breast cancer risk.

Ribosomes can slide over and beyond “hungry” codons, resuming protein chain elongation many nucleotides downstream

In cells subjected to moderate aminoacyl-tRNA limitation, the peptidyl-tRNA–ribosome complex stalled at the “hungry” codon can slide well beyond it on the messenger RNA and resume translation further downstream. This behavior is proved by unequivocal amino acid sequence data, showing a protein that lacks the bypassed sequence encoded between the hungry codon and specific landing sites. The landing sites are codons cognate to the anticodon of the peptidyl-tRNA. The efficiency of this behavior can be as high as 10–20% but declines with the length of the slide. Interposition of “trap” sites (nonproductive landing sites) in the bypassed region reduces the frequency of successful slides, confirming that the ribosome–peptidyl-tRNA complex passes through the untranslated region of the message. This behavior appears to be quite general: it can occur at the two kinds of hungry codons tested, AUA and AAG; the sliding peptidyl-tRNA can be any of three species tested, phenylalanine, tyrosine, or leucine tRNA; the peptidyl component can be either of two very different peptide sequences; and translation can resume at any of the three codons tested.

Statistical prediction of single-stranded regions in RNA secondary structure and application to predicting effective antisense target sites and beyond

Single-stranded regions in RNA secondary structure are important for RNA–RNA and RNA–protein interactions. We present a probability profile approach for the prediction of these regions based on a statistical algorithm for sampling RNA secondary structures. For the prediction of phylogenetically-determined single-stranded regions in secondary structures of representative RNA sequences, the probability profile offers substantial improvement over the minimum free energy structure. In designing antisense oligonucleotides, a practical problem is how to select a secondary structure for the target mRNA from the optimal structure(s) and many suboptimal structures with similar free energies. By summarizing the information from a statistical sample of probable secondary structures in a single plot, the probability profile not only presents a solution to this dilemma, but also reveals ‘well-determined’ single-stranded regions through the assignment of probabilities as measures of confidence in predictions. In antisense application to the rabbit β-globin mRNA, a significant correlation between hybridization potential predicted by the probability profile and the degree of inhibition of in vitro translation suggests that the probability profile approach is valuable for the identification of effective antisense target sites. Coupling computational design with DNA–RNA array technique provides a rational, efficient framework for antisense oligonucleotide screening. This framework has the potential for high-throughput applications to functional genomics and drug target validation.

Radial and longitudinal diffusion of myoglobin in single living heart and skeletal muscle cells

We have used a fluorescence recovery after photobleaching (FRAP) technique to measure radial diffusion of myoglobin and other proteins in single skeletal and cardiac muscle cells. We compare the radial diffusivities, Dr (i.e., diffusion perpendicular to the long fiber axis), with longitudinal ones, Dl (i.e., parallel to the long fiber axis), both measured by the same technique, for myoglobin (17 kDa), lactalbumin (14 kDa), and ovalbumin (45 kDa). At 22°C, Dl for myoglobin is 1.2 × 10−7 cm2/s in soleus fibers and 1.1 × 10−7 cm2/s in cardiomyocytes. Dl for lactalbumin is similar in both cell types. Dr for myoglobin is 1.2 × 10−7 cm2/s in soleus fibers and 1.1 × 10−7 cm2/s in cardiomyocytes and, again, similar for lactalbumin. Dl and Dr for ovalbumin are 0.5 × 10−7 cm2/s. In the case of myoglobin, both Dl and Dr at 37°C are about 80% higher than at 22°C. We conclude that intracellular diffusivity of myoglobin and other proteins (i) is very low in striated muscle cells, ≈1/10 of the value in dilute protein solution, (ii) is not markedly different in longitudinal and radial direction, and (iii) is identical in heart and skeletal muscle. A Krogh cylinder model calculation holding for steady-state tissue oxygenation predicts that, based on these myoglobin diffusivities, myoglobin-facilitated oxygen diffusion contributes 4% to the overall intracellular oxygen transport of maximally exercising skeletal muscle and less than 2% to that of heart under conditions of high work load.

Maspin acts at the cell membrane to inhibit invasion and motility of mammary and prostatic cancer cells.

Maspin, a novel serine protease inhibitor (serpin), inhibits tumor invasion and metastasis of mammary carcinoma. We show here that recombinant maspin protein blocks the motility of these carcinoma cells in culture over 12 h, as demonstrated by time-lapse video microscopy. Lamellopodia are withdrawn but ruffling continues. Both exogenous recombinant maspin and maspin expressed by tumor transfectants exhibit inhibitory effects on cell motility and cell invasion as shown in modified Boyden chamber assays. In addition, three prostatic cancer cell lines treated with recombinant maspin exhibited similar inhibition of both invasion and motility, suggesting a similar mode of maspin action in these two glandular epithelial cancers. When mammary carcinoma cells were treated with recombinant maspin, the protein was shown by immunostaining to bind specifically to the cell surface, suggesting that maspin activity is membrane associated. When pretreated with antimaspin antibody, maspin loses its inhibitory effects on both invasion and motility. However, when maspin is added to these cells preceding antibody treatment, the activity of maspin is no longer inhibited by subsequent addition of the antibody. It is concluded therefore that the inhibition of invasion and motility by maspin is initially localized to the cell surface.