Sequential resonance assignments in protein 1H nuclear magnetic resonance spectra: Computation of sterically allowed proton-proton distances and statistical analysis of proton-proton distances in single crystal protein conformations

Two different, theoretical studies of intramolecular proton-proton distances in polypeptide chains are described. Firstly, the distances between amide, C^α and C^β protons of neighbouring residues in the amino acid sequence, which correspond to the sterically allowed values for the dihedral angles φ_i, ψ_i and χ_i¹, were computed. Secondly, the frequency with which short distances occur between amide, C^α and C^β protons of neighbouring and distant residues in the amino acid sequence were statistically evaluated in a representative sample of globular protein crystal structures. Both approaches imply that semi-quantitative measurements of short, non-bonding proton-proton distances, e.g. by nuclear Overhauser experiments, should present a reliable and generally applicable method for sequential, individual resonance assignments in protein ¹H nuclear magnetic resonance spectra. Similar calculations imply that corresponding distance measurements can be used for resonance assignments in the side-chains of the aromatic amino acid residues, asparagine and glutamine, where the complete spin systems cannot usually be identified from through-bond spin-spin coupling connectivities.     

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.     

Application of Elastic Network Models to Proteins in the Crystalline State

Normal mode analysis using elastic network models has grown popular for probing the low-frequency collective dynamics of proteins and other biomolecular assemblies. In most previous studies, these models were validated by comparing calculated atomic fluctuations for isolated proteins with experimental temperature factors determined in the crystalline state, although there were also hints that including crystal contacts in the calculations has an impact on the comparison. In this study, a set of 83 ultra-high resolution crystal structures with experimentally determined anisotropic displacement parameters is used to evaluate several C_α-based elastic network models that either ignore or treat the crystal environment in different ways; the latter include using periodic boundary conditions defined with respect to the asymmetric unit or the primitive unit cell as well as using the Born-von Kármán boundary condition that accounts for lattice vibrations. For all elastic network models, treating the crystal environment leads to better agreement with experimental anisotropic displacement parameters with the Born-von Kármán boundary condition giving the best agreement. Atomic correlations over the entire protein are clearly affected by the presence of the crystal contacts and fairly sensitive to the way that the crystal environment is treated. These observations highlight the importance of properly treating the protein system in an environment consistent with experiment when either evaluating approximate protein models or using approximate dynamic models in structural refinement application types. Finally, investigation of the scaling behaviors of the cumulative density of states and the heat capacity indicates that there are still gaps between simplified elastic models and all-atom models for proteins.     

Spatial structure of dimeric a genetically engineered variant of green fluorescent protein EGFP-K162Q in the P61 crystal space group

The spatial structure of dimeric green fluorescent protein EGFP-K162Q with MDELYK (EGFPv) C-terminal deletion has been assigned in the P6₁ space group with resolution 1.34 Å by X-ray diffraction analysis. The results have been compared with X-ray diffraction data of monomeric EGFP (green biomarker with enhanced photophysical properties) assigned in another crystal space group, P2₁2₁2₁, with resolution 1.50 and 1.35 Å. Subunits in the EGFPv dimeric structure are located at 75° angle with the contact area ～800 Ų. The dimeric framework is stabilized by the six hydrogen bonds and central hydrophobic core of six residues. The root-mean-square deviation value for C^α atoms in 3–230 residues of the P6₁ and P2₁2₁2₁ crystal structures is 0.55 Å. The differential characteristics of EGFPv-P6₁ structure, compared to that of P2₁2₁2₁, is a noticeably different orientation of the Glu222 side chain, and a new conformation of the 155–159 loop fragment, characterized by deviations among the C^α atoms of superimposed structures reaching 4.6 Å for Lys156 and 5.5 Å for Lys158.     

Variability of the canonical loop conformations in serine proteinases inhibitors and other proteins

