Side-chain Torsion Angles Research Articles

Fine grained sampling of residue characteristics using molecular dynamics simulation

In a fine-grained computational analysis of protein structure, we investigated the relationships between a residue's backbone conformations and its side-chain packing as well as conformations. To produce continuous distributions in high resolution, we ran molecular dynamics simulations over a set of protein folds (dynameome). In effect, the dynameome dataset samples not only the states well represented in the PDB but also the known states that are not well represented in the structural database. In our analysis, we characterized the mutual influence among the backbone ϕ, ψ angles with the first side-chain torsion angles ( χ 1) and the volumes occupied by the side-chains. The dependencies of these relationships on side-chain environment and amino acids are further explored. We found that residue volumes exhibit dependency on backbone 2° structure conformation: side-chains pack more densely in extended β-sheet than in α-helical structures. As expected, residue volumes on the protein surface were larger than those in the interior. The first side-chain torsion angles are found to be dependent on the backbone conformations in agreement with previous studies, but the dynameome dataset provides higher resolution of rotamer preferences based on the backbone conformation. All three gauche −, gauche +, and trans rotamers show different patterns of ϕ, ψ dependency, and variations in χ 1 value are skewed from their canonical values to relieve the steric strains. By demonstrating the utility of dynameomic modeling on the native state ensemble, this study reveals details of the interplay among backbone conformations, residue volumes and side-chain conformations.

Computational Design of a PAK1 Binding Protein

We describe a computational protocol, called DDMI, for redesigning scaffold proteins to bind to a specified region on a target protein. The DDMI protocol is implemented within the Rosetta molecular modeling program and uses rigid-body docking, sequence design, and gradient-based minimization of backbone and side-chain torsion angles to design low-energy interfaces between the scaffold and target protein. Iterative rounds of sequence design and conformational optimization were needed to produce models that have calculated binding energies that are similar to binding energies calculated for native complexes. We also show that additional conformation sampling with molecular dynamics can be iterated with sequence design to further lower the computed energy of the designed complexes. To experimentally test the DDMI protocol, we redesigned the human hyperplastic discs protein to bind to the kinase domain of p21-activated kinase 1 (PAK1). Six designs were experimentally characterized. Two of the designs aggregated and were not characterized further. Of the remaining four designs, three bound to the PAK1 with affinities tighter than 350 μM. The tightest binding design, named Spider Roll, bound with an affinity of 100 μM. NMR-based structure prediction of Spider Roll based on backbone and 13C β chemical shifts using the program CS-ROSETTA indicated that the architecture of human hyperplastic discs protein is preserved. Mutagenesis studies confirmed that Spider Roll binds the target patch on PAK1. Additionally, Spider Roll binds to full-length PAK1 in its activated state but does not bind PAK1 when it forms an auto-inhibited conformation that blocks the Spider Roll target site. Subsequent NMR characterization of the binding of Spider Roll to PAK1 revealed a comparably small binding ‘on-rate’ constant (≪ 10 5 M − 1 s − 1 ). The ability to rationally design the site of novel protein–protein interactions is an important step towards creating new proteins that are useful as therapeutics or molecular probes.

Conformational analysis and vibrational spectroscopic investigation of L-alanyl-L-glutamine dipeptide

In this study conformational behavior of anticancer chemotherapy dipeptide Ala-Gln and its dimers have been investigated by molecular mechanic and ab-initio calculations. The calculations on Ala-Gln dipeptide as a function of side chain torsion angles, enable us to determine their energetically preferred conformations. The relative positions of the side chain residues of the stable conformations of dipeptide were obtained, depending on the obtained conformational analysis results. The lowest energy conformation of the dipeptide has been determined by using the Ramachandran maps (Biopolymers6(1963), 1494;J. Mol. Biol.7(1963), 95) and compared with the quantum chemical ab-initio results. The geometry optimization, vibrational wavenumbers and intensity calculations of Ala-Gln dipeptide were carried out with theGaussian03program by using DFT with B3LYP functional and 6-31++G(d,p) basis set. The IR (4000–400 cm−1) and Raman spectra of the Ala-Gln dipeptide have been reported in solid phase, and compared with the theoretical vibrational data.

