The accuracy of protein structures in solution determined by AlphaFold and NMR

1. Download : Download high-res image (139KB)
2. Download : Download full-size image

     

Systematic solution to homo‐oligomeric structures determined by NMR

Protein structure determination by NMR has predominantly relied on simulated annealing‐based conformational search for a converged fold using primarily distance constraints, including constraints derived from nuclear Overhauser effects, paramagnetic relaxation enhancement, and cysteine crosslinkings. Although there is no guarantee that the converged fold represents the global minimum of the conformational space, it is generally accepted that good convergence is synonymous to the global minimum. Here, we show such a criterion breaks down in the presence of large numbers of ambiguous constraints from NMR experiments on homo‐oligomeric protein complexes. A systematic evaluation of the conformational solutions that satisfy the NMR constraints of a trimeric membrane protein, DAGK, reveals 9 distinct folds, including the reported NMR and crystal structures. This result highlights the fundamental limitation of global fold determination for homo‐oligomeric proteins using ambiguous distance constraints and provides a systematic solution for exhaustive enumeration of all satisfying solutions. Proteins 2015; 83:651–661. © 2015 Wiley Periodicals, Inc.     

Structure of epidermal growth factor bound to perdeuterated dodecylphosphocholine micelles determined by two-dimensional NMR and simulated annealing calculations.

The interaction of mouse epidermal growth factor (mEGF) with micelles of a phospholipid analogue, perdeuterated dodecylphosphocholine (DPC), was investigated by two-dimensional 1H NMR. Sequence-specific resonance assignments of the micelle-bound mEGF have been made, and the chemical shifts were compared with those in the absence of DPC. DPC induced large chemical shift changes of the resonances from the residues in the C-terminal tail (residues 46-53) but little perturbation on the residues in the main core (residues 1-45). Starting from the three-dimensional structure in the absence of DPC, micelle-bound structures were calculated using the program XPLOR with interproton distance data obtained from NOESY spectra recorded in the presence of DPC. The C-terminal tail of mEGF was found to change conformation to form an amphiphilic structure when bound to the micelles. It is possible that induced fit in the C-terminal tail of mEGF occurs upon binding to a putative hydrophobic pocket of the EGF receptor.     

The disordered conformation of kappa-carrageenan in solution as determined by NMR experiments and molecular modeling

The conformation of kappa-carrageenan in solution was studied combining 1H and 13C NMR with molecular mechanics. The experimental conditions were chosen to characterize the disordered conformation of the polymer. Particular attention has been given to explore a wide range of experimental conditions as to the dependence on solvent (water and Me2SO), polymer concentration, temperature, pH, presence of a denaturing agent (guanidinium chloride), and of ions otherwise able to induce conformational order of the carrageenan chains, either in solution (I-) or in the gel state (Rb+). Two-dimensional NOE experiments were analyzed to obtain information on internuclear distances, and molecular mechanics provided the range of energetically accessible conformations. Two inter-residue topological constraints were clearly identified: their combination is rather restricting for the chain and suggests that the disordered conformation of kappa-carrageenan is characterized by an intrinsic stiffness with high values of persistent length and characteristic ratio. They also rule out any postulated interchain hydrogen bonds. In contrast, experiments on the temperature dependence of the chemical shift in Me2SO reveal the existence of two inter-residue intramolecular H-bonds which might contribute positively to the rigidity of the polymer chain. The overall picture emerging from the present results is that of a locally elongated 'loose single helix'.     

Optimal molecular structures of prion AGAAAAGA amyloid fibrils formatted by simulated annealing

To date, there is little structural data available on the AGAAAAGA palindrome in the hydrophobic region (113–120) of prion proteins, although many experimental studies have shown that this region has amyloid fibril forming properties. This region belongs to the N-terminal unstructured region (1–123) of prions, the structure of which has proved hard to determine using NMR or X-ray crystallography. This paper reports the successful construction of three amyloid fibril models for this region. The models were formatted by standard simulated annealing using suitable templates from the Protein Data Bank, and were refined using several traditional optimization methods within AMBER. Because the NMR or X-ray structure of the hydrophobic region AGAAAAGA of prion proteins has not yet been determined, these models can be used as a reference for experimental studies on this region. The results presented here confirm standard simulated annealing as an effective tool in molecular modeling. The three constructed models for amyloid fibrils may be useful in furthering the goals of medicinal chemistry in this field.     

NMR and simulated annealing investigations of bradykinin in presence of polyphenols

