3d interaction homology: The structurally known rotamers of tyrosine derive from a surprisingly limited set of information‐rich hydropathic interaction environments described by maps

Sidechain rotamer libraries are obtained through exhaustive statistical analysis of existing crystallographic structures of proteins and have been applied in multiple aspects of structural biology, for example, crystallography of relatively low‐resolution structures, in homology model building and in biomolecular NMR. Little is known, however, about the driving forces that lead to the preference or suitability of one rotamer over another. Construction of 3D hydropathic interaction maps for nearly 30,000 tyrosines reveals the environment around each, in terms of hydrophobic (π–π stacking, etc.) and polar (hydrogen bonding, etc.) interactions. After partitioning the tyrosines into backbone‐dependent (?, ψ) bins, a map similarity metric based on the correlation coefficient was applied to each map‐map pair to build matrices suitable for clustering with k‐means. The first bin (?200° ≤ ? < –155°; ?205° ≤ ψ < –160°), representing 631 tyrosines, reduced to 14 unique hydropathic environments, with most diversity arising from favorable hydrophobic interactions with many different residue partner types. Polar interactions for tyrosine include surprisingly ubiquitous hydrogen bonding with the phenolic OH and a handful of unique environments surrounding the tyrosine backbone. The memberships of all but one of the 14 environments are dominated (>50%) by a single χ₁/χ₂ rotamer. The last environment has weak or no interactions with the tyrosine ring and its χ₁/χ₂ rotamer is indeterminate, which is consistent with it being composed of mostly surface residues. Each tyrosine residue attempts to fulfill its hydropathic valence and thus, structural water molecules are seen in a variety of roles throughout protein structure. Proteins 2015; 83:1118–1136. © 2015 Wiley Periodicals, Inc.     

ω‐Turn: A novel β‐turn mimic in globular proteins stabilized by main‐chain to side‐chain CH···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10‐membered hydrogen‐bonded ring involving the main‐chain C?O in a β‐turn can be formed using a side‐chain carbonyl group leading to Asx‐turn. We show that the N? H component of hydrogen bond can be replaced by a C^γ‐H group in the side chain, culminating in a nonconventional C? H···O interaction. Because of its shape this β‐turn mimic is designated as ω‐turn, which is found to occur ～three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C? H···O interaction occurring between the terminal residues, constraining the torsion angles ?_{i + 1}, ψ_{i + 1}, ?_{i + 2} and χ′_{1(i + 2)} (using the interacting C^γ atom). Based on these angles there are two types of ω‐turns, each of which can be further divided into two groups. C^β‐branched side‐chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal‐binding sites. N‐linked glycosylation occurs at the consensus pattern Asn‐Xaa‐Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω‐turn, which may be the recognition site for protein modification. Location between two β‐strands is the most common occurrence in protein tertiary structure, and being generally exposed ω‐turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design. Proteins 2015; 83:203–214. © 2014 Wiley Periodicals, Inc.     

Protein structure determination by conformational space annealing using NMR geometric restraints

We have carried out numerical experiments to investigate the applicability of the global optimization method of conformational space annealing (CSA) to the enhanced NMR protein structure determination over existing PDB structures. The NMR protein structure determination is driven by the optimization of collective multiple restraints arising from experimental data and the basic stereochemical properties of a protein‐like molecule. By rigorous and straightforward application of CSA to the identical NMR experimental data used to generate existing PDB structures, we redetermined 56 recent PDB protein structures starting from fully randomized structures. The quality of CSA‐generated structures and existing PDB structures were assessed by multiobjective functions in terms of their consistencies with experimental data and the requirements of protein‐like stereochemistry. In 54 out of 56 cases, CSA‐generated structures were better than existing PDB structures in the Pareto‐dominant manner, while in the remaining two cases, it was a tie with mixed results. As a whole, all structural features tested improved in a statistically meaningful manner. The most improved feature was the Ramachandran favored portion of backbone torsion angles with about 8.6% improvement from 88.9% to 97.5% (P‐value <10^?17). We show that by straightforward application of CSA to the efficient global optimization of an energy function, NMR structures will be of better quality than existing PDB structures. Proteins 2015; 83:2251–2262. © 2015 Wiley Periodicals, Inc.     

