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.

Top-down free-energy minimization on protein potential energy landscapes

The hierarchical properties of potential energy landscapes have been used to gain insight into thermodynamic and kinetic properties of protein ensembles. It also may be possible to use them to direct computational searches for thermodynamically stable macroscopic states, i.e., computational protein folding. To this end, we have developed a top-down search procedure in which conformation space is recursively dissected according to the intrinsic hierarchical structure of a landscape's effective-energy barriers. This procedure generates an inverted tree similar to the disconnectivity graphs generated by local minima-clustering methods, but it fundamentally differs in the manner in which the portion of the tree that is to be computationally explored is selected. A key ingredient is a branch-selection algorithm that takes advantage of statistically predictive properties of the landscape to guide searches down the tree branches that are most likely to lead to the physically relevant macroscopic states. Using the computational folding of a β-hairpin-forming peptide as an example, we show that such predictive properties indeed exist and can be used for structure prediction by free-energy global minimization.

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.

Linking topography of its potential surface with the dynamics of folding of a protein model

The “3-color, 46-bead” model of a folding polypeptide is the vehicle for adapting to proteins a mode of analysis used heretofore for atomic clusters, to relate the topography of the potential surface to the dynamics that lead to formation of selected structures. The analysis is based on sequences of stationary points—successive minima, joined by saddles—that rise monotonically in energy from basin bottoms. Like structure-seeking clusters, the potential surface of the model studied here is staircase-like, rather than sawtooth-like, with highly collective motions required for passage from one minimum to the next. The surface has several deep basins whose minima correspond to very similar structures, but which are separated by high energy barriers.

Acute decrease in net glutamate uptake during energy deprivation

The extracellular glutamate concentration ([glu]o) rises during cerebral ischemia, reaching levels capable of inducing delayed neuronal death. The mechanisms underlying this glutamate accumulation remain controversial. We used N-methyl-d-aspartate receptors on CA3 pyramidal neurons as a real-time, on-site, glutamate sensor to identify the source of glutamate release in an in vitro model of ischemia. Using glutamate and l-trans-pyrrolidine-2,4-dicarboxylic acid (tPDC) as substrates and dl-threo-β-benzyloxyaspartate (TBOA) as an inhibitor of glutamate transporters, we demonstrate that energy deprivation decreases net glutamate uptake within 2–3 min and later promotes reverse glutamate transport. This process accounts for up to 50% of the glutamate accumulation during energy deprivation. Enhanced action potential-independent vesicular release also contributes to the increase in [glu]o, by ≈50%, but only once glutamate uptake is inhibited. These results indicate that a significant rise in [glu]o already occurs during the first minutes of energy deprivation and is the consequence of reduced uptake and increased vesicular and nonvesicular release of glutamate.

Energy-based de novo protein folding by conformational space annealing and an off-lattice united-residue force field: Application to the 10-55 fragment of staphylococcal protein A and to apo calbindin D9K

The conformational space annealing (CSA) method for global optimization has been applied to the 10-55 fragment of the B-domain of staphylococcal protein A (protein A) and to a 75-residue protein, apo calbindin D9K (PDB ID code 1CLB), by using the UNRES off-lattice united-residue force field. Although the potential was not calibrated with these two proteins, the native-like structures were found among the low-energy conformations, without the use of threading or secondary-structure predictions. This is because the CSA method can find many distinct families of low-energy conformations. Starting from random conformations, the CSA method found that there are two families of low-energy conformations for each of the two proteins, the native-like fold and its mirror image. The CSA method converged to the same low-energy folds in all cases studied, as opposed to other optimization methods. It appears that the CSA method with the UNRES force field, which is based on the thermodynamic hypothesis, can be used in prediction of protein structures in real time.

Relationship between the oxidation potential and electron spin density of the primary electron donor in reaction centers from Rhodobacter sphaeroides

The primary electron donor in bacterial reaction centers is a dimer of bacteriochlorophyll a molecules, labeled L or M based on their proximity to the symmetry-related protein subunits. The electronic structure of the bacteriochlorophyll dimer was probed by introducing small systematic variations in the bacteriochlorophyll–protein interactions by a series of site-directed mutations that replaced residue Leu M160 with histidine, tyrosine, glutamic acid, glutamine, aspartic acid, asparagine, lysine, and serine. The midpoint potentials for oxidation of the dimer in the mutants showed an almost continuous increase up to ≈60 mV compared with wild type. The spin density distribution of the unpaired electron in the cation radical state of the dimer was determined by electron–nuclear–nuclear triple resonance spectroscopy in solution. The ratio of the spin density on the L side of the dimer to the M side varied from ≈2:1 to ≈5:1 in the mutants compared with ≈2:1 for wild type. The correlation between the midpoint potential and spin density distribution was described using a simple molecular orbital model, in which the major effect of the mutations is assumed to be a change in the energy of the M half of the dimer, providing estimates for the coupling and energy levels of the orbitals in the dimer. These results demonstrate that the midpoint potential can be fine-tuned by electrostatic interactions with amino acids near the dimer and show that the properties of the electronic structure of a donor or acceptor in a protein complex can be directly related to functional properties such as the oxidation–reduction midpoint potential.