Epidemiological studies have shown that the incidence of some cardiovascular degenerative diseases appears to be lower in populations with regular but moderate drinking of red wine rich in polyphenols. One of the most important properties of polyphenols is to form complexes with proteins. The linear nonapeptide hormone bradykinin (H-Arg1-Pro2-Pro3-Gly4-Phe5-Ser6-Pro7-Phe8-Arg9-OH) is involved in a variety of physiological processes such as the cardiovascular processes. Thus, the goal of this work was to study the effects of tannins on the peptide structure by NMR investigations and molecular modeling. The results of these investigations show that in the presence of catechin, the peptide conformation is not affected and is in a random coil structure. On the contrary, the peptide structure is modified by the addition of dimeric proanthocyanidin B3 (catechin 4alpha-->8 catechin). The dimer leads to the formation of a large flexible turn between the 6-9 residues. Thus, the biological activities of bradykinin in the presence of polyphenols could be affected.     

Determination of solution structures of paramagnetic proteins by NMR

Standard procedures for using nuclear Overhauser enhancements (NOE) between protons to generate structures for diamagnetic proteins in solution from NMR data may be supplemented by using dipolar shifts if the protein is paramagnetic. This is advantageous since the electron-nuclear dipolar coupling provides relatively long-range geometric information with respect to the paramagnetic centre which complements the short-range distance constraints from NOEs. Several different strategies have been developed to date, but none of these attempts to combine data from NOEs and dipolar shifts in the initial stages of structure calculation or to determine three dimensional protein structures together with their magnetic properties. This work shows that the magnetic and atomic structures are highly correlated and that it is important to have additional constraints both to provide starting parameters for the magnetic properties and to improve the definition of the best fit. Useful parameters can be obtained for haem proteins from Fermi contact shifts; this approach is compared with a new method based on the analysis of dipolar shifts in haem methyl groups with respect to data from horse and tuna ferricytochromes c. The methods developed for using data from NOEs and dipolar shifts have been incorporated in a new computer program, PARADYANA, which is demonstrated in application to a model data set for the sequence of the haem octapeptide known as microperoxidase-8. Received: 13 October 1997 / Accepted: 19 December 1997     

Solution structures of alpha-conotoxin G1 determined by two-dimensional NMR spectroscopy

Two-dimensional NMR data have been used to generate solution structures of alpha-conotoxin G1, a potent peptide antagonist of the acetylcholine receptor. Structural information was obtained in the form of proton-proton internuclear distance constraints, and initial structures were produced with a distance geometry algorithm. Energetically more favorable structures were generated by using the distance geometry structures as input for a constrained energy minimization program. The results of both of these calculations indicate that the overall backbone conformation of the molecule is well-defined by the NMR data whereas the side-chain conformations are generally less well-defined. The main structural features derived from the NMR data were the presence of tight turns centered on residues Pro5 and Arg9. The solution structures are compared with previous proposed models of conotoxin G1, and the NMR data are interpreted in conjunction with chemical modification studies and structural properties of other antagonists of the acetylcholine receptor to gain insight into structure-activity relationships in these peptide toxins.     

Three-dimensional solution structure of the calcium-signaling protein apo-S100A1 as determined by NMR.

S100A1, a member of the S100 protein family, is an EF-hand containing Ca(2+)-binding protein (93 residues per subunit) with noncovalent interactions at its dimer interface. Each subunit of S100A1 has four alpha-helices and a small antiparallel beta-sheet consistent with two helix-loop-helix calcium-binding domains [Baldiserri et al. (1999) J. Biomol. NMR 14, 87-88]. In this study, the three-dimensional structure of reduced apo-S100A1 was determined by NMR spectroscopy using a total of 2220 NOE distance constraints, 258 dihedral angle constraints, and 168 backbone hydrogen bond constraints derived from a series of 2D, 3D, and 4D NMR experiments. The final structure was found to be globular and compact with the four helices in each subunit aligning to form a unicornate-type four-helix bundle. Intermolecular NOE correlations were observed between residues in helices 1 and 4 from one subunit to residues in helices 1' and 4' of the other subunit, respectively, consistent with the antiparallel alignment of the two subunits to form a symmetric X-type four-helix bundle as found for other members of the S100 protein family. Because of the similarity of the S100A1 dimer interface to that found for S100B, it was possible to calculate a model of the S100A1/B heterodimer. This model is consistent with a number of NMR chemical shift changes observed when S100A1 is titrated into a sample of (15)N-labeled S100B. Helix 3 (and 3') of S100A1 was found to have an interhelical angle of -150 degrees with helix 4 (and 4') in the apo state. This crossing angle is quite different (>50 degrees ) from that typically found in other EF-hand containing proteins such as apocalmodulin and apotroponin C but more similar to apo-S100B, which has an interhelical angle of -166 degrees. As with S100B, it is likely that the second EF-hand of apo-S100A1 reorients dramatically upon the addition of Ca(2+), which can explain the Ca(2+) dependence that S100A1 has for binding several of its biological targets.     

