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.

The RIIβ regulatory subunit of protein kinase A binds to cAMP response element: An alternative cAMP signaling pathway

cAMP, through the activation of cAMP-dependent protein kinase (PKA), is involved in transcriptional regulation. In eukaryotic cells, cAMP is not considered to alter the binding affinity of CREB/ATF to cAMP-responsive element (CRE) but to induce serine phosphorylation and consequent increase in transcriptional activity. In contrast, in prokaryotic cells, cAMP enhances the DNA binding of the catabolite repressor protein to regulate the transcription of several operons. The structural similarity of the cAMP binding sites in catabolite repressor protein and regulatory subunit of PKA type II (RII) suggested the possibility of a similar role for RII in eukaryotic gene regulation. Herein we report that RIIβ subunit of PKA is a transcription factor capable of interacting physically and functionally with a CRE. In contrast to CREB/ATF, the binding of RIIβ to a CRE was enhanced by cAMP, and in addition, RIIβ exhibited transcriptional activity as a Gal4-RIIβ fusion protein. These experiments identify RIIβ as a component of an alternative pathway for regulation of CRE-directed transcription in eukaryotic cells.

The Ti plasmid increases the efficiency of Agrobacterium tumefaciens as a recipient in virB-mediated conjugal transfer of an IncQ plasmid

The T-DNA transfer apparatus of Agrobacterium tumefaciens mediates the delivery of the T-DNA into plant cells, the transfer of the IncQ plasmid RSF1010 into plant cells, and the conjugal transfer of RSF1010 between Agrobacteria. We show in this report that the Agrobacterium-to-Agrobacterium conjugal transfer efficiencies of RSF1010 increase dramatically if the recipient strain, as well as the donor strain, carries a wild-type Ti plasmid and is capable of vir gene expression. Investigation of possible mechanisms that could account for this increased efficiency revealed that the VirB proteins encoded by the Ti plasmid were required. Although, with the exception of VirB1, all of the proteins that form the putative T-DNA transfer apparatus (VirB1–11, VirD4) are required for an Agrobacterium strain to serve as an RSF1010 donor, expression of only a subset of these proteins is required for the increase in conjugal transfer mediated by the recipient. Specifically, VirB5, 6, 11, and VirD4 are essential donor components but are dispensable for the increased recipient capacity. Defined point mutations in virB9 affected donor and recipient capacities to the same relative extent, suggesting that similar functions of VirB9 are important in both of these contexts.

Exploring the origins of topological frustration: Design of a minimally frustrated model of fragment B of protein A

Topological frustration in an energetically unfrustrated off-lattice model of the helical protein fragment B of protein A from Staphylococcus aureus was investigated. This Gō-type model exhibited thermodynamic and kinetic signatures of a well-designed two-state folder with concurrent collapse and folding transitions and single exponential kinetics at the transition temperature. Topological frustration is determined in the absence of energetic frustration by the distribution of Fersht φ values. Topologically unfrustrated systems present a unimodal distribution sharply peaked at intermediate φ, whereas highly frustrated systems display a bimodal distribution peaked at low and high φ values. The distribution of φ values in protein A was determined both thermodynamically and kinetically. Both methods yielded a unimodal distribution centered at φ = 0.3 with tails extending to low and high φ values, indicating the presence of a small amount of topological frustration. The contacts with high φ values were located in the turn regions between helices I and II and II and III, intimating that these hairpins are in large part required in the transition state. Our results are in good agreement with all-atom simulations of protein A, as well as lattice simulations of a three- letter code 27-mer (which can be compared with a 60-residue helical protein). The relatively broad unimodal distribution of φ values obtained from the all-atom simulations and that from the minimalist model for the same native fold suggest that the structure of the transition state ensemble is determined mostly by the protein topology and not energetic frustration.

Targeted overexpression of protein kinase C β2 isoform in myocardium causes cardiomyopathy

