A local electrostatic change is the cause of the large-scale protein conformation shift in bacteriorhodopsin

During light-driven proton transport bacteriorhodopsin shuttles between two protein conformations. A large-scale structural change similar to that in the photochemical cycle is produced in the D85N mutant upon raising the pH, even without illumination. We report here that (i) the pKa values for the change in crystallographic parameters and for deprotonation of the retinal Schiff base are the same, (ii) the retinal isomeric configuration is nearly unaffected by the protein conformation, and (iii) preventing rotation of the C13—C14 double bond by replacing the retinal with an all-trans locked analogue makes little difference to the Schiff base pKa. We conclude that the direct cause of the conformational shift is destabilization of the structure upon loss of interaction of the positively charged Schiff base with anionic residues that form its counter-ion.

Structure and function in rhodopsin: Rhodopsin mutants with a neutral amino acid at E134 have a partially activated conformation in the dark state*

The Glu-134–Arg-135 residues in rhodopsin, located near the cytoplasmic end of the C helix, are involved in G protein binding, or activation, or both. Furthermore, the charge-neutralizing mutation Glu-134 to Gln-134 produces hyperactivity in the activated state and produces constitutive activity in opsin. The Glu/Asp-Arg charge pair is highly conserved in equivalent positions in other G protein-coupled receptors. To investigate the structural consequences of charge-neutralizing mutations at Glu-134 and Arg-135 in rhodopsin, single spin-labeled side chains were introduced at sites in the cytoplasmic domains of helices C (140), E (227), F (250), or G (316) to serve as “molecular sensors” of the local helix bundle conformation. In each of the spin-labeled rhodopsins, a Gln substitution was introduced at either Glu-134 or Arg-135, and the electron paramagnetic resonance spectrum of the spin label was used to monitor the structural response of the helix bundle. The results indicate that a Gln substitution at Glu-134 induces a photoactivated conformation around helices C and G even in the dark state, an observation of potential relevance to the hyperactivity and constitutive activity of the mutant. In contrast, little change is induced in helix F, which has been shown to undergo a dominant motion upon photoactivation. This result implies that the multiple helix motions accompanying photoactivation are not strongly coupled and can be induced to take place independently. Gln substitution at Arg-135 produces only minor structural changes in the dark- or light-activated conformation, suggesting that this residue is not a determinant of structure in the regions investigated, although it may be functionally important.

Transition-state structure as a unifying basis in protein-folding mechanisms: Contact order, chain topology, stability, and the extended nucleus mechanism

I attempt to reconcile apparently conflicting factors and mechanisms that have been proposed to determine the rate constant for two-state folding of small proteins, on the basis of general features of the structures of transition states. Φ-Value analysis implies a transition state for folding that resembles an expanded and distorted native structure, which is built around an extended nucleus. The nucleus is composed predominantly of elements of partly or well-formed native secondary structure that are stabilized by local and long-range tertiary interactions. These long-range interactions give rise to connecting loops, frequently containing the native loops that are poorly structured. I derive an equation that relates differences in the contact order of a protein to changes in the length of linking loops, which, in turn, is directly related to the unfavorable free energy of the loops in the transition state. Kinetic data on loop extension mutants of CI2 and α-spectrin SH3 domain fit the equation qualitatively. The rate of folding depends primarily on the interactions that directly stabilize the nucleus, especially those in native-like secondary structure and those resulting from the entropy loss from the connecting loops, which vary with contact order. This partitioning of energy accounts for the success of some algorithms that predict folding rates, because they use these principles either explicitly or implicitly. The extended nucleus model thus unifies the observations of rate depending on both stability and topology.

Biodiversity of Costa Rican salamanders: Implications of high levels of genetic differentiation and phylogeographic structure for species formation

