Accuracy and reliability of the scaling-relaxation method for loop closure: An evaluation based on extensive and multiple copy conformational samplings

The accuracy and reliability of the recently proposed scaling-relaxation method for loop closure were examined by using extensive conformational sampling. For each of the eight heptapeptides chosen to represent a variety of protein conformations, 1,000–2,000 conformations were sampled. Each segment contained 14 rotatable backbone dihedral angles. The average root mean square deviations (RMSDs) between the predicted and the native conformations were 0.7 Å for the backbone and 1.2 Å for the side chain atoms. These predictions were substantially more accurate than the previous predictions (1.1 Å for the backbone and 2.2 Å for the side chain atoms) of the same eight protein segments based on limited conformational sampling (100 conformations for each segment). Large prediction errors mostly occurred at polar and surface side chains that are unlikely to have any meaningful conformation. Moreover, the reliability of seven of the eight predictions was demonstrated with their energy-RMSD and stability-RMSD correlations of the low-energy conformations, where the conformational stability was estimated by using the multiple copy simultaneous sampling method.     

A method to search for similar protein local structures at ligand-binding sites and its application to adenine recognition

We have developed a method of searching for similar spatial arrangements of atoms around a given chemical moiety in proteins that bind a common ligand. The first step in this method is to consider a set of atoms that closely surround a given chemical moiety. Then, to compare the spatial arrangements of such surrounding atoms in different proteins, they are translated and rotated so that the chemical moieties are superposed on each other. Spatial arrangements of surrounding atoms in a pair of proteins are judged to be similar, when there are many corresponding atoms occupying similar spatial positions. Because the method focuses on the arrangements of surrounding atoms, it can detect structural similarities of binding sites in proteins that are dissimilar in their amino acid sequences or in their chain folds. We have applied this method to identify modes of nucleotide base recognition by proteins. An all-against-all comparison of the arrangements of atoms surrounding adenine moieties revealed an unexpected structural similarity between protein kinases, cAMP-dependent protein kinase (cAPK), and casein kinase-1 (CK1), and D-Ala:D-Ala ligase (DD-ligase) at their adenine-binding sites, despite a lack of similarity in their chain folds. The similar local structure consists of a four-residue segment and three sequentially separated residues. In particular the four-residue segments of these enzymes were found to have nearly identical conformations in their backbone parts, which are involved in the recognition of adenine. This common local structure was also found in substrate-free three-dimensional structures of other proteins that are similar to DD-ligase in the chain fold and of other protein kinases. As the proteins with different folds were found to share a common local structure, these proteins seem to constitute a remarkable example of convergent evolution for the same recognition mechanism. Received: 9 December 1996 / Accepted: 7 February 1997     

Constructing optimal backbone segments for joining fixed DNA base pairs.

A method is presented to link a sequence of space-fixed base pairs by the sugar-phosphate segments of single nucleotides and to evaluate the effects in the backbone caused by this positioning of the bases. The entire computational unit comprises several nucleotides that are energy-minimized, subject to constraints imposed by the sugar-phosphate backbone segments being anchored to space-fixed base pairs. The minimization schemes are based on two stages, a conjugate gradient method followed by a Newton-Raphson algorithm. Because our purpose is to examine the response, or relaxation, of an artificially stressed backbone, it is essential to be able to obtain, as closely as possible, a lowest minimum energy conformation of the backbone segment in conformational space. For this purpose, an algorithm is developed that leads to the generation of an assembly of many local energy minima. From these sets of local minima, one conformation corresponding to the one with the lowest minimum is then selected and designated to represent the backbone segment at its minimum. The effective electrostatic potential of mean force is expressed in terms of adjustable parameters that incorporate solvent screening action in the Coulombic interactions between charged backbone atoms; these parameters are adjusted to obtain the best fit of the nearest-neighbor phosphorous atoms in an x-ray structure.     

SCWRL and MolIDE: computer programs for side-chain conformation prediction and homology modeling