Three-dimensional structure of acyl carrier protein in solution determined by nuclear magnetic resonance and the combined use of dynamical simulated annealing and distance geometry

The solution conformation of acyl carrier protein from Escherichia coli (77 residues) has been determined on the basis of 423 interproton-distance restraints and 32 hydrogen-bonding restraints derived from NMR measurements. A total of nine structures were computed using a hybrid approach combining metric matrix distance geometry and dynamic simulated annealing. The polypeptide fold is well defined with an average backbone atomic root-mean-square difference of 0.20 +/- 0.03 nm between the final nine converged structures and the mean structure obtained by averaging their coordinates. The principal structural motif is composed of three helices: 1 (residues 3-12), 2 (residues 37-47) and 4 (residues 65-75) which line a hydrophobic cavity. Helices 2 and 4 are approximately parallel to each other and anti-parallel at an angle of approximately equal to 150 degrees to helix 1. The smaller helix 3 (residues 56-63) is at an angle of approximately equal to 100 degrees to helix 4.     

Comparison of experimentally determined protein structures by solution of Bloch equations

Nuclear Overhauser enhancement (NOESY) spectra were theoretically generated by solving the generalized Bloch equations with the appropriate initial conditions. The input to the equations were the coordinates of the protons of two similar crystal structures of basic pancreatic trypsin inhibitor. The two NOESY spectra obtained were compared to published experimental spectra of the protein in solution. It was found that the two crystal structures of basic pancreatic trypsin inhibitor give different theoretical spectra. The solution of the Bloch equations is very sensitive to small variations in the distance between protons (approx. 0.2 A), and to differences in the surrounding configurations. The method allows a detailed comparison of the crystal and solution structures of proteins. The structure of the trypsin inhibitor in solution was found to be similar to either one or the other crystal forms in different regions of the molecule.     

The PP-fold solution structure of human polypeptide YY and human PYY3-36 as determined by NMR

PYY3-36 is a biopharmaceutical antiobesity agent under development as well as an endogenous satiety hormone, which is generated by dipeptidyl peptidase-IV digestion of polypetide YY (PYY), and in contrast to the parent hormone, PYY is highly selective for the Y2 versus the Y1 receptor. NMR analysis revealed a highly ordered, back-folded structure for human PYY in aqueous solution similar to the classical PP-fold structure of pancreatic polypeptide. The NMR analysis of PYY3-36 also showed a folded structure resembling a PP-fold, which however was characterized by far fewer long distance NOEs than the PP-fold observed in the full-length peptide. This suggests that either a conformational change has occurred in the N-terminal segment of PYY3-36 or that this segments is characterized by larger dynamics. The study supports the notion that the PP-fold is crucial for establishing simultaneous interactions with two subsites in the receptor for binding of, respectively, the N- and C-terminal ends of PYY. The Y2 receptor only requires recognition of the C-terminal segment of the molecule as displayed by the Y2 selective PYY3-36.     

The solution structures of calyculin A and dephosphonocalyculin A by NMR

The NMR solution structure of calyculin A (1) in chloroform exhibits intramolecular interactions, resembling the original crystal structure. In methanol, calyculin A has the hydrogen bonding moieties solvent exposed. Dephosphonocalyculin A in chloroform resembles calyculin A in chloroform and the crystal structure of calyculin A. Dephosphonocalyculin A in methanol resembles calyculin A in methanol.     

Differences in the solution structures of the parallel beta-helical pectate lyases as determined by limited proteolysis

The pectate lyase family of proteins has been shown to fold into a novel domain motif, the right-handed parallel beta-helix. As a means of gaining insight to the solution structure of the pectate lyases, the enzymes were subjected to limited proteolytic digestion by the endoproteases AspN, GluC and trypsin. The effects of proteolytic cleavage on enzymatic activity were determined, and the early products of proteolysis were identified by capillary electrophoresis, MALDI-TOF mass spectrometry and HPLC. A single peptide bond between Lys158 and Asp159 in pectate lyase B (PLb) was cleaved by both AspN and trypsin, with no detectable hydrolysis of PLb by GluC. Pectate lyase E (PLe) was hydrolyzed by trypsin between Lys164 and Asp165, a bond on an analogous loop structure found to be susceptible to proteolytic attack in PLb. AspN and GluC preferentially hydrolyzed peptide bonds (at Asp127 and Glu124, respectively) on another loop extending from the central beta-helical core of PLe. A single beta-strand of the central cylinder of the pectate lyase C (PLc) molecule was susceptible to all three proteases used. These data demonstrate that the most susceptible peptide bonds to proteolytic scission within the native enzymes lie on or near one of the three parallel beta-sheets that compose the core domain motif Despite the proximity of the proteolytic cleavages to the catalytic sites of the enzymes, significant retention of lyase activity was observed after partial proteolysis, indicating preservation of functional tertiary structure in the proteolytic products.     