Optimized distance‐dependent atom‐pair‐based potential DOOP for protein structure prediction

The DOcking decoy‐based Optimized Potential (DOOP) energy function for protein structure prediction is based on empirical distance‐dependent atom‐pair interactions. To optimize the atom‐pair interactions, native protein structures are decomposed into polypeptide chain segments that correspond to structural motives involving complete secondary structure elements. They constitute near native ligand–receptor systems (or just pairs). Thus, a total of 8609 ligand–receptor systems were prepared from 954 selected proteins. For each of these hypothetical ligand–receptor systems, 1000 evenly sampled docking decoys with 0–10 Å interface root‐mean‐square‐deviation (iRMSD) were generated with a method used before for protein–protein docking. A neural network‐based optimization method was applied to derive the optimized energy parameters using these decoys so that the energy function mimics the funnel‐like energy landscape for the interaction between these hypothetical ligand–receptor systems. Thus, our method hierarchically models the overall funnel‐like energy landscape of native protein structures. The resulting energy function was tested on several commonly used decoy sets for native protein structure recognition and compared with other statistical potentials. In combination with a torsion potential term which describes the local conformational preference, the atom‐pair‐based potential outperforms other reported statistical energy functions in correct ranking of native protein structures for a variety of decoy sets. This is especially the case for the most challenging ROSETTA decoy set, although it does not take into account side chain orientation‐dependence explicitly. The DOOP energy function for protein structure prediction, the underlying database of protein structures with hypothetical ligand–receptor systems and their decoys are freely available at http://agknapp.chemie.fu‐berlin.de/doop/ . Proteins 2015; 83:881–890. © 2015 Wiley Periodicals, Inc.     

Flap flexibility amongst plasmepsins I,II, III,IV, and V: Sequence,structural, and molecular dynamics analyses

Herein, for the first time, we comparatively report the opening and closing of apo plasmepsin I – V. Plasmepsins belong the aspartic protease family of enzymes, and are expressed during the various stages of the P. falciparum lifecycle, the species responsible for the most lethal and virulent malaria to infect humans. Plasmepsin I, II, IV and HAP degrade hemoglobin from infected red blood cells, whereas plasmepsin V transport proteins crucial to the survival of the malaria parasite across the endoplasmic reticulum. Flap‐structures covering the active site of aspartic proteases (such as HIV protease) are crucial to the conformational flexibility and dynamics of the protein, and ultimately control the binding landscape. The flap‐structure in plasmepsins is made up of a flip tip in the N‐terminal lying perpendicular to the active site, adjacent to the flexible loop region in the C‐terminal. Using molecular dynamics, we propose three parameters to better describe the opening and closing of the flap‐structure in apo plasmepsins. Namely, the distance, d₁, between the flap tip and the flexible region; the dihedral angle, ?, to account for the twisting motion; and the TriCα angle, θ₁. Simulations have shown that as the flap‐structure twists, the flap and flexible region move apart opening the active site, or move toward each other closing the active site. The data from our study indicate that of all the plasmepsins investigated in the present study, Plm IV and V display the highest conformational flexibility and are more dynamic structures versus Plm I, II, and HAP. Proteins 2015; 83:1693–1705. © 2015 Wiley Periodicals, Inc.     

Toward optimal fragment generations for ab initio protein structure assembly

Fragment assembly using structural motifs excised from other solved proteins has shown to be an efficient method for ab initio protein‐structure prediction. However, how to construct accurate fragments, how to derive optimal restraints from fragments, and what the best fragment length is are the basic issues yet to be systematically examined. In this work, we developed a gapless‐threading method to generate position‐specific structure fragments. Distance profiles and torsion angle pairs are then derived from the fragments by statistical consistency analysis, which achieved comparable accuracy with the machine‐learning‐based methods although the fragments were taken from unrelated proteins. When measured by both accuracies of the derived distance profiles and torsion angle pairs, we come to a consistent conclusion that the optimal fragment length for structural assembly is around 10, and at least 100 fragments at each location are needed to achieve optimal structure assembly. The distant profiles and torsion angle pairs as derived by the fragments have been successfully used in QUARK for ab initio protein structure assembly and are provided by the QUARK online server at http://zhanglab.ccmb. med.umich.edu/QUARK/ . Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Predicting kinetic constants of protein-protein interactions based on structural properties

