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.

Comparison of structures and energies of CH52+• with CH4+• and their possible role in superacidic methane activation

Contrary to previous theoretical studies at the UHF/6-31G* level, the methonium radical dication CH52+ is not a Cs symmetrical structure with a 2e—3c bond but a C2v symmetrical structure 1 with two 2e—3c bonds (at the UHF/6-31G**, UMP2/6-31G**, and UQCISD(T)/6-311G** levels). The Cs symmetrical structure is not even a minimum at the higher level of calculations. The four hydrogen atoms in 1 are bonded to the carbon atom by two 2e—3c bonds and the fifth hydrogen atom by a 2e—2c bond. The unpaired electron of 1 is located in a formal p-orbital (of the sp2-hybridized carbon atom) perpendicular to the plane of the molecule. Hydrogen scrambling in 1 is however extremely facile, as is in other C1 cations. It is found that the protonation of methane to CH5+ decreases the energy for subsequent homolytic cleavage resulting in the exothermic (24.1 kcal/mol) formation of CH4+•. Subsequent reaction with neutral methane while reforming CH5+ gives the methyl radical enabling reaction with excess methane to ethane and H2. The overall reaction is endothermic by 11.4 kcal/mol, but offers under conditions of oxidative removal of H2 an alternative to the more energetic carbocationic conversion of methane.

Fluorescence energy transfer detection as a homogeneous DNA diagnostic method

A homogeneous DNA diagnostic assay based on template-directed primer extension detected by fluorescence resonance energy transfer, named template-directed dye-terminator incorporation (TDI) assay, has been developed for mutation detection and high throughput genome analysis. Here, we report the successful application of the TDI assay to detect mutations in the cystic fibrosis transmembrane conductance regulator (CFTR) gene, the human leukocyte antigen H (HLA-H) gene, and the receptor tyrosin kinase (RET) protooncogene that are associated with cystic fibrosis, hemochromatosis, and multiple endocrine neoplasia, type 2, respectively. Starting with total human DNA, the samples are amplified by the PCR followed by enzymatic degradation of excess primers and deoxyribonucleoside triphosphates before the primer extension reaction is performed. All these standardized steps are performed in the same tube, and the fluorescence changes are monitored in real time, making it a useful clinical DNA diagnostic method.

Molecular dynamics and free-energy calculations applied to affinity maturation in antibody 48G7

We investigated the relative free energies of hapten binding to the germ line and mature forms of the 48G7 antibody Fab fragments by applying a continuum model to structures sampled from molecular dynamics simulations in explicit solvent. Reasonable absolute and very good relative free energies were obtained. As a result of nine somatic mutations that do not contact the hapten, the affinity-matured antibody binds the hapten >104 tighter than the germ line antibody. Energetic analysis reveals that van der Waals interactions and nonpolar contributions to solvation are similar and drive the formations of both the germ line and mature antibody–hapten complexes. Affinity maturation of the 48G7 antibody therefore appears to occur through reorganization of the combining site geometry in a manner that optimizes the balance of gaining favorable electrostatic interactions with the hapten and losing those with solvent during the binding process. As reflected by lower rms fluctuations in the antibody–hapten complex, the mature complex undergoes more restricted fluctuations than the germ line complex. The dramatically increased affinity of the 48G7 antibody over its germ line precursor is thus made possible by electrostatic optimization.

Uracil-DNA glycosylase–DNA substrate and product structures: Conformational strain promotes catalytic efficiency by coupled stereoelectronic effects

Enzymatic transformations of macromolecular substrates such as DNA repair enzyme/DNA transformations are commonly interpreted primarily by active-site functional-group chemistry that ignores their extensive interfaces. Yet human uracil–DNA glycosylase (UDG), an archetypical enzyme that initiates DNA base-excision repair, efficiently excises the damaged base uracil resulting from cytosine deamination even when active-site functional groups are deleted by mutagenesis. The 1.8-Å resolution substrate analogue and 2.0-Å resolution cleaved product cocrystal structures of UDG bound to double-stranded DNA suggest enzyme–DNA substrate-binding energy from the macromolecular interface is funneled into catalytic power at the active site. The architecturally stabilized closing of UDG enforces distortions of the uracil and deoxyribose in the flipped-out nucleotide substrate that are relieved by glycosylic bond cleavage in the product complex. This experimentally defined substrate stereochemistry implies the enzyme alters the orientation of three orthogonal electron orbitals to favor electron transpositions for glycosylic bond cleavage. By revealing the coupling of this anomeric effect to a delocalization of the glycosylic bond electrons into the uracil aromatic system, this structurally implicated mechanism resolves apparent paradoxes concerning the transpositions of electrons among orthogonal orbitals and the retention of catalytic efficiency despite mutational removal of active-site functional groups. These UDG/DNA structures and their implied dissociative excision chemistry suggest biology favors a chemistry for base-excision repair initiation that optimizes pathway coordination by product binding to avoid the release of cytotoxic and mutagenic intermediates. Similar excision chemistry may apply to other biological reaction pathways requiring the coordination of complex multistep chemical transformations.

