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.

cDNA cloning of Batis maritima methyl chloride transferase and purification of the enzyme

Methyl chloride transferase catalyzes the synthesis of methyl chloride from S-adenosine-l-methionine and chloride ion. This enzyme has been purified 2,700-fold to homogeneity from Batis maritima, a halophytic plant that grows abundantly in salt marshes. The purification of the enzyme was accomplished by a combination of ammonium sulfate fractionation, column chromatography on Sephadex G100 and adenosine-agarose, and TSK-250 size-exclusion HPLC. The purified enzyme exhibits a single band on SDS/PAGE with a molecular mass of approximately 22.5 kDa. The molecular mass of the purified enzyme was 22,474 Da as determined by matrix-associated laser desorption ionization mass spectrometry. The methylase can function in either a monomeric or oligomeric form. A 32-aa sequence of an internal fragment of the methylase was determined (GLVPGCGGGYDVVAMANPER FMVGLDIXENAL, where X represents unknown residue) by Edman degradation, and a full-length cDNA of the enzyme was obtained by rapid amplification of cDNA ends–PCR amplification of cDNA oligonucleotides. The cDNA gene contains an ORF of 690 bp encoding an enzyme of 230 aa residues having a predicted molecular mass of 25,761 Da. The disparity between the observed and calculated molecular mass suggests that the methylase undergoes posttranslational cleavage, possibly during purification. Sequence homologies suggest that the B. maritima methylase defines a new family of plant methyl transferases. A possible function for this novel methylase in halophytic plants is discussed.

Binocular vision in insects: How mantids solve the correspondence problem

Praying mantids use binocular cues to judge whether their prey is in striking distance. When there are several moving targets within their binocular visual field, mantids need to solve the correspondence problem. They must select between the possible pairings of retinal images in the two eyes so that they can strike at a single real target. In this study, mantids were presented with two targets in various configurations, and the resulting fixating saccades that precede the strike were analyzed. The distributions of saccades show that mantids consistently prefer one out of several possible matches. Selection is in part guided by the position and the spatiotemporal features of the target image in each eye. Selection also depends upon the binocular disparity of the images, suggesting that insects can perform local binocular computations. The pairing rules ensure that mantids tend to aim at real targets and not at “ghost” targets arising from false matches.

Neural mechanisms underlying binocular fusion and stereopsis: Position vs. phase

The visual system utilizes binocular disparity to discriminate the relative depth of objects in space. Since the striate cortex is the first site along the central visual pathways at which signals from the left and right eyes converge onto a single neuron, encoding of binocular disparity is thought to begin in this region. There are two possible mechanisms for encoding binocular disparity through simple cells in the striate cortex: a difference in receptive field (RF) position between the two eyes (RF position disparity) and a difference in RF profile between the two eyes (RF phase disparity). Although there have been studies supporting each of the two encoding mechanisms, both mechanisms have not been examined in a single study. Therefore, the relative roles of the two mechanisms have not been determined. To address this issue, we have mapped left and right eye RFs of simple cells in the cat’s striate cortex using binary m-sequence noise, and then we have estimated RF position and phase disparities. We find that RF position disparities are generally limited to small values that are not sufficient to encode large binocular disparities. In contrast, RF phase disparities cover a wide range of binocular disparities and exhibit dependencies on orientation and spatial frequency in a manner expected for a mechanism that encodes binocular disparity. These results indicate that binocular disparity is mainly encoded through RF phase disparity. However, RF position disparity may play a significant role for cells with high spatial frequency selectivity, which are constrained to small RF phase disparities.

Definition of HLA-DQ as a transplantation antigen

Recent studies have demonstrated the importance of recipient HLA-DRB1 allele disparity in the development of acute graft-versus-host disease (GVHD) after unrelated donor marrow transplantation. The role of HLA-DQB1 allele disparity in this clinical setting is unknown. To elucidate the biological importance of HLA-DQB1, we conducted a retrospective analysis of 449 HLA-A, -B, and -DR serologically matched unrelated donor transplants. Molecular typing of HLA-DRB1 and HLA-DQB1 alleles revealed 335 DRB1 and DQB1 matched pairs; 41 DRB1 matched and DQB1 mismatched pairs; 48 DRB1 mismatched and DQB1 matched pairs; and 25 DRB1 and DQB1 mismatched pairs. The conditional probabilities of grades III-IV acute GVHD were 0.42, 0.61, 0.55, and 0.71, respectively. The relative risk of acute GVHD associated with a single locus HLA-DQB1 mismatch was 1.8 (1.1, 2.7; P = 0.01), and the risk associated with any HLA-DQB1 and/or HLA-DRB1 mismatch was 1.6 (1.2, 2.2; P = 0.003). These results provide evidence that HLA-DQ is a transplant antigen and suggest that evaluation of both HLA-DQB1 and HLA-DRB1 is necessary in selecting potential donors.

Explaining the moon illusion

