Biomedical subjects

T Head-Gordon

Publications and source records attributed to T Head-Gordon.

17 recordsLinked to original sources

PubMed2000

A global optimization strategy for predicting alpha-helical protein tertiary structure.

We present a global optimization strategy that incorporates predicted restraints in both a local optimization context and as directives for global optimization approaches, to predict protein tertiary structure for alpha-helical proteins. Specifically, neural networks are used to predict the secondary structure of a protein, restraints are defined as manifestations of the network with a predicted secondary structure and the secondary structure is formed using local minimizations on a protein energy surface, in the presence of the restraints. Those residues predicted to be coil, by the network, define a conformational sub-space that is subject to optimization using a global approach known as stochastic perturbation that has been found to be effective for Lennard-Jones clusters and homo-polypeptides. Our energy surface is an all-atom 'gas phase' molecular mechanics force field, that is combined with a new solvation energy function that penalizes hydrophobic group exposure. This energy function gives the crystal structure of four different alpha-helical proteins as the lowest energy structure relative to other conformations, with correct secondary structure but incorrect tertiary structure. We demonstrate this global optimization strategy by determining the tertiary structure of the A-chain of the alpha-helical protein, uteroglobin and of a four-helix bundle, DNA binding protein.

Algorithms↗

PubMed2000

Matching simulation and experiment: a new simplified model for simulating protein folding.

Simulations of simplified protein folding models have provided much insight into solving the protein folding problem. We propose here a new off-lattice bead model, capable of simulating several different fold classes of small proteins. We present the sequence for an alpha/beta protein resembling the IgG-binding proteins L and G. The thermodynamics of the folding process for this model are characterized using the multiple multihistogram method combined with constant-temperature Langevin simulations. The folding is shown to be highly cooperative, with chain collapse nearly accompanying folding. Two parallel folding pathways are shown to exist on the folding free energy landscape. One pathway contains an intermediate--similar to experiments on protein G, and one pathway contains no intermediates-similar to experiments on protein L. The folding kinetics are characterized by tabulating mean-first passage times, and we show that the onset of glasslike kinetics occurs at much lower temperatures than the folding temperature. This model is expected to be useful in many future contexts: investigating questions of the role of local versus nonlocal interactions in various fold classes, addressing the effect of sequence mutations affecting secondary structure propensities, and providing a computationally feasible model for studying the role of solvation forces in protein folding.

Computational Biology↗

PubMed1999

Redesigning the hydrophobic core of a model beta-sheet protein: destabilizing traps through a threading approach.

An off-lattice 46-bead model of a small all-beta protein has been recently criticized for possessing too many traps and long-lived intermediates compared with the folding energy landscape predicted for real proteins and models using the principle of minimal frustration. Using a novel sequence design approach based on threading for finding beneficial mutations for destabilizing traps, we proposed three new sequences for folding in the beta-sheet model. Simulated annealing on these sequences found the global minimum more reliably, indicative of a smoother energy landscape, and simulated thermodynamic variables found evidence for a more cooperative collapse transition, lowering of the collapse temperature, and higher folding temperatures. Folding and unfolding kinetics were acquired by calculating first-passage times, and the new sequences were found to fold significantly faster than the original sequence, with a concomitant lowering of the glass temperature, although none of the sequences have highly stable native structures. The new sequences found here are more representative of real proteins and are good folders in the T(f) > T(g) sense, and they should prove useful in future studies of the details of transition states and the nature of folding intermediates in the context of simplified folding models. These results show that our sequence design approach using threading can improve models possessing glasslike folding dynamics.

Drug Design↗

PubMed1999

Chemical bonding effects in the determination of protein structures by electron crystallography.

Scattering of electrons is affected by the distribution of valence electrons that participate in chemical bonding and thus change the electrostatic shielding of the nucleus. This effect is particularly significant for low-angle scattering. Thus, while chemical bonding effects are difficult to measure with small-unit cell materials, they can be substantial in the study of proteins by electron crystallography. This work investigates the magnitude of chemical bonding effects for a representative collection of protein fragments and a model ligand for nucleotide-binding proteins within the resolution range generally used in determining protein structures by electron crystallography. Electrostatic potentials were calculated by ab initio methods for both the test molecules and for superpositions of their free atoms. Differences in scattering amplitudes can be well over 10% in the resolution range below 5 A and are especially large in the case of ionized side chains and ligands. We conclude that the use of molecule-based scattering factors can provide a much more accurate representation of the low-resolution data obtained in electron crystallographic studies. The comparison of neutral and ionic structure factors at resolutions below 5 A can also provide a sensitive determination of charge states, important for biological function, that is not accessible from X-ray crystallographic measurements.

Crystallography↗

PubMed1999

Evidence for microscopic, long-range hydration forces for a hydrophobic amino acid.