Increased cardiovascular mortality occurs in diabetic patients with or without coronary artery disease and is attributed to the presence of diabetic cardiomyopathy. One potential mechanism is hyperglycemia that has been reported to activate protein kinase C (PKC), preferentially the β isoform, which has been associated with the development of micro- and macrovascular pathologies in diabetes mellitus. To establish that the activation of the PKCβ isoform can cause cardiac dysfunctions, we have established lines of transgenic mice with the specific overexpression of PKCβ2 isoform in the myocardium. These mice overexpressed the PKCβ2 isoform transgene by 2- to 10-fold as measured by mRNA, and proteins exhibited left ventricular hypertrophy, cardiac myocyte necrosis, multifocal fibrosis, and decreased left ventricular performance without vascular lesions. The severity of the phenotypes exhibited gene dose-dependence. Up-regulation of mRNAs for fetal type myosin heavy chain, atrial natriuretic factor, c-fos, transforming growth factor, and collagens was also observed. Moreover, treatment with a PKCβ-specific inhibitor resulted in functional and histological improvement. These findings have firmly established that the activation of the PKCβ2 isoform can cause specific cardiac cellular and functional changes leading to cardiomyopathy of diabetic or nondiabetic etiology.

Loss of protein kinase C inhibition in the β-T594M variant of the amiloride-sensitive Na+ channel

We previously reported the presence of a novel variant (β-T594M) of the amiloride-sensitive Na+ channel (ASSC) in which the threonine residue at position 594 in the β-subunit has been replaced by a methionine residue. Electrophysiological studies of the ASSC on Epstein–Barr virus (EBV)-transformed lymphocytes carrying this variant showed that the 8-(4-chlorophenylthio) adenosine 3′:5′-cyclic monophosphate (8cpt-cAMP)-induced responses were enhanced when compared to wild-type EBV-transformed lymphocytes. Furthermore, in wild-type EBV-transformed cells, the 8cpt-cAMP-induced response was totally blocked by the phorbol ester, phorbol 12-myristate 13-acetate (PMA). This inhibitory effect of PMA was blocked by a protein kinase C inhibitor, chelerythrine. We now have identified individuals who are homozygous for this variant, and showed that PMA had no effect on the 8cpt-cAMP-induced responses in the EBV-transformed lymphocytes from such individuals. Cells heterozygous for this variant showed mixed responses to PMA, with the majority of cells partially inhibited by PMA. Our results demonstrate that an alteration in a single amino acid residue in the β-subunit of the ASSC can lead to a total loss of inhibition to PMA, and establish the β-subunit as having an important role in conferring a regulatory effect on the ASSC of lymphocytes.

Activation of protein kinase C by tyrosine phosphorylation in response to H2O2

Protein kinase C (PKC) isoforms, α, βI, and γ of cPKC subgroup, δ and ɛ of nPKC subgroup, and ζ of aPKC subgroup, were tyrosine phosphorylated in COS-7 cells in response to H2O2. These isoforms isolated from the H2O2-treated cells showed enhanced enzyme activity to various extents. The enzymes, PKC α and δ, recovered from the cells were independent of lipid cofactors for their catalytic activity. Analysis of mutated molecules of PKC δ showed that tyrosine residues, which are conserved in the catalytic domain of the PKC family, are critical for PKC activation induced by H2O2. These results suggest that PKC isoforms can be activated through tyrosine phosphorylation in a manner unrelated to receptor-coupled hydrolysis of inositol phospholipids.

B cell receptor-associated protein α4 displays rapamycin-sensitive binding directly to the catalytic subunit of protein phosphatase 2A