An old explanation of the moon illusion holds that various cues place the horizon moon at an effectively greater distance than the elevated moon. Although both moons have the same angular size, the horizon moon must be perceived as larger. More recent explanations hold that differences in accommodation or other factors cause the elevated moon to appear smaller. As a result of this illusory difference in size, the elevated moon appears to be more distant than the horizon moon. These two explanations, both based on the geometry of stereopsis, lead to two diametrically opposed hypotheses. That is, a depth interval at a long distance is associated with a smaller binocular disparity, whereas an equal depth interval at a smaller distance is associated with a larger disparity. We conducted experiments involving artificial moons and confirmed the hypothesis that the horizon moon is at a greater perceptual distance. Moreover, when a moon of constant angular size was moved closer it was also perceived as growing smaller, which is consistent with the older explanation. Although Emmert's law does not predict the size–distance relationship over long distances, we conclude that the horizon moon is perceived as larger because the perceptual system treats it as though it is much farther away. Finally, we observe that recent explanations substitute perceived size for angular size as a cue to distance. Thus, they imply that perceptions cause perceptions.

How the owl resolves auditory coding ambiguity

The barn owl (Tyto alba) uses interaural time difference (ITD) cues to localize sounds in the horizontal plane. Low-order binaural auditory neurons with sharp frequency tuning act as narrow-band coincidence detectors; such neurons respond equally well to sounds with a particular ITD and its phase equivalents and are said to be phase ambiguous. Higher-order neurons with broad frequency tuning are unambiguously selective for single ITDs in response to broad-band sounds and show little or no response to phase equivalents. Selectivity for single ITDs is thought to arise from the convergence of parallel, narrow-band frequency channels that originate in the cochlea. ITD tuning to variable bandwidth stimuli was measured in higher-order neurons of the owl’s inferior colliculus to examine the rules that govern the relationship between frequency channel convergence and the resolution of phase ambiguity. Ambiguity decreased as stimulus bandwidth increased, reaching a minimum at 2–3 kHz. Two independent mechanisms appear to contribute to the elimination of ambiguity: one suppressive and one facilitative. The integration of information carried by parallel, distributed processing channels is a common theme of sensory processing that spans both modality and species boundaries. The principles underlying the resolution of phase ambiguity and frequency channel convergence in the owl may have implications for other sensory systems, such as electrolocation in electric fish and the computation of binocular disparity in the avian and mammalian visual systems.

Efficient pyrophosphorolysis by a hepatitis B virus polymerase may be a primer-unblocking mechanism

Effective antiviral agents are thought to inhibit hepatitis B virus (HBV) DNA synthesis irreversibly by chain termination because reverse transcriptases (RT) lack an exonucleolytic activity that can remove incorporated nucleotides. However, since the parameters governing this inhibition are poorly defined, fully delineating the catalytic mechanism of the HBV-RT promises to facilitate the development of antiviral drugs for treating chronic HBV infection. To this end, pyrophosphorolysis and pyrophosphate exchange, two nonhydrolytic RT activities that result in the removal of newly incorporated nucleotides, were characterized by using endogenous avian HBV replication complexes assembled in vivo. Although these activities are presumed to be physiologically irrelevant for every polymerase examined, the efficiency with which they are catalyzed by the avian HBV-RT strongly suggests that it is the first known polymerase to catalyze these reactions under replicative conditions. The ability to remove newly incorporated nucleotides during replication has important biological and clinical implications: these activities may serve a primer-unblocking function in vivo. Analysis of pyrophosphorolysis on chain-terminated DNA revealed that the potent anti-HBV drug β-l-(−)-2′,3′-dideoxy-3′-thiacytidine (3TC) was difficult to remove by pyrophosphorolysis, in contrast to ineffective chain terminators such as ddC. This disparity may account for the strong antiviral efficacy of 3TC versus that of ddC. The HBV-RT pyrophosphorolytic activity may therefore be a novel determinant of antiviral drug efficacy, and could serve as a target for future antiviral drug therapy. The strong inhibitory effect of cytoplasmic pyrophosphate concentrations on viral DNA synthesis may also partly account for the apparent slow rate of HBV genome replication.

Library outreach: addressing Utah's “Digital Divide”

A “Digital Divide” in information and technological literacy exists in Utah between small hospitals and clinics in rural areas and the larger health care institutions in the major urban area of the state. The goals of the outreach program of the Spencer S. Eccles Health Sciences Library at the University of Utah address solutions to this disparity in partnership with the National Network of Libraries of Medicine—Midcontinental Region, the Utah Department of Health, and the Utah Area Health Education Centers. In a circuit-rider approach, an outreach librarian offers classes and demonstrations throughout the state that teach information-access skills to health professionals. Provision of traditional library services to unaffiliated health professionals is integrated into the library's daily workload as a component of the outreach program. The paper describes the history, methodology, administration, funding, impact, and results of the program.

NMR evidence for the participation of a low-barrier hydrogen bond in the mechanism of delta 5-3-ketosteroid isomerase.

