The photoisomerization of retinal in bacteriorhodopsin: Experimental evidence for a three-state model

The primary events in the all-trans to 13-cis photoisomerization of retinal in bacteriorhodopsin have been investigated with femtosecond time-resolved absorbance spectroscopy. Spectra measured over a broad range extending from 7000 to 22,400 cm−1 reveal features whose dynamics are inconsistent with a model proposed earlier to account for the highly efficient photoisomerization process. Emerging from this work is a new three-state model. Photoexcitation of retinal with visible light accesses a shallow well on the excited state potential energy surface. This well is bounded by a small barrier, arising from an avoided crossing that separates the Franck–Condon region from the nearby reactive region of the photoisomerization coordinate. At ambient temperatures, the reactive region is accessed with a time constant of ≈500 fs, whereupon the retinal rapidly twists and encounters a second avoided crossing region. The protein mediates the passage into the second avoided crossing region and thereby exerts control over the quantum yield for forming 13-cis retinal. The driving force for photoisomerization resides in the retinal, not in the surrounding protein. This view contrasts with an earlier model where photoexcitation was thought to access directly a reactive region of the excited-state potential and thereby drive the retinal to a twisted conformation within 100–200 fs.

A model for amplification of hair-bundle motion by cyclical binding of Ca2+ to mechanoelectrical-transduction channels

Amplification of auditory stimuli by hair cells augments the sensitivity of the vertebrate inner ear. Cell-body contractions of outer hair cells are thought to mediate amplification in the mammalian cochlea. In vertebrates that lack these cells, and perhaps in mammals as well, active movements of hair bundles may underlie amplification. We have evaluated a mathematical model in which amplification stems from the activity of mechanoelectrical-transduction channels. The intracellular binding of Ca2+ to channels is posited to promote their closure, which increases the tension in gating springs and exerts a negative force on the hair bundle. By enhancing bundle motion, this force partially compensates for viscous damping by cochlear fluids. Linear stability analysis of a six-state kinetic model reveals Hopf bifurcations for parameter values in the physiological range. These bifurcations signal conditions under which the system’s behavior changes from a damped oscillatory response to spontaneous limit-cycle oscillation. By varying the number of stereocilia in a bundle and the rate constant for Ca2+ binding, we calculate bifurcation frequencies spanning the observed range of auditory sensitivity for a representative receptor organ, the chicken’s cochlea. Simulations using prebifurcation parameter values demonstrate frequency-selective amplification with a striking compressive nonlinearity. Because transduction channels occur universally in hair cells, this active-channel model describes a mechanism of auditory amplification potentially applicable across species and hair-cell types.

Simple model of protein folding kinetics.

A simple model of the kinetics of protein folding is presented. The reaction coordinate is the "correctness" of a configuration compared with the native state. The model has a gap in the energy spectrum, a large configurational entropy, a free energy barrier between folded and partially folded states, and a good thermodynamic folding transition. Folding kinetics is described by a master equation. The folding time is estimated by means of a local thermodynamic equilibrium assumption and then is calculated both numerically and analytically by solving the master equation. The folding time has a maximum near the folding transition temperature and can have a minimum at a lower temperature.

Population regulation in snowshoe hare and Canadian lynx: Asymmetric food web configurations between hare and lynx

The snowshoe hare and the Canadian lynx in the boreal forests of North America show 9- to 11-year density cycles. These are generally assumed to be linked to each other because lynx are specialist predators on hares. Based on time series data for hare and lynx, we show that the dominant dimensional structure of the hare series appears to be three whereas that of the lynx is two. The three-dimensional structure of the hare time series is hypothesized to be due to a three-trophic level model in which the hare may be seen as simultaneously regulated from below and above. The plant species in the hare diet appear compensatory to one another, and the predator species may, likewise, be seen as an internally compensatory guild. The lynx time series are, in contrast, consistent with a model of donor control in which their populations are regulated from below by prey availability. Thus our analysis suggests that the classic view of a symmetric hare–lynx interaction is too simplistic. Specifically, we argue that the classic food chain structure is inappropriate: the hare is influenced by many predators other than the lynx, and the lynx is primarily influenced by the snowshoe hare.

