Systematic derivation of partition functions for ligand binding to two-dimensional lattices

The Ising problem consists in finding the analytical solution of the partition function of a lattice once the interaction geometry among its elements is specified. No general analytical solution is available for this problem, except for the one-dimensional case. Using site-specific thermodynamics, it is shown that the partition function for ligand binding to a two-dimensional lattice can be obtained from those of one-dimensional lattices with known solution. The complexity of the lattice is reduced recursively by application of a contact transformation that involves a relatively small number of steps. The transformation implemented in a computer code solves the partition function of the lattice by operating on the connectivity matrix of the graph associated with it. This provides a powerful new approach to the Ising problem, and enables a systematic analysis of two-dimensional lattices that model many biologically relevant phenomena. Application of this approach to finite two-dimensional lattices with positive cooperativity indicates that the binding capacity per site diverges as Na (N = number of sites in the lattice) and experiences a phase-transition-like discontinuity in the thermodynamic limit N → ∞. The zeroes of the partition function tend to distribute on a slightly distorted unit circle in complex plane and approach the positive real axis already for a 5×5 square lattice. When the lattice has negative cooperativity, its properties mimic those of a system composed of two classes of independent sites with the apparent population of low-affinity binding sites increasing with the size of the lattice, thereby accounting for a phenomenon encountered in many ligand-receptor interactions.

Membrane-bound state of the colicin E1 channel domain as an extended two-dimensional helical array

Atomic level structures have been determined for the soluble forms of several colicins and toxins, but the structural changes that occur after membrane binding have not been well characterized. Changes occurring in the transition from the soluble to membrane-bound state of the C-terminal 190-residue channel polypeptide of colicin E1 (P190) bound to anionic membranes are described. In the membrane-bound state, the α-helical content increases from 60–64% to 80–90%, with a concomitant increase in the average length of the helical segments from 12 to 16 or 17 residues, close to the length required to span the membrane bilayer in the open channel state. The average distance between helical segments is increased and interhelix interactions are weakened, as shown by a major loss of tertiary structure interactions, decreased efficiency of fluorescence resonance energy transfer from an energy donor on helix V of P190 to an acceptor on helix IX, and decreased resonance energy transfer at higher temperatures, not observed in soluble P190, implying freedom of motion of helical segments. Weaker interactions are also shown by a calorimetric thermal transition of low cooperativity, and the extended nature of the helical array is shown by a 3- to 4-fold increase in the average area subtended per molecule to 4,200 Å2 on the membrane surface. The latter, with analysis of the heat capacity changes, implies the absence of a developed hydrophobic core in the membrane-bound P190. The membrane interfacial layer thus serves to promote formation of a highly helical extended two-dimensional flexible net. The properties of the membrane-bound state of the colicin channel domain (i.e., hydrophobic anchor, lengthened and loosely coupled α-helices, and close association with the membrane interfacial layer) are plausible structural features for the state that is a prerequisite for voltage gating, formation of transmembrane helices, and channel opening.

Two-dimensional freezing in the liquid–vapor interface of a dilute Pb:Ga alloy

We report the results of x-ray reflectivity and grazing incidence x-ray diffraction studies of the liquid–vapor interface of a dilute alloy of Pb in Ga over the temperature range of 23–76°C. Our data show that the liquid–vapor interface of this alloy is stratified for several atomic diameters into the bulk liquid and that a monolayer of Pb forms the outermost stratum of the interface. Over the temperature range of 23–56°C, the monolayer of Pb is in an ordered hexagonal phase. At about 58°C, this monolayer undergoes a first-order transition to a hexatic phase, which remains stable to 76°C. An analogy between the observed transition and the first-order melting transition in a one-component classical plasma is suggested.

Two-dimensional IR correlation spectroscopy: Sequential events in the unfolding process of the λ Cro-V55C repressor protein

A question often posed in protein folding/unfolding studies is whether the process is fully cooperative or whether it contains sequential elements. To address this question, one needs tools capable of resolving different events. It seems that, at least in certain cases, two-dimensional (2D) IR correlation spectroscopy can provide answers to this question. To illustrate this point, we have turned to the Cro-V55C dimer of the λ Cro repressor, a protein known to undergo thermal unfolding in two discrete steps through a stable equilibrium intermediate. The secondary structure of this intermediate is compatible with that of a partially unfolded protein and involves a reorganization of the N terminus, whereas the antiparallel β-ribbon formed by the C-terminal part of each subunit remains largely intact. To establish whether the unfolding process involves sequential events, we have performed a 2D correlation analysis of IR spectra recorded over the temperature range of 20–95°C. The 2D IR correlation analysis indeed provides evidence for a sequential formation of the stable intermediate, which is created in three (closely related) steps. A first step entails the unfolding of the short N-terminal β-strand, followed by the unfolding of the α-helices in a second step, and the third step comprises the reorganization of the remaining β-sheet and of some unordered segments in the protein. The complete unfolding of the stable intermediate at higher temperatures also undergoes sequential events that ultimately end with the breaking of the H bonds between the two β-strands at the dimer interface.