Conformational analysis and vibrational spectroscopic investigation of L-alanyl-L-glutamine dipeptide

In this study conformational behavior of anticancer chemotherapy dipeptide Ala-Gln and its dimers have been investigated by molecular mechanic and ab-initio calculations. The calculations on Ala-Gln dipeptide as a function of side chain torsion angles, enable us to determine their energetically preferred conformations. The relative positions of the side chain residues of the stable conformations of dipeptide were obtained, depending on the obtained conformational analysis results. The lowest energy conformation of the dipeptide has been determined by using the Ramachandran maps (Biopolymers 6 (1963), 1494; J. Mol. Biol.7 (1963), 95) and compared with the quantum chemical ab-initio results. The geometry optimization, vibrational wavenumbers and intensity calculations of Ala-Gln dipeptide were carried out with the Gaussian03 program by using DFT with B3LYP functional and 6-31

A Webserver for Generating Stereochemically-Acceptable Protein Pathways and Movies

We introduce a new, quick method for generating stereochemically-acceptable pathways in proteins, The method, called geometric targeting, is an alternative to the computationally-intensive targeted molecular dynamics approach. Geometric targeting takes as input two distinct protein conformations and produces an all atom pathway between the two states, guided by geometric considerations that will be described. We also present our new webserver for protein pathways. The user submits two protein structures to the webserver, and the geometric targeting method is run automatically to generate a pathway. The webserver also includes tools for visualization of the pathway and downloading of pathway movie files for use in presentations. The strategy behind the geometric targeting method is to take random steps while gradually decreasing the RMSD to the target, and while imposing various geometric constraints to make sure that each snapshot has good stereochemistry. The pathways maintain good covalent bond distances and angles, keep backbone dihedral angles in allowed Ramachandran regions, avoid eclipsed side-chain torsion angles, avoid non-bonded overlap, and maintain a set of hydrogen bonds and hydrophobic contacts. The method does not necessarily produce the optimal pathway, but rather a stereochemically-acceptable pathway. By running multiple times, a collection of random pathways between the two states can be generated. These pathways will be useful for further quantitative analysis, such as to study free energy changes or search for transition states.

Zero-quantum frequency-selective recoupling of homonuclear dipole-dipole interactions in solid state nuclear magnetic resonance

We describe a method for measuring magnetic dipole-dipole interactions, and hence distances, between pairs of like nuclear spins in a many-spin system under magic-angle spinning (MAS). This method employs a homonuclear dipolar recoupling sequence that creates an average dipole-dipole coupling Hamiltonian under MAS with full zero-quantum symmetry, including both secular and flip-flop terms. Flip-flop terms are then attenuated by inserting rotor-synchronized periods of chemical shift evolution between recoupling blocks, leaving an effective Hamiltonian that contains only secular terms to a good approximation. Couplings between specific pairs of nuclear spins can then be selected with frequency-selective pi pulses. We demonstrate this technique, which we call zero-quantum shift evolution assisted homonuclear recoupling, in a series of one-dimensional and two-dimensional (13)C NMR experiments at 17.6 T and 40.00 kHz MAS frequency on uniformly (13)C-labeled L-threonine powder and on the helix-forming peptide MB(i+4)EK, synthesized with a pair of uniformly (13)C-labeled L-alanine residues. Experimental demonstrations include measurements of distances between (13)C sites that are separated by three bonds, placing quantitative constraints on both sidechain and backbone torsion angles in polypeptides.

Assessing the accuracy of protein structures by quantum mechanical computations of 13C(alpha) chemical shifts.

