On the evolutionary conservation of hydrogen bonds made by buried polar amino acids: the hidden joists,braces and trusses of protein architecture

Background  

The hydrogen bond patterns between mainchain atoms in protein structures not only give rise to regular secondary structures but also satisfy mainchain hydrogen bond potential. However, not all mainchain atoms can be satisfied through hydrogen bond interactions that arise in regular secondary structures; in some locations sidechain-to-mainchain hydrogen bonds are required to provide polar group satisfaction. Buried polar residues that are hydrogen-bonded to mainchain amide atoms tend to be highly conserved within protein families, confirming that mainchain architecture is a critical restraint on the evolution of proteins. We have investigated the stabilizing roles of buried polar sidechains on the backbones of protein structures by performing an analysis of solvent inaccessible residues that are entirely conserved within protein families and superfamilies and hydrogen bonded to an equivalent mainchain atom in each family member.     

D-Amino acids in protein de novo design. II. Protein-diastereomerism versus protein-enantiomerism

Summary With a few exceptions, proteins in our biosphere are based exclusively onl-amino acids. The inversion of configuration of all the stereogenic centers in a protein leads to anall-d compound with ‘mirror image’ properties and ‘mirror image’ structure. We propose to use the termprotein-enantiomerism to describe the relationship between two proteins that have the same sequence but whose amino acids have opposite configuration. We will use the termprotein-diastereomerism to define the relationship between two proteins that have the same sequence in which some amino acids have opposite configurations. A classification of type I, II, III, and IV protein-diastereomerism is proposed. By extension, a diastereoprotein is a protein where some amino acids have the same configuration (l ord) while others have the opposite one (d orl). A particular case of diastereoproteins aremesoproteins, also analyzed in this article. In addition to the goal of making proteins resistant to protease degradation, the use ofd-amino acids in protein de novo design may give rise to proteins with structures, and perhaps properties, very different to those of nativeall-l-proteins.     

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.     

Predicting the Structure of the Flavodoxin from Eschericia coli by Homology Modeling,Distance Geometry and Molecular Dynamics

Abstract

As part of our on-going development of a method, based upon distance geometry calculations, for predicting the structures of proteins from the known structures of their homologues, we have predicted the structure of the 176 residue Flavodoxin from Escherichia coli. This prediction was based upon the crystal structures of the homologous Flavodoxins from Anacystis nidulans, Chondrus crispus, Desulfovibrio vulgaris and Clostridium beijerinckii, whose sequence identities with Escherichia coli were 44%, 33%, 23% and 16%, respectively. A total of 13,043 distance constraints among the alpha-carbons of the Escherichia coli structure were derived from the sequence alignments with the known structures, together with 8,893 distance constraints among backbone and sidechain atoms of adjacent residues, 978 between the alpha-carbons and selected atoms of the flavin mononucleotide cofactor, 116 constraints to enforce conserved hydrogen bonds, and 452 constraints on the torsion angles in conserved residues. An ensemble of ten random Escherichia coli structures was computed from these constraints, with an average root mean square coordinate deviation (RMSD) among the alpha carbons of 0.85 Ångstroms (excluding the first 1 and last 6 residues, which have no corresponding residues in any of the homologues and hence were unconstrained); the corresponding average heavy-atom RMSD was 1.60 Å.

Since the distance geometry calculations were performed without hydrogen atoms, protons were added to the resulting structures and these structures embedded in a 50 × 50 × 40 Å solvent box with periodic boundary conditions. They were then subjected to a 20 picosecond dynamical simulated annealing procedure, starting at 300 K and gradually reduced to 10K, in which all the distance and torsion angle constraints were maintained by means of harmonic restraint functions. This was followed up by 1000 iterations of unrestrained conjugate gradients minimization. The goal of this energy refinement procedure was not to drastically modify the structures in an attempt at a priori prediction, but merely to improve upon the predictions obtained from the geometric constraints, particularly with regard to their local backbone and sidechain conformations and their hydrogen bonds. The resulting structures differed from the respective starting structures by an average of 1.52 Å in their heavy atom RMSD's, while the average RMSD among the heavy atoms of residues 2-170 increased slightly to 1.66 Å. We hope these structures will be good enough to enable the phase problem to be solved for the crystallographic data that is now being collected on this protein.     

Prediction of structures of zinc‐binding proteins through explicit modeling of metal coordination geometry