Elucidating kinetic processes of protein–protein interactions (PPI) helps to understand how basic building blocks affect overall behavior of living systems. In this study, we used structure‐based properties to build predictive models for kinetic constants of PPI. A highly diverse PPI dataset, protein–protein kinetic interaction data and structures (PPKIDS), was built. PPKIDS contains 62 PPI with complex structures and kinetic constants measured experimentally. The influence of structural properties on kinetics of PPI was studied using 35 structure‐based features, describing different aspects of complex structures. Linear models for the prediction of kinetic constants were built by fitting with selected subsets of structure‐based features. The models gave correlation coefficients of 0.801, 0.732, and 0.770 for k_off, k_on, and K_d, respectively, in leave‐one‐out cross validations. The predictive models reported here use only protein complex structures as input and can be generally applied in PPI studies as well as systems biology modeling. Our study confirmed that different properties play different roles in the kinetic process of PPI. For example, k_on was affected by overall structural features of complexes, such as the composition of secondary structures, the change of translational and rotational entropy, and the electrostatic interaction; while k_off was determined by interfacial properties, such as number of contacted atom pairs per 100 Ų. This information provides useful hints for PPI design. Proteins 2010;79:720–734. © 2010 Wiley‐Liss, Inc.     

Shape evolution with temperature of a thermotolerant protein (PeaT1) in solution detected by small angle X‐ray scattering

The protein elicitor from Alternaria tenuissima (PeaT1) presented excellent thermotolerance and potential application in agriculture as a pesticide. Previous synchrotron radiation circular dichroism study demonstrated that the secondary structures in PeaT1 protein are reversible with temperature change. To further clarify the mechanism of its thermotolerance, synchrotron radiation small angle x‐ray scattering (SAXS) technique was used to study the shape change of PeaT1 protein with temperature in this article. Ab initio structure restorations based on the SAXS data revealed that PeaT1 protein has a prolate shape with a P₂ symmetry axis along the prolate anisometric direction. With temperature increase, a gooseneck vase‐like (25°C), to jug‐like (55°C), then to oval (85°C) shape change can be found, and these shape changes are also approximately reversible with temperature decrease. PeaT1 protein contains two homogenous molecules, and each of them consists of F, NAC, T, and UBA domains. The structures of the four domains were predicted. Simulated annealing algorithm was used to superimpose the domain structures onto the SAXS shapes. It was found that all the structural domains have position rotation and translation with temperature change, but the NAC domains are relatively stable, playing a role of frame. This shape change information provides clues for further exploring its biological function and application. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

Structural characterizations of the chloroplast translocon protein Tic110

Tic110 is a major component of the chloroplast protein import translocon. Two functions with mutually exclusive structures have been proposed for Tic110: a protein‐conducting channel with six transmembrane domains and a scaffold with two N‐terminal transmembrane domains followed by a large soluble domain for binding transit peptides and other stromal translocon components. To investigate the structure of Tic110, Tic110 from Cyanidioschyzon merolae (CmTic110) was characterized. We constructed three fragments, CmTic110_A, CmTic110_B and CmTic110_C, with increasing N‐terminal truncations, to perform small‐angle X‐ray scattering (SAXS) and X‐ray crystallography analyses and Dali structural comparison. Here we report the molecular envelope of CmTic110_B and CmTic110_C determined by SAXS, and the crystal structure of CmTic110_C at 4.2 Å. Our data indicate that the C‐terminal half of CmTic110 possesses a rod‐shaped helix‐repeat structure that is too flattened and elongated to be a channel. The structure is most similar to the HEAT‐repeat motif that functions as scaffolds for protein–protein interactions.     

A proteomic Ramachandran plot (PRplot)