We have combined neutron solution scattering experiments with molecular dynamics simulation to isolate an excess experimental signal that is caused solely by N-acetyl-leucine-amide (NALA) correlations in aqueous solution. This excess signal contains information about how NALA molecule centers are correlated in water, and we show how these solute-solute correlations might be determined at dilute concentrations in the small angle region. We have tested qualitatively different pair distribution functions for NALA molecule centers-gas, cluster, and aqueous forms of gc(r)-and have found that the excess experimental signal is adequate enough to rule out gas and cluster pair distribution functions. The aqueous form of gc(r) that exhibits a solvent-separated minimum, and possibly longer-ranged correlations as well, is not only physically sound but reproduces the experimental data reasonably well. This work demonstrates that important information in the small angle region can be mined to resolve solute-solute correlations, their lengthscales, and thermodynamic consequences even at dilute concentrations. The hydration forces that operate on the microscopic scale of individual amino acid side chains, implied by the small angle scattering data, could have significant effects on the early stages of protein folding, on ligand binding, and on other intermolecular interactions.

Amides↗

PubMed1998

The importance of hydration for the kinetics and thermodynamics of protein folding: simplified lattice models.

BACKGROUND: Recent studies have proposed various sources for the origin of cooperativity in simplified protein folding models. Important contributions to cooperativity that have been discussed include backbone hydrogen bonding, sidechain packing and hydrophobic interactions. Related work has also focused on which interactions are responsible for making the free energy of the native structure a pronounced global minimum in the free energy landscape. In addition, two-flavor bead models have been found to exhibit poor folding cooperativity and often lack unique native structures. We propose a simple multibody description of hydration with expectations that it might modify the free energy surface in such a way as to increase the cooperativity of folding and improve the performance of two-flavor models. RESULTS: We study the thermodynamics and kinetics of folding for designed 36-mer sequences on a cubic lattice using both our solvation model and the corresponding model without solvation terms. Degeneracies of the native states are studied by enumerating the maximally compact states. The histogram Monte Carlo method is used to obtain folding temperatures, densities of states and heat capacity curves. Folding kinetics are examined by accumulating mean first-passage times versus temperature. Sequences in the proposed solvation model are found to have more unique ground states, fold faster and fold with more cooperativity than sequences in the nonsolvation model. CONCLUSIONS: We find that the addition of a multibody description of solvation can improve the poor performance of two-flavor lattice models and provide an additional source for more cooperative folding. Our results suggest that a better description of solvation will be important for future theoretical protein folding studies.

Animals↗

PubMed1997

Differences in hydration structure near hydrophobic and hydrophilic amino acids.

We use molecular dynamics to simulate recent neutron scattering experiments on aqueous solutions of N-acetyl-leucine-amide and N-acetyl-glutamine-amide, and break down the total scattering function into contributions from solute-solute, solute-water, water-water, and intramolecular correlations. We show that the shift of the main diffraction peak to smaller angle that is observed for leucine, but not for glutamine, is attributable primarily to alterations in water-water correlations relative to bulk. The perturbation of the water hydrogen-bonded network extends roughly two solvation layers from the hydrophobic side chain surface, and is characterized by a distribution of hydrogen bonded ring sizes that are more planar and are dominated by pentagons in particular than those near the hydrophilic side chain. The different structural organization of water near the hydrophobic solute that gives rise to the inward shift in the main neutron diffraction peak under ambient conditions may also provide insight into the same directional shift for pure liquid water as it is cooled and supercooled.

Amino Acids↗

PubMed1996

Direct evidence for modified solvent structure within the hydration shell of a hydrophobic amino acid.

Neutron scattering experiments are used to determine scattering profiles for aqueous solutions of hydrophobic and hydrophilic amino acid analogs. Solutions of hydrophobic solutes show a shift in the main diffraction peak to smaller angle as compared with pure water, whereas solutions of hydrophilic solutes do not. The same difference for solutions of hydrophobic and hydrophilic side chains is also predicted by molecular dynamics simulations. The neutron scattering curves of aqueous solutions of hydrophobic amino acids at room temperature are qualitatively similar to differences between the liquid molecular structure functions measured for ambient and supercooled water. The nonpolar solute-induced expansion of water structure reported here is also complementary to recent neutron experiments where compression of aqueous solvent structure has been observed at high salt concentration.

Amides↗

PubMed1995

Is water structure around hydrophobic groups clathrate-like?

The term "clathrate structure" is quantified for solvation of nonpolar groups by enumerating hydrogen-bonded ring sizes both in the solvation shell and through the shell-bulk interface and comparing it to a bulk control using the ST4 water model. For clathrate-like structure to be evident, the distributions along the hydrophobic surface are expected to be dominated by pentagons, with significant depletion of hexagons and larger polygons. While the distribution in this region is indeed distinguished by a large number of pentagons, there are significant contributions from hexagons and larger rings as well. Calculated polygon distributions through the shell-bulk interface indicate that when water structure is highly cooperative along the hydrophobic surface, hydrogen-bonded pathways leading back into bulk are then reduced. These results are qualitatively consistent with the observation that hydrophobicity is proportional to the nonpolar solute surface area.

Journal Article↗

PubMed1995