Metal ions play an essential role in stabilizing protein structures and contributing to protein function. Ions such as zinc have well‐defined coordination geometries, but it has not been easy to take advantage of this knowledge in protein structure prediction efforts. Here, we present a computational method to predict structures of zinc‐binding proteins given knowledge of the positions of zinc‐coordinating residues in the amino acid sequence. The method takes advantage of the “atom‐tree” representation of molecular systems and modular architecture of the Rosetta3 software suite to incorporate explicit metal ion coordination geometry into previously developed de novo prediction and loop modeling protocols. Zinc cofactors are tethered to their interacting residues based on coordination geometries observed in natural zinc‐binding proteins. The incorporation of explicit zinc atoms and their coordination geometry in both de novo structure prediction and loop modeling significantly improves sampling near the native conformation. The method can be readily extended to predict protein structures bound to other metal and/or small chemical cofactors with well‐defined coordination or ligation geometry.     

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.     

Reduced representation model of protein structure prediction: statistical potential and genetic algorithms.

A reduced representation model, which has been described in previous reports, was used to predict the folded structures of proteins from their primary sequences and random starting conformations. The molecular structure of each protein has been reduced to its backbone atoms (with ideal fixed bond lengths and valence angles) and each side chain approximated by a single virtual united-atom. The coordinate variables were the backbone dihedral angles phi and psi. A statistical potential function, which included local and nonlocal interactions and was computed from known protein structures, was used in the structure minimization. A novel approach, employing the concepts of genetic algorithms, has been developed to simultaneously optimize a population of conformations. With the information of primary sequence and the radius of gyration of the crystal structure only, and starting from randomly generated initial conformations, I have been able to fold melittin, a protein of 26 residues, with high computational convergence. The computed structures have a root mean square error of 1.66 A (distance matrix error = 0.99 A) on average to the crystal structure. Similar results for avian pancreatic polypeptide inhibitor, a protein of 36 residues, are obtained. Application of the method to apamin, an 18-residue polypeptide with two disulfide bonds, shows that it folds apamin to native-like conformations with the correct disulfide bonds formed.     

On the Relation of Attempting to Lose Weight,Restraint, and Binge Eating in Outpatients with Binge Eating Disorder

Objective: To examine the relationship among attempts to lose weight, restraint, and eating behavior in outpatients with binge eating disorder (BED). Research Methods and Procedures: Participants were 93 consecutive outpatients evaluated for a clinical trial who met Diagnostic and Statistical Manual, Fourth edition criteria for BED. The Eating Disorder Examination Interview was administered to assess attempts at weight loss, restraint, different forms of overeating, and the attitudinal psychopathology of eating disorders (i.e., concerns regarding eating, shape, and weight). In addition, the Three‐Factor Eating Questionnaire was used to assess cognitive restraint, hunger, and disinhibition. Psychometrically established measures were given to assess body dissatisfaction, depression, and self‐esteem. Results: The majority of participants (75.3%; N = 70) reported attempting to lose weight, but only 37.6% (N = 35) reported dietary restraint on at least half the days of the month. Dietary restraint and cognitive restraint were not associated with any form of binge eating or overeating. Dietary restraint and cognitive restraint were positively correlated with weight concern, shape concern, and body dissatisfaction, and negatively correlated with body mass index. To further examine the interplay between attempting to lose weight and restraint, three study groups were created: unrestrained nonattempters (21.5%, N = 20), unrestrained attempters (40.9%; N = 38), and restrained attempters (34.4%; N = 32). The three groups did not differ significantly on binge eating or other eating behaviors; however, significant differences were observed for weight concern, shape concern, and body dissatisfaction. Discussion: Attempts to lose weight and restraint are not synonymous for patients with BED. Although 75.3% of BED patients reported that they were attempting to lose weight, only 37.6% reported dietary restraint on at least half the days of the previous month. While restraint was negatively associated with body mass index, it was not related to binge eating or overeating. Our findings raise questions about prevailing models that posit restraint as a predominant factor in the maintenance of binge eating in BED.     

1H and 15N NMR assignment and solution structure of the SH3 domain of spectrin: Comparison of unrefined and refined structure sets with the crystal structure

The assignment of the ¹H and ¹⁵Nnuclear magnetic resonance spectra of the Src-homology region 3 domain ofchicken brain -spectrin has been obtained. A set of solutionstructures has been determined from distance and dihedral angle restraints,which provide a reasonable representation of the protein structure insolution, as evaluated by a principal component analysis of the globalpairwise root-mean-square deviation (rmsd) in a large set of structuresconsisting of the refined and unrefined solution structures and the crystalstructure. The solution structure is well defined, with a lower degree ofconvergence between the structures in the loop regions than in the secondarystructure elements. The average pairwise rmsd between the 15 refinedsolution structures is 0.71 ± 0.13 Å for the backbone atoms and1.43 ± 0.14 Å for all heavy atoms. The solution structure isbasically the same as the crystal structure. The average rmsd between the 15refined solution structures and the crystal structure is 0.76 Å forthe backbone atoms and 1.45 ± 0.09 Å for all heavy atoms. Thereare, however, small differences probably caused by intermolecular contactsin the crystal structure.     

