Folding protein models with a simple hydrophobic energy function: The fundamental importance of monomer inside/outside segregation

The present study explores a “hydrophobic” energy function for folding simulations of the protein lattice model. The contribution of each monomer to conformational energy is the product of its “hydrophobicity” and the number of contacts it makes, i.e., E(h⃗, c⃗) = −Σi=1N cihi = −(h⃗.c⃗) is the negative scalar product between two vectors in N-dimensional cartesian space: h⃗ = (h1, … , hN), which represents monomer hydrophobicities and is sequence-dependent; and c⃗ = (c1, … , cN), which represents the number of contacts made by each monomer and is conformation-dependent. A simple theoretical analysis shows that restrictions are imposed concomitantly on both sequences and native structures if the stability criterion for protein-like behavior is to be satisfied. Given a conformation with vector c⃗, the best sequence is a vector h⃗ on the direction upon which the projection of c⃗ − c̄⃗ is maximal, where c̄⃗ is the diagonal vector with components equal to c̄, the average number of contacts per monomer in the unfolded state. Best native conformations are suggested to be not maximally compact, as assumed in many studies, but the ones with largest variance of contacts among its monomers, i.e., with monomers tending to occupy completely buried or completely exposed positions. This inside/outside segregation is reflected on an apolar/polar distribution on the corresponding sequence. Monte Carlo simulations in two dimensions corroborate this general scheme. Sequences targeted to conformations with large contact variances folded cooperatively with thermodynamics of a two-state transition. Sequences targeted to maximally compact conformations, which have lower contact variance, were either found to have degenerate ground state or to fold with much lower cooperativity.

Recent improvements in prediction of protein structure by global optimization of a potential energy function

Recent improvements of a hierarchical ab initio or de novo approach for predicting both α and β structures of proteins are described. The united-residue energy function used in this procedure includes multibody interactions from a cumulant expansion of the free energy of polypeptide chains, with their relative weights determined by Z-score optimization. The critical initial stage of the hierarchical procedure involves a search of conformational space by the conformational space annealing (CSA) method, followed by optimization of an all-atom model. The procedure was assessed in a recent blind test of protein structure prediction (CASP4). The resulting lowest-energy structures of the target proteins (ranging in size from 70 to 244 residues) agreed with the experimental structures in many respects. The entire experimental structure of a cyclic α-helical protein of 70 residues was predicted to within 4.3 Å α-carbon (Cα) rms deviation (rmsd) whereas, for other α-helical proteins, fragments of roughly 60 residues were predicted to within 6.0 Å Cα rmsd. Whereas β structures can now be predicted with the new procedure, the success rate for α/β- and β-proteins is lower than that for α-proteins at present. For the β portions of α/β structures, the Cα rmsd's are less than 6.0 Å for contiguous fragments of 30–40 residues; for one target, three fragments (of length 10, 23, and 28 residues, respectively) formed a compact part of the tertiary structure with a Cα rmsd less than 6.0 Å. Overall, these results constitute an important step toward the ab initio prediction of protein structure solely from the amino acid sequence.

The role of structure, energy landscape, dynamics, and allostery in the enzymatic function of myoglobin

The grail of protein science is the connection between structure and function. For myoglobin (Mb) this goal is close. Described as only a passive dioxygen storage protein in texts, we argue here that Mb is actually an allosteric enzyme that can catalyze reactions among small molecules. Studies of the structural, spectroscopic, and kinetic properties of Mb lead to a model that relates structure, energy landscape, dynamics, and function. Mb functions as a miniature chemical reactor, concentrating and orienting diatomic molecules such as NO, CO, O2, and H2O2 in highly conserved internal cavities. Reactions can be controlled because Mb exists in distinct taxonomic substates with different catalytic properties and connectivities of internal cavities.

Single-pair fluorescence resonance energy transfer on freely diffusing molecules: Observation of Förster distance dependence and subpopulations

