Cα and Cβ Carbon-13 Chemical Shifts in Proteins From an Empirical Database

We have constructed an extensive database of 13C Cα and Cβ chemical shifts in proteins of solution, for proteins of which a high-resolution crystal structure exists, and for which the crystal structure has been shown to be essentially identical to the solution structure. There is no systematic effect of temperature, reference compound, or pH on reported shifts, but there appear to be differences in reported shifts arising from referencing differences of up to 4.2 ppm. The major factor affecting chemical shifts is the backbone geometry, which causes differences of ca. 4 ppm between typical α- helix and β-sheet geometries for Cα, and of ca. 2 ppm for Cβ. The side-chain dihedral angle χ1 has an effect of up to 0.5 ppm on the Cα shift, particularly for amino acids with branched side-chains at Cβ. Hydrogen bonding to main-chain atoms has an effect of up to 0.9 ppm, which depends on the main- chain conformation. The sequence of the protein and ring-current shifts from aromatic rings have an insignificant effect (except for residues following proline). There are significant differences between different amino acid types in the backbone geometry dependence; the amino acids can be grouped together into five different groups with different φ,ψ shielding surfaces. The overall fit of individual residues to a single non-residue-specific surface, incorporating the effects of hydrogen bonding and χ1 angle, is 0.96 ppm for both Cα and Cβ. The results from this study are broadly similar to those from ab initio studies, but there are some differences which could merit further attention.     

基于图论模型的蛋白质结构类型的研究

提出了基于图论模型的H系数分类蛋白质结构为H结型和NH结型的方法.论述了蛋白质结构中序列不相邻的C_α原子之间的空间距离与序列相邻的C_α原子之间空间距离的关系.用此方法对PDB的66个单链蛋白质结构进行分类,结果显示H结型占18.2%.H结在全α型中出现比例较高,在全β型中出现比例较小,所以H结倾向出现在含有α螺旋的蛋白质结构中.     

Refined solution structure of the anti-mammal and anti-insect LqqIII scorpion toxin: Comparison with other scorpion toxins

The solution structure of the anti-mammal and anti-insect LqqIII toxin from the scorpion Leiurus quinquestriatus quinquestriatuswas refined and compared with other long-chain scorpion toxins. This structure, determined by ¹H-NMR and molecular modeling, involves an α-helix (18–29) linked to a three-stranded β-sheet (2–6, 33–39, and 43–51) by two disulfide bridges. The average RMSD between the 15 best structures and the mean structure is 0.71 Å for Cα atoms. Comparison between LqqIII, the potent anti-mammal AaHII, and the weakly active variant-3 toxins revealed that the LqqIII three-dimensional structure is closer to that of AaHII than to the variant-3 structure. Moreover, striking analogies were observed between the electrostatic and hydrophobic potentials of LqqIII and AaHII. Several residues are well conserved in long-chain scorpion toxin sequences and seem to be important in protein structure stability and function. Some of them are involved in the CSαβ (Cysteine Stabilized α-helix β-sheet) motif. A comparison between the sequences of the RII rat brain and the Drosophila extracellular loops forming scorpion toxin binding-sites of Na⁺ channels displays differences in the subsites interacting with anti-mammal or anti-insect toxins. This suggests that hydrophobic as well as electrostatic interactions are essential for the binding and specificity of long-chain scorpion toxins. Proteins 28:360–374, 1997 © 1997 Wiley-Liss, Inc.     

On the satisfaction of backbone‐carbonyl lone pairs of electrons in protein structures