Linking genome and proteome by mass spectrometry: Large-scale identification of yeast proteins from two dimensional gels

The function of many of the uncharacterized open reading frames discovered by genomic sequencing can be determined at the level of expressed gene products, the proteome. However, identifying the cognate gene from minute amounts of protein has been one of the major problems in molecular biology. Using yeast as an example, we demonstrate here that mass spectrometric protein identification is a general solution to this problem given a completely sequenced genome. As a first screen, our strategy uses automated laser desorption ionization mass spectrometry of the peptide mixtures produced by in-gel tryptic digestion of a protein. Up to 90% of proteins are identified by searching sequence data bases by lists of peptide masses obtained with high accuracy. The remaining proteins are identified by partially sequencing several peptides of the unseparated mixture by nanoelectrospray tandem mass spectrometry followed by data base searching with multiple peptide sequence tags. In blind trials, the method led to unambiguous identification in all cases. In the largest individual protein identification project to date, a total of 150 gel spots—many of them at subpicomole amounts—were successfully analyzed, greatly enlarging a yeast two-dimensional gel data base. More than 32 proteins were novel and matched to previously uncharacterized open reading frames in the yeast genome. This study establishes that mass spectrometry provides the required throughput, the certainty of identification, and the general applicability to serve as the method of choice to connect genome and proteome.

The two-dimensional IR nonlinear spectroscopy of a cyclic penta-peptide in relation to its three-dimensional structure

A form of two-dimensional (2D) vibrational spectroscopy, which uses two ultrafast IR laser pulses, is used to examine the structure of a cyclic penta-peptide in solution. Spectrally resolved cross peaks occur in the off-diagonal region of the 2D IR spectrum of the amide I region, analogous to those in 2D NMR spectroscopy. These cross peaks measure the coupling between the different amide groups in the structure. Their intensities and polarizations relate directly to the three-dimensional structure of the peptide. With the help of a model coupling Hamiltonian, supplemented by density functional calculations, the spectra of this penta-peptide can be regenerated from the known solution phase structure. This 2D-IR measurement, with an intrinsic time resolution of less than 1 ps, could be used in all time regimes of interest in biology.

Two-dimensional infrared spectroscopy of peptides by phase-controlled femtosecond vibrational photon echoes

Two-dimensional infrared spectra of peptides are introduced that are the direct analogues of two- and three-pulse multiple quantum NMR. Phase matching and heterodyning are used to isolate the phase and amplitudes of the electric fields of vibrational photon echoes as a function of multiple pulse delays. Structural information is made available on the time scale of a few picoseconds. Line narrowed spectra of acyl-proline-NH2 and cross peaks implying the coupling between its amide-I modes are obtained, as are the phases of the various contributions to the signals. Solvent-sensitive structural differences are seen for the dipeptide. The methods show great promise to measure structure changes in biology on a wide range of time scales.

Two-dimensional confocal images of organization, density, and gating of focal Ca2+ release sites in rat cardiac myocytes

In cardiac myocytes Ca2+ cross-signaling between Ca2+ channels and ryanodine receptors takes place by exchange of Ca2+ signals in microdomains surrounding dyadic junctions, allowing first the activation and then the inactivation of the two Ca2+-transporting proteins. To explore the details of Ca2+ signaling between the two sets of receptors we measured the two-dimensional cellular distribution of Ca2+ at 240 Hz by using a novel confocal imaging technique. Ca2+ channel-triggered Ca2+ transients could be resolved into dynamic “Ca2+ stripes” composed of hundreds of discrete focal Ca2+ releases, appearing as bright fluorescence spots (radius ≅ 0.5 μm) at reproducible sites, which often coincided with t-tubules as visualized with fluorescent staining of the cell membrane. Focal Ca2+ releases triggered stochastically by Ca2+ current (ICa) changed little in duration (≅7 ms) and size (≅100,000 Ca ions) between −40 and +60 mV, but their frequency of activation and first latency mirrored the kinetics and voltage dependence of ICa. The resolution of 0.95 ± 0.13 reproducible focal Ca2+ release sites per μm3 in highly Ca2+-buffered cells, where diffusion of Ca2+ is limited to 50 nm, suggests the presence of about one independent, functional Ca2+ release site per half sarcomere. The density and distribution of Ca2+ release sites suggest they correspond to dyadic junctions. The abrupt onset and termination of focal Ca2+ releases indicate that the cluster of ryanodine receptors in individual dyadic junctions may operate in a coordinated fashion.

Behavior of several two-dimensional fluid equations in singular scenarios

We give conditions that rule out formation of sharp fronts for certain two-dimensional incompressible flows. We show that a necessary condition of having a sharp front is that the flow has to have uncontrolled velocity growth. In the case of the quasi-geostrophic equation and two-dimensional Euler equation, we obtain estimates on the formation of semi-uniform fronts.

Two-dimensional protein crystallization via metal-ion coordination by naturally occurring surface histidines.