The elasticity and failure of fluid-filled cellular solids: Theory and experiment

We extend and apply theories of filled foam elasticity and failure to recently available data on foods. The predictions of elastic modulus and failure mode dependence on internal pressure and on wall integrity are borne out by photographic evidence of distortion and failure under compressive loading and under the localized stress applied by a knife blade, and by mechanical data on vegetables differing only in their turgor pressure. We calculate the dry modulus of plate-like cellular solids and the cross over between dry-like and fully fluid-filled elastic response. The bulk elastic properties of limp and aging cellular solids are calculated for model systems and compared with our mechanical data, which also show two regimes of response. The mechanics of an aged, limp beam is calculated, thus offering a practical procedure for comparing experiment and theory. This investigation also thereby offers explanations of the connection between turgor pressure and crispness and limpness of cellular materials.

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.

Native protein sequences are close to optimal for their structures

How large is the volume of sequence space that is compatible with a given protein structure? Starting from random sequences, low free energy sequences were generated for 108 protein backbone structures by using a Monte Carlo optimization procedure and a free energy function based primarily on Lennard–Jones packing interactions and the Lazaridis–Karplus implicit solvation model. Remarkably, in the designed sequences 51% of the core residues and 27% of all residues were identical to the amino acids in the corresponding positions in the native sequences. The lowest free energy sequences obtained for ensembles of native-like backbone structures were also similar to the native sequence. Furthermore, both the individual residue frequencies and the covariances between pairs of positions observed in the very large SH3 domain family were recapitulated in core sequences designed for SH3 domain structures. Taken together, these results suggest that the volume of sequence space optimal for a protein structure is surprisingly restricted to a region around the native sequence.

Sulfur K-edge x-ray absorption spectroscopy: A spectroscopic tool to examine the redox state of S-containing metabolites in vivo

The sulfur K-edge x-ray absorption spectra for the amino acids cysteine and methionine and their corresponding oxidized forms cystine and methionine sulfoxide are presented. Distinct differences in the shape of the edge and the inflection point energy for cysteine and cystine are observed. For methionine sulfoxide the inflection point energy is 2.8 eV higher compared with methionine. Glutathione, the most abundant thiol in animal cells, also has been investigated. The x-ray absorption near-edge structure spectrum of reduced glutathione resembles that of cysteine, whereas the spectrum of oxidized glutathione resembles that of cystine. The characteristic differences between the thiol and disulfide spectra enable one to determine the redox status (thiol to disulfide ratio) in intact biological systems, such as unbroken cells, where glutathione and cyst(e)ine are the two major sulfur-containing components. The sulfur K-edge spectra for whole human blood, plasma, and erythrocytes are shown. The erythrocyte sulfur K-edge spectrum is similar to that of fully reduced glutathione. Simulation of the plasma spectrum indicated 32% thiol and 68% disulfide sulfur. The whole blood spectrum can be simulated by a combination of 46% disulfide and 54% thiol sulfur.

Mode matches and their locations in the hydrophobic free energy sequences of peptide ligands and their receptor eigenfunctions

Patterns in sequences of amino acid hydrophobic free energies predict secondary structures in proteins. In protein folding, matches in hydrophobic free energy statistical wavelengths appear to contribute to selective aggregation of secondary structures in “hydrophobic zippers.” In a similar setting, the use of Fourier analysis to characterize the dominant statistical wavelengths of peptide ligands’ and receptor proteins’ hydrophobic modes to predict such matches has been limited by the aliasing and end effects of short peptide lengths, as well as the broad-band, mode multiplicity of many of their frequency (power) spectra. In addition, the sequence locations of the matching modes are lost in this transformation. We make new use of three techniques to address these difficulties: (i) eigenfunction construction from the linear decomposition of the lagged covariance matrices of the ligands and receptors as hydrophobic free energy sequences; (ii) maximum entropy, complex poles power spectra, which select the dominant modes of the hydrophobic free energy sequences or their eigenfunctions; and (iii) discrete, best bases, trigonometric wavelet transformations, which confirm the dominant spectral frequencies of the eigenfunctions and locate them as (absolute valued) moduli in the peptide or receptor sequence. The leading eigenfunction of the covariance matrix of a transmembrane receptor sequence locates the same transmembrane segments seen in n-block-averaged hydropathy plots while leaving the remaining hydrophobic modes unsmoothed and available for further analyses as secondary eigenfunctions. In these receptor eigenfunctions, we find a set of statistical wavelength matches between peptide ligands and their G-protein and tyrosine kinase coupled receptors, ranging across examples from 13.10 amino acids in acid fibroblast growth factor to 2.18 residues in corticotropin releasing factor. We find that the wavelet-located receptor modes in the extracellular loops are compatible with studies of receptor chimeric exchanges and point mutations. A nonbinding corticotropin-releasing factor receptor mutant is shown to have lost the signatory mode common to the normal receptor and its ligand. Hydrophobic free energy eigenfunctions and their transformations offer new quantitative physical homologies in database searches for peptide-receptor matches.