Canonical loops of protein inhibitors of serine proteinases occur in proteins having completely different folds. In this article, conformations of the loops have been analyzed for inhibitors belonging to 10 structurally different families. Using deviation in C_α-C_α distances as a criterion for loop similarity, we found that the P3-P3′ segment defines most properly the length of the loop. When conformational differences among loops of individual inhibitors were compared using root mean square deviation (rmsd) in atomic coordinates for all main chain atoms (Δr method) and rmsd operating in main chain torsion angles (Δt method), differences of up to 2.1 Å and 72.3°, respectively, were observed. Such large values indicate significant conformational differences among individual loops. Nevertheless, the overall geometry of the inhibitor-proteinase interaction is very well preserved, as judged from the similarity of C_α-C_α distances between C_α of catalytic Ser and C_α of P3-P3′ residues in various enzyme-inhibitor complexes. The mode of interaction is very well preserved both in the chymotrypsin and subtilisin families, as distances calculated for subtilisin-inhibitor complexes are almost always within the range of those for chymotrypsin-inhibitor complexes. Complex formation leads to conformational changes of up to 160° for χ₁ angle. Side chains of residue P1 and P2′ adopt in different complexes a similar orientation (χ₁ angle = −60° and −180°, respectively). To check whether the canonical conformation can be found among non–proteinase-inhibitor Brookhaven Protein Data Bank structures, two selection criteria—the allowed main chain dihedral angles and C_α-C_α distances for the P3-P3′ segment—were applied to all these structures. This procedure detected 132 unique hexapeptide segments in 121 structurally and functionally unrelated proteins. Partial preferences for certain amino acids occurring at particular positions in these hexapeptides could be noted. Further restriction of this set to hexapeptides with a highly exposed P1 residue side chain resulted in 40 segments. The possibility of complexes formation between these segments and serine proteinases was ruled out in molecular modeling due to steric clashes. Several structural features that determine the canonical conformation of the loop both in inhibitors and in other proteins can be distinguished. They include main chain hydrogen bonds both within the P3-P3′ segment and with the scaffold region, P3-P4 and P3′-P4′ hydrophobic interactions, and finally either hydrophobic or polar interactions involving the P1′ residue. Proteins 32:459–474, 1998. © 1998 Wiley-Liss, Inc.     

A three dimensional visualisation approach to protein heavy-atom structure reconstruction

Background

A commonly recurring problem in structural protein studies, is the determination of all heavy atom positions from the knowledge of the central α-carbon coordinates.

Results

We employ advances in virtual reality to address the problem. The outcome is a 3D visualisation based technique where all the heavy backbone and side chain atoms are treated on equal footing, in terms of the C_α coordinates. Each heavy atom is visualised on the surfaces of a different two-sphere, that is centered at another heavy backbone and side chain atoms. In particular, the rotamers are visible as clusters, that display a clear and strong dependence on the underlying backbone secondary structure.

Conclusions

We demonstrate that there is a clear interdependence between rotameric states and secondary structure. Our method easily detects those atoms in a crystallographic protein structure which are either outliers or have been likely misplaced, possibly due to radiation damage. Our approach forms a basis for the development of a new generation, visualization based side chain construction, validation and refinement tools. The heavy atom positions are identified in a manner which accounts for the secondary structure environment, leading to improved accuracy.

     

Distance restraints from crosslinking mass spectrometry: Mining a molecular dynamics simulation database to evaluate lysine–lysine distances

Integrative structural biology attempts to model the structures of protein complexes that are challenging or intractable by classical structural methods (due to size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such experimental method is chemical crosslinking mass spectrometry (XL‐MS), in which protein complexes are crosslinked and characterized using liquid chromatography‐mass spectrometry to pinpoint specific amino acid residues in close structural proximity. The commonly used lysine‐reactive N‐hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS³) have a linker arm that is 11.4 Å long when fully extended, allowing C_α (alpha carbon of protein backbone) atoms of crosslinked lysine residues to be up to ～24 Å apart. However, XL‐MS studies on proteins of known structure frequently report crosslinks that exceed this distance. Typically, a tolerance of ～3 Å is added to the theoretical maximum to account for this observation, with limited justification for the chosen value. We used the Dynameomics database, a repository of high‐quality molecular dynamics simulations of 807 proteins representative of diverse protein folds, to investigate the relationship between lysine–lysine distances in experimental starting structures and in simulation ensembles. We conclude that for DSS/BS³, a distance constraint of 26–30 Å between C_α atoms is appropriate. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL‐MS results to structures or in modeling. We also discuss the comparison of XL‐MS results to MD simulations and known structures as a means to test and validate experimental XL‐MS methods.     