Recently, TAP42 was isolated as a high copy suppressor of sit4−, a yeast phosphatase related to protein phosphatase 2A (PP2A). TAP42 is related to the murine α4 protein, which was discovered independently by its association with Ig-α in the B cell receptor complex. Herein we show that a glutathione S-transferase (GST)–α4 fusion protein bound the catalytic subunit (C) of human PP2A from monomeric or multimeric preparations of PP2A in a “pull-down” assay. In an overlay assay, the GST–α4 protein bound to the phosphorylated and unphosphorylated forms of C that were separated in two-dimensional gels and immobilized on filters. The results show direct and exclusive binding of α4 to C. This is unusual because all known regulatory B subunits, or tumor virus antigens, bind stably only to the AC dimer of PP2A. The α4–C form of PP2A had an increased activity ratio compared with the AC form of PP2A when myelin basic protein phosphorylated by mitogen-activated protein kinase and phosphorylase a were used as substrates. Recombinant α4 cleaved from GST was phosphorylated by p56lck tyrosine kinase and protein kinase C. A FLAG-tagged α4 expressed in COS7 cells was recovered as a protein containing phosphoserine and coimmunoprecipitated with the C but not the A subunit of PP2A. Treatment of cells with rapamycin prevented the association of PP2A with FLAG-α4. The results reveal a novel heterodimer α4–C form of PP2A that may be involved in rapamycin-sensitive signaling pathways in mammalian cells.

Protein folding kinetics exhibit an Arrhenius temperature dependence when corrected for the temperature dependence of protein stability

The anomalous temperature dependence of protein folding has received considerable attention. Here we show that the temperature dependence of the folding of protein L becomes extremely simple when the effects of temperature on protein stability are corrected for; the logarithm of the folding rate is a linear function of 1/T on constant stability contours in the temperature–denaturant plane. This convincingly demonstrates that the anomalous temperature dependence of folding derives from the temperature dependence of the interactions that stabilize proteins, rather than from the super Arrhenius temperature dependence predicted for the configurational diffusion constant on a rough energy landscape. However, because of the limited temperature range accessible to experiment, the results do not rule out models with higher order temperature dependences. The significance of the slope of the stability-corrected Arrhenius plots is discussed.

The efficiency of propulsion by a rotating flagellum

[At very low Reynolds number, the regime in which fluid dynamics is governed by Stokes equations, a helix that translates along its axis under an external force but without an external torque will necessarily rotate. By the linearity of the Stokes equations, the same helix that is caused to rotate due to an external torque will necessarily translate. This is the physics that underlies the mechanism of flagellar propulsion employed by many microorganisms. Here, I examine the linear relationships between forces and torques and translational and angular velocities of helical objects to understand the nature of flagellar propulsion.]

Identification and characterization of a conserved family of protein serine/threonine phosphatases homologous to Drosophila retinal degeneration C (rdgC)

The Drosophila retinal degeneration C (rdgC) gene encodes an unusual protein serine/threonine phosphatase in that it contains at least two EF-hand motifs at its carboxy terminus. By a combination of large-scale sequencing of human retina cDNA clones and searches of expressed sequence tag and genomic DNA databases, we have identified two sequences in mammals [Protein Phosphatase with EF-hands-1 and 2 (PPEF-1 and PPEF-2)] and one in Caenorhabditis elegans (PPEF) that closely resemble rdgC. In the adult, PPEF-2 is expressed specifically in retinal rod photoreceptors and the pineal. In the retina, several isoforms of PPEF-2 are predicted to arise from differential splicing. The isoform that most closely resembles rdgC is localized to rod inner segments. Together with the recently described localization of PPEF-1 transcripts to primary somatosensory neurons and inner ear cells in the developing mouse, these data suggest that the PPEF family of protein serine/threonine phosphatases plays a specific and conserved role in diverse sensory neurons.

On roads not taken in the evolution of protein catalysts: Antibody steroid isomerases that use an enamine mechanism