Each protein structure can be characterized by the average values of the main chain torsion angles ? and ψ and, as a consequence, be plotted on a bidimensional diagram, which resembles the Ramachandran plot. Here, we describe a proteomic ?–ψ plot (PRplot) where each protein structure is associated with one point, allowing in this way to represent the entire protein structure universe. It was verified that the PRplot is a robust tool since it does not depend on the dimension of the proteins, on the crystallographic resolution of the structures, nor on the biological source; moreover, it is little affected by disordered and structurally uncharacterized residues. The proteins mapped on the PRplot tend to cluster in three regions that correspond to the structures rich in alpha-helices, in beta-strands, and in both helices and strands, and are distributed along a sigmoidal curve that connect these three highly populated regions. PRplots are a unique instrument to project all protein structures on a single bidimensional plane where the entire structural complexity is reduced to a striking simplicity, with the sigmoid curve clearly delineating the space fraction accessible to a stable protein.     

One-bond and two-bond J couplings help annotate protein secondary-structure motifs: J-coupling indexing applied to human endoplasmic reticulum protein ERp18

NMR coupling constants, both direct one‐bond (¹J) and geminal two‐bond (²J), are employed to analyze the protein secondary structure of human oxidized ERp18. Coupling constants collected and evaluated for the 18 kDa protein comprise 1268 values of ¹J_CαHα, ¹J_CαCβ, ¹J_CαC′, ¹J_C′N′, ¹J_N′Cα, ¹J_N′HN, ²J_CαN′, ²J_HNCα, ²J_C′HN, and ²J_HαC′. Comparison with ¹J and ²J data from reference proteins and pattern analysis on a per‐residue basis permitted main‐chain φ,ψ torsion‐angle combinations of many of the 149 amino‐acid residues in ERp18 to be narrowed to particular secondary‐structure motifs. J‐coupling indexing is here being developed on statistical criteria and used to devise a ternary grid for interpreting patterns of relative values of J. To account for the influence of the varying substituent pattern in different amino‐acid sidechains, a table of residue‐type specific threshold values was compiled for discriminating small, medium, and large categories of J. For the 15‐residue insertion that distinguishes the ERp18 fold from that of thioredoxin, the J‐coupling data hint at a succession of five isolated Type‐I β turns at progressively shorter sequence intervals, in agreement with the crystal structure. Proteins 2011. © 2010 Wiley‐Liss, Inc.     

The structure and NO binding properties of the nitrophorin‐like heme‐binding protein from Arabidopsis thaliana gene locus At1g79260.1

The protein from Arabidopsis thaliana gene locus At1g79260.1 is comprised of 166‐residues and is of previously unknown function. Initial structural studies by the Center for Eukaryotic Structural Genomics (CESG) suggested that this protein might bind heme, and consequently, the crystal structures of apo and heme‐bound forms were solved to near atomic resolution of 1.32 Å and 1.36 Å, respectively. The rate of hemin loss from the protein was measured to be 3.6 × 10^?5 s^?1, demonstrating that it binds heme specifically and with high affinity. The protein forms a compact 10‐stranded β‐barrel that is structurally similar to the lipocalins and fatty acid binding proteins (FABPs). One group of lipocalins, the nitrophorins (NP), are heme proteins involved in nitric oxide (NO) transport and show both sequence and structural similarity to the protein from At1g79260.1 and two human homologues, all of which contain a proximal histidine capable of coordinating a heme iron. Rapid‐mixing and laser photolysis techniques were used to determine the rate constants for carbon monoxide (CO) binding to the ferrous form of the protein (k′_CO = 0.23 μM^?1 s^?1, k_CO = 0.050 s^?1) and NO binding to the ferric form (k′_NO = 1.2 μM^–1 s^–1, k_NO = 73 s^?1). Based on both structural and functional similarity to the nitrophorins, we have named the protein nitrobindin and hypothesized that it plays a role in NO transport. However, one of the two human homologs of nitrobindin contains a THAP domain, implying a possible role in apoptosis. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

The influence of Mg2+ coordination on 13C and 15N chemical shifts in CKI1RD protein domain from experiment and molecular dynamics/density functional theory calculations

