Problems In Structural Biology Research Articles

CHOP: parsing proteins into structural domains.

Sequence-based domain assignment is one of the most important and challenging problems in structural biology. We have developed a method, CHOP, that chops proteins into domain-like fragments. The basic idea is to cut proteins from entirely sequenced organisms beginning from very reliable experimental information (Protein Data Bank), proceeding to expert annotations of domain-like regions (Pfam-A) and completing through cuts based on termini of native protein ends. The CHOP server takes protein sequences as input and returns the dissections supported by homology transfer. CHOP results are precompiled for many entirely sequenced proteomes. The service is available at http://www.rostlab.org/services/CHOP/.

Identification of Protein–Protein Interaction Sites from Docking Energy Landscapes

Protein recognition is one of the most challenging and intriguing problems in structural biology. Despite all the available structural, sequence and biophysical information about protein–protein complexes, the physico-chemical patterns, if any, that make a protein surface likely to be involved in protein–protein interactions, remain elusive. Here, we apply protein docking simulations and analysis of the interaction energy landscapes to identify protein–protein interaction sites. The new protocol for global docking based on multi-start global energy optimization of an all-atom model of the ligand, with detailed receptor potentials and atomic solvation parameters optimized in a training set of 24 complexes, explores the conformational space around the whole receptor without restrictions. The ensembles of the rigid-body docking solutions generated by the simulations were subsequently used to project the docking energy landscapes onto the protein surfaces. We found that highly populated low-energy regions consistently corresponded to actual binding sites. The procedure was validated on a test set of 21 known protein–protein complexes not used in the training set. As much as 81% of the predicted high-propensity patch residues were located correctly in the native interfaces. This approach can guide the design of mutations on the surfaces of proteins, provide geometrical details of a possible interaction, and help to annotate protein surfaces in structural proteomics.

A model for the cooperative free energy transduction and kinetics of ATP hydrolysis by F1-ATPase.

Although the binding change mechanism of rotary catalysis by which F1-ATPase hydrolyzes ATP has been supported by equilibrium, kinetic, and structural observations, many questions concerning the function remain unanswered. Because of the importance of this enzyme, the search for a full understanding of its mechanism is a key problem in structural biology. Making use of the results of free energy simulations and experimental binding constant measurements, a model is developed for the free energy change during the hydrolysis cycle. This model makes possible the development of a kinetic scheme for ATP hydrolysis by F1-ATPase, in which the rate constants are associated with specific configurations of the beta subunits. An essential new element is that the strong binding site for ADP,Pi is shown to be the betaDP site, in contrast to the strong binding site for ATP, which is betaTP. This result provides a rationale for the rotation of the gamma subunit, which induces the cooperativity required for a tri-site binding change mechanism. The model explains a series of experimental data, including the ATP concentration dependence of the rate of hydrolysis and catalytic site occupation for both the Escherichia coli F1-ATPase (EcF1) and Thermophilic Bacillus PS3 F1-ATPase (TF1), which have different behavior.

Rapid protein domain assignment from amino acid sequence using predicted secondary structure.

The elucidation of the domain content of a given protein sequence in the absence of determined structure or significant sequence homology to known domains is an important problem in structural biology. Here we address how successfully the delineation of continuous domains can be accomplished in the absence of sequence homology using simple baseline methods, an existing prediction algorithm (Domain Guess by Size), and a newly developed method (DomSSEA). The study was undertaken with a view to measuring the usefulness of these prediction methods in terms of their application to fully automatic domain assignment. Thus, the sensitivity of each domain assignment method was measured by calculating the number of correctly assigned top scoring predictions. We have implemented a new continuous domain identification method using the alignment of predicted secondary structures of target sequences against observed secondary structures of chains with known domain boundaries as assigned by Class Architecture Topology Homology (CATH). Taking top predictions only, the success rate of the method in correctly assigning domain number to the representative chain set is 73.3%. The top prediction for domain number and location of domain boundaries was correct for 24% of the multidomain set (+/-20 residues). These results have been put into context in relation to the results obtained from the other prediction methods assessed.

Three key residues form a critical contact network in a protein folding transition state.