Protein-protein docking with a reduced protein model accounting for side-chain flexibility

A protein-protein docking approach has been developed based on a reduced protein representation with up to three pseudo atoms per amino acid residue. Docking is performed by energy minimization in rotational and translational degrees of freedom. The reduced protein representation allows an efficient search for docking minima on the protein surfaces within. During docking, an effective energy function between pseudo atoms has been used based on amino acid size and physico-chemical character. Energy minimization of protein test complexes in the reduced representation results in geometries close to experiment with backbone root mean square deviations (RMSDs) of approximately 1 to 3 A for the mobile protein partner from the experimental geometry. For most test cases, the energy-minimized experimental structure scores among the top five energy minima in systematic docking studies when using both partners in their bound conformations. To account for side-chain conformational changes in case of using unbound protein conformations, a multicopy approach has been used to select the most favorable side-chain conformation during the docking process. The multicopy approach significantly improves the docking performance, using unbound (apo) binding partners without a significant increase in computer time. For most docking test systems using unbound partners, and without accounting for any information about the known binding geometry, a solution within approximately 2 to 3.5 A RMSD of the full mobile partner from the experimental geometry was found among the 40 top-scoring complexes. The approach could be extended to include protein loop flexibility, and might also be useful for docking of modeled protein structures.     

Improving the physical realism and structural accuracy of protein models by a two-step atomic-level energy minimization

Most protein structural prediction algorithms assemble structures as reduced models that represent amino acids by a reduced number of atoms to speed up the conformational search. Building accurate full-atom models from these reduced models is a necessary step toward a detailed function analysis. However, it is difficult to ensure that the atomic models retain the desired global topology while maintaining a sound local atomic geometry because the reduced models often have unphysical local distortions. To address this issue, we developed a new program, called ModRefiner, to construct and refine protein structures from Cα traces based on a two-step, atomic-level energy minimization. The main-chain structures are first constructed from initial Cα traces and the side-chain rotamers are then refined together with the backbone atoms with the use of a composite physics- and knowledge-based force field. We tested the method by performing an atomic structure refinement of 261 proteins with the initial models constructed from both ab initio and template-based structure assemblies. Compared with other state-of-art programs, ModRefiner shows improvements in both global and local structures, which have more accurate side-chain positions, better hydrogen-bonding networks, and fewer atomic overlaps. ModRefiner is freely available at http://zhanglab.ccmb.med.umich.edu/ModRefiner.     

Design,synthesis, cytotoxic evaluation,and QSAR study of some 6H-indolo[2,3-b]quinoxaline derivatives

In the pathway of anticancer drug development, we designed and synthesized some 6H-indolo[2,3-b]quinoxaline derivatives (which act as DNA intercalators) by structural modification. The structure of the 6H-indolo[2,3-b]quinoxaline derivatives was confirmed by IR, NMR, Mass and elemental analysis. The compounds (IDQ-5, IDQ-10, IDQ-11, IDQ-13, and IDQ-14) exhibited significant in vitro activity against a human leukemia (HL-60) cell line. The QSAR derived for modeling the cytotoxic activity of 6H-indolo[2,3-b]quinoxaline derivatives suggests that candidate structures for increased cytotoxic potency should incorporate cyclic substituents or substituents with primary carbon atoms.     

Stereochemical studies on cyclic peptides. IX. Conformational studies on cyclic tetrapeptides containing alternating cis and trans peptide units

Conformational analyses of cyclic tetrapeptides consisting of alternating cis and trans peptide units have been made using contact criteria and energy calculations. This study has been restricted to those structures having a symmetry element in the backbone ring, such as a twofold axis (d) or a center of inversion (i). There are five main results. (1) There are two distinct types of conformations, which are stereochemically favorable corresponding to each of twofold and inversion-symmetrical structures, designated as d₁, d₂ (for twofold symmetrical) and i₁, i₂ (for inversion-symmetrical). Among these, the i₁ type has the lowest energy when glycyl residues occur at all four α-carbon atoms. (2) With the glycyl residue at all four α-carbon atoms, methyl substitution at the cis peptide nitrogen atoms is possible in all the four types, whereas the substitution at trans peptide nitrogen atoms is possible only for the i₁ type. Thus only in the i₁ type can all the nitrogen atoms be methylated simultaneously. The conformation of the molecule in the crystal structure of cyclotetrasarcosyl belongs to the i₁ type. (3) When alanyl residues occur at all four α-carbon atoms, the possible symmetrical type is dependent on the enantiomorphic form and the actual sequence of the alanyl residues. (4) The methyl substitution at peptide nitrogen atoms for cyclic tetrapeptides having alanyl residues causes more stereochemical restriction in the allowed conformations than with glycyl residues. (5) The prolyl residue can be incorporated favorably at the cis-trans junction of both d and i types of structures. The results of the present study are compared with the data on cyclic tetrapeptides available from the crystal structure and nmr studies. The results show an overall agreement both regarding the type of symmetry and the conformational parameters.     