Protein structures are stabilized by a variety of noncovalent interactions (NCIs), including the hydrophobic effect, hydrogen bonds, electrostatic forces and van der Waals’ interactions. Our knowledge of the contributions of NCIs, and the interplay between them remains incomplete. This has implications for computational modeling of NCIs, and our ability to understand and predict protein structure, stability, and function. One consideration is the satisfaction of the full potential for NCIs made by backbone atoms. Most commonly, backbone‐carbonyl oxygen atoms located within α‐helices and β‐sheets are depicted as making a single hydrogen bond. However, there are two lone pairs of electrons to be satisfied for each of these atoms. To explore this, we used operational geometric definitions to generate an inventory of NCIs for backbone‐carbonyl oxygen atoms from a set of high‐resolution protein structures and associated molecular‐dynamics simulations in water. We included more‐recently appreciated, but weaker NCIs in our analysis, such as n→π* interactions, Cα‐H bonds and methyl‐H bonds. The data demonstrate balanced, dynamic systems for all proteins, with most backbone‐carbonyl oxygen atoms being satisfied by two NCIs most of the time. Combinations of NCIs made may correlate with secondary structure type, though in subtly different ways from traditional models of α‐ and β‐structure. In addition, we find examples of under‐ and over‐satisfied carbonyl‐oxygen atoms, and we identify both sequence‐dependent and sequence‐independent secondary‐structural motifs in which these reside. Our analysis provides a more‐detailed understanding of these contributors to protein structure and stability, which will be of use in protein modeling, engineering and design.     

Refinement of human lysozyme at 1.5 Å resolution analysis of non-bonded and hydrogen-bond interactions

The structure of human lysozyme has been crystallographically refined at 1.5 Å resolution by difference map and restrained least-squares procedures to an R factor of 0.187. A comprehensive analysis of the non-bonded and hydrogen-bonded contacts in the lysozyme molecule, which were not restrained, revealed by the refinement has been carried out. The non-bonded CC contacts begin at ~3.45 Å, and the shorter contacts are dominated, as expected, by interactions between trigonal and tetrahedral carbon atoms. The CO contact distances have a “foot” at 3.05 Å. The CN distance plot shows a significant peak at 3.25 Å, which results from close contact between peptide NHs and carbonyl carbons involved in N_iC′_{i ? 2} interactions in α-helices and reverse turns. The distances involving sulphur atoms discriminate SC trigonal interactions at 3.4 to 3.6 Å from SC tetrahedral interactions at 3.7 Å. All these types of non-bonded interactions show minimum distances close to standard van der Waals' separations.Analysis of hydrogen-bond distances has been carried out by using standard geometry to place hydrogen atoms and measuring the XHO distances. On this basis, there are 130 intramolecular hydrogens: 111 NHO bonds, of which 69 are between main-chain atoms, 13 between side-chain atoms and 29 between mainchain and side-chain atoms. If a cluster of four well-defined internal water molecules is included in the protein structure, there is a total of 19 OHO hydrogen bonds. The mean NO, NHO distances and HN?O angles are 2.96 ± 0.17 Å, 2.05 ± 0.18 Å and 18.5 ± 9.6 °, and the mean OO, OHO distances and HÔO angles are 2.83 ± 0.19 Å, 1.98 ± 0.26 Å and 23.8 ± 13.4 °. The distances agree well with standard values, although the hydrogen bonds are consistently more non-linear than in equivalent small molecules. An analysis of the hydrogen-bond angles at the receptor atom indicates that the α-helix, β-sheet and reverse turn have characteristic angular values. A detailed analysis of the regularity of the α-helices and reverse turns shows small but consistent differences between the α-helices in lysozyme and the current standard model, which may now need revision. Of the 21 reverse turns that include a hydrogen bond, the conformations of 19 agree very closely with four of the five standard types. We conclude that the restrained least-squares method of refinement has been validated by these analyses.     

Combining optimal control theory and molecular dynamics for protein folding

A new method to develop low-energy folding routes for proteins is presented. The novel aspect of the proposed approach is the synergistic use of optimal control theory with Molecular Dynamics (MD). In the first step of the method, optimal control theory is employed to compute the force field and the optimal folding trajectory for the Cα atoms of a Coarse-Grained (CG) protein model. The solution of this CG optimization provides an harmonic approximation of the true potential energy surface around the native state. In the next step CG optimization guides the MD simulation by specifying the optimal target positions for the Cα atoms. In turn, MD simulation provides an all-atom conformation whose Cα positions match closely the reference target positions determined by CG optimization. This is accomplished by Targeted Molecular Dynamics (TMD) which uses a bias potential or harmonic restraint in addition to the usual MD potential. Folding is a dynamical process and as such residues make different contacts during the course of folding. Therefore CG optimization has to be reinitialized and repeated over time to accomodate these important changes. At each sampled folding time, the active contacts among the residues are recalculated based on the all-atom conformation obtained from MD. Using the new set of contacts, the CG potential is updated and the CG optimal trajectory for the Cα atoms is recomputed. This is followed by MD. Implementation of this repetitive CG optimization-MD simulation cycle generates the folding trajectory. Simulations on a model protein Villin demonstrate the utility of the method. Since the method is founded on the general tools of optimal control theory and MD without any restrictions, it is widely applicable to other systems. It can be easily implemented with available MD software packages.     