Ecological intensification of cereal production systems: Yield potential, soil quality, and precision agriculture

Wheat (Triticum aestivum L.), rice (Oryza sativa L.), and maize (Zea mays L.) provide about two-thirds of all energy in human diets, and four major cropping systems in which these cereals are grown represent the foundation of human food supply. Yield per unit time and land has increased markedly during the past 30 years in these systems, a result of intensified crop management involving improved germplasm, greater inputs of fertilizer, production of two or more crops per year on the same piece of land, and irrigation. Meeting future food demand while minimizing expansion of cultivated area primarily will depend on continued intensification of these same four systems. The manner in which further intensification is achieved, however, will differ markedly from the past because the exploitable gap between average farm yields and genetic yield potential is closing. At present, the rate of increase in yield potential is much less than the expected increase in demand. Hence, average farm yields must reach 70–80% of the yield potential ceiling within 30 years in each of these major cereal systems. Achieving consistent production at these high levels without causing environmental damage requires improvements in soil quality and precise management of all production factors in time and space. The scope of the scientific challenge related to these objectives is discussed. It is concluded that major scientific breakthroughs must occur in basic plant physiology, ecophysiology, agroecology, and soil science to achieve the ecological intensification that is needed to meet the expected increase in food demand.

Potential effects of gas hydrate on human welfare

For almost 30 years. serious interest has been directed toward natural gas hydrate, a crystalline solid composed of water and methane, as a potential (i) energy resource, (ii) factor in global climate change, and (iii) submarine geohazard. Although each of these issues can affect human welfare, only (iii) is considered to be of immediate importance. Assessments of gas hydrate as an energy resource have often been overly optimistic, based in part on its very high methane content and on its worldwide occurrence in continental margins. Although these attributes are attractive, geologic settings, reservoir properties, and phase-equilibria considerations diminish the energy resource potential of natural gas hydrate. The possible role of gas hydrate in global climate change has been often overstated. Although methane is a “greenhouse” gas in the atmosphere, much methane from dissociated gas hydrate may never reach the atmosphere, but rather may be converted to carbon dioxide and sequestered by the hydrosphere/biosphere before reaching the atmosphere. Thus, methane from gas hydrate may have little opportunity to affect global climate change. However, submarine geohazards (such as sediment instabilities and slope failures on local and regional scales, leading to debris flows, slumps, slides, and possible tsunamis) caused by gas-hydrate dissociation are of immediate and increasing importance as humankind moves to exploit seabed resources in ever-deepening waters of coastal oceans. The vulnerability of gas hydrate to temperature and sea level changes enhances the instability of deep-water oceanic sediments, and thus human activities and installations in this setting can be affected.

Kinetics of photo-induced electron transfer from high-potential iron-sulfur protein to the photosynthetic reaction center of the purple phototroph Rhodoferax fermentans.

The kinetics of photo-induced electrontransfer from high-potential iron-sulfur protein (HiPIP) to the photosynthetic reaction center (RC) of the purple phototroph Rhodoferarfermentans were studied. The rapid photooxidation of heme c-556 belonging to RC is followed, in the presence of HiPIP, by a slower reduction having a second-order rate constant of 4.8 x 10(7) M(-1) x s(-1). The limiting value of kobs at high HiPIP concentration is 95 s(-1). The amplitude of this slow process decreases with increasing HiPIP concentration. The amplitude of a faster phase, observed at 556 and 425 nm and involving heme c-556 reduction, increases proportionately. The rate constant of this fast phase, determined at 425 and 556 nm, is approximately 3 x 10(5) s(-1). This value is not dependent on HiPIP concentration, indicating that it is related to a first-order process. These observations are interpreted as evidence for the formation of a HiPIP-RC complex prior to the excitation flash, having a dissociation constant of -2.5 microM. The fast phase is absent at high ionic strength, indicating that the complex involves mainly electrostatic interactions. The ionic strength dependence of kobs for the slow phase yields a second-order rate constant at infinite ionic strength of 5.4 x 10(6) M(-1) x s(-1) and an electrostatic interaction energy of -2.1 kcal/mol (1 cal = 4.184 J). We conclude that Rhodoferar fermentans HiPIP is a very effective electron donor to the photosynthetic RC.