Determining how a protein folds is a central problem in structural biology. The rate of folding of many proteins is determined by the transition state, so that a knowledge of its structure is essential for understanding the protein folding reaction. Here we use mutation measurements--which determine the role of individual residues in stabilizing the transition state--as restraints in a Monte Carlo sampling procedure to determine the ensemble of structures that make up the transition state. We apply this approach to the experimental data for the 98-residue protein acylphosphatase, and obtain a transition-state ensemble with the native-state topology and an average root-mean-square deviation of 6 A from the native structure. Although about 20 residues with small positional fluctuations form the structural core of this transition state, the native-like contact network of only three of these residues is sufficient to determine the overall fold of the protein. This result reveals how a nucleation mechanism involving a small number of key residues can lead to folding of a polypeptide chain to its unique native-state structure.

Deciphering the folding kinetics of transmembrane helical proteins.

Nearly a quarter of genomic sequences and almost half of all receptors that are likely to be targets for drug design are integral membrane proteins. Understanding the detailed mechanisms of the folding of membrane proteins is a largely unsolved, key problem in structural biology. Here, we introduce a general model and use computer simulations to study the equilibrium properties and the folding kinetics of a C(alpha)-based two-helix bundle fragment (comprised of 66 aa) of bacteriorhodopsin. Various intermediates are identified and their free energy are calculated together with the free energy barrier between them. In 40% of folding trajectories, the folding rate is considerably increased by the presence of nonobligatory intermediates acting as traps. In all cases, a substantial portion of the helices is rapidly formed. This initial stage is followed by a long period of consolidation of the helices accompanied by their correct packing within the membrane. Our results provide the framework for understanding the variety of folding pathways of helical transmembrane proteins.

Multiple quantum solid-state NMR indicates a parallel, not antiparallel, organization of beta-sheets in Alzheimer's beta-amyloid fibrils.

Senile plaques associated with Alzheimer's disease contain deposits of fibrils formed by 39- to 43-residue beta-amyloid peptides with possible neurotoxic effects. X-ray diffraction measurements on oriented fibril bundles have indicated an extended beta-sheet structure for Alzheimer's beta-amyloid fibrils and other amyloid fibrils, but the supramolecular organization of the beta-sheets and other structural details are not well established because of the intrinsically noncrystalline, insoluble nature of amyloid fibrils. Here we report solid-state NMR measurements, using a multiple quantum (MQ) (13)C NMR technique, that probe the beta-sheet organization in fibrils formed by the full-length, 40-residue beta-amyloid peptide (Abeta(1-40)). Although an antiparallel beta-sheet organization often is assumed and is invoked in recent structural models for full-length beta-amyloid fibrils, the MQNMR data indicate an in-register, parallel organization. This work provides site-specific, atomic-level structural constraints on full-length beta-amyloid fibrils and applies MQNMR to a significant problem in structural biology.

Interpreting the folding kinetics of helical proteins.

The detailed mechanism of protein folding is one of the major problems in structural biology. Its solution is of practical as well as fundamental interest because of its possible role in utilizing the many sequences becoming available from genomic analysis. Although the Levinthal paradox (namely, that a polypeptide chain can find its unique native state in spite of the astronomical number of configurations in the denatured state) has been resolved, the reasons for the differences in the folding behaviour of individual proteins remain to be elucidated. Here a Calpha-based three-helix-bundle-like protein model with a realistic thermodynamic phase diagram is used to calculate several hundred folding trajectories. By varying a single parameter, the difference between the strength of native and non-native contacts, folding is changed from a diffusion-collision mechanism to one that involves simultaneous collapse and partial secondary-structure formation, followed by reorganization to the native structure. Non-obligatory intermediates are important in the former, whereas there is an obligatory on-pathway intermediate in the latter. Our results provide a basis for understanding the range of folding behaviour that is observed in helical proteins.

Hartree−Fock and Møller−Plesset (MP2) Treatment of Oxygen-Containing Phosphorus Compounds