QSE: A new 3-D solvent exposure measure for the analysis of protein structure

Solvent exposure of amino acids measures how deep residues are buried in tertiary structure of proteins, and hence it provides important information for analyzing and predicting protein structure and functions. Existing methods of calculating solvent exposure such as accessible surface area, relative accessible surface area, residue depth, contact number, and half-sphere exposure still have some limitations. In this article, we propose a novel solvent exposure measure named quadrant-sphere exposure (QSE) based on eight quadrants derived from spherical neighborhood. The proposed measure forms a microenvironment around Cα atom as a sphere with a radius of 13??, and subdivides it into eight quadrants according to a rectangular coordinate system constructed based on geometric relationships of backbone atoms. The number of neighboring Cα atoms whose labels are the same is given as the QSE value of the center Cα atom at hand. As evidenced by histograms that show very different distributions for different structure configurations, the proposed measure captures local properties that are characteristic for a residue's eight-directional neighborhood within a sphere. Compared with other measures, QSE provides a different view of solvent exposure, and provides information that is specific for different tertiary structure. As the experimental results show, QSE measure can potentially be used in protein structure analysis and predictions.     

Evaluation of local structure alphabets based on residue burial

Residue burial, which describes a protein residue's exposure to solvent and neighboring atoms, is key to protein structure prediction, modeling, and analysis. We assessed 21 alphabets representing residue burial, according to their predictability from amino acid sequence, conservation in structural alignments, and utility in one fold-recognition scenario. This follows upon our previous work in assessing nine representations of backbone geometry.1 The alphabet found to be most effective overall has seven states and is based on a count of C(beta) atoms within a 14 A-radius sphere centered at the C(beta) of a residue of interest. When incorporated into a hidden Markov model (HMM), this alphabet gave us a 38% performance boost in fold recognition and 23% in alignment quality.     

Patterns in ionizable side chain interactions in protein structures

In a selected set of 44 high-resolution, non-homologous protein structures, the intramolecular hydrogen bonds or salt bridges formed by ionizable amino acid side chains were identified and analyzed. The analysis was based on the investigation of several properties of the involved residues such as their solvent exposure, their belonging to a certain secondary structural element, and their position relative to the N- and C-termini of their respective structural element. It was observed that two-thirds of the interactions made by basic or acidic side chains are hydrogen bonds to polar uncharged groups. In particular, the majority (78%) of the hydrogen bonds between ionizable side chains and main chain polar groups (sch:mch bonds) involved at least one buried atom, and in 42% of the cases both interacting atoms were buried. In α-helices, the sch:mch bonds observed in the proximity of the C- and N-termini show a clear preference for acidic and basic side chains, respectively. This appears to be due to the partial charges of peptide group atoms at the termini of α-helices, which establish energetically favorable electrostatic interactions with side chain carrying opposite charge, at distances even greater than 4.5 Å. The sch:mch interactions involving ionizable side chains that belong either to β-strands or to the central part of α-helices are based almost exclusively on basic residues. This results from the presence of main chain carbonyl oxygen atoms in the protein core which have unsatisfied hydrogen bonding capabilities.     

The thiol proteinase inhibitor E-64-d ameliorates amyloid-β-induced reduction of sAPPα secretion by reversing ceramide-induced protein kinase C down-regulation in SH-SY5Y neuroblastoma cells