High resolution ultrastructural mapping of total calcium: electron spectroscopic imaging/electron energy loss spectroscopy analysis of a physically/chemically processed nerve-muscle preparation.

We report on a procedure for tissue preparation that combines thoroughly controlled physical and chemical treatments: quick-freezing and freeze-drying followed by fixation with OsO4 vapors and embedding by direct resin infiltration. Specimens of frog cutaneous pectoris muscle thus prepared were analyzed for total calcium using electron spectroscopic imaging/electron energy loss spectroscopy (ESI/EELS) approach. The preservation of the ultrastructure was excellent, with positive K/Na ratios revealed in the fibers by x-ray microanalysis. Clear, high-resolution EELS/ESI calcium signals were recorded from the lumen of terminal cisternae of the sarcoplasmic reticulum but not from longitudinal cisternae, as expected from previous studies carried out with different techniques. In many mitochondria, calcium was below detection whereas in others it was appreciable although at variable level. Within the motor nerve terminals, synaptic vesicles as well as some cisternae of the smooth endoplasmic reticulum yielded positive signals at variance with mitochondria, that were most often below detection. Taken as a whole, the present study reveals the potential of our experimental approach to map with high spatial resolution the total calcium within individual intracellular organelles identified by their established ultrastructure, but only where the element is present at high levels.

How optimization of potential functions affects protein folding.

The relationship between the optimization of the potential function and the foldability of theoretical protein models is studied based on investigations of a 27-mer cubic-lattice protein model and a more realistic lattice model for the protein crambin. In both the simple and the more complicated systems, optimization of the energy parameters achieves significant improvements in the statistical-mechanical characteristics of the systems and leads to foldable protein models in simulation experiments. The foldability of the protein models is characterized by their statistical-mechanical properties--e.g., by the density of states and by Monte Carlo folding simulations of the models. With optimized energy parameters, a high level of consistency exists among different interactions in the native structures of the protein models, as revealed by a correlation function between the optimized energy parameters and the native structure of the model proteins. The results of this work are relevant to the design of a general potential function for folding proteins by theoretical simulations.

Purified inositol hexakisphosphate kinase is an ATP synthase: diphosphoinositol pentakisphosphate as a high-energy phosphate donor.

Diphosphoinositol pentakisphosphate (PP-IP5) and bis(diphospho)inositol tetrakisphosphate (bis-PP-IP4) are recently identified inositol phosphates that possess pyrophosphate bonds. We have purified an inositol hexakisphosphate (IP6) kinase from rat brain supernatants. The pure protein, a monomer of 54 kDa, displays high affinity (Km = 0.7 microM) and selectivity for inositol hexakisphosphate as substrate. It can be dissociated from bis(diphospho)inositol tetrakisphosphate synthetic activity. The purified enzyme transfers a phosphate from PP-IP5 to ADP to form ATP. This ATP synthase activity indicates the high phosphate group transfer potential of PP-IP5 and may represent a physiological role for PP-IP5.

Chronic mitochondrial energy impairment produces selective striatal degeneration and abnormal choreiform movements in primates.

Although the gene defect responsible for Huntington disease (HD) has recently been identified, the pathogenesis of the disease remains obscure. One potential mechanism is that the gene defect may lead to an impairment of energy metabolism followed by slow excitotoxic neuronal injury. In the present study we examined whether chronic administration of 3-nitropropionic acid (3-NP), an irreversible inhibitor of succinate dehydrogenase, can replicate the neuropathologic and clinical features of HD in nonhuman primates. After 3-6 weeks of 3-NP administration, apomorphine treatment induced a significant increase in motor activity as compared with saline-treated controls. Animals showed both choreiform movements, as well as foot and limb dystonia, which are characteristic of HD. More prolonged 3-NP treatment in two additional primates resulted in spontaneous dystonia and dyskinesia accompanied by lesions in the caudate and putamen seen by magnetic resonance imaging. Histologic evaluation showed that there was a depletion of calbindin neurons, astrogliosis, sparing of NADPH-diaphorase neurons, and growth-related proliferative changes in dendrites of spiny neurons similar to changes in HD. The striosomal organization of the striatum and the nucleus accumbens were spared. These findings show that chronic administration of 3-NP to nonhuman primates can replicate many of the characteristic motor and histologic features of HD, further strengthening the possibility that a subtle impairment of energy metabolism may play a role in its pathogenesis.