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.

Neuronal nicotinic threonine-for-leucine 247 α7 mutant receptors show different gating kinetics when activated by acetylcholine or by the noncompetitive agonist 5-hydroxytryptamine

Mutation of the highly conserved leucine residue (Leu-247) converts 5-hydroxytryptamine (5HT) from an antagonist into an agonist of neuronal homomeric α7 nicotinic acetylcholine receptor expressed in Xenopus oocytes. We show here that acetylcholine (AcCho) activates two classes of single channels with conductances of 44 pS and 58 pS, similar to those activated by 5HT. However, the mean open time of AcCho-gated ion channels (11 ms) is briefer than that of 5HT-gated ion channels (18 ms). Furthermore, whereas the open time of AcCho channels lengthens with hyperpolarization, that of 5HT channels is decreased. In voltage-clamped oocytes, the apparent affinity of the α7 mutant receptor for 5HT is not modified by the presence of dihydro-β-erythroidine, which acts on the AcCho binding site in a competitive manner. This indicates a noncompetitive action of 5HT on nicotinic acetylcholine receptors. Considered together, our findings show that AcCho gates α7 mutant channels with similar conductance but with different kinetic profile than the channels gated by 5HT, suggesting that the two agonists act on different docking sites. These results will help to understand the crosstalk between cholinergic and serotonergic systems in the central nervous system.

Voltage dependence of the pattern and frequency of discrete Ca2+ release events after brief repriming in frog skeletal muscle

Applying a brief repolarizing pre-pulse to a depolarized frog skeletal muscle fiber restores a small fraction of the transverse tubule membrane voltage sensors from the inactivated state. During a subsequent depolarizing test pulse we detected brief, highly localized elevations of myoplasmic Ca2+ concentration (Ca2+ “sparks”) initiated by restored voltage sensors in individual triads at all test pulse voltages. The latency histogram of these events gives the gating pattern of the sarcoplasmic reticulum (SR) calcium release channels controlled by the restored voltage sensors. Both event frequency and clustering of events near the start of the test pulse increase with test pulse depolarization. The macroscopic SR calcium release waveform, obtained from the spark latency histogram and the estimated open time of the channel or channels underlying a spark, exhibits an early peak and rapid marked decline during large depolarizations. For smaller depolarizations, the release waveform exhibits a smaller peak and a slower decline. However, the mean use time and mean amplitude of the individual sparks are quite similar at all test depolarizations and at all times during a given depolarization, indicating that the channel open times and conductances underlying sparks are essentially independent of voltage. Thus, the voltage dependence of SR Ca2+ release is due to changes in the frequency and pattern of occurrence of individual, voltage-independent, discrete release events.

Empirical laws of survival and evolution: Their universality and implications

Presented analysis of human and fly life tables proves that with the specified accuracy their entire survival and mortality curves are uniquely determined by a single point (e.g., by the birth mortality q0), according to the law, which is universal for species as remote as humans and flies. Mortality at any age decreases with the birth mortality q0. According to life tables, in the narrow vicinity of a certain q0 value (which is the same for all animals of a given species, independent of their living conditions), the curves change very rapidly and nearly simultaneously for an entire population of different ages. The change is the largest in old age. Because probability to survive to the mean reproductive age quantifies biological fitness and evolution, its universal rapid change with q0 (which changes with living conditions) manifests a new kind of an evolutionary spurt of an entire population. Agreement between theoretical and life table data is explicitly seen in the figures. Analysis of the data on basic metabolism reduces it to the maximal mean lifespan (for animals from invertebrates to mammals), or to the maximal mean fission time (for bacteria), and universally scales them with the total number of body atoms only. Phenomenological origin of this unification and universality of metabolism, survival, and evolution is suggested. Their implications and challenges are discussed.

Ca2+-induced inhibition of the cardiac Ca2+ channel depends on calmodulin

Ca2+-induced inhibition of α1C voltage-gated Ca2+ channels is a physiologically important regulatory mechanism that shortens the mean open time of these otherwise long-lasting high-voltage-activated channels. The mechanism of action of Ca2+ has been a matter of some controversy, as previous studies have proposed the involvement of a putative Ca2+-binding EF hand in the C terminus of α1C and/or a sequence downstream from this EF-hand motif containing a putative calmodulin (CaM)-binding IQ motif. Previously, using site directed mutagenesis, we have shown that disruption of the EF-hand motif does not remove Ca2+ inhibition. We now show that the IQ motif binds CaM and that disruption of this binding activity prevents Ca2+ inhibition. We propose that Ca2+ entering through the voltage-gated pore binds to CaM and that the Ca/CaM complex is the mediator of Ca2+ inhibition.

Nociception in cyclooxygenase isozyme-deficient mice