Although salamanders are characteristic amphibians in Holarctic temperate habitats, in tropical regions they have diversified evolutionarily only in tropical America. An adaptive radiation centered in Middle America occurred late in the history of a single clade, the supergenus Bolitoglossa (Plethodontidae), and large numbers of species now occur in diverse habitats. Sublineages within this clade decrease in number from the northern to southern parts of Middle America, and in Costa Rica, there are but three. Despite this phylogenetic constraint, Costa Rica has many species; the number of salamander species on one local elevational transect in the Cordillera de Talamanca may be the largest for any such transect in the world. Extraordinary variation in sequences of the mitochondrial gene cytochrome b within a clade of the genus Bolitoglossa in Costa Rica reveals strong phylogeographic structure within a single species, Bolitoglossa pesrubra. Allozymic variation in 19 proteins reveals a pattern largely concordant with the mitochondrial DNA phylogeography. More species exist than are currently recognized. Diversification occurs in restricted geographic areas and involves sharp geographic and elevational differentiation and zonation. In their degree of genetic differentiation at a local scale, these species of the deep tropics exceed the known variation of extratropical salamanders, which also differ in being less restricted in elevational range. Salamanders display “tropicality” in that although speciose, they are usually local in distribution and rare. They display strong ecological and physiological differentiation that may contribute importantly to morphological divergence and species formation.

Crystal structure of the Holliday junction DNA in complex with a single RuvA tetramer

In the major pathway of homologous DNA recombination in prokaryotic cells, the Holliday junction intermediate is processed through its association with RuvA, RuvB, and RuvC proteins. Specific binding of the RuvA tetramer to the Holliday junction is required for the RuvB motor protein to be loaded onto the junction DNA, and the RuvAB complex drives the ATP-dependent branch migration. We solved the crystal structure of the Holliday junction bound to a single Escherichia coli RuvA tetramer at 3.1-Å resolution. In this complex, one side of DNA is accessible for cleavage by RuvC resolvase at the junction center. The refined junction DNA structure revealed an open concave architecture with a four-fold symmetry. Each arm, with B-form DNA, in the Holliday junction is predominantly recognized in the minor groove through hydrogen bonds with two repeated helix-hairpin-helix motifs of each RuvA subunit. The local conformation near the crossover point, where two base pairs are disrupted, suggests a possible scheme for successive base pair rearrangements, which may account for smooth Holliday junction movement without segmental unwinding.

NMR structure of the bovine prion protein

The NMR structures of the recombinant 217-residue polypeptide chain of the mature bovine prion protein, bPrP(23–230), and a C-terminal fragment, bPrP(121–230), include a globular domain extending from residue 125 to residue 227, a short flexible chain end of residues 228–230, and an N-terminal flexibly disordered “tail” comprising 108 residues for the intact protein and 4 residues for bPrP(121–230), respectively. The globular domain contains three α-helices comprising the residues 144–154, 173–194, and 200–226, and a short antiparallel β-sheet comprising the residues 128–131 and 161–164. The best-defined parts of the globular domain are the central portions of the helices 2 and 3, which are linked by the only disulfide bond in bPrP. Significantly increased disorder and mobility is observed for helix 1, the loop 166–172 leading from the β-strand 2 to helix 2, the end of helix 2 and the following loop, and the last turn of helix 3. Although there are characteristic local differences relative to the conformations of the murine and Syrian hamster prion proteins, the bPrP structure is essentially identical to that of the human prion protein. On the other hand, there are differences between bovine and human PrP in the surface distribution of electrostatic charges, which then appears to be the principal structural feature of the “healthy” PrP form that might affect the stringency of the species barrier for transmission of prion diseases between humans and cattle.

Photochemically induced nuclear spin polarization in reaction centers of photosystem II observed by 13C-solid-state NMR reveals a strongly asymmetric electronic structure of the P680.+ primary donor chlorophyll

We report 13C magic angle spinning NMR observation of photochemically induced dynamic nuclear spin polarization (photo- CIDNP) in the reaction center (RC) of photosystem II (PS2). The light-enhanced NMR signals of the natural abundance 13C provide information on the electronic structure of the primary electron donor P680 (chlorophyll a molecules absorbing around 680 nm) and on the pz spin density pattern in its oxidized form, P680⨥. Most centerband signals can be attributed to a single chlorophyll a (Chl a) cofactor that has little interaction with other pigments. The chemical shift anisotropy of the most intense signals is characteristic for aromatic carbon atoms. The data reveal a pronounced asymmetry of the electronic spin density distribution within the P680⨥. PS2 shows only a single broad and intense emissive signal, which is assigned to both the C-10 and C-15 methine carbon atoms. The spin density appears shifted toward ring III. This shift is remarkable, because, for monomeric Chl a radical cations in solution, the region of highest spin density is around ring II. It leads to a first hypothesis as to how the planet can provide itself with the chemical potential to split water and generate an oxygen atmosphere using the Chl a macroaromatic cycle. A local electrostatic field close to ring III can polarize the electronic charge and associated spin density and increase the redox potential of P680 by stabilizing the highest occupied molecular orbital, without a major change of color. This field could be produced, e.g., by protonation of the keto group of ring V. Finally, the radical cation electronic structure in PS2 is different from that in the bacterial RC, which shows at least four emissive centerbands, indicating a symmetric spin density distribution over the entire bacteriochlorophyll macrocycle.

