Cross-validated maximum likelihood enhances crystallographic simulated annealing refinement

Recently, the target function for crystallographic refinement has been improved through a maximum likelihood analysis, which makes proper allowance for the effects of data quality, model errors, and incompleteness. The maximum likelihood target reduces the significance of false local minima during the refinement process, but it does not completely eliminate them, necessitating the use of stochastic optimization methods such as simulated annealing for poor initial models. It is shown that the combination of maximum likelihood with cross-validation, which reduces overfitting, and simulated annealing by torsion angle molecular dynamics, which simplifies the conformational search problem, results in a major improvement of the radius of convergence of refinement and the accuracy of the refined structure. Torsion angle molecular dynamics and the maximum likelihood target function interact synergistically, the combination of both methods being significantly more powerful than each method individually. This is demonstrated in realistic test cases at two typical minimum Bragg spacings (dmin = 2.0 and 2.8 Å, respectively), illustrating the broad applicability of the combined method. In an application to the refinement of a new crystal structure, the combined method automatically corrected a mistraced loop in a poor initial model, moving the backbone by 4 Å.

Structure and function in rhodopsin: Topology of the C-terminal polypeptide chain in relation to the cytoplasmic loops*

Cysteine mutagenesis and site-directed spin labeling in the C-terminal region of rhodopsin have been used to probe the local structure and proximity of that region to the cytoplasmic loops. Each of the native amino acids in the sequence T335–T340 was replaced with Cys, one at a time. The sulfhydryl groups of all mutants reacted rapidly with the sulfhydryl reagent 4,4′-dithiodipyridine, which indicated a high degree of solvent accessibility. Furthermore, to probe the proximity relationships, a series of double Cys mutants was constructed. One Cys in all sets was at position 338 and the other was at a position in the sequence S240–V250 in the EF interhelical loop, at position 65 in the AB interhelical loop, or at position 140 in the CD interhelical loop. In the dark state, no significant disulfide formation was observed between C338 and C65 or C140 under the conditions used, whereas a relatively rapid disulfide formation was observed between C338 and C242 or C245. Spin labels in the double Cys mutants showed the strongest magnetic interactions between the nitroxides attached to C338 and C245 or C246. Light activation of the double mutant T242C/S338C resulted in slower disulfide formation, whereas interactions between nitroxides at C338 and C245 or C246 decreased. These results suggest the proximity of the C-terminal residue C338 to residues located on the outer face of a cytoplasmic helical extension of the F helix with an apparent increase of distance upon photoactivation.

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.

Interaction of batrachotoxin with the local anesthetic receptor site in transmembrane segment IVS6 of the voltage-gated sodium channel

The voltage-gated sodium channel is the site of action of more than six classes of neurotoxins and drugs that alter its function by interaction with distinct, allosterically coupled receptor sites. Batrachotoxin (BTX) is a steroidal alkaloid that binds to neurotoxin receptor site 2 and causes persistent activation. BTX binding is inhibited allosterically by local anesthetics. We have investigated the interaction of BTX with amino acid residues I1760, F1764, and Y1771, which form part of local anesthetic receptor site in transmembrane segment IVS6 of type IIA sodium channels. Alanine substitution for F1764 (mutant F1764A) reduces tritiated BTX-A-20-α-benzoate binding affinity, causing a 60-fold increase in Kd. Alanine substitution for I1760, which is adjacent to F1764 in the predicted IVS6 transmembrane alpha helix, causes only a 4-fold increase in Kd. In contrast, mutant Y1771A shows no change in BTX binding affinity. For wild-type and mutant Y1771A, BTX shifted the voltage for half-maximal activation ≈40 mV in the hyperpolarizing direction and increased the percentage of noninactivating sodium current to ≈60%. In contrast, these BTX effects were eliminated completely for the F1764A mutant and were reduced substantially for mutant I1760A. Our data suggest that the BTX receptor site shares overlapping but nonidentical molecular determinants with the local anesthetic receptor site in transmembrane segment IVS6 as well as having unique molecular determinants in transmembrane segment IS6, as demonstrated in previous work. Evidently, BTX conforms to a domain–interface allosteric model of ligand binding and action, as previously proposed for calcium agonist and antagonist drugs acting on l-type calcium channels.

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.