Methyl 6-deoxy-6-iodo-α-d-glucopyranoside: crystal structure and high-resolution NMR spectroscopy

Single-crystal X-ray diffraction and high-resolution ¹H and ¹³C NMR spectral data for methyl 6-deoxy-6-iodo-α-d-glucopyranoside are reported. The ⁴C₁ conformation was found to be the preferred form for this compound, both in the crystal lattice and in solution. The rotational preferences of all the groups bound to the pyranose ring are presented. The stabilization of the crystal structure by a network of O-H···O intra- and intermolecular interactions as well as the short contacts of the iodine atoms is discussed.     

Crystal structure of bis(N-methylethylenediammonium) hexabromodicuprate(I)

The crystal structure of the title compound (C₃H₁₂N₂)₂Cu₂Br₆, is monoclinic, space group P2₁/n, with lattice constants a=8.222(2), b=11.214(2), c=10.646(2) Å, β=91.97(1)° and V=981.0(3) ų. The structure contains anionic Cu₂Br₆²⁻ dimers and N-methylethylenediammonium (henceforth MEDA²⁺) cations. The centrosymmetric dimer units are composed of two edge-shared CuBr₄ tetrahedra, with bridging Cu-Br distances 0.08 Å longer than the terminal distances. A pair of MEDA²⁺ cations hydrogen bond to each dimer via two of the -NH₃⁺ hydrogen atoms and one of the -NH₂⁺- hydrogen atoms. Additional hydrogen bonding between the cations and anions tie the structure together in a complex supramolecular network.     

Characterization of a mononuclear copper carboxylate complex: Bis(acetylsalicylato)bis(pyridine)copper(II)

The preparation, spectral properties, and crystal structure of a mononuclear copper(II) complex of acetylsalicylate and pyridine are reported. The complex exists as bis(acetylsalicylato)bis(pyridine)copper(II) both in the solid state and in chloroform solution. The crystal is monoclinic, space group P2₁/n, with a = 17.823(5), b = 10.903(4), c = 6.598(2) Å, β = 95.74(2)°. The final refinement used 1472 observed reflections and gave an R of 0.046. The copper atom is surrounded by four atoms in a trans square planar arrangement with two short CuO distances of 1.949(3) Å and two CuN distances of 2.003(4) Å. Two longer CuO distances of 2.623(3) Å are made with the remaining oxygen atoms of the aspirin carboxylate groups.     

Detecting protein atom correlations using correlation of probability of recurrence

The dynamic cross‐correlation Map(DCCM) technique has been used extensively to study protein dynamics. In this work, we introduce the use of the method of correlation of probability of recurrence (CPR) as a complementary method to detect correlations between protein residue atoms. Time series of the distances of the C_α atoms of the β‐lactamase inhibitory protein (BLIP) from a reference position are analyzed using CPR and mutual information (MI). The results are compared to those provided by DCCM. In comparison to MI, CPR is found to detect more of the correlations present in DCCM. It is also able to detect a small number of significant correlations between distant residues that are not detected by DCCM. Proteins 2014; 82:2180–2189. © 2014 Wiley Periodicals, Inc.     

Coarse-Grained Simulation of Full-Length Integrin Activation

Integrin conformational dynamics are critical to their receptor and signaling functions in many cellular processes, including spreading, adhesion, and migration. However, assessing integrin conformations is both experimentally and computationally challenging because of limitations in resolution and dynamic sampling. Thus, structural changes that underlie transitions between conformations are largely unknown. Here, focusing on integrin αvβ₃, we developed a modified form of the coarse-grained heterogeneous elastic network model (hENM), which allows sampling conformations at the onset of activation by formally separating local fluctuations from global motions. Both local fluctuations and global motions are extracted from all-atom molecular dynamics simulations of the full-length αvβ₃ bent integrin conformer, but whereas the former are incorporated in the hENM as effective harmonic interactions between groups of residues, the latter emerge by systematically identifying and treating weak interactions between long-distance domains with flexible and anharmonic connections. The new hENM model allows integrins and single-point mutant integrins to explore various conformational states, including the initiation of separation between α- and β-subunit cytoplasmic regions, headpiece extension, and legs opening.     