A prediction of tertiary structures of peptide by the Monte Carlo simulated annealing method

The Monte Carlo simulated annealing method has been applied to the prediction of three-dimensional structures of enkephalin. The low-energy conformations obtained were classified into a few groups of similar structures, which indicates that our method is effective. New low-energy structures were identified together with previously proposed structures.     

The solvent protection of alzheimer amyloid-beta-(1-42) fibrils as determined by solution NMR spectroscopy

Alzheimer disease is a neurodegenerative disorder that is tightly linked to the self-assembly and amyloid formation of the 39-43-residue-long amyloid-beta (Abeta) peptide. Considerable evidence suggests a correlation between Alzheimer disease development and the longer variants of the peptide, Abeta-(1-42/43). Currently, a molecular understanding for this behavior is lacking. In the present study, we have investigated the hydrogen/deuterium exchange of Abeta-(1-42) fibrils under physiological conditions, using solution NMR spectroscopy. The obtained residue-specific and quantitative map of the solvent protection within the Abeta-(1-42) fibril shows that there are two protected core regions, Glu11-Gly25 and Lys28-Ala42, and that the residues in between, Ser26 and Asn27, as well as those in the N terminus, Asp1-Tyr10, are solvent-accessible. This result reveals considerable discrepancies when compared with a previous investigation on Abeta-(1-40) fibrils and suggests that the additional residues in Abeta-(1-42), Ile41 and Ala42, significantly increase the solvent protection and stability of the C-terminal region Lys28-Ala42. Consequently, our findings provide a molecular explanation for the increased amyloidogenicity and toxicity of Abeta-(1-42) compared with shorter Abeta variants found in vivo.     

Three-dimensional solution structure of apo-neocarzinostatin from Streptomyces carzinostaticus determined by NMR spectroscopy.

The three-dimensional solution structure of apo-neocarzinostatin has been resolved from nuclear magnetic resonance spectroscopy data. Up to 1034 constraints were used to generate an initial set of 45 structures using a distance geometry algorithm (DSPACE). From this set, ten structures were subjected to refinement by restrained energy minimization and molecular dynamics. The average atomic root mean square deviations between the final ten structures and the mean structure obtained by averaging their coordinates run from 0.085 nm for the best defined beta-sheet regions of the protein to 0.227 nm for the side chains of the most flexible loops. The solution structure of apo-neocarzinostatin is closely similar to that of the related proteins, macromomycin and actinoxanthin. It contains a seven-stranded antiparallel beta-barrel which forms, together with two external loops, a deep cavity that is the chromophore binding site. It is noteworthy that aromatic side chains extend into the binding cleft. They may be responsible for the stabilization of the holo-protein complex and for the chromophore specificity within the antitumoral family.     

Solution structure of a sweet protein single-chain monellin determined by nuclear magnetic resonance and dynamical simulated annealing calculations

Single-chain monellin (SCM), which is an engineered 94-residue polypeptide, has proven to be as sweet as native two-chain monellin. SCM is more stable than the native monellin for both heat and acidic environments. Data from gel filtration HPLC and NMR indicate that the SCM exists as a monomer in aqueous solution. The solution structure of SCM has been determined by nuclear magnetic resonance (NMR) spectroscopy and dynamical simulated annealing calculations. A stable alpha-helix spanning residues Phe11-Ile26 and an antiparallel beta-sheet formed by residues 2-5, 36-38, 41-47, 54-64, 69-75, and 83-88 have been identified. The sheet was well defined by backbone-backbone NOEs, and the corresponding beta-strands were further confirmed by hydrogen bond networks based on amide hydrogen exchange data. Strands beta2 and beta3 are connected by a small bulge comprising residues Ile38-Cys41. A total of 993 distance and 56 dihedral angle restraints were used for simulated annealing calculations. The final simulated annealing structures (k) converged well with a root-mean-square deviation (rmsd) between backbone atoms of 0.49 A for secondary structural regions and 0.70 A for backbone atoms excluding two loop regions. The average restraint energy-minimized (REM) structure exhibited root-mean-square deviations of 1.19 A for backbone atoms and 0.85 A for backbone atoms excluding two loop regions with respect to 20 k structures. The solution structure of SCM revealed that the long alpha-helix was folded into the concave side of a six-stranded antiparallel beta-sheet. The side chains of Tyr63 and Asp66 which are common to all sweet peptides showed an opposite orientation relative to H1 helix, and they were all solvent-exposed. Residues at the proposed dimeric interface in the X-ray structure were observed to be mostly solvent-exposed and demonstrated high degrees of flexibility.