Two major techniques have been used to determine the three-dimensional structures of proteins: X-ray diffraction and NMR spectroscopy. In particular, the validation of NMR-derived protein structures is one of the most challenging problems in NMR spectroscopy. Therefore, researchers have proposed a plethora of methods to determine the accuracy and reliability of protein structures. Despite these proposals, there is a growing need for more sophisticated, physics-based structure validation methods. This approach will enable us to (a) characterize the "quality" of the NMR-derived ensemble as a whole by a single parameter, (b) unambiguously identify flaws in the sequence at a residue level, and (c) provide precise information, such as sets of backbone and side-chain torsional angles, that we can use to detect local flaws. Rather than reviewing all of the existing validation methods, this Account describes the contributions of our research group toward a solution of the long-standing problem of both global and local structure validation of NMR-derived protein structures. We emphasize a recently introduced physics-based methodology that makes use of observed and computed (13)C(alpha) chemical shifts (at the density functional theory (DFT) level of theory) for an accurate validation of protein structures in solution and in crystals. By assessing the ability of computed (13)C(alpha) chemical shifts to reproduce observed (13)C(alpha) chemical shifts of a single structure or ensemble of structures in solution and in crystals, we accomplish a global validation by using the conformationally averaged root-mean-square deviation, ca-rmsd, as a scoring function. In addition, the method enables us to provide local validation by identifying a set of individual amino acid conformations for which the computed and observed (13)C(alpha) chemical shifts do not agree within a certain error range and may represent a nonreliable fold of the protein model. Although it is computationally intensive, our validation method has several advantages, which we illustrate through a series of applications. This method makes use of the (13)C(alpha) chemical shifts, not shielding, that are ubiquitous to proteins and can be computed precisely from the phi, psi, and chi torsional angles. There is no need for a priori knowledge of the oligomeric state of the protein, and no knowledge-based information or additional NMR data are required. The primary limitation at this point is the computational cost of such calculations. However, we anticipate that enhancements in the speed of calculating these chemical shifts coupled with the ever-increasing computational power should soon make this a standard method accessible to the general NMR community.

Prediction of Protein Loop Structures Using a Local Move Monte Carlo Approach and a Grid-Based Force Field

We have developed an improved local move Monte Carlo (LMMC) loop sampling approach for loop predictions. The method generates loop conformations based on simple moves of the torsion angles of side chains and local moves of backbone of loops. To reduce the computational costs for energy evaluations, we developed a grid-based force field to represent the protein environment and solvation effect. Simulated annealing has been used to enhance the efficiency of the LMMC loop sampling and identify low-energy loop conformations. The prediction quality is evaluated on a set of protein loops with known crystal structure that has been previously used by others to test different loop prediction methods. The results show that this approach can reproduce the experimental results with

Structure and conformation of linear peptides. I. Structure of L-tyrosyl-L-tyrosine.

L-tyrosyl-L-tyrosine crystallizes as a dihydrate in the orthorhombic system, space group C222(1), with a = 12.105(2), b = 12.789(2), c = 24.492(3) A, Z = 8. The structure was solved by direct methods and refined to a final R-value of 0.059 for 1740 observed reflections. The molecule exists as a zwitterion, the peptide unit is trans planar, and the backbone torsion angles correspond to an extended conformation, with psi 1 = 149.4 degrees, phi 2 = -161.2 degrees, psi 2 = 158.3 degrees. The values of the side-chain torsion angles (chi 1, chi 2) are (-58.8 degrees, -63.1 degrees) for the first tyrosine and (-171.7 degrees, -116.5 degrees) for the second. The planes of the aromatic rings are nearly parallel (dihedral angle of 6.1 degrees), and their centers are separated by 10.9 A. The carboxyl plane forms a dihedral angle of 23.8 degrees with the plane of the peptide bond.

Prediction of protein loop structures using a local move Monte Carlo approach and a grid-based force field.