Prostaglandins formed by cyclooxygenase-1 (COX-1) or COX-2 produce hyperalgesia in sensory nerve endings. To assess the relative roles of the two enzymes in pain processing, we compared responses of COX-1- or COX-2-deficient homozygous and heterozygous mice with wild-type controls in the hot plate and stretching tests for analgesia. Preliminary observational studies determined that there were no differences in gross parameters of behavior between the different groups. Surprisingly, on the hot plate (55°C), the COX-1-deficient heterozygous groups showed less nociception, because mean reaction time was longer than that for controls. All other groups showed similar reaction times. In the stretching test, there was less nociception in COX-1-null and COX-1-deficient heterozygotes and also, unexpectedly, in female COX-2-deficient heterozygotes, as shown by a decreased number of writhes. Measurements of mRNA levels by reverse transcription–PCR demonstrated a compensatory increase of COX-1 mRNA in spinal cords of COX-2-null mice but no increase in COX-2 mRNA in spinal cords of COX-1-null animals. Thus, compensation for the absence of COX-1 may not involve increased expression of COX-2, whereas up-regulation of COX-1 in the spinal cord may compensate for the absence of COX-2. The longer reaction times on the hot plate of COX-1-deficient heterozygotes are difficult to explain, because nonsteroid anti-inflammatory drugs have no analgesic action in this test. Reduction in the number of writhes of the COX-1-null and COX-1-deficient heterozygotes may be due to low levels of COX-1 at the site of stimulation with acetic acid. Thus, prostaglandins made by COX-1 mainly are involved in pain transmission in the stretching test in both male and female mice, whereas those made by COX-2 also may play a role in the stretching response in female mice.

Indirect recognition of allopeptides promotes the development of cardiac allograft vasculopathy

Graft loss from chronic rejection has become the major obstacle to the long-term success of whole organ transplantation. In cardiac allografts, chronic rejection is manifested as a diffuse and accelerated form of arteriosclerosis, termed cardiac allograft vasculopathy. It has been suggested that T-cell recognition of processed alloantigens (allopeptides) presented by recipient antigen-presenting cells through the indirect pathway of allorecognition plays a critical role in the development and progression of chronic rejection. However, definitive preclinical evidence to support this hypothesis is lacking. To examine the role of indirect allorecognition in a clinically relevant large animal model of cardiac allograft vasculopathy, we immunized MHC inbred miniature swine with synthetic polymorphic peptides spanning the α1 domain of an allogeneic donor-derived swine leukocyte antigen class I gene. Pigs immunized with swine leukocyte antigen class I allopeptides showed in vitro proliferative responses and in vivo delayed-type hypersensitivity responses to the allogeneic peptides. Donor MHC class I disparate hearts transplanted into peptide-immunized cyclosporine-treated pigs not only rejected faster than unimmunized cyclosporine-treated controls (mean survival time = 5.5 +/−1.7 vs. 54.7 +/−3.8 days, P < 0.001), but they also developed obstructive fibroproliferative coronary artery lesions much earlier than unimmunized controls (<9 vs. >30 days). These results definitively link indirect allorecognition and cardiac allograft vasculopathy.

Insulin promotes rapid delivery of N-methyl-d- aspartate receptors to the cell surface by exocytosis

Insulin potentiates N-methyl-d-aspartate receptors (NMDARs) in neurons and Xenopus oocytes expressing recombinant NMDARs. The present study shows that insulin induced (i) an increase in channel number times open probability (nPo) in outside-out patches excised from Xenopus oocytes, with no change in mean open time, unitary conductance, or reversal potential, indicating an increase in n and/or Po; (ii) an increase in charge transfer during block of NMDA-elicited currents by the open channel blocker MK-801, indicating increased number of functional NMDARs in the cell membrane with no change in Po; and (iii) increased NR1 surface expression, as indicated by Western blot analysis of surface proteins. Botulinum neurotoxin A greatly reduced insulin potentiation, indicating that insertion of new receptors occurs via SNARE-dependent exocytosis. Thus, insulin potentiation occurs via delivery of new channels to the plasma membrane. NMDARs assembled from mutant subunits lacking all known sites of tyrosine and serine/threonine phosphorylation in their carboxyl-terminal tails exhibited robust insulin potentiation, suggesting that insulin potentiation does not require direct phosphorylation of NMDAR subunits. Because insulin and insulin receptors are localized to glutamatergic synapses in the hippocampus, insulin-regulated trafficking of NMDARs may play a role in synaptic transmission and plasticity, including long-term potentiation.

A Calcium-Selective Channel from Root-Tip Endomembranes of Garden Cress1

A Ca2+ channel from root-tip endomembranes of garden cress (Lepidium sativum L.) (LCC1) was characterized using the planar lipid-bilayer technique. Investigation of single-channel recordings revealed that LCC1 is voltage gated and strongly rectifying. In symmetrical 50 mm CaCl2 solutions, the single-channel conductance was 24 picosiemens. LCC1 showed a moderate selectivity for Ca2+ over K+ (9.4:1) and was permeable for a range of divalent cations (Ca2+, Ba2+, and Sr2+). In contrast to Bryonia dioica Ca2+ channel 1, a Ca2+-selective channel from the endoplasmic reticulum of touch-sensitive tendrils, LCC1 showed no bursting channel activity and had a low open probability and mean open time (2.83 ms at 50 mV). Inhibitor studies demonstrated that LCC1 is blocked by micromolar concentrations of erythrosin B (inhibitor concentration for 50% inhibition [IC50] = 1.8 μm) and the trivalent cations La3+ (IC50 = 5 μm) and Gd3+ (IC50 = 10 μm), whereas verapamil showed no blocking effect. LCC1 may play an important role in the regulation of the cytoplasmic free Ca2+ concentration in root-tip and/or root-cap cells. The question of whether this ion channel is part of the gravitropic signal transduction pathway deserves further investigation.