Mechanistic Studies of ABCR, the ABC Transporter in Photoreceptor Outer Segments Responsible for Autosomal Recessive Stargardt Disease

ABCR is an ABC transporter that is found exclusively in vertebrate photoreceptor outer segments. Mutations in the human ABCR gene are responsible for autosomal recessive Stargardt disease, the most common cause of early onset macular degeneration. In this paper we review our recent work with purified and reconstituted ABCR derived from bovine retina and from cultured cells expressing wild type or site-directed mutants of human ABCR. These experiments implicate all-trans-retinal (or Schiff base adducts between all-trans-retinal and phosphatidylethanolamine) as the transport substrate, and they reveal asymmetric roles for the two nucleotide binding domains in the transport reaction. A model for the retinal transport reaction is presented which accounts for these experimental observations.     

Condensed Representations of Protein Secondary Structure Sequences and Their Application

Abstract

In this paper, we propose a nongraphical representation for protein secondary structures. By counting the frequency of occurrence of all possible four-tuples (i.e., four-letter words) of a protein secondary structure sequence, we construct a set of 3 × 3 matrices for the corresponding protein secondary structure sequence. Furthermore, the leading eigenvalues of these matrices are computed and considered as invariants for the protein secondary structure sequences. To illustrate the utility of our approach, we apply it to a set of real data to distinguish protein structural classes. The result indicates that it can be used to complement the classification of protein secondary structures.     

Macrophage Recognition of Thiol-Group Oxidized Cells: Recognition of Carbohydrate Chains by Macrophage Surface Nucleolin as Apoptotic Cells

Three kinds of xylo-oligosaccharides having structures of 3²-β-xylosylxylobiose, 3²-β-xylobiosylxylobiose, and 2²-β-xylobiosylxylobiose were isolated from an enzymatic hydrolysate of hardwood xylan with Streptomyces β-xylanase. The structures suggest that the hardwood xylan has both (1 → 2)- and (1 → 3)-β-d-xylopyranosyl linkages in the structure, and the specificity of Streptomyces β-xylanase toward the stubs is similar to that toward glucuronic acid stubs, but is somewhat different from that toward arabinose and xylosylarabinose stubs.     

The Eating Inventory and Body Adiposity from Leanness to Massive Obesity: a Study of 2509 Adults

Objective: To examine how chronic dietary restraint, disinhibition, and hunger, as assessed by the Eating Inventory, vary over a broad range of BMI values from leanness to massive obesity, in subjects with family obesity. Eating Inventory factors were also studied as a function of personal weight history. Research Methods and Procedures: Subjects were 2509 participants in a genetic study of obesity. BMIs ranged from 15 to 87 kg/m². Six BMI groups were formed (<27, 27 to 30, 30 to 35, 35 to 40, 40 to 45, and >45). Results: Multivariate analyses showed that restraint and disinhibition were significantly associated with BMI in men, whereas only disinhibition was in women. Disinhibition scores correlated strongly with hunger scores in both genders in all BMI categories; dietary restraint tended to correlate with the other two factors positively in leaner subjects and negatively in the highest BMI categories. Highly restrained normal‐weight subjects were likely to exhibit disinhibition and hunger, whereas massively obese persons with very high disinhibition scores showed high hunger but little restraint. The highest restraint scores were observed in nonobese adult women with previous obesity in childhood and/or adolescence. Discussion: The factor most strongly associated with BMI in this large population was disinhibition, suggesting that obesity treatment should target behaviors associated with disinhibition, especially in individuals showing a low level of dietary restraint. High restraint scores in formerly obese normal‐weight women suggest that dietary restraint may exert a beneficial influence on body weight control under conditions that deserve further investigation.     

Can molecular dynamics simulations help in discriminating correct from erroneous protein 3D models?

Background  

Recent approaches for predicting the three-dimensional (3D) structure of proteins such asde novoor fold recognition methods mostly rely on simplified energy potential functions and a reduced representation of the polypeptide chain. These simplifications facilitate the exploration of the protein conformational space but do not permit to capture entirely the subtle relationship that exists between the amino acid sequence and its native structure. It has been proposed that physics-based energy functions together with techniques for sampling the conformational space, e.g., Monte Carlo or molecular dynamics (MD) simulations, are better suited to the task of modelling proteins at higher resolutions than those of models obtained with the former type of methods. In this study we monitor different protein structural properties along MD trajectories to discriminate correct from erroneous models. These models are based on the sequence-structure alignments provided by our fold recognition method, FROST. We define correct models as being built from alignments of sequences with structures similar to their native structures and erroneous models from alignments of sequences with structures unrelated to their native structures.