Folding thermodynamics of a model three-helix-bundle protein

The calculated folding thermodynamics of a simple off-lattice three-helix-bundle protein model under equilibrium conditions shows the experimentally observed protein transitions: a collapse transition, a disordered-to-ordered globule transition, a globule to native-state transition, and the transition from the active native state to a frozen inactive state. The cooperativity and physical origin of the various transitions are explored with a single “optimization” parameter and characterized with the Lindemann criterion for liquid versus solid-state dynamics. Below the folding temperature, the model has a simple free energy surface with a single basin near the native state; the surface is similar to that calculated from a simulation of the same three-helix-bundle protein with an all-atom representation [Boczko, E. M. & Brooks III, C. L. (1995) Science 269, 393–396].

The "allosteric three-site model" of elongation cannot be confirmed in a well-defined ribosome system from Escherichia coli.

For the functional role of the ribosomal tRNA exit (E) site, two different models have been proposed. It has been suggested that transient E-site binding of the tRNA leaving the peptidyl (P) site promotes elongation factor G (EF-G)-dependent translocation by lowering the energetic barrier of tRNA release [Lill, R., Robertson, J. M. & Wintermeyer, W. (1989) EMBO J. 8, 3933-3938]. The alternative "allosteric three-site model" [Nierhaus, K.H. (1990) Biochemistry 29, 4997-5008] features stable, codon-dependent tRNA binding to the E site and postulates a coupling between E and aminoacyl (A) sites that regulates the tRNA binding affinity of the two sites in an anticooperative manner. Extending our testing of the two conflicting models, we have performed translocation experiments with fully active ribosomes programmed with heteropolymeric mRNA. The results confirm that the deacylated tRNA released from the P site is bound to the E site in a kinetically labile fashion, and that the affinity of binding, i.e., the occupancy of the E site, is increased by Mg2+ or polyamines. At conditions of high E-site occupancy in the posttranslocation complex, filling the A site with aminoacyl-tRNA had no influence on the E site, i.e., there was no detectable anticooperative coupling between the two sites, provided that second-round translocation was avoided by removing EF-G. On the basis of these results, which are entirely consistent with our previous results, we consider the allosteric three-site model of elongation untenable. Rather, as proposed earlier, the E site-bound state of the leaving tRNA is a transient intermediate and, as such, is a mechanistic feature of the classic two-state model of the elongating ribosome.

Comparing in vitro, in situ, and in vivo experimental data in a three-dimensional model of mammalian cochlear mechanics

Normal mammalian hearing is refined by amplification of the motion of the cochlear partition. This partition, comprising the organ of Corti sandwiched between the basilar and tectorial membranes, contains the outer hair cells that are thought to drive this amplification process. Force generation by outer hair cells has been studied extensively in vitro and in situ, but, to understand cochlear amplification fully, it is necessary to characterize the role played by each of the components of the cochlear partition in vivo. Observations of cochlear partition motion in vivo are severely restricted by its inaccessibility and sensitivity to surgical trauma, so, for the present study, a computer model has been used to simulate the operation of the cochlea under different experimental conditions. In this model, which uniquely retains much of the three-dimensional complexity of the real cochlea, the motions of the basilar and tectorial membranes are fundamentally different during in situ- and in vivo-like conditions. Furthermore, enhanced outer hair cell force generation in vitro leads paradoxically to a decrease in the gain of the cochlear amplifier during sound stimulation to the model in vivo. These results suggest that it is not possible to extrapolate directly from experimental observations made in vitro and in situ to the normal operation of the intact organ in vivo.

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.

Exploring the folding free energy surface of a three-helix bundle protein