The area code hypothesis revisited: Olfactory receptors and other related transmembrane receptors may function as the last digits in a cell surface code for assembling embryos

Recent evidence emerging from several laboratories, integrated with new data obtained by searching the genome databases, suggests that the area code hypothesis provides a good heuristic model for explaining the remarkable specificity of cell migration and tissue assembly that occurs throughout embryogenesis. The area code hypothesis proposes that cells assemble organisms, including their brains and nervous systems, with the aid of a molecular-addressing code that functions much like the country, area, regional, and local portions of the telephone dialing system. The complexity of the information required to code cells for the construction of entire organisms is so enormous that we assume that the code must make combinatorial use of members of large multigene families. Such a system would reuse the same receptors as molecular digits in various regions of the embryo, thus greatly reducing the total number of genes required. We present the hypothesis that members of the very large families of olfactory receptors and vomeronasal receptors fulfill the criteria proposed for area code molecules and could serve as the last digits in such a code. We discuss our evidence indicating that receptors of these families are expressed in many parts of developing embryos and suggest that they play a key functional role in cell recognition and targeting not only in the olfactory system but also throughout the brain and numerous other organs as they are assembled.

Local translation of classes of mRNAs that are targeted to neuronal dendrites

The functioning of the neuronal dendrite results from a variety of biological processes including mRNA transport to and protein translation in the dendrite. The complexity of the mRNA population in dendrites suggests that specific biological processes are modulated through the regulation of dendritic biology. There are various classes of mRNAs in dendrites whose translation modulates the ability of the dendrite to receive and integrate presynaptic information. Among these mRNAs are those encoding selective transcription factors that function in the neuronal soma and ionotropic glutamate receptors that function on the neuronal membrane. Conclusive evidence that these mRNAs can be translated is reviewed, and identification of the endogenous sites of translation in living dendrites is presented. These data, as well as those described in the other articles resulting from this colloquium, highlight the complexity of dendritic molecular biology and the exquisitely selective and sensitive modulatory role played by the dendrite in facilitating intracellular and intercellular communication.

Function of the oxidative burst in hypersensitive disease resistance.

Microbial elicitors or attempted infection with an avirulent pathogen strain causes the rapid production of reactive oxygen intermediates. Recent findings indicate that H2O2 from this oxidative burst plays a central role in the orchestration of the hypersensitive response: (i) as the substrate driving the cross-linking of cell wall structural proteins to slow microbial ingress prior to the deployment of transcription-dependent defenses and to trap pathogens in cells destined to undergo hypersensitive cell death, (ii) as a local threshold trigger of this programmed death in challenged cells, and (iii) as a diffusible signal for the induction in adjacent cells of genes encoding cellular protectants such as glutathione S-transferase and glutathione peroxidase. These findings provide the basis for an integrated model for the orchestration of the localized hypersensitive resistance response to attack by an avirulent pathogen.

Bipartite function of a small RNA hairpin in transcription antitermination in bacteriophage lambda.

Transcription of downstream genes in the early operons of phage lambda requires a promoter-proximal element known as nut. This site acts in cis in the form of RNA to assemble a transcription antitermination complex which is composed of lambda N protein and at least four host factors. The nut-site RNA contains a small stem-loop structure called boxB. Here, we show that boxB RNA binds to N protein with high affinity and specificity. While N binding is confined to the 5' subdomain of the stem-loop, specific N recognition relies on both an intact stem-loop structure and two critical nucleotides in the pentamer loop. Substitutions of these nucleotides affect both N binding and antitermination. Remarkably, substitutions of other loop nucleotides also diminish antitermination in vivo, yet they have no detectable effect on N binding in vitro. These 3' loop mutants fail to support antitermination in a minimal system with RNA polymerase (RNAP), N, and the host factor NusA. Furthermore, the ability of NusA to stimulate the formation of the RNAP-boxB-N complex is diminished with these mutants. Hence, we suggest that boxB RNA performs two critical functions in antitermination. First, boxB binds to N and secures it near RNAP to enhance their interaction, presumably by increasing the local concentration of N. Second, boxB cooperates with NusA, most likely to bring N and RNAP in close contact and transform RNAP to the termination-resistant state.