In Alzheimer’s disease (AD), enhancing α-secretase processing of amyloid precursor protein (APP) is an important pathway to decrease neurotoxic amyloid β (Aβ) secretion. The α-secretase is reported to be regulated by protein kinase C (PKC) and various endogenous proteins or cell surface receptors. In this report, we first examined whether Aβ reduces α-secretase activity, and showed that Aβ peptide 1–40 (0.001 and 0.01 μM) reduced the secretion of soluble amyloid precursor protein α (sAPPα) in carbachol-stimulated SH-SY5Y neuroblastoma cells. E-64-d (3 μM), which is a potent calpain inhibitor that prevents PKC degradation, ameliorated the Aβ-induced reduction of sAPPα secretion. In addition, we observed that Aβ significantly enhanced ceramide production by activating neutral sphingomyelinase. The cell-permeable ceramide analog, C₂-ceramide (1 μg/mL), also reduced sAPPα secretion, and in addition, E-64-d eliminated the observed decrease of sAPPα secretion. C₂-ceramide induced down-regulation of PKC-α, -β₁, and -β₂ isozymes in SH-SY5Y cells. These findings suggest that ceramide may play an important role in sAPPα processing by modulating PKC activity.     

Accurate geometries for “Mountain pass” regions of the Ramachandran plot using quantum chemical calculations

Unusual local arrangements of protein in Ramachandran space are not well represented by standard geometry tools used in either protein structure refinement using simple harmonic geometry restraints or in protein simulations using molecular mechanics force fields. In contrast, quantum chemical computations using small poly‐peptide molecular models can predict accurate geometries for any well‐defined backbone Ramachandran orientation. For conformations along transition regions—ϕ from −60 to 60°—a very good agreement with representative high‐resolution experimental X‐ray (≤1.5 Å) protein structures is obtained for both backbone C⁻¹‐N‐C_α angle and the nonbonded O⁻¹…C distance, while “standard geometry” leads to the “clashing” of O…C atoms and Amber FF99SB predicts distances too large by about 0.15 Å. These results confirm that quantum chemistry computations add valuable support for detailed analysis of local structural arrangements in proteins, providing improved or missing data for less understood high‐energy or unusual regions.     

Hydrogen bonds involving sulfur atoms in proteins.

Intrachain hydrogen bonds are a hallmark of globular proteins. Traditionally, these involve oxygen and nitrogen atoms. The electronic structure of sulfur is compatible with hydrogen bond formation as well. We surveyed a set of 85 high-resolution protein structures in order to evaluate the prevalence and geometry of sulfur-containing hydrogen bonds. This information should be of interest to experimentalists and theoreticians interested in protein structure and protein engineering.     

Analysis of six protein structures predicted by comparative modeling techniques

The protein structures of six comparative modeling targets were predicted in a procedure that relied on improved energy minimization, without empirical rules, to position all new atoms. The structures of human nucleoside diphosphate kinase NM23-H2, HPr from Mycoplasma capricolum, 2Fe-2S ferredoxin from Haloarcula marismortui, eosinophil-derived neurotoxin (EDN), mouse cellular retinoic acid protein I (CRABP1), and P450eryf were predicted with root mean square deviations on Cα atoms of 0.69, 0.73, 1.11, 1.48, 1.69, and 1.73 Å, respectively, compared to the target crystal structures. These differences increased as the sequence similarity between the target and parent proteins decreased from about 60 to 20% identity. More residues were predicted than form the common region shared by the two crystal structures. In most cases insertions or deletions between the target and the related protein of known structure were not correctly positioned. One two residue insertion in CRABP1 was predicted in the correct conformation, while a nine residue insertion in EDN was predicted in the correct spatial region, although not in the correct conformation. The positions of common cofactors and their binding sites were predicted correctly, even when overall sequence similarity was low. © 1995 Wiley-Liss, Inc.     

Interaction between bound water molecules and local protein structures: A statistical analysis of the hydrogen bond structures around bound water molecules

Water molecules play an important role in protein folding and protein interactions through their structural association with proteins. Examples of such structural association can be found in protein crystal structures, and can often explain protein functionality in the context of structure. We herein report the systematic analysis of the local structures of proteins interacting with water molecules, and the characterization of their geometric features. We first examined the interaction of water molecules with a large local interaction environment by comparing the preference of water molecules in three regions, namely, the protein–protein interaction (PPI) interfaces, the crystal contact (CC) interfaces, and the non‐interfacial regions. High preference of water molecules to the PPI and CC interfaces was found. In addition, the bound water on the PPI interface was more favorably associated with the complex interaction structure, implying that such water‐mediated structures may participate in the shaping of the PPI interface. The pairwise water‐mediated interaction was then investigated, and the water‐mediated residue–residue interaction potential was derived. Subsequently, the types of polar atoms surrounding the water molecules were analyzed, and the preference of the hydrogen bond acceptor was observed. Furthermore, the geometries of the structures interacting with water were analyzed, and it was found that the major structure on the protein surface exhibited planar geometry rather than tetrahedral geometry. Several previously undiscovered characteristics of water–protein interactions were unfolded in this study, and are expected to lead to a better understanding of protein structure and function. Proteins 2016; 84:43–51. © 2015 Wiley Periodicals, Inc.     