A database and tools for 3-D protein structure comparison and alignment using the Combinatorial Extension (CE) algorithm

The database reported here is derived using the Combinatorial Extension (CE) algorithm which compares pairs of protein polypeptide chains and provides a list of structurally similar proteins along with their structure alignments. Using CE, structure–structure alignments can provide insights into biological function. When a protein of known function is shown to be structurally similar to a protein of unknown function, a relationship might be inferred; a relationship not necessarily detectable from sequence comparison alone. Establishing structure–structure relationships in this way is of great importance as we enter an era of structural genomics where there is a likelihood of an increasing number of structures with unknown functions being determined. Thus the CE database is an example of a useful tool in the annotation of protein structures of unknown function. Comparisons can be performed on the complete PDB or on a structurally representative subset of proteins. The source protein(s) can be from the PDB (updated monthly) or uploaded by the user. CE provides sequence alignments resulting from structural alignments and Cartesian coordinates for the aligned structures, which may be analyzed using the supplied Compare3D Java applet, or downloaded for further local analysis. Searches can be run from the CE web site, http://cl.sdsc.edu/ce.html, or the database and software downloaded from the site for local use.

RNA–protein hybrid ribozymes that efficiently cleave any mRNA independently of the structure of the target RNA

Ribozyme activity in vivo depends on achieving high-level expression, intracellular stability, target colocalization, and cleavage site access. At present, target site selection is problematic because of unforeseeable secondary and tertiary RNA structures that prevent cleavage. To overcome this design obstacle, we wished to engineer a ribozyme that could access any chosen site. To create this ribozyme, the constitutive transport element (CTE), an RNA motif that has the ability to interact with intracellular RNA helicases, was attached to our ribozymes so that the helicase-bound, hybrid ribozymes would be produced in cells. This modification significantly enhanced ribozyme activity in vivo, permitting cleavage of sites previously found to be inaccessible. To confer cleavage enhancement, the CTE must retain helicase-binding activity. Binding experiments demonstrated the likely involvement of RNA helicase(s). We found that attachment of the RNA motif to our tRNA ribozymes leads to cleavage in vivo at the chosen target site regardless of the local RNA secondary or tertiary structure.

A model of self-thinning through local competition.

Explanations of self-thinning in plant populations have focused on plant shape and packing. A dynamic model based on the structure of local interactions successfully reproduces the pattern and can be approximated to identify key parameters and relationships. The approach generates testable new explanations for differences between species and populations, unifies self-thinning with other patterns in plant population dynamics, and indicates why organisms other than plants can follow the law.

Crystal structure at 2.6-A resolution of human macrophage migration inhibitory factor.

Macrophage migration inhibitory factor (MIF) was the first cytokine to be described, but for 30 years its role in the immune response remained enigmatic. In recent studies, MIF has been found to be a novel pituitary hormone and the first protein identified to be released from immune cells on glucocorticoid stimulation. Once secreted, MIF counterregulates the immunosuppressive effects of steroids and thus acts as a critical component of the immune system to control both local and systemic immune responses. We report herein the x-ray crystal structure of human MIF to 2.6 angstrom resolution. The protein is a trimer of identical subunits. Each monomer contains two antiparallel alpha-helices that pack against a four-stranded beta-sheet. The monomer has an additional two beta-strands that interact with the beta-sheets of adjacent subunits to form the interface between monomers. The three beta-sheets are arranged to form a barrel containing a solvent-accessible channel that runs through the center of the protein along a molecular 3-fold axis. Electrostatic potential maps reveal that the channel has a positive potential, suggesting that it binds negatively charged molecules. The elucidated structure for MIF is unique among cytokines or hormonal mediators, and suggests that this counterregulator of glucocorticoid action participates in novel ligand-receptor interactions.