Ab initio calculations were carried out in order to obtain conformational, structural, and vibrational data for a set of small oxygen-containing phosphorus compounds. Many of these molecules are fundamentally important since they are analogous to macrostructures that play a central role in biological processes. It is anticipated that the chemical and physical properties of these phosphorus compounds will provide valuable insight into problems in structural biology. These calculations involved the use of the standard 6-31G** basis set at both the Hartree−Fock and Møller−Plesset levels of theory for the determination of rotational profiles, molecular geometries, and vibrational frequencies. Comparisons to experimental results are made when possible.

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS)

Matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) has been responsible for solving many problems in structural biology. Mass analysis is now used routinely to confirm proper expression and processing of proteins, and to locate and identify post-translational modifications. Innovative advances in instrumentation have led to higher mass resolution and mass accuracy. New sample preparation methods are likewise yielding higher sensitivity plus greater tolerance for buffer components that have in the past suppressed signals at higher concentrations. Advancements in the technique have also led to new or improved applications in many areas, including peptide sequencing and the identification of proteins by database searching with peptide masses. Instruments with lower cost, smaller size, and higher performance are making mass measurements available to an increasing number of laboratories. MALDI-MS is poised to continue to improve in performance and in its usefulness for current and new applications.

Toward the Accurate Modeling of DNA: The Importance of Long-Range Electrostatics

The relationship between DNA structure and function is fundamental to the understanding of biological processes. Currently, the most reliable source of biomolecular structural information comes from X-ray crystallographic data. Recent advances in theoretical modeling techniques have allowed molecular simulations to approach the accuracy obtainable from X-ray crystallography for proteins. We report the results of a 2.2 ns simulation of the B-DNA dodecamer d[CGCGAATTCGCG12 in a crystal unit cell, and demonstrate that, with rigorous accommodation of long-range forces, molecular simulation may be extended to provide atomic level accuracy of polynucleotide structures. Most of what is known about DNA structure has been derived from crystal structure analyses. A concern that has been expressed brings into question the degree to which the intrinsic structure of DNA is disrupted by the crystal lattice. It has been recently pointed out by Dickerson’ that lattice forces in DNA crystals are relatively weak, and that in fact these forces can be exploited to obtain useful information about DNA deformability. This notion strikes at the heart of a fundamental problem in structural biology: to understand the nature of DNA structure in the presence of environmental forces in vivo. A powerful theoretical technique for studying biomolecular structure is molecular simulation. The reliability of molecular simulation to the study of polynucleotide structure, however, has lagged behind that of proteins. This is largely because of the difficulty associated with accommodation of long-range electrostatic forces. In conventional macromolecular simulations, pairwise Coulomb interactions are evaluated using a finite cutoff to decrease computational effort. Ultimately, this truncation of the Coulomb potential leads to unrealistic behaviors2 The adverse effects of the cutoff methods are amplified in highly ionic systems such as DNA. In an attempt to circumvent these problems, artificial constructs are frequently introduced such as reduced phosphate charges on the nucleotide ba~kbone,~ modified Coulomb potentials in the form of distance dependent dielectric function^^,^ or switching function^,^,^ and distance restraints to enforce Watson-Crick hydrogen Perhaps the most studied DNA sequence both experimentally and theoretically is the dodecamer sequence d[CGCGAATTCGCG]2.8,9 Recent attempts to simulate the dodecamer in solution using full charges on the phosphates and explicit counterions to balance the charge led to structures that were

15] Quantitation of water in membranes by neutron diffraction and x-ray techniques

The general principle of placing neutron and X-ray scattering density profiles on an absolute scale is being applied to an increasing number of problems in structural biology. This maximizes the information from the experiments by facilitating the identification of various molecular species. The greater detail available on the membrane water distribution has been highlighted in this chapter. The quantitative analysis of water in the headgroup region and the intermembrane water layer provides valuable information on membrane structure and function. The single most important limitation of the method is the lack of resolution. Improvements in experimental techniques will improve the resolution in a number of situations.

Computer analysis and structure prediction of nucleic acids and proteins.

We have developed an integrated computer system for analysis of nucleic acid and protein sequences, which consists of sequence and structure databases, a relational database, and software for structural analysis. The system is potentially applicable to a number of problems in structural biology including predictive classification of the function and location of oncogene products.