Differential Isotope Labeling Strategy for Determining the Structure of Myristoylated Recoverin by NMR Spectroscopy

The three-dimensional solution structure of recombinant bovine myristoylated recoverin in the Ca2+-free state has been refined using an array of isotope-assisted multidimensional heteronuclear NMR techniques. In some experiments, the myristoyl group covalently attached to the protein N-terminus was labeled with 13C and the protein was unlabeled or vice versa; in others, both were 13C-labeled. This differential labeling strategy was essential for structural refinement and can be applied to other acylated proteins. Stereospecific assignments of 41 pairs of -methylene protons and 48 methyl groups of valine and leucine were included in the structure refinement. The refined structure was constructed using a total of 3679 experimental NMR restraints, comprising 3242 approximate interproton distance restraints (including 153 between the myristoyl group and the polypeptide), 140 distance restraints for 70 backbone hydrogen bonds, and 297 torsion angle restraints. The atomic rms deviations about the averaged minimized coordinate positions for the secondary structure region of the N-terminal and C-terminal domains are 0.44 ± 0.07 and 0.55 ± 0.18 Å for backbone atoms, and 1.09 ± 0.07 and 1.10 ± 0.15 Å for all heavy atoms, respectively. The refined structure allows for a detailed analysis of the myristoyl binding pocket. The myristoyl group is in a slightly bent conformation: the average distance between C1 and C14 atoms of the myristoyl group is 14.6 Å. Hydrophobic residues Leu28, Trp31, and Tyr32 form a cluster that interacts with the front end of the myristoyl group (C1-C8), whereas residues Phe49, Phe56, Tyr86, Val87, and Leu90 interact with the tail end (C9-C14). The relatively deep hydrophobic pocket that binds the myristoyl group (C14:0) could also accommodate other naturally occurring acyl groups such as C12:0, C14:1, and C14:2 chains.     

Comprehensive mutagenesis of HIV-1 protease: a computational geometry approach

A computational geometry technique based on Delaunay tessellation of protein structure, represented by C(alpha) atoms, is used to study effects of single residue mutations on sequence-structure compatibility in HIV-1 protease. Profiles of residue scores derived from the four-body statistical potential are constructed for all 1881 mutants of the HIV-1 protease monomer and compared with the profile of the wild-type protein. The profiles for an isolated monomer of HIV-1 protease and the identical monomer in a dimeric state with an inhibitor are analyzed to elucidate changes to structural stability. Protease residues shown to undergo the greatest impact are those forming the dimer interface and flap region, as well as those known to be involved in inhibitor binding.     

The surface of β-sheet proteins contains amphiphilic regions which may provide clues about protein folding

A major bottleneck in the field of biochemistry is our limited understanding of the processes by which a protein folds into its native conformation. Much of the work on this issue has focused on the conserved core of the folded protein. However, one might imagine that a ubiquitous motif for unaided folding or for the recognition of chaperones may involve regions on the surface of the native structure. We explore this possibility by an analysis of the spatial distribution of regions with amphiphilic α-helical potential on the surface of β-sheet proteins. All proteins, Including β-sheet proteins, contain regions with amphiphilic α-helical potential. That is, any α-helix formed by that region would be amphiphilic, having both hydrophobic and hydrophilic surfaces. In the three-dimensional structure of all β-sheet proteins analyzed, we have found a distinct pattern in the spatial distribution of sequences with amphiphilic α-helical potential. The amphiphilic regions occur in ring shaped clusters approximately 20 to 30 Å in diameter on the surface of the protein. In addition, these regions have a strong preference for positively charged amino acids and a lower preference for residues not favorable to α-helix formation. Although the purpose of these amphiphilic regions which are not associated with naturally occurring α-helix is unknown, they may play a critical role in highly conserved processes such as protein folding. © 1996 Wiley-Liss, Inc.