The multidimensional free energy surface for a small fast folding helical protein is explored based on first-principle calculations. The model represents the 46-residue segment from fragment B of staphylococcal protein A. The relationship between collapse and tertiary structure formation, and the order of collapse and secondary structure formation, are investigated. We find that the initial collapse process gives rise to a transition state with about 30% of the native tertiary structure and 50–70% of the native helix content. We also observe two distinct distributions of native helix in this collapsed state (Rg ≈ 12 Å), one with about 20% of the native helical hydrogen bonds, the other with near 70%. The former corresponds to a local minimum. The barrier from this metastable state to the native state is about 2 kBT. In the latter case, folding is essentially a downhill process involving topological assembly. In addition, the order of formation of secondary structure among the three helices is examined. We observe cooperative formation of the secondary structure in helix I and helix II. Secondary structure in helix III starts to form following the formation of certain secondary structure in both helix I and helix II. Comparisons of our results with those from theory and experiment are made.

Three-dimensional domain swapping in p13suc1 occurs in the unfolded state and is controlled by conserved proline residues

p13suc1 has two native states, a monomer and a domain-swapped dimer. We show that their folding pathways are connected by the denatured state, which introduces a kinetic barrier between monomer and dimer under native conditions. The barrier is lowered under conditions that speed up unfolding, thereby allowing, to our knowledge for the first time, a quantitative dissection of the energetics of domain swapping. The monomer–dimer equilibrium is controlled by two conserved prolines in the hinge loop that connects the exchanging domains. These two residues exploit backbone strain to specifically direct dimer formation while preventing higher-order oligomerization. Thus, the loop acts as a loaded molecular spring that releases tension in the monomer by adopting its alternative conformation in the dimer. There is an excellent correlation between domain swapping and aggregation, suggesting they share a common mechanism. These insights have allowed us to redesign the domain-swapping propensity of suc1 from a fully monomeric to a fully dimeric protein.

Photoactive yellow protein: A structural prototype for the three-dimensional fold of the PAS domain superfamily

PAS domains are found in diverse proteins throughout all three kingdoms of life, where they apparently function in sensing and signal transduction. Although a wealth of useful sequence and functional information has become recently available, these data have not been integrated into a three-dimensional (3D) framework. The very early evolutionary development and diverse functions of PAS domains have made sequence analysis and modeling of this protein superfamily challenging. Limited sequence similarities between the ∼50-residue PAS repeats and one region of the bacterial blue-light photosensor photoactive yellow protein (PYP), for which ground-state and light-activated crystallographic structures have been determined to high resolution, originally were identified in sequence searches using consensus sequence probes from PAS-containing proteins. Here, we found that by changing a few residues particular to PYP function, the modified PYP sequence probe also could select PAS protein sequences. By mapping a typical ∼150-residue PAS domain sequence onto the entire crystallographic structure of PYP, we show that the PAS sequence similarities and differences are consistent with a shared 3D fold (the PAS/PYP module) with obvious potential for a ligand-binding cavity. Thus, PYP appears to prototypically exhibit all the major structural and functional features characteristic of the PAS domain superfamily: the shared PAS/PYP modular domain fold of ∼125–150 residues, a sensor function often linked to ligand or cofactor (chromophore) binding, and signal transduction capability governed by heterodimeric assembly (to the downstream partner of PYP). This 3D PAS/PYP module provides a structural model to guide experimental testing of hypotheses regarding ligand-binding, dimerization, and signal transduction.

Three-dimensional model of the potyviral genome-linked protein.

The full sequence of the genome-linked viral protein (VPg) cistron located in the central part of potato virus Y (common strain) genome has been identified. The VPg gene codes for a protein of 188 amino acids, with significant homology to other known potyviral VPg polypeptides. A three-dimensional model structure of VPg is proposed on the basis of similarity of hydrophobic-hydrophilic residue distribution to the sequence of malate dehydrogenase of known crystal structure. The 5' end of the viral RNA can be fitted to interact with the protein through the exposed hydroxyl group of Tyr-64, in agreement with experimental data. The complex favors stereochemically the formation of a phosphodiester bond [5'-(O4-tyrosylphospho)adenylate] typical for representatives of picornavirus-like viruses. The chemical mechanisms of viral RNA binding to VPg are discussed on the basis of the model structure of protein-RNA complex.

DNA-mediated immunization in a transgenic mouse model of the hepatitis B surface antigen chronic carrier state.