Delta 5-3-Ketosteroid isomerase (EC 5.3.3.1) promotes an allylic rearrangement involving intramolecular proton transfer via a dienolic intermediate. This enzyme enhances the catalytic rate by a factor of 10(10). Two residues, Tyr-14, the general acid that polarizes the steroid 3-carbonyl group and facilitates enolization, and Asp-38 the general base that abstracts and transfers the 4 beta-proton to the 6 beta-position, contribute 10(4.7) and 10(5.6) to the rate increase, respectively. A major mechanistic enigma is the huge disparity between the pKa values of the catalytic groups and their targets. Upon binding of an analog of the dienolate intermediate to isomerase, proton NMR detects a highly deshielded resonance at 18.15 ppm in proximity to aromatic protons, and with a 3-fold preference for protium over deuterium (fractionation factor, phi = 0.34), consistent with formation of a short, strong (low-barrier) hydrogen bond to Tyr-14. The strength of this hydrogen bond is estimated to be at least 7.1 kcal/mol. This bond is relatively inaccessible to bulk solvent and is pH insensitive. Low-barrier hydrogen bonding of Tyr-14 to the intermediate, in conjunction with the previously demonstrated tunneling contribution to the proton transfer by Asp-38, provide a plausible and quantitative explanation for the high catalytic power of this isomerase.

Radiation target analysis of RNA.

Ribozymes are polynucleotide molecules with intrinsic catalytic activity, capable of cleaving nucleic acid substrates. Large RNA molecules were synthesized containing a hammerhead ribozyme moiety of 52 nucleotides linked to an inactive leader sequence, for total lengths of either 262 or 1226 nucleotides. Frozen RNAs were irradiated with high energy electrons. Surviving ribozyme activity was determined using the ability of the irradiated ribozymes to cleave a labeled substrate. The amount of intact RNA remaining was determined from the same irradiated samples by scanning the RNA band following denaturing gel electrophoresis. Radiation target analyses of these data revealed a structural target size of 80 kDa and a ribozyme activity target size of 15 kDa for the smaller ribozyme, and 319 kDa and 16 kDa, respectively, for the larger ribozyme. The disparity in target size for activity versus structure indicates that, in contrast to proteins, there is no spread of radiation damage far from the primary site of ionization in RNA molecules. The smaller target size for activity indicates that only primary ionizations occurring in the specific active region are effective. This is similar to the case for oligosaccharides. We concluded that the presence of the ribose sugar in the polymer chain restricts radiation damage to a small region and prevents major energy transfer throughout the molecule. Radiation target analysis should be a useful technique for evaluating local RNA:RNA and RNA:protein interactions in vitro.

Binocular visual surface perception.

Binocular disparity, the differential angular separation between pairs of image points in the two eyes, is the well-recognized basis for binocular distance perception. Without denying disparity's role in perceiving depth, we describe two perceptual phenomena, which indicate that a wider view of binocular vision is warranted. First, we show that disparity can play a critical role in two-dimensional perception by determining whether separate image fragments should be grouped as part of a single surface or segregated as parts of separate surfaces. Second, we show that stereoscopic vision is not limited to the registration and interpretation of binocular disparity but that it relies on half-occluded points, visible to one eye and not the other, to determine the layout and transparency of surfaces. Because these half-visible points are coded by neurons carrying eye-of-origin information, we suggest that the perception of these surface properties depends on neural activity available at visual cortical area V1.

Deletion of high-avidity T cells by thymic epithelium.

Tolerance induction by thymic epithelium induces a state of so-called "split tolerance," characterized in vivo by tolerance and in vitro by reactivity to a given thymically expressed antigen. Using a model major histocompatibility complex class I antigen, H-2Kb (Kb), three mechanisms of thymic epithelium-induced tolerance were tested: induction of tolerance of tissue-specific antigens exclusively, selective inactivation of T helper cell-independent cytotoxic T lymphocytes, and deletion of high-avidity T cells. To this end, thymic anlagen from Kb-transgenic embryonic day 10 mouse embryos, taken before colonization by cells of hemopoietic origin, were grafted to nude mice. Tolerance by thymic epithelium was not tissue-specific, since Kb-bearing skin and spleen grafts were maintained indefinitely. Only strong priming in vivo could partially overcome the tolerant state and induce rejection of some skin grafts overexpressing transgenic Kb. Furthermore, the hypothesis that thymic epithelium selectively inactivates those T cells that reject skin grafts in a T helper-independent fashion could not be supported. Thus, when T-cell help was provided by a second skin graft bearing an additional major histocompatibility complex class II disparity, tolerance to the Kb skin graft was not broken. Finally, direct evidence could be obtained for the avidity model of thymic epithelium-induced negative selection, using Kb-specific T-cell receptor (TCR) transgenic mice. Thymic epithelium-grafted TCR transgenic mice showed a selective deletion of those CD8+ T cells with the highest density of the clonotypic TCR. These cells presumably represent the T cells with the highest avidity for Kb. We conclude that split tolerance induced by thymic epithelium was mediated by the deletion of those CD8+ T lymphocytes that have the highest avidity for antigen.