Photon bursts from single diffusing donor-acceptor labeled macromolecules were used to measure intramolecular distances and identify subpopulations of freely diffusing macromolecules in a heterogeneous ensemble. By using DNA as a rigid spacer, a series of constructs with varying intramolecular donor-acceptor spacings were used to measure the mean and distribution width of fluorescence resonance energy transfer (FRET) efficiencies as a function of distance. The mean single-pair FRET efficiencies qualitatively follow the distance dependence predicted by Förster theory. Possible contributions to the widths of the FRET efficiency distributions are discussed, and potential applications in the study of biopolymer conformational dynamics are suggested. The ability to measure intramolecular (and intermolecular) distances for single molecules implies the ability to distinguish and monitor subpopulations of molecules in a mixture with different distances or conformational states. This is demonstrated by monitoring substrate and product subpopulations before and after a restriction endonuclease cleavage reaction. Distance measurements at single-molecule resolution also should facilitate the study of complex reactions such as biopolymer folding. To this end, the denaturation of a DNA hairpin was examined by using single-pair FRET.

Factors controlling the competition among rotational and vibrational energy transfer channels in glyoxal

The state-to-state transfer of rotational and vibrational energy has been studied for S1 glyoxal (CHOCHO) in collisions with D2, N2, CO and C2H4 using crossed molecular beams. A laser is used to pump glyoxal seeded in He to its S1 zero point level with zero angular momentum about its top axis (K′ = 0). The inelastic scattering to each of at least 26 S1 glyoxal rotational and rovibrational levels is monitored by dispersed S1–S0 fluorescence. Various collision partners are chosen to investigate the relative influences of reduced mass and the collision pair interaction potential on the competition among the energy transfer channels. When the data are combined with that obtained previously from other collision partners whose masses range from 2 to 84 amu, it is seen that the channel competition is controlled primarily by the kinematics of the collisional interaction. Variations in the intermolecular potential play strictly a secondary role.

Single-molecule protein folding: Diffusion fluorescence resonance energy transfer studies of the denaturation of chymotrypsin inhibitor 2

We report single-molecule folding studies of a small, single-domain protein, chymotrypsin inhibitor 2 (CI2). CI2 is an excellent model system for protein folding studies and has been extensively studied, both experimentally (at the ensemble level) and theoretically. Conformationally assisted ligation methodology was used to synthesize the proteins and site-specifically label them with donor and acceptor dyes. Folded and denatured subpopulations were observed by fluorescence resonance energy transfer (FRET) measurements on freely diffusing single protein molecules. Properties of these subpopulations were directly monitored as a function of guanidinium chloride concentration. It is shown that new information about different aspects of the protein folding reaction can be extracted from such subpopulation properties. Shifts in the mean transfer efficiencies are discussed, FRET efficiency distributions are translated into potentials, and denaturation curves are directly plotted from the areas of the FRET peaks. Changes in stability caused by mutation also are measured by comparing pseudo wild-type CI2 with a destabilized mutant (K17G). Current limitations and future possibilities and prospects for single-pair FRET protein folding investigations are discussed.

Dynamics and folding of single two-stranded coiled-coil peptides studied by fluorescent energy transfer confocal microscopy

We report single-molecule measurements on the folding and unfolding conformational equilibrium distributions and dynamics of a disulfide crosslinked version of the two-stranded coiled coil from GCN4. The peptide has a fluorescent donor and acceptor at the N termini of its two chains and a Cys disulfide near its C terminus. Thus, folding brings the two N termini of the two chains close together, resulting in an enhancement of fluorescent resonant energy transfer. End-to-end distance distributions have thus been characterized under conditions where the peptide is nearly fully folded (0 M urea), unfolded (7.4 M urea), and in dynamic exchange between folded and unfolded states (3.0 M urea). The distributions have been compared for the peptide freely diffusing in solution and deposited onto aminopropyl silanized glass. As the urea concentration is increased, the mean end-to-end distance shifts to longer distances both in free solution and on the modified surface. The widths of these distributions indicate that the molecules are undergoing millisecond conformational fluctuations. Under all three conditions, these fluctuations gave nonexponential correlations on 1- to 100-ms time scale. A component of the correlation decay that was sensitive to the concentration of urea corresponded to that measured by bulk relaxation kinetics. The trajectories provided effective intramolecular diffusion coefficients as a function of the end-to-end distances for the folded and unfolded states. Single-molecule folding studies provide information concerning the distributions of conformational states in the folded, unfolded, and dynamically interconverting states.