Transgenic mice expressing the sequences coding for the envelope proteins of the hepatitis B virus (HBV) in the liver have been used as a model of the HBV chronic carrier state. We evaluated the possibility of inducing a specific immune response to the viral envelope antigens and thus potentially controlling chronic HBV infection. Using HBV-specific DNA-mediated immunization in this transgenic model, we show that the immune response induced after a single intramuscular injection of DNA resulted in the complete clearance of circulating hepatitis B surface antigen and in the long-term control of transgene expression in hepatocytes. This response does not involve a detectable cytopathic effect in the liver. Adoptive transfer of fractionated primed spleen cells from DNA-immunized mice shows that T cells are responsible for the down-regulation of HBV mRNA in the liver of transgenic mice. To our knowledge, this is the first demonstration of a potential immunotherapeutic application of DNA-mediated immunization against an infectious disease and raises the possibility of designing more effective ways of treating HBV chronic carriers.

A model for the transcriptional regulation of the CYP2B1/B2 gene in rat liver.

The phenobarbitone-responsive minimal promoter has been shown to lie between nt -179 and nt + 1 in the 5' (upstream) region of the CYP2B1/B2 gene in rat liver, on the basis of the drug responsiveness of the sequence linked to human growth hormone gene as reporter and targeted to liver as an asialoglycoprotein-DNA complex in vivo. Competition analyses of the nuclear protein-DNA complexes formed in gel shift assays with the positive (nt -69 to -98) and negative (nt -126 to -160) cis elements (PE and NE, respectively) identified within this region earlier indicate that the same protein may be binding to both the elements. The protein species purified on PE and NE affinity columns appear to be identical based on SDS/PAGE analysis, where it migrates as a protein of 26-28 kDa. Traces of a high molecular weight protein (94-100 kDa) are also seen in the preparation obtained after one round of affinity chromatography. The purified protein stimulates transcription of a minigene construct containing the 179 nt on the 5' side of the CYP2B1/B2 gene linked to the I exon in a cell-free system from liver nuclei. The purified protein can give rise to all the three complexes (I, II, and III) with the PE, just as the crude nuclear extract, under appropriate conditions. Manipulations in vitro indicate that the NE has a significantly higher affinity for the dephosphorylated form than for the phosphorylated form of the protein. The PE binds both forms. Phenobarbitone treatment of the animal leads to a significant increase in the phosphorylation of the 26- to 28-kDa and 94-kDa proteins in nuclear labeling experiments followed by isolation on a PE affinity column. We propose that the protein binding predominantly to the NE in the dephosphorylated state characterizes the basal level of transcription of the CYP2B1/B2 gene. Phenobarbitone treatment leads to phosphorylation of the protein, shifting the equilibrium toward binding to the PE. This can promote interaction with an upstream enhancer through other proteins such as the 94-kDa protein and leads to a significant activation of transcription.

Genetic and redox determinants of nitric oxide cytotoxicity in a Salmonella typhimurium model.

Paradoxically, nitric oxide (NO) has been found to exhibit cytotoxic, antiproliferative, or cytoprotective activity under different conditions. We have utilized Salmonella mutants deficient in antioxidant defenses or peptide transport to gain insights into NO actions. Comparison of three NO donor compounds reveals distinct and independent cellular responses associated with specific redox forms of NO. The peroxynitrite (OONO-) generator 3-morpholinosydnonimine hydrochloride mediates oxygen-dependent Salmonella killing, whereas S-nitrosoglutathione (GSNO) causes oxygen-independent cytostasis, and the NO. donor diethylenetriamine-nitric oxide adduct has no antibacterial activity. GSNO has the greatest activity for stationary cells, a characteristic relevant to latent or intracellular pathogens. Moreover, the cytostatic activity of GSNO may best correlate with antiproliferative or antimicrobial effects of NO, which are unassociated with overt cell injury. dpp mutants defective in active dipeptide transport are resistant to GSNO, implicating heterolytic NO+ transfer rather than homolytic NO. release in the mechanism of cytostasis. This transport system may provide a specific pathway for GSNO-mediated signaling in biological systems. The redox state and associated carrier molecules are critical determinants of NO activity.