A powerful and potentially general approach to the targeting and crystallization of proteins on lipid interfaces through coordination of surface histidine residues to lipid-chelated divalent metal ions is presented. This approach, which should be applicable to the crystallization of a wide range of naturally occurring or engineered proteins, is illustrated here by the crystallization of streptavidin on a monolayer of an iminodiacetate-Cu(II) lipid spread at the air-water interface. This method allows control of the protein orientation at interfaces, which is significant for the facile production of highly ordered protein arrays and for electron density mapping in structural analysis of two-dimensional crystals. Binding of native streptavidin to the iminodiacetate-Cu lipids occurs via His-87, located on the protein surface near the biotin binding pocket. The two-dimensional streptavidin crystals show a previously undescribed microscopic shape that differs from that of crystals formed beneath biotinylated lipids.

Two Drosophila nervous system antigens, Nervana 1 and 2, are homologous to the beta subunit of Na+,K(+)-ATPase.

A nervous system-specific glycoprotein antigen from adult Drosophila heads, designated Nervana (Nrv), has been purified on the basis of reactivity of its carbohydrate epitope(s) with anti-horseradish peroxidase (HRP) antibodies that are specific markers for Drosophila neurons. Anti-Nrv monoclonal antibodies (mAbs), specific for the protein moiety of Nrv, were used to screen a Drosophila embryo cDNA expression library. Three cDNA clones (designated Nrv1, Nrv2.1, and Nrv2.2) were isolated that code for proteins recognized by anti-Nrv mAbs on Western blots. DNA sequencing and Southern blot analyses established that the cDNA clones are derived from two different genes. In situ hybridization to Drosophila polytene chromosomes showed that the cDNA clones map to the third chromosome near 92C-D. Nrv1 and Nrv2.1/2.2 have open reading frames of 309 and 322/323 amino acids, respectively, and they are 43.4% identical at the amino acid level. The proteins deduced from these clones exhibit significant homology in both primary sequence and predicted topology to the beta subunit of Na+,K(+)-ATPase. Immunoaffinity-purified Nrv is associated with a protein (M(r) 100,000) recognized on Western blots by anti-ATPase alpha-subunit mAb. Our results suggest that the Drosophila nervous system-specific antigens Nrv1 and -2 are neuronal forms of the beta subunit of Na+,K(+)-ATPase.

Rapid mass spectrometric peptide sequencing and mass matching for characterization of human melanoma proteins isolated by two-dimensional PAGE.

We report a general mass spectrometric approach for the rapid identification and characterization of proteins isolated by preparative two-dimensional polyacrylamide gel electrophoresis. This method possesses the inherent power to detect and structurally characterize covalent modifications. Absolute sensitivities of matrix-assisted laser desorption ionization and high-energy collision-induced dissociation tandem mass spectrometry are exploited to determine the mass and sequence of subpicomole sample quantities of tryptic peptides. These data permit mass matching and sequence homology searching of computerized peptide mass and protein sequence data bases for known proteins and design of oligonucleotide probes for cloning unknown proteins. We have identified 11 proteins in lysates of human A375 melanoma cells, including: alpha-enolase, cytokeratin, stathmin, protein disulfide isomerase, tropomyosin, Cu/Zn superoxide dismutase, nucleoside diphosphate kinase A, galaptin, and triosephosphate isomerase. We have characterized several posttranslational modifications and chemical modifications that may result from electrophoresis or subsequent sample processing steps. Detection of comigrating and covalently modified proteins illustrates the necessity of peptide sequencing and the advantages of tandem mass spectrometry to reliably and unambiguously establish the identity of each protein. This technology paves the way for studies of cell-type dependent gene expression and studies of large suites of cellular proteins with unprecedented speed and rigor to provide information complementary to the ongoing Human Genome Project.

The ethylene hormone response in Arabidopsis: a eukaryotic two-component signaling system.

The simple gas ethylene affects numerous physiological processes in the growth and development of higher plants. With the use of molecular genetic approaches, we are beginning to learn how plants perceive ethylene and how this signal is transduced. Components of ethylene signal transduction are defined by ethylene response mutants in Arabidopsis thaliana. The genes corresponding to two of these mutants, etr1 and etr1, have been cloned. The ETR1 gene encodes a homolog of two-component regulators that are known almost exclusively in prokaryotes. The two-component regulators in prokaryotes are involved in the perception and transduction of a wide range of environmental signals leading to adaptive responses. The CTR1 gene encodes a homolog of the Raf family of serine/threonine protein kinases. Raf is part of a mitogen-activated protein kinase cascade known to regulate cell growth and development in mammals, worms, and flies. The ethylene response pathway may, therefore, exemplify a conserved protein kinase cascade regulated by a two-component system. The dominance of all known mutant alleles of ETR1 may be due to either constitutive activation of the ETR1 protein or dominant interference of wild-type activity. The discovery of Arabidopsis genes encoding proteins related to ETR1 suggests that the failure to recover recessive etr1 mutant alleles may be due to the presence of redundant genes.