Local densities orthogonal to beta-sheet amide planes: patterns of packing in globular proteins.

We have investigated the efficiency of packing by calculating intramolecular packing density above and below peptide planes of internal beta-pleated sheet residues in five globular proteins. The orientation of interest was chosen to allow study of regions that are approximately perpendicular to the faces of beta-pleated sheets. In these locations, nonbonded van der Waals packing interactions predominate over hydrogen bonding and solvent interactions. We observed considerable variability in packing densities within these regions, confirming that the interior packing of a protein does not result in uniform occupation of the available space. Patterns of fluctuation in packing density suggest that the regular backbone-to-backbone network of hydrogen bonds is not likely to be interrupted to maximize van der Waals interactions. However, high-density packing tends to occur toward the ends of beta-structure strands where hydrogen bonds are more likely to involve nonpolar side-chain groups or solvent molecules. These features result in internal protein folding with a central low-density core surrounded by a higher-density subsurface shell, consistent with our previous calculations regarding overall protein packing density.

From local to global in quasiconformal structures.

We exhibit a large class of metric spaces whose infinitesimal quasiconformal structure is strong enough to capture the global quasiconformal structure. A sufficient condition for this to happen is described in terms of a Poincaré-type inequality.

Tertiary structure of an amyloid immunoglobulin light chain protein: a proposed model for amyloid fibril formation.

An immunoglobulin light chain protein was isolated from the urine of an individual (BRE) with systemic amyloidosis. Complete amino acid sequence of the variable region of the light chain (VL) protein established it as a kappa I, which when compared with other kappa I amyloid associated proteins had unique residues, including Ile-34, Leu-40, and Tyr-71. To study the tertiary structure, BRE VL was expressed in Escherichia coli by using a PCR product amplified from the patient BRE's bone marrow DNA. The PCR product was ligated into pCZ11, a thermal-inducible replication vector. Recombinant BRE VL was isolated, purified to homogeneity, and crystallized by using ammonium sulfate as the precipitant. Two crystal forms were obtained. In crystal form I the BRE VL kappa domain crystallizes as a dimer with unit cell constants isomorphous to previously published kappa protein structures. Comparison with a nonamyloid VL kappa domain from patient REI, identified significant differences in position of residues in the hypervariable segments plus variations in framework region (FR) segments 40-46 (FR2) and 66-67 (FR3). In addition, positional differences can be seen along the two types of local diads, corresponding to the monomer-monomer and dimer-dimer interfaces. From the packing diagram, a model for the amyloid light chain (AL) fibril is proposed based on a pseudohexagonal spiral structure with a rise of approximately the width of two dimers per 360 degree turn. This spiral structure could be consistent with the dimensions of amyloid fibrils as determined by electron microscopy.

Periodicity of polar and nonpolar amino acids is the major determinant of secondary structure in self-assembling oligomeric peptides.

The tendency of a polypeptide chain to form alpha-helical or beta-strand secondary structure depends upon local and nonlocal effects. Local effects reflect the intrinsic propensities of the amino acid residues for particular secondary structures, while nonlocal effects reflect the positioning of the individual residues in the context of the entire amino acid sequence. In particular, the periodicity of polar and nonpolar residues specifies whether a given sequence is consistent with amphiphilic alpha-helices or beta-strands. The importance of intrinsic propensities was compared to that of polar/nonpolar periodicity by a direct competition. Synthetic peptides were designed using residues with intrinsic propensities that favored one or the other type of secondary structure. The polar/nonpolar periodicities of the peptides were designed either to be consistent with the secondary structure favored by the intrinsic propensities of the component residues or in other cases to oppose these intrinsic propensities. Characterization of the synthetic peptides demonstrated that in all cases the observed secondary structure correlates with the periodicity of the peptide sequence--even when this secondary structure differs from that predicted from the intrinsic propensities of the component amino acids. The observed secondary structures are concentration dependent, indicating that oligomerization of the amphiphilic peptides is responsible for the observed secondary structures. Thus, for self-assembling oligomeric peptides, the polar/nonpolar periodicity can overwhelm the intrinsic propensities of the amino acid residues and serves as the major determinant of peptide secondary structure.