Reactive immunization has emerged as a new tool for the study of biological catalysis. A powerful application resulted in catalytic antibodies that use an enamine mechanism akin to that used by the class I aldolases. With regard to the evolution of enzyme mechanisms, we investigated the utility of an enamine pathway for the allylic rearrangement exemplified by Δ5-3-ketosteroid isomerase (KSI; EC 5.3.3.1). Our aldolase antibodies were found to catalyze the isomerization of both steroid model compounds and steroids. The kinetic and chemical studies showed that the antibodies afforded rate accelerations up to a factor of 104 by means of an enamine mechanism in which imine formation was the rate-determining step. In light of our observations and the enzyme studies by other workers, we suggest that an enamine pathway could have been an early, viable KSI mechanism. Although this pathway is amenable to optimization for increased catalytic power, it appears that certain factors precluded its evolution in known KSI enzymes.

A structural census of the current population of protein sequences

We examine the occurrence of the ≈300 known protein folds in different groups of organisms. To do this, we characterize a large fraction of the currently known protein sequences (≈140,000) in structural terms, by matching them to known structures via sequence comparison (or by secondary-structure class prediction for those without structural homologues). Overall, we find that an appreciable fraction of the known folds are present in each of the major groups of organisms (e.g., bacteria and eukaryotes share 156 of 275 folds), and most of the common folds are associated with many families of nonhomologous sequences (i.e., >10 sequence families for each common fold). However, different groups of organisms have characteristically distinct distributions of folds. So, for instance, some of the most common folds in vertebrates, such as globins or zinc fingers, are rare or absent in bacteria. Many of these differences in fold usage are biologically reasonable, such as the folds of metabolic enzymes being common in bacteria and those associated with extracellular transport and communication being common in animals. They also have important implications for database-based methods for fold recognition, suggesting that an unknown sequence from a plant is more likely to have a certain fold (e.g., a TIM barrel) than an unknown sequence from an animal.

Long-term expression of protein kinase C in adult mouse hearts improves postischemic recovery

Activation of protein kinase C (PKC) protects the heart from ischemic injury; however, its mechanism of action is unknown, in part because no model for chronic activation of PKC has been available. To test whether chronic, mild elevation of PKC activity in adult mouse hearts results in myocardial protection during ischemia or reperfusion, hearts isolated from transgenic mice expressing a low level of activated PKCβ throughout adulthood (β-Tx) were compared with control hearts before ischemia, during 12 or 28 min of no-flow ischemia, and during reperfusion. Left-ventricular-developed pressure in isolated isovolumic hearts, normalized to heart weight, was similar in the two groups at baseline. However, recovery of contractile function was markedly improved in β-Tx hearts after either 12 (97 ± 3% vs. 69 ± 4%) or 28 min of ischemia (76 ± 8% vs. 48 ± 3%). Chelerythrine, a PKC inhibitor, abolished the difference between the two groups, indicating that the beneficial effect was PKC-mediated. 31P NMR spectroscopy was used to test whether modification of intracellular pH and/or preservation of high-energy phosphate levels during ischemia contributed to the cardioprotection in β-Tx hearts. No difference in intracellular pH or high-energy phosphate levels was found between the β-Tx and control hearts at baseline or during ischemia. Thus, long-term modest increase in PKC activity in adult mouse hearts did not alter baseline function but did lead to improved postischemic recovery. Furthermore, our results suggest that mechanisms other than reduced acidification and preservation of high-energy phosphate levels during ischemia contribute to the improved recovery.

Conformational influences of glycosylation of a peptide: A possible model for the effect of glycosylation on the rate of protein folding

Improved strategies for synthesis make it possible to expand the range of glycopeptides available for detailed conformational studies. The glycopeptide 1 was synthesized using a new solid phase synthesis of carbohydrates and a convergent coupling to peptide followed by deprotection. Its conformational properties were subjected to NMR analysis and compared with a control peptide 2 prepared by conventional solid phase methods. Whereas peptide 2 fails to manifest any appreciable secondary structure, the glycopeptide 1 does show considerable conformational bias suggestive of an equilibrium between an ordered and a random state. The implications of this ordering effect for the larger issue of protein folding are considered.