Predictions of protein backbone bond distances and angles from first principles

A procedure is described, based on a spline-function representation of ab initio peptide conformational geometry maps, that allows one to predict backbone bond distances and angles of proteins as functions of the peptide ?(N-C^α)/Ψ(C^α-C′) torsions with an accuracy comparable to that of high-resolution protein crystallography. For example, for the more than 40 residues of crambin, the rms deviation between predicted and crystallographic values of N-C^α-C′ is 1.9° for the 1.5 Å resolution structure and 1.9° for the 0.83 Å resolution structure, compared with angle variations of < 10°. Accurate information on protein backbone geometries is important for establishing dictionaries of flexible geometry functions for use in empirical peptide and protein modeling. © 1995 John Wiley & Sons, Inc.     

The characterization of pc-polylines representing protein backbones

The backbone of a protein is typically represented as either a C _α-polyline, a three-dimensional (3D) polyline that passes through the C _α atoms, or a tuple of ϕ,ψ pairs while its fold is usually assigned using the 3D topological arrangement of the secondary structure elements (SSEs). It is tricky to obtain the SSE composition for a protein from the C _α-polyline representation while its 3D SSE arrangement is not apparent in the two-dimensional (2D) ϕ,ψ representation. In this article, we first represent the backbone of a protein as a pc-polyline that passes through the centers of its peptide planes. We then analyze the pc-polylines for six different sets of proteins with high quality crystal structures. The results show that SSE composition becomes recognizable in pc-polyline presentation and consequently the geometrical property of the pc-polyline of a protein could be used to assign its secondary structure. Furthermore, our analysis finds that for each of the six sets the total length of a pc-polyline increases linearly with the number of the peptide planes. Interestingly a comparison of the six regression lines shows that they have almost identical slopes but different intercepts. Most interestingly there exist decent linear correlations between the intercepts of the six lines and either the average helix contents or the average sheet contents and between the intercepts and the average backbone hydrogen bonding energetics. Finally, we discuss the implications of the identified correlations for structure classification and protein folding, and the potential applications of pc-polyline representation to structure prediction and protein design.     

Crystal structures of diketopiperazines containing α-aminoisobutyric acid: Cyclo(Aib-Aib) and cyclo(Aib-L-Ile)

The crystal and molecular structures of two α-aminoisobutyric acid (Aib)-containing diketopiperazines, cyclo(Aib-Aib) 1 and cyclo(Aib-L -Ile) 2 , are reported. Cyclo(Aib-Aib) crystallizes in the space group P1 with a = 5.649(3), b = 5.865(2), c = 8.363(1), α = 69.89(6), β = 113.04(8), γ = 116.0(3), and Z = 1, while 2 occurs in the space group P2₁2₁2₁ with a = 6.177(1), b = 10.791(1), c = 16.676(1), and Z = 4. The structures of 1 and 2 have been refined to final R factors of 0.085 and 0.086, respectively. In both structures the diketopiperazine ring shows small but significant deviation from planarity. A very flat chair conformation is adopted by 1, in which the C^α atoms are displaced by 0.07 Å on each side of the mean plane, passing through the other four atoms of the ring. Cyclo(Aib-Ile) favors a slight boat conformation, with Aib C^α and Ile C^α atoms displaced by 0.11 and 0.05 Å on the same side of the mean plane formed by the other ring atoms. Structural features in these two molecules are compared with other related diketopiperazines.     

Backbone chemical shift assignments of the acyl-acyl carrier protein intermediates of the fatty acid biosynthesis pathway of Plasmodium falciparum

We report the backbone chemical shift assignments of the acyl-acyl carrier protein (ACP) intermediates of the fatty acid biosynthesis pathway of Plasmodium falciparum. The acyl-ACP intermediates butyryl (C₄), -octanoyl (C₈), -decanoyl (C₁₀), -dodecanoyl (C₁₂) and -tetradecanoyl (C₁₄)-ACPs display marked changes in backbone HN, C_α and C_β chemical shifts as a result of acyl chain insertion into the hydrophobic core. Chemical shift changes cast light on the mechanism of expansion of the acyl carrier protein core.