Solvent effects on the energy landscapes and folding kinetics of polyalanine

The effect of a solvation on the thermodynamics and kinetics of polyalanine (Ala12) is explored on the basis of its energy landscapes in vacuum and in an aqueous solution. Both energy landscapes are characterized by two basins, one associated with α-helical structures and the other with coil and β-structures of the peptide. In both environments, the basin that corresponds to the α-helical structure is considerably narrower than the basin corresponding to the β-state, reflecting their different contributions to the entropy of the peptide. In vacuum, the α-helical state of Ala12 constitutes the native state, in agreement with common helical propensity scales, whereas in the aqueous medium, the α-helical state is destabilized, and the β-state becomes the native state. Thus solvation has a dramatic effect on the energy landscape of this peptide, resulting in an inverted stability of the two states. Different folding and unfolding time scales for Ala12 in hydrophilic and hydrophobic chemical environments are caused by the higher entropy of the native state in water relative to vacuum. The concept of a helical propensity has to be extended to incorporate environmental solvent effects.

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.

Crystal structures of the vitamin D receptor complexed to superagonist 20-epi ligands

The crystal structures of the ligand-binding domain (LBD) of the vitamin D receptor complexed to 1α,25(OH)2D3 and the 20-epi analogs, MC1288 and KH1060, show that the protein conformation is identical, conferring a general character to the observation first made for retinoic acid receptor (RAR) that, for a given LBD, the agonist conformation is unique, the ligands adapting to the binding pocket. In all complexes, the A- to D-ring moieties of the ligands adopt the same conformation and form identical contacts with the protein. Differences are observed only for the 17β-aliphatic chains that adapt their conformation to anchor the 25-hydroxyl group to His-305 and His-397. The inverted geometry of the C20 methyl group induces different paths of the aliphatic chains. The ligands exhibit a low-energy conformation for MC1288 and a more strained conformation for the two others. KH1060 compensates this energy cost by additional contacts. Based on the present data, the explanation of the superagonist effect is to be found in higher stability and longer half-life of the active complex, thereby excluding different conformations of the ligand binding domain.

Evidence for nonrandom hydrophobicity structures in protein chains.

The question of whether proteins originate from random sequences of amino acids is addressed. A statistical analysis is performed in terms of blocked and random walk values formed by binary hydrophobic assignments of the amino acids along the protein chains. Theoretical expectations of these variables from random distributions of hydrophobicities are compared with those obtained from functional proteins. The results, which are based upon proteins in the SWISS-PROT data base, convincingly show that the amino acid sequences in proteins differ from what is expected from random sequences in a statistically significant way. By performing Fourier transforms on the random walks, one obtains additional evidence for nonrandomness of the distributions. We have also analyzed results from a synthetic model containing only two amino acid types, hydrophobic and hydrophilic. With reasonable criteria on good folding properties in terms of thermodynamical and kinetic behavior, sequences that fold well are isolated. Performing the same statistical analysis on the sequences that fold well indicates similar deviations from randomness as for the functional proteins. The deviations from randomness can be interpreted as originating from anticorrelations in terms of an Ising spin model for the hydrophobicities. Our results, which differ from some previous investigations using other methods, might have impact on how permissive with respect to sequence specificity protein folding process is-only sequences with nonrandom hydrophobicity distributions fold well. Other distributions give rise to energy landscapes with poor folding properties and hence did not survive the evolution.

Temporally and spectrally resolved subpicosecond energy transfer within the peripheral antenna complex (LH2) and from LH2 to the core antenna complex in photosynthetic purple bacteria.

We report studies of energy transfer from the 800-nm absorbing pigment (B800) to the 850-nm absorbing pigment (B850) of the LH2 peripheral antenna complex and from LH2 to the core antenna complex (LH1) in Rhodobacter (Rb.) sphaeroides. The B800 to B850 process was studied in membranes from a LH2-reaction center (no LH1) mutant of Rb. sphaeroides and the LH2 to LH1 transfer was studied in both the wild-type species and in LH2 mutants with blue-shifted B850. The measurements were performed by using approximately 100-fs pulses to probe the formation of acceptor excitations in a two-color pump-probe measurement. Our experiments reveal a B800 to B850 transfer time of approximately 0.7 ps at 296 K and energy transfer from LH2 to LH1 is characterized by a time constant of approximately 3 ps at 296 K and approximately 5 ps at 77 K. In the blue-shifted B850 mutants, the transfer time from B850 to LH1 becomes gradually longer with increasing blue-shift of the B850 band as a result of the decreasing spectral overlap between the antennae. The results have been used to produce a model for the association between the ring-like structures that are characteristic of both the LH2 and LH1 antennae.