An unusual pathway of excitation energy deactivation in carotenoids: Singlet-to-triplet conversion on an ultrafast timescale in a photosynthetic antenna

Carotenoids are important biomolecules that are ubiquitous in nature and find widespread application in medicine. In photosynthesis, they have a large role in light harvesting (LH) and photoprotection. They exert their LH function by donating their excited singlet state to nearby (bacterio)chlorophyll molecules. In photosynthetic bacteria, the efficiency of this energy transfer process can be as low as 30%. Here, we present evidence that an unusual pathway of excited state relaxation in carotenoids underlies this poor LH function, by which carotenoid triplet states are generated directly from carotenoid singlet states. This pathway, operative on a femtosecond and picosecond timescale, involves an intermediate state, which we identify as a new, hitherto uncharacterized carotenoid singlet excited state. In LH complex-bound carotenoids, this state is the precursor on the reaction pathway to the triplet state, whereas in extracted carotenoids in solution, this state returns to the singlet ground state without forming any triplets. We discuss the possible identity of this excited state and argue that fission of the singlet state into a pair of triplet states on individual carotenoid molecules constitutes the mechanism by which the triplets are generated. This is, to our knowledge, the first ever direct observation of a singlet-to-triplet conversion process on an ultrafast timescale in a photosynthetic antenna.

Overexpression of Arabidopsis Phytochrome B Inhibits Phytochrome A Function in the Presence of Sucrose1

Overexpression of phytochrome B (phyB) in Arabidopsis has previously been demonstrated to result in dominant negative interference of phytochrome A (phyA)-mediated hypocotyl growth inhibition in far-red (FR) light. This phenomenon has been examined further in this study and has been found to be dependent on the FR fluence rate and on the availability of metabolizable sugars in the growth medium. Poorly metabolized sugars capable of activating the putative hexokinase sensory function were not effective in eliciting the phytochrome interference response. Overexpressed phyB lacking the chromophore-binding site was also effective at inhibiting the phyA response, especially at higher fluence rates of FR. Overexpressed phyB produces the dominant negative phenotype without any apparent effect on phyA abundance or degradation. It is possible that phyA and phyB interact with a common reaction partner but that either the energy state of the cell or a separate sugar-signaling mechanism modulates the phytochrome-signaling interactions.

The Xanthophyll Cycle Modulates the Kinetics of Nonphotochemical Energy Dissipation in Isolated Light-Harvesting Complexes, Intact Chloroplasts, and Leaves of Spinach1

We analyzed the kinetics of nonphotochemical quenching of chlorophyll fluorescence (qN) in spinach (Spinacia oleracea) leaves, chloroplasts, and purified light-harvesting complexes. The characteristic biphasic pattern of fluorescence quenching in dark-adapted leaves, which was removed by preillumination, was evidence of light activation of qN, a process correlated with the de-epoxidation state of the xanthophyll cycle carotenoids. Chloroplasts isolated from dark-adapted and light-activated leaves confirmed the nature of light activation: faster and greater quenching at a subsaturating transthylakoid pH gradient. The light-harvesting chlorophyll a/b-binding complexes of photosystem II were isolated from dark-adapted and light-activated leaves. When isolated from light-activated leaves, these complexes showed an increase in the rate of quenching in vitro compared with samples prepared from dark-adapted leaves. In all cases, the quenching kinetics were fitted to a single component hyperbolic function. For leaves, chloroplasts, and light-harvesting complexes, the presence of zeaxanthin was associated with an increased rate constant for the induction of quenching. We discuss the significance of these observations in terms of the mechanism and control of qN.