We have developed an improved local move Monte Carlo (LMMC) loop sampling approach for loop predictions. The method generates loop conformations based on simple moves of the torsion angles of side chains and local moves of backbone of loops. To reduce the computational costs for energy evaluations, we developed a grid-based force field to represent the protein environment and solvation effect. Simulated annealing has been used to enhance the efficiency of the LMMC loop sampling and identify low-energy loop conformations. The prediction quality is evaluated on a set of protein loops with known crystal structure that has been previously used by others to test different loop prediction methods. The results show that this approach can reproduce the experimental results with the root mean square deviation within 1.8 A for all the test cases. The LMMC loop prediction approach developed here could be useful for improvement in the quality the loop regions in homology models, flexible protein-ligand and protein-protein docking studies.

The NMR Structures of the Major Intermediates of the Two-domain Tick Carboxypeptidase Inhibitor Reveal Symmetry in Its Folding and Unfolding Pathways

There is a lack of experimental structural information about folding intermediates of multidomain proteins. Tick carboxypeptidase inhibitor (TCI) is a small, disulfide-rich protein consisting of two domains that fold and unfold autonomously through the formation of two major intermediates, IIIa and IIIb. Each intermediate contains three native disulfide bonds in one domain and six free cysteines in the other domain. Here we have determined the NMR structures of these two intermediates trapped and isolated at acidic pH in which they are stable and compared their structures with that of the native protein analyzed under the same conditions. Both IIIa and IIIb were found to contain a folded region that corresponds to the N- and C-terminal domains of TCI, respectively, with structures very similar to the corresponding regions of the native protein. The remainder of the polypeptide chains of the intermediates was shown to be unfolded in a random coil conformation. Solvent exchange measurements further indicated that the two protein domains are not completely independent, but affect each other in terms of dynamics and stability, in agreement with reported inhibitory activity data. The derived results provide structural evidence for symmetric TCI folding and unfolding mechanisms that converge in IIIa and IIIb and reveal the structural basis that accounts for the strong and simultaneous accumulation of both intermediates. Altogether, this work has important implications for a better understanding of the folding mechanisms of multidomain, disulfide-rich proteins.

Factors affecting the use of13Cα chemical shifts to determine, refine, and validate protein structures