Sequence dependence of ¹³C and ¹⁵N chemical shifts in the receiver domain of CKI1 protein from Arabidopsis thaliana, CKI1_RD, and its complexed form, CKI1_RD?Mg²⁺, was studied by means of MD/DFT calculations. MD simulations of a 20–ns production run length were performed. Nine explicitly hydrated structures of increasing complexity were explored, up to a 40‐amino‐acid structure. The size of the model necessary depended on the type of nucleus, the type of amino acid and its sequence neighbors, other spatially close amino acids, and the orientation of amino acid NH groups and their surface/interior position. Using models covering a 10 and a 15 Å environment of Mg²⁺, a semi‐quantitative agreement has been obtained between experiment and theory for the V67?I73 sequence. The influence of Mg²⁺ binding was described better by the 15 Å as compared to the 10 Å model. Thirteen chemical shifts were analyzed in terms of the effect of Mg²⁺ insertion and geometry preparation. The effect of geometry was significant and opposite in sign to the effect of Mg²⁺ binding. The strongest individual effects were found for ¹⁵N of D70, S74, and V68, where the electrostatics dominated; for ¹³Cβ of D69 and ¹⁵N of K76, where the influences were equal, and for ¹³Cα of F72 and ¹³Cβ of K76, where the geometry adjustment dominated. A partial correlation between dominant geometry influence and torsion angle shifts upon the coordination has been observed. Proteins 2016; 84:686–699. © 2016 Wiley Periodicals, Inc.     

Experimental conformational energy maps of proteins and peptides

We have presented an extensive analysis of the peptide backbone dihedral angles in the PDB structures and computed experimental Ramachandran plots for their distributions seen under a various constraints on X‐ray resolution, representativeness at different sequence identity percentages, and hydrogen bonding distances. These experimental distributions have been converted into isoenergy contour plots using the approach employed previously by F. M. Pohl. This has led to the identification of energetically favored minima in the Ramachandran (? , ψ ) plots in which global minima are predominantly observed either in the right‐handed α‐helical or the polyproline II regions. Further, we have identified low energy pathways for transitions between various minima in the (? ,ψ ) plots. We have compared and presented the experimental plots with published theoretical plots obtained from both molecular mechanics and quantum mechanical approaches. In addition, we have developed and employed a root mean square deviation (RMSD) metric for isoenergy contours in various ranges, as a measure (in kcal.mol^?1) to compare any two plots and determine the extent of correlation and similarity between their isoenergy contours. In general, we observe a greater degree of compatibility with experimental plots for energy maps obtained from molecular mechanics methods compared to most quantum mechanical methods. The experimental energy plots we have investigated could be helpful in refining protein structures obtained from X‐ray, NMR, and electron microscopy and in refining force field parameters to enable simulations of peptide and protein structures that have higher degree of consistency with experiments. Proteins 2017; 85:979–1001. © 2017 Wiley Periodicals, Inc.     

Toward engineering efficient peptidomimetics. Screening conformational landscape of two modified dehydroaminoacids

Effective peptidomimetics should posses structural rigidity and appropriate interaction pattern leading to potential spatial and electronic matching to the target receptor site. Rational design of such small bioactive molecules could push chemical synthesis and molecular modeling toward faster progress in medicinal chemistry. Conformational properties of N‐t‐butoxycarbonyl‐glycine‐(E/Z)‐dehydrophenylalanine N′,N′‐dimethylamides (Boc‐Gly‐(E/Z)‐ΔPhe‐NMe₂) in chloroform were studied by NMR and IR spectroscopy. The experimental findings were supported by extensive calculations at DFT(B3LYP, M06‐2X) and MP2 levels of theory and the β‐turn tendency for both isomers of the studied dipeptide were determined in vacuum and in solution. The theoretical data and experimental IR results were used as an additional information for the NMR‐based determination of the detailed solution conformations of the peptides. The obtained results reveal that N‐methylation of C‐terminal amide group changes dramatically the conformational properties of studied dehydropeptides. Theoretical conformational analysis reveals that the tendency to adopt β‐turn conformations is much weaker for the N‐methylated Z isomer (Boc‐Gly‐(Z)‐ΔPhe‐NMe₂), both in vacuum and in polar environment. On the contrary, N‐methylated E isomer (Boc‐Gly‐(E)‐ΔPhe‐NMe₂) can easily adopt β‐turn conformation, but the backbone torsion angles (φ₁, ψ₁, φ₂, ψ₂) are off the limits for common β‐turn types. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 28–40, 2014.     