SCWRL and MolIDE are software applications for prediction of protein structures. SCWRL is designed specifically for the task of prediction of side-chain conformations given a fixed backbone usually obtained from an experimental structure determined by X-ray crystallography or NMR. SCWRL is a command-line program that typically runs in a few seconds. MolIDE provides a graphical interface for basic comparative (homology) modeling using SCWRL and other programs. MolIDE takes an input target sequence and uses PSI-BLAST to identify and align templates for comparative modeling of the target. The sequence alignment to any template can be manually modified within a graphical window of the target-template alignment and visualization of the alignment on the template structure. MolIDE builds the model of the target structure on the basis of the template backbone, predicted side-chain conformations with SCWRL and a loop-modeling program for insertion-deletion regions with user-selected sequence segments. SCWRL and MolIDE can be obtained at (http://dunbrack.fccc.edu/Software.php).     

A 175-psec molecular dynamics simulation of camphor-bound cytochrome P-450cam.

The structure and internal motions of cytochrome P-450cam, a monooxygenase heme enzyme with 414 amino acid residues, with camphor bound at the active site have been evaluated on the basis of a 175-psec molecular dynamics simulation carried out at 300 K. All hydrogen atoms were explicitly modeled, and 204 crystallographic waters were included in the simulation. Based on an analysis of the time course of the trajectory versus potential energy, root mean square deviation, radius of gyration, and hydrogen bonding, the simulation was judged to be stable and representative of the average experimental structure. The averaged structural properties of the enzyme were evaluated from the final 135 psec of the simulation. The average atomic displacement from the X-ray structure was 1.39 A for all heavy atoms and 1.17 A for just C-alpha atoms. The average root-mean-square (rms) fluctuations of all heavy atoms and backbone atoms were 0.42 and 0.37 A, respectively. The computed rms fluctuations were in reasonable agreement with the experimentally determined temperature factors. All 13 segments of alpha-helix and 5 segments of beta-sheet were well preserved with the exception of the N-terminal half of helix F which alternated between an alpha-helix and a 310-helix. In addition there were in general only small variations in the relative orientation of adjacent alpha-helices. The rms fluctuations of the backbone dihedral angles in the secondary structure elements were almost uniformly smaller, with the fluctuation in alpha-helices and beta-sheets, 31 and 10% less, respectively, than those in nonsecondary structure regions. The reported crystal structure contains kinks in both helices C and I. In the simulation, both of these regions showed high mobility and large deviations from their starting positions. Since the kink in the I helix is at the oxygen binding site, these motions may have mechanistic implications.     

Backbone and side-chain assignments of an effector membrane localization domain from Vibrio vulnificus MARTX toxin

¹H, ¹³C, and ¹⁵N chemical shift assignments are presented for the isolated four-helical bundle membrane localization domain from the domain of unknown function 5 (DUF5) effector (MLD_VvDUF5) of the MARTX toxin from Vibrio vulnificus in its solution state. We have assigned 97 % of all backbone and side-chain carbon atoms, including 96 % of all backbone residues. Secondary chemical shift analysis using TALOS+ demonstrates four helices that align with those predicted by structure homology modeling using the MLDs of Pasteurella multocida toxin (PMT) and the clostridial TcdB and TcsL toxins as templates. Future studies will be towards solving the structure and determining the dynamics in the solution state.     

The interaction of cholesterol with bilayers of phosphatidylethanolamine

The interaction of cholesterol with the glycerol backbone segments of phospholipids was studied in bilayers of phosphatidylethanolamine containing equimolar amounts of cholesterol. Glycerol selectively deuterated at various positions was supplied to the growth medium of Escherichia coli strain 131 GP which is defective in endogeneous glycerol synthesis. The procedure enables the stereospecific labeling of the three glycerol backbone segments of the membrane phospholipids. Phosphatidylethanolamine with wild-type fatty acid composition was purified from E. coli cells and deuterium magnetic resonance spectra were obtained either from dispersions of pure phosphatidylethanolamine or from equimolar mixtures of phosphatidylethanolamine with cholesterol. For comparative purposes 1,2-di[9,10-²H₂]elaidoyl-sn-glycero-3-phosphoethanolamine and [3-α-²H]cholesterol were synthesized in order to monitor the behavior of the fatty acyl chains and of the cholesterol molecule itself. For all deuterated segments the deuterium quadrupole splittings as well as the deuterium spin-lattice (T₁) relaxation times were measured as a function of temperature. The glycerol backbone was found to be a remarkably stable structural element of the phospholipid molecule. The quadrupole splittings of the backbone segments changed only by at most 2 kHz upon incorporation of 50 mol % cholesterol. This was in contrast to the fatty acyl chains where the same amount of cholesterol increased the quadrupole splitting by more than 20 kHz. The glycerol segments exhibited the shortest T₁ relaxation times of all CH₂ segments indicating that the glycerol backbone is the slowest motional moiety of the lipid molecule. Addition of cholesterol has no effect on the backbone motion but the fast reorientation rate of the trans-double bonds in 1,2-dielaidoyl-sn-glycero-3-phosphoethanolamine is increased dramatically.     

On the relationship between residue structural environment and sequence conservation in proteins

Residues that are crucial to protein function or structure are usually evolutionarily conserved. To identify the important residues in protein, sequence conservation is estimated, and current methods rely upon the unbiased collection of homologous sequences. Surprisingly, our previous studies have shown that the sequence conservation is closely correlated with the weighted contact number (WCN), a measure of packing density for residue's structural environment, calculated only based on the C_α positions of a protein structure. Moreover, studies have shown that sequence conservation is correlated with environment‐related structural properties calculated based on different protein substructures, such as a protein's all atoms, backbone atoms, side‐chain atoms, or side‐chain centroid. To know whether the C_α atomic positions are adequate to show the relationship between residue environment and sequence conservation or not, here we compared C_α atoms with other substructures in their contributions to the sequence conservation. Our results show that C_α positions are substantially equivalent to the other substructures in calculations of various measures of residue environment. As a result, the overlapping contributions between C_α atoms and the other substructures are high, yielding similar structure–conservation relationship. Take the WCN as an example, the average overlapping contribution to sequence conservation is 87% between C_α and all‐atom substructures. These results indicate that only C_α atoms of a protein structure could reflect sequence conservation at the residue level.     

Predicting the helix packing of globular proteins by self-correcting distance geometry.

A new self-correcting distance geometry method for predicting the three-dimensional structure of small globular proteins was assessed with a test set of 8 helical proteins. With the knowledge of the amino acid sequence and the helical segments, our completely automated method calculated the correct backbone topology of six proteins. The accuracy of the predicted structures ranged from 2.3 A to 3.1 A for the helical segments compared to the experimentally determined structures. For two proteins, the predicted constraints were not restrictive enough to yield a conclusive prediction. The method can be applied to all small globular proteins, provided the secondary structure is known from NMR analysis or can be predicted with high reliability.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

In silico predictions of 3D structures of linear and cyclic peptides with natural and non‐proteinogenic residues

We extended the use of Peplook, an in silico procedure for the prediction of three‐dimensional (3D) models of linear peptides to the prediction of 3D models of cyclic peptides and thanks to the ab initio calculation procedure, to the calculation of peptides with non‐proteinogenic amino acids. Indeed, such peptides cannot be predicted by homology or threading. We compare the calculated models with NMR and X‐ray models and for the cyclic peptides, with models predicted by other in silico procedures (Pep‐Fold and I‐Tasser). For cyclic peptides, on a set of 38 peptides, average root mean square deviation of backbone atoms (BB‐RMSD) was 3.8 and 4.1 Å for Peplook and Pep‐Fold, respectively. The best results are obtained with I‐Tasser (2.5 Å) although evaluations were biased by the fact that the resolved Protein Data Bank models could be used as template by the server. Peplook and Pep‐Fold give similar results, better for short (up to 20 residues) than for longer peptides. For peptides with non‐proteinogenic residues, performances of Peplook are sound with an average BB‐RMSD of 3.6 Å for ‘non‐natural peptides’ and 3.4 Å for peptides combining non‐proteinogenic residues and cyclic structure. These results open interesting possibilities for the design of peptidic drugs. Copyright © 2011 European Peptide Society and John Wiley & Sons, Ltd.     

Three-dimensional solution structure of the E3-binding domain of the dihydrolipoamide succinyltransferase core from the 2-oxoglutarate dehydrogenase multienzyme complex of Escherichia coli.

The three-dimensional solution structure of a 51-residue synthetic peptide comprising the dihydrolipoamide dehydrogenase (E3)-binding domain of the dihydrolipoamide succinyltransferase (E2) core of the 2-oxoglutarate dehydrogenase multienzyme complex of Escherichia coli has been determined by nuclear magnetic resonance spectroscopy and hybrid distance geometry-dynamical simulated annealing calculations. The structure is based on 630 approximate interproton distance and 101 torsion angle (phi, psi, chi 1) restraints. A total of 56 simulated annealing structures were calculated, and the atomic rms distribution about the mean coordinate positions for residues 12-48 of the synthetic peptide is 1.24 A for the backbone atoms, 1.68 A for all atoms, and 1.33 A for all atoms excluding the six side chains which are disordered at chi 1 and the seven which are disordered at chi 2; when the irregular partially disordered loop from residues 31 to 39 is excluded, the rms distribution drops to 0.77 A for the backbone atoms, 1.55 A for all atoms, and 0.89 A for ordered side chains. Although proton resonance assignments for the N-terminal 11 residues and the C-terminal 3 residues were obtained, these two segments of the polypeptide are disordered in solution as evidenced by the absence of nonsequential nuclear Overhauser effects. The solution structure of the E3-binding domain consists of two parallel helices (residues 14-23 and 40-48), a short extended strand (24-26), a five-residue helical-like turn, and an irregular (and more disordered) loop (residues 31-39). This report presents the first structure of an E3-binding domain from a 2-oxo acid dehydrogenase complex.(ABSTRACT TRUNCATED AT 250 WORDS)     

Letter to the Editor: Backbone 1H, 15N and 13C Assignments for the Subunit a of the E. Coli ATP Synthase

The structure of the 30 KDa subunit a of the membrane component (F(0)) of E. coli ATP synthase is investigated in a mixture of chloroform, methanol and water, a solvent previously used for solving the structure of another integral membrane protein, subunit c. Near complete backbone chemical shift assignments were made from a set of TROSY experiments including HNCO, HNCA, HN(CA)CB, HN(CO)CACB and 4D HNCOCA and HNCACO. Secondary structure of subunit a was predicted from the backbone chemical shifts using TALOS program. The protein was found to consist of multiple elongated alpha-helical segments. This finding is generally consistent with previous predictions of multiple transmembrane alpha-helices in this polytopic protein.     

NMR solution structure of plastocyanin from the photosynthetic prokaryote, Prochlorothrix hollandica

The solution structure of a divergent plastocyanin (PC) from the photosynthetic prokaryote Prochlorothrix hollandica was determined by homonuclear 1H NMR spectroscopy. Nineteen structures were calculated from 1222 distance restraints, yielding a family of structures having an average rmsd of 0.42 +/- 0.08 A for backbone atoms and 0.71 +/- 0.07 A for heavy atoms to the mean structure. No distance constraint was violated by more than 0.26 A in the structure family. Despite the low number of conserved residues shared with other PC homologues, the overall folding pattern of P. hollandica PC is similar to other PCs, in that the protein forms a two-sheet beta-barrel tertiary structure. The greatest variability among the backbone structures is seen in the loop region from residues 47-60. The differences seen in the P. hollandica PC homologue likely arise due to a small deletion of 2-4 residues compared to the PC consensus; this yields a less extended loop containing a short alpha-helix from residues Ala52-Leu55. Additionally, the protein has an altered hydrophobic patch thought to be important in binding reaction partners. Whereas the backbone structure is very similar within the loops of the hydrophobic region, the presence of two unique residues (Tyr12 and Pro14) yields a structurally different hydrophobic surface likely important in binding P. hollandica Photosystem I.     

The NMR structure of the activation domain isolated from porcine procarboxypeptidase B.

The three-dimensional structure of the activation domain isolated from porcine pancreatic procarboxypeptidase B was determined using 1H NMR spectroscopy. A group of 20 conformers is used to describe the solution structure of this 81 residue polypeptide chain, which has a well-defined backbone fold from residues 11-76 with an average root mean square distance for the backbone atoms of 1.0 +/- 0.1 A relative to the mean of the 20 conformers. The molecular architecture contains a four-stranded beta-sheet with the polypeptide segments 11-17, 36-39, 50-56 and 75-76, two well defined alpha-helices from residues 20-30 and 60-70, and a 3(10) helix from residues 43-46. The three helices are oriented almost exactly antiparallel to each other, are all on the same side of the beta-sheet, and the helix axes from an angle of approximately 45 degrees relative to the direction of the beta-strands. Three segments linking beta-strands and helical secondary structures, with residues 32-35, 39-43 and 56-61, are significantly less well ordered than the rest of the molecule. In the three-dimensional structure two of these loops (residues 32-35 and 56-61) are located close to each other near the protein surface, forming a continuous region of increased mobility, and the third disordered loop is separated from this region only by the peripheral beta-strand 36-39 and precedes the short 3(10) helix.     

Human epidermal growth factor. High resolution solution structure and comparison with human transforming growth factor alpha.

The solution structure of the 53 amino acid peptide hormone, human epidermal growth factor (hEGF), has been determined to high resolution from nuclear magnetic resonance (n.m.r.) data. A large number of internuclear distance and dihedral restraints was obtained, including data from uniformly 15N-labelled hEGF. Dynamical simulated annealing methods using the program XPLOR were used for structure calculation. An improved protocol was developed combining efficient conformational searching at a reduced computational cost. The general fold of the calculated structures compared well with that of a derivative of the carboxy-terminally truncated hEGF determined previously. A group of 44 structures were calculated with no violations greater than 0.3 A and 3 degrees for distance and dihedral restraints, respectively. The average pairwise root mean square (r.m.s.) deviation of all backbone atoms for these structures was 2.25 A for all 53 residues, 0.92 A for the bulk of the protein, and 0.23 A for the functionally important carboxy-terminal domain. Two new helical segments containing highly conserved amino acids have been identified; one between cysteines 6 and 14 and a second at the end of the carboxy-terminal domain. New insight into the molecular architecture of the site of putative receptor binding was provided by comparing the structure of hEGF with its biologically equipotent analogue, human transforming growth factor alpha. This comparison revealed a close structural relationship between the two growth factors and provides an improved understanding of the structure/function relationships in EGF.     

Packing of transmembrane helices in bacteriorhodopsin folding: structure and thermodynamics

We propose a coarse-grained (CG) model to study the native structure and physical properties of helical membrane proteins (HMPs) using off-lattice computer simulations. Instead of considering sequence heterogeneity explicitly, we model its effect on the packing of helices by employing a mean packing parameter r(0), which is calculated from an all-atom (AA) model. Specifically, this CG model is applied to investigate the packing of helices in bacteriorhodopsin (BR), and predicts the seven helix bundle structure of BR with a root mean square deviation (RMSD) in coordinates of helix backbone atoms (N, C, C(alpha)) of 3.99 A from its crystal structure. This predicted structure is further refined in an AA model by Amber and the refined structure has a RMSD (in coordinates of helix backbone atoms) of 2.64 A. The predicted packing position, tilting angle, and orientation angle of each helix in the refined structure are consistent with experimental data and their physical origins can be well understood in our model. Our results show that a reasonably good structure of BR can be predicted by using such a dual-scale approach, provided that its secondary structure is known. Starting from a random initial configuration, the folded structure can be obtained in days using a regular desktop computer. Various thermodynamic properties of helix packing of BR are also investigated in this CG model.     

Extracting information on folding from the amino acid sequence: accurate predictions for protein regions with preferred conformation in the absence of tertiary interactions.

A recently developed procedure to predict backbone structure from the amino acid sequence [Rooman, M., Kocher, J. P., & Wodak, S. (1991) J. Mol. Biol, 221, 961-979] is fine tuned to identify protein segments, of length 5-15 residues, that adopt well-defined conformations in the absence of tertiary interactions. These segments are obtained by requiring that their predicted lowest energy structures have a sizable energy gap relative to other computed conformations. Applying this procedure to 69 proteins of known structure, we find that regions with largest energy gaps--those having highly preferred conformations--are also the most accurately predicted ones. On the basis of previous findings that such regions correlate well with sites that become structured early during folding, our approach provides the means of identifying such sites in proteins without prior knowledge of the tertiary structure. Furthermore, when predictions are performed so as to ignore the influence of residues flanking each segment along the sequence, a situation akin to excising the considered peptide from the rest of the chain, they offer the possibility of identifying protein segments liable to adopt well-defined conformations on their own. The described approach should have useful applications in experimental and theoretical investigations of protein folding and stability, and aid in designing peptide drugs and vaccines.