Interest centers here on the analysis of two different, but related, phenomena that affect side-chain conformations and consequently 13C(alpha) chemical shifts and their applications to determine, refine, and validate protein structures. The first is whether 13C(alpha) chemical shifts, computed at the DFT level of approximation with charged residues is a better approximation of observed 13C(alpha) chemical shifts than those computed with neutral residues for proteins in solution. Accurate computation of 13C(alpha) chemical shifts requires a proper representation of the charges, which might not take on integral values. For this analysis, the charges for 139 conformations of the protein ubiquitin were determined by explicit consideration of protein binding equilibria, at a given pH, that is, by exploring the 2(xi) possible ionization states of the whole molecule, with xi being the number of ionizable groups. The results of this analysis, as revealed by the shielding/deshielding of the 13C(alpha) nucleus, indicated that: (i) there is a significant difference in the computed 13C(alpha) chemical shifts, between basic and acidic groups, as a function of the degree of charge of the side chain; (ii) this difference is attributed to the distance between the ionizable groups and the 13C(alpha) nucleus, which is shorter for the acidic Asp and Glu groups as compared with that for the basic Lys and Arg groups; and (iii) the use of neutral, rather than charged, basic and acidic groups is a better approximation of the observed 13C(alpha) chemical shifts of a protein in solution. The second is how side-chain flexibility influences computed 13C(alpha) chemical shifts in an additional set of ubiquitin conformations, in which the side chains are generated from an NMR-derived structure with the backbone conformation assumed to be fixed. The 13C(alpha) chemical shift of a given amino acid residue in a protein is determined, mainly, by its own backbone and side-chain torsional angles, independent of the neighboring residues; the conformation of a given residue itself, however, depends on the environment of this residue and, hence, on the whole protein structure. As a consequence, this analysis reveals the role and impact of an accurate side-chain computation in the determination and refinement of protein conformation. The results of this analysis are: (i) a lower error between computed and observed 13C(alpha) chemical shifts (by up to 3.7 ppm), was found for approximately 68% and approximately 63% of all ionizable residues and all non-Ala/Pro/Gly residues, respectively, in the additional set of conformations, compared with results for the model from which the set was derived; and (ii) all the additional conformations exhibit a lower root-mean-square-deviation (1.97 ppm < or = rmsd < or = 2.13 ppm), between computed and observed 13C(alpha) chemical shifts, than the rmsd (2.32 ppm) computed for the starting conformation from which this additional set was derived. As a validation test, an analysis of the additional set of ubiquitin conformations, comparing computed and observed values of both 13C(alpha) chemical shifts and chi(1) torsional angles (given by the vicinal coupling constants, 3J(N-Cgamma) and 3J(C'-Cgamma), is discussed.

Solid-State 19F NMR Spectroscopy Reveals That Trp41 Participates in the Gating Mechanism of the M2 Proton Channel of Influenza A Virus

The integral membrane protein M2 of influenza A virus assembles as a tetrameric bundle to form a proton-conducting channel that is activated by low pH. The side chain of His37 in the transmembrane alpha-helix is known to play an important role in the pH activation of the proton channel. It has also been suggested that Trp41, which is located in an adjacent turn of the helix, forms part of the gating mechanism. Here, a synthetic 25-residue peptide containing the M2 transmembrane domain was labeled with 6F-Trp41 and studied in lipid membranes by solid-state 19F NMR. We monitored the pH-dependent differences in the 19F dipolar couplings and motionally narrowed chemical shift anisotropies of this 6F-Trp41 residue, and we discuss the pH activation mechanism of the H+ channel. At pH 8.0, the structural parameters implicate an inactivated state, while at pH 5.3 the tryptophan conformation represents the activated state. With the aid of COSMOS force field simulations, we have obtained new side-chain torsion angles for Trp41 in the inactivated state (chi1 = -100 degrees +/- 10 degrees , chi2 = +110 degrees +/- 10 degrees ), and we predict a most probable activated state with chi1 = -50 degrees +/- 10 degrees and chi2 = +115 degrees +/- 10 degrees . We have also validated the torsion angles of His37 in the inactivated state as chi1 = -175 degrees +/- 10 degrees and chi2 = -170 degrees +/- 10 degrees .

Efficiency of Rotational Operators for Geometric Manipulation of Chain Molecules

Geometric manipulation of molecules is an essential elementary component in computational modeling programs for molecular structure, stability, dynamics, and design. The computational complexity of transformation of internal coordinates to Cartesian coordinates was discussed before.1 The use of rotation matrices was found to be slightly more efficient than that of quaternion although quaternion operators have been widely advertised for rotational operations, especially in molecular dynamics simulations of liquids where the orientation is a dynamical variable.2 The discussion on computational efficiency is extended here to a more general case in which bond angles and sidechain torsion angles are allowed to vary. The algorithm of Thompson3 is derived again in terms of quaternions as well as rotation matrices, and an algorithm with optimal efficiency is described. The algorithm based on rotation matrices is again found to be slightly more efficient than that based on quaternions.

Intrinsic Differences in the Response of the Human Lutropin Receptor Versus the Human Follitropin Receptor to Activating Mutations

In contrast to the human lutropin receptor (hLHR), very few naturally occurring activating mutations of the structurally related human follitropin receptor (hFSHR) have been identified. The present study was undertaken to determine if one aspect underlying this discrepancy might be a general resistance of the hFSHR to mutation-induced constitutive activity. Five different mutations were introduced into both the hLHR and hFSHR (four based on activating mutations of the hLHR gene, one based on an activating mutation of the hFSHR gene). Our results demonstrate that hFSHR constitutively activating mutants (CAMs) were not as active as hLHR CAMs containing the comparable mutation. Furthermore, although all hFSHR CAMs exhibited strong promiscuous activation by high concentrations of the other glycoprotein hormone receptors, hLHR CAMs showed little or no promiscuous activation. Our in vitro findings are consistent with in vivo observations of known pathophysiological conditions associated with hLHR CAMs, but not hFSHR CAMs, and with promiscuous activation of hFSHR CAMs, but not hLHR CAMs. Computational experiments suggest that the mechanisms through which homologous mutations increase the basal activity of the hLHR and the hFSHR are similar. This is particularly true for the strongest CAMs like L460(3.43)R. Disparate properties of the hLHR versus hFSHR CAMs may, therefore, be due to differences in shape and electrostatics features of the solvent-exposed cytosolic receptor domains involved in the receptor-G protein interface rather than to differences in the nature of local perturbation at the mutation site or in the way local perturbation is transferred to the putative G protein binding domains.

FlexE ensemble docking approach to virtual screening for CDK2 inhibitors

In spite of a proven potential and effectiveness of FlexE in docking flexible ligands into an ensemble of protein structures, FlexE has rarely been successful in virtual screening situations. In this study, we constructed cyclin-dependent kinase 2 (CDK2) ensemble structures which have exactly the same backbone conformations as 1AQ1 but differ only at the side chain torsion angles of the key amino acid (Lys33, Phe80, Lys89 and Asp145) residues: the torsion angles observed in the 17 CDK2 crystal structures were adapted to represent conformational flexibility. FlexE ensemble docking protocol then completely samples the full conformational fields generated by combination of torsions of the four amino acids in the ATP binding site. Virtual screening for CDK2 inhibitors by using the FlexE ensemble docking of a database composed of 48,703 inactives and 82 actives showed significant enrichment factor (EF = 18.5), and successfully identified 71 actives among the top 132 ligands (53.8%) ranked by total energy scores. Moreover, total energy scoring followed by visual inspection filtered-off non-specific binders among the highly-ranked ligands to increase the ratio of actives-to-inactives to 71:13 at the top 5% of the virtual screening solutions.

Use of 13Cα Chemical Shifts in Protein Structure Determination

A physics-based method aimed at determining protein structures by using NOE-derived distances together with observed and computed 13C chemical shifts is proposed. The approach makes use of 13Calpha chemical shifts, computed at the density functional level of theory, to obtain torsional constraints for all backbone and side-chain torsional angles without making a priori use of the occupancy of any region of the Ramachandran map by the amino acid residues. The torsional constraints are not fixed but are changed dynamically in each step of the procedure, following an iterative self-consistent approach intended to identify a set of conformations for which the computed 13Calpha chemical shifts match the experimental ones. A test is carried out on a 76-amino acid, all-alpha-helical protein; namely, the Bacillus subtilis acyl carrier protein. It is shown that, starting from randomly generated conformations, the final protein models are more accurate than an existing NMR-derived structure model of this protein, in terms of both the agreement between predicted and observed 13Calpha chemical shifts and some stereochemical quality indicators, and of similar accuracy as one of the protein models solved at a high level of resolution. The results provide evidence that this methodology can be used not only for structure determination but also for additional protein structure refinement of NMR-derived models deposited in the Protein Data Bank.

Steiner minimal trees, twist angles, and the protein folding problem

The Steiner Minimal Tree (SMT) problem determines the minimal length network for connecting a given set of vertices in three-dimensional space. SMTs have been shown to be useful in the geometric modeling and characterization of proteins. Even though the SMT problem is an NP-Hard Optimization problem, one can define planes within the amino acids that have a surprising regularity property for the twist angles of the planes. This angular property is quantified for all amino acids through the Steiner tree topology structure. The twist angle properties and other associated geometric properties unique for the remaining amino acids are documented in this paper. We also examine the relationship between the Steiner ratio rho and the torsion energy in amino acids with respect to the side chain torsion angle chi(1). The rho value is shown to be inversely proportional to the torsion energy. Hence, it should be a useful approximation to the potential energy function. Finally, the Steiner ratio is used to evaluate folded and misfolded protein structures. We examine all the native proteins and their decoys at http://dd.stanford.edu. and compare their Steiner ratio values. Because these decoy structures have been delicately misfolded, they look even more favorable than the native proteins from the potential energy viewpoint. However, the rho value of a decoy folded protein is shown to be much closer to the average value of an empirical Steiner ratio for each residue involved than that of the corresponding native one, so that we recognize the native folded structure more easily. The inverse relationship between the Steiner ratio and the energy level in the protein is shown to be a significant measure to distinguish native and decoy structures. These properties should be ultimately useful in the ab initio protein folding prediction.

Solution Structure of a Post-transition State Analog of the Phosphotransfer Reaction between the A and B Cytoplasmic Domains of the Mannitol Transporter IIMannitol of the Escherichia coli Phosphotransferase System

The solution structure of the post-transition state complex between the isolated cytoplasmic A (IIAMtl) and phosphorylated B (phospho-IIBMtl) domains of the mannitol transporter of the Escherichia coli phosphotransferase system has been solved by NMR. The active site His-554 of IIAMtl was mutated to glutamine to block phosphoryl transfer activity, and the active site Cys-384 of IIBMtl (residues of IIBMtl are denoted in italic type) was substituted by serine to permit the formation of a stable phosphorylated form of IIBMtl. The two complementary interaction surfaces are predominantly hydrophobic, and two methionines on IIBMtl, Met-388 and Met-393, serve as anchors by interacting with two deep pockets on the surface of IIAMtl. With the exception of a salt bridge between the conserved Arg-538 of IIAMtl and the phosphoryl group of phospho-IIBMtl, electrostatic interactions between the two proteins are limited to the outer edges of the interface, are few in number, and appear to be weak. This accounts for the low affinity of the complex (Kd approximately 3.7 mm), which is optimally tuned to the intact biological system in which the A and B domains are expressed as a single polypeptide connected by a flexible 21-residue linker. The phosphoryl transition state can readily be modeled with no change in protein-protein orientation and minimal perturbations in both the backbone immediately adjacent to His-554 and Cys-384 and the side chains in close proximity to the phosphoryl group. Comparison with the previously solved structure of the IIAMtl-HPr complex reveals how IIAMtl uses the same interaction surface to recognize two structurally unrelated proteins and explains the much higher affinity of IIAMtl for HPr than IIBMtl.

Molecular structure of amyloid fibrils: insights from solid-state NMR

Solid-state nuclear magnetic resonance (NMR) measurements have made major contributions to our understanding of the molecular structures of amyloid fibrils, including fibrils formed by the beta-amyloid peptide associated with Alzheimer's disease, by proteins associated with fungal prions, and by a variety of other polypeptides. Because solid-state NMR techniques can be used to determine interatomic distances (both intramolecular and intermolecular), place constraints on backbone and side-chain torsion angles, and identify tertiary and quaternary contacts, full molecular models for amyloid fibrils can be developed from solid-state NMR data, especially when supplemented by lower-resolution structural constraints from electron microscopy and other sources. In addition, solid-state NMR data can be used as experimental tests of various proposals and hypotheses regarding the mechanisms of amyloid formation, the nature of intermediate structures, and the common structural features within amyloid fibrils. This review introduces the basic experimental and conceptual principles behind solid-state NMR methods that are applicable to amyloid fibrils, reviews the information about amyloid structures that has been obtained to date with these methods, and discusses how solid-state NMR data provide insights into the molecular interactions that stabilize amyloid structures, the generic propensity of polypeptide chains to form amyloid fibrils, and a number of related issues that are of current interest in the amyloid field.