The dataset for protein–RNA binding affinity

We have developed a non‐redundant protein–RNA binding benchmark dataset derived from the available protein–RNA structures in the Protein Database Bank. It consists of 73 complexes with measured binding affinity. The experimental conditions (pH and temperature) for binding affinity measurements are also listed in our dataset. This binding affinity dataset can be used to compare and develop protein–RNA scoring functions. The predicted binding free energy of the 73 complexes from three available scoring functions for protein–RNA docking has a low correlation with the binding Gibbs free energy calculated from K_d. © 2013 The Protein Society     

On interpretation of protein X‐ray structures: Planarity of the peptide unit

Pauling's mastery of peptide stereochemistry—based on small molecule crystal structures and the theory of chemical bonding—led to his realization that the peptide unit is planar and then to the Pauling–Corey–Branson model of the α‐helix. Similarly, contemporary protein structure refinement is based on experimentally determined diffraction data together with stereochemical restraints. However, even an X‐ray structure at ultra‐high resolution is still an under‐determined model in which the linkage among refinement parameters is complex. Consequently, restrictions imposed on any given parameter can affect the entire structure. Here, we examine recent studies of high resolution protein X‐ray structures, where substantial distortions of the peptide plane are found to be commonplace. Planarity is assessed by the ω‐angle, a dihedral angle determined by the peptide bond (C? N) and its flanking covalent neighbors; for an ideally planar trans peptide, ω = 180°. By using a freely available refinement package, Phenix [Afonine et al. (2012) Acta Cryst. D, 68:352–367], we demonstrate that tightening default restrictions on the ω‐angle can significantly reduce apparent deviations from peptide unit planarity without consequent reduction in reported evaluation metrics (e.g., R‐factors). To be clear, our result does not show that substantial non‐planarity is absent, only that an equivalent alternative model is possible. Resolving this disparity will ultimately require improved understanding of the deformation energy. Meanwhile, we urge inclusion of ω‐angle statistics in new structure reports in order to focus critical attention on the usual practice of assigning default values to ω‐angle constraints during structure refinement. Proteins 2015; 83:1687–1692. © 2015 Wiley Periodicals, Inc.     

Bridging of partially negative atoms by hydrogen bonds from main‐chain NH groups in proteins: The crown motif

The backbone NH groups of proteins can form N1N3‐bridges to δ‐ve or anionic acceptor atoms when the tripeptide in which they occur orients them appropriately, as in the RL and LR nest motifs, which have dihedral angles 1,2‐α_Rα_L and 1,2‐α_Lα_R, respectively. We searched a protein database for structures with backbone N1N3‐bridging to anionic atoms of the polypeptide chain and found that RL and LR nests together accounted for 92% of examples found (88% RL nests, 4% LR nests). Almost all the remaining 8% of N1N3‐bridges were found within a third tripeptide motif which has not been described previously. We term this a “crown,” because of the disposition of the tripeptide CO groups relative to the three NH groups and the acceptor oxygen anion, and the crown together with its bridged anion we term a “crown bridge.” At position 2 of these structures the dihedral angles have a tight α_R distribution, but at position 1 they have a wider distribution, with ? and ψ values generally being lower than those at position 1. Over half of crown bridges involve the backbone CO group three residues N‐terminal to the tripeptide, the remainder being to other main‐chain or side‐chain carbonyl groups. As with nests, bridging of crowns to oxygen atoms within ligands was observed, as was bridging to the sulfur atom of an iron‐sulfur cluster. This latter property may be of significance for protein evolution. Proteins 2015; 83:2067–2076. © 2015 Wiley Periodicals, Inc.