A cardiac myocyte vascular endothelial growth factor paracrine pathway is required to maintain cardiac function

The role of the cardiac myocyte as a mediator of paracrine signaling in the heart has remained unclear. To address this issue, we generated mice with cardiac myocyte-specific deletion of the vascular endothelial growth factor gene, thereby producing a cardiomyocyte-specific knockout of a secreted factor. The hearts of these mice had fewer coronary microvessels, thinned ventricular walls, depressed basal contractile function, induction of hypoxia-responsive genes involved in energy metabolism, and an abnormal response to β-adrenergic stimulation. These findings establish the critical importance of cardiac myocyte-derived vascular endothelial growth factor in cardiac morphogenesis and determination of heart function. Further, they establish an adult murine model of hypovascular nonnecrotic cardiac contractile dysfunction.

VanX, a bacterial d-alanyl-d-alanine dipeptidase: Resistance, immunity, or survival function?

The zinc-containing d-alanyl-d-alanine (d-Ala-d-Ala) dipeptidase VanX has been detected in both Gram-positive and Gram-negative bacteria, where it appears to have adapted to at least three distinct physiological roles. In pathogenic vancomycin-resistant enterococci, vanX is part of a five-gene cluster that is switched on to reprogram cell-wall biosynthesis to produce peptidoglycan chain precursors terminating in d-alanyl-d-lactate (d-Ala-d-lactate) rather than d-Ala-d-Ala. The modified peptidoglycan exhibits a 1,000-fold decrease in affinity for vancomycin, accounting for the observed phenotypic resistance. In the glycopeptide antibiotic producers Streptomyces toyocaensis and Amylocatopsis orientalis, a vanHAX operon may have coevolved with antibiotic biosynthesis genes to provide immunity by reprogramming cell-wall termini to d-Ala-d-lactate as antibiotic biosynthesis is initiated. In the Gram-negative bacterium Escherichia coli, which is never challenged by the glycopeptide antibiotics because they cannot penetrate the outer membrane permeability barrier, the vanX homologue (ddpX) is cotranscribed with a putative dipeptide transport system (ddpABCDF) in stationary phase by the transcription factor RpoS (σs). The combined action of DdpX and the permease would permit hydrolysis of d-Ala-d-Ala transported back into the cytoplasm from the periplasm as cell-wall crosslinks are refashioned. The d-Ala product could then be oxidized as an energy source for cell survival under starvation conditions.

Sequence analysis of a mannitol dehydrogenase cDNA from plants reveals a function for the pathogenesis-related protein ELI3.

Mannitol is the most abundant sugar alcohol in nature, occurring in bacteria, fungi, lichens, and many species of vascular plants. Celery (Apium graveolens L.), a plant that forms mannitol photosynthetically, has high photosynthetic rates thought to results from intrinsic differences in the biosynthesis of hexitols vs. sugars. Celery also exhibits high salt tolerance due to the function of mannitol as an osmoprotectant. A mannitol catabolic enzyme that oxidizes mannitol to mannose (mannitol dehydrogenase, MTD) has been identified. In celery plants, MTD activity and tissue mannitol concentration are inversely related. MTD provides the initial step by which translocated mannitol is committed to central metabolism and, by regulating mannitol pool size, is important in regulating salt tolerance at the cellular level. We have now isolated, sequenced, and characterized a Mtd cDNA from celery. Analyses showed that Mtd RNA was more abundant in cells grown on mannitol and less abundant in salt-stressed cells. A protein database search revealed that the previously described ELI3 pathogenesis-related proteins from parsley and Arabidopsis are MTDs. Treatment of celery cells with salicylic acid resulted in increased MTD activity and RNA. Increased MTD activity results in an increased ability to utilize mannitol. Among other effects, this may provide an additional source of carbon and energy for response to pathogen attack. These responses of the primary enzyme controlling mannitol pool size reflect the importance of mannitol metabolism in plant responses to divergent types of environmental stress.