An improved capping unit for stabilizing the ends of associated β-strands

Understanding protein beta structures has been hindered by the challenge of designing small, well-folded β-sheet systems. A β-capping motif was previously designed to help solve this problem, but not without limitations, as the termini of this β-cap were not fully available for chain extension. Combining Coulombic side chain attractions with a Trp/Trp edge-to-face interaction we produced a new capping motif that provided greater β-sheet stability. This stability was maintained even in systems lacking a turn locus with a high propensity for chain direction reversal. The Coulombic cap was shown to improve β-sheet stability in a number of difficult systems, hence providing an additional tool for protein structure and folding studies.     

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.     

High Molecular Weight Forms of Mammalian Respiratory Chain Complex II

Mitochondrial respiratory chain is organised into supramolecular structures that can be preserved in mild detergent solubilisates and resolved by native electrophoretic systems. Supercomplexes of respiratory complexes I, III and IV as well as multimeric forms of ATP synthase are well established. However, the involvement of complex II, linking respiratory chain with tricarboxylic acid cycle, in mitochondrial supercomplexes is questionable. Here we show that digitonin-solubilised complex II quantitatively forms high molecular weight structures (CII_hmw) that can be resolved by clear native electrophoresis. CII_hmw structures are enzymatically active and differ in electrophoretic mobility between tissues (500 – over 1000 kDa) and cultured cells (400–670 kDa). While their formation is unaffected by isolated defects in other respiratory chain complexes, they are destabilised in mtDNA-depleted, rho0 cells. Molecular interactions responsible for the assembly of CII_hmw are rather weak with the complexes being more stable in tissues than in cultured cells. While electrophoretic studies and immunoprecipitation experiments of CII_hmw do not indicate specific interactions with the respiratory chain complexes I, III or IV or enzymes of the tricarboxylic acid cycle, they point out to a specific interaction between CII and ATP synthase.     

Nanotubules formed by highly hydrophobic amphiphilic alpha-helical peptides and natural phospholipids

We previously reported that the 18-mer amphiphilic alpha-helical peptide, Hel 13-5, consisting of 13 hydrophobic residues and five hydrophilic amino acid residues, can induce neutral liposomes (egg yolk phosphatidylcholine) to adopt long nanotubular structures and that the interaction of specific peptides with specific phospholipid mixtures induces the formation of membrane structures resembling cellular organelles such as the Golgi apparatus. In the present study we focused our attention on the effects of peptide sequence and chain length on the nanotubule formation occurring in mixture systems of Hel 13-5 and various neutral and acidic lipid species by means of turbidity measurements, dynamic light scattering measurements, and electron microscopy. We designed and synthesized two sets of Hel 13-5 related peptides: 1) Five peptides to examine the role of hydrophobic or hydrophilic residues in amphiphilic alpha-helical structures, and 2) Six peptides to examine the role of peptide length, having even number residues from 12 to 24. Conformational, solution, and morphological studies showed that the amphiphilic alpha-helical structure and the peptide chain length (especially 18 amino acid residues) are critical determinants of very long tubular structures. A mixture of alpha-helix and beta-structures determines the tubular shapes and assemblies. However, we found that the charged Lys residues comprising the hydrophilic regions of amphiphilic structures can be replaced by Arg or Glu residues without a loss of tubular structures. This suggests that the mechanism of microtubule formation does not involve the charge interaction. The immersion of the hydrophobic part of the amphiphilic peptides into liposomes initially forms elliptic-like structures due to the fusion of small liposomes, which is followed by a transformation into tubular structures of various sizes and shapes.     

Spectroscopic studies and PM5 semiempirical calculations of new Schiff bases of gossypol with polyoxaalkylamines

Three Schiff bases of racemic gossypol with polyoxaalkylamines were synthesized and studied by FTIR and (1)H-NMR spectroscopy, and their structures were calculated by the PM5 semiempirical method. These Schiff bases exist in the solid state and in solutions as enamine forms. An increasing length of the polyoxaalkyl chain causes the increase of the interaction of the oxygen atoms of this chain with the OH groups in the 6,6' positions. This interaction is very well evidenced in the FTIR and (1)H-NMR spectra. The structures of the Schiff bases and the hydrogen bonds within these structures are discussed.     

Structural basis of fullerene derivatives as novel potent inhibitors of protein acetylcholinesterase without catalytic active site interaction: insight into the inhibitory mechanism through molecular modeling studies

Abstract

Acetylcholinesterase (AChE) is an important kind of esterase that plays a key biological role in the central and peripheral nervous systems. Recent research has demonstrated that some fullerene derivatives serve as a new nanoscale class of potent inhibitors of AChE, but the specific inhibition mechanism remains unclear. In the present article, several molecular modeling methods, such as molecular docking, molecular dynamics simulations and molecular mechanics/generalized Born surface area calculations, were used for the investigation of the binding mode and inhibition mechanism of fullerene inhibitions for AChE. Results revealed that fullerene inhibitors block the entrance of substrates by binding with the peripheral anionic site (PAS) region. Thus, fullerene derivatives might mainly serve as competitive inhibitors. The interactions of a fullerene backbone with AChE are at the same level in different single side chain systems and seem to be related to the length or aromaticity of the side chain. The inhibitor with multihydroxyl side chains shows an obviously large electrostatic interaction as it forms additional hydrogen bonds with AChE. Moreover, fullerene derivatives might exhibit noncompetitive inhibition partly by affecting some secondary structures around them. Thus, the destructions of these secondary structures can lead to conformational changes in some important regions, such as the catalytic triad and acyl pocket. The investigation is of great importance to the discovery of good fullerene inhibitors.

Communicated by Ramaswamy H. Sarma     

Persistence in food webs—I Lotka-Volterra food chains

Persistence-extinction in simple food chains modelled by Lotka-Volterra dynamics is governed by a single parameter which depends upon the interspecific interaction coefficients, the intraspecific interaction coefficients, and the length of the food chain. In persistent systems with nonzero carrying capacity, two new features predominate. Trophic level influence factors relate persistence on different trophic levels and determine, in conjunction with the persistence parameter, the magnitude of persistence. Equilibrium component ordering, which results in persistent systems, mandates once again that systems need to be studied on the complete ecosystem level; static field measurements reflect species location in the food chain, the total length of the food chain and assume characteristics according to these factors.     

EmPLiCS: an empirical approach for structure-based design of natural peptide drugs

The computer implementation of a peptide drug-design strategy has been developed. The system is named EmPLiCS (Empirical Peptide Ligand Construction System) according to the strategy of the system, which searches for peptide-ligand structures by referring to empirical rules that are derived from known protein 3D structures. The system was tested on several known peptide-protein complexes. The results demonstrated the ability of this system to detect key residues of peptides that are crucial for interaction with their specific proteins. The system also showed the ability to detect the main chain trace of these peptides. Some of the main chain atoms were detected even though the complete primary structures were not reproduced, suggesting that main chain structure is important in peptide-protein recognition. The results of the present study demonstrated that the empirical rules-based system can generate significant information for use in the design of natural peptide drugs.     

An algorithm for predicting protein-protein interaction sites: Abnormally exposed amino acid residues and secondary structure elements

Multiprotein systems mediate most regulatory processes in living organisms. Although the structures of the individual proteins are often defined, less is known of the structures of multiprotein systems. Computational methods for predicting interfaces, using evolutionary conservation and/or physicochemical data, have been developed. Here we consider the use of solvent accessibility, residue propensity, and hydrophobicity, in conjunction with secondary structure data, as prediction parameters. We analyze the influence of residue type and secondary structure on solvent accessibility and define a measure of "relative exposedness." Clustering abnormally high scoring residues provides a basis for predicting interaction sites. The analysis is extended to investigate abnormally exposed secondary structure elements, particularly beta-sheet strands. We show that surface-exposed beta-strands lacking protective features are more likely to be found at protein-protein interfaces, allowing us to create an algorithm with approximately 68% and approximately 75% accuracy in differentiating between interacting and edge strands in isolated beta-strands and beta-sheet strands, respectively. These methods of identifying abnormally exposed surface regions are combined in an algorithm, which, on a data set of 77 unbound and disjoint (single chain extracted from complex) structures, predicts 79% of the protein-protein interfaces correctly. If enzyme-inhibitor complexes, where the inhibitor mimics a nonprotein substrate, are excluded, the accuracy increases to 85%.     

Electron transfer and energy transduction in the terminal part of the respiratory chain — Lessons from bacterial model systems

This review focuses on the terminal part of the respiratory chain where, macroscopically speaking, electron transfer (ET) switches from the two-electron donor, ubiquinol, to the single-electron carrier, cytochrome c, to finally reduce the four-electron acceptor dioxygen. With 3-D structures of prominent representatives of such multi-subunit membrane complexes known for some time, this section of the ET chain still leaves a number of key questions unanswered. The two relevant enzymes, ubiquinol:cytochrome c oxidoreductase and cytochrome c oxidase, appear as rather diverse modules, differing largely in their design for substrate interaction, internal ET, and moreover, in their mechanisms of energy transduction. While the canonical mitochondrial complexes have been investigated for almost five decades, the corresponding bacterial enzymes have been established only recently as attractive model systems to address basic reactions in ET and energy transduction. Lacking the intricate coding background and mitochondrial assembly pathways, bacterial respiratory enzymes typically offer a much simpler subunit composition, while maintaining all fundamental functions established for their complex “relatives”. Moreover, related issues ranging from primary steps in cofactor insertion to supramolecular architecture of ET complexes, can also be favourably addressed in prokaryotic systems to hone our views on prototypic structures and mechanisms common to all family members.     

Conformational correlation and coupled motion between residue A21 and B25 side chain observed in crystal structures of insulin mutants at position A21

The C-terminal residue of the insulin A chain is invariant and kept as asparagine in all known insulin molecules from hagfish through birds to mammals. To get information on the role of this conserved residue, which is still unclear, the three-dimensional structures of four human insulin mutants, A21 Asn-->Gly, A21 Asn-->Asp, A21 Asn-->Ala, and A21 Asn-->Gln DesB30, were determined by X-ray crystallography. The four mutants crystallize separately into two kinds (rhombohedral and cubic) of crystals. In the refined structures, conformational correlation and coupled motion between the A chain C-terminal residue A21 and the B25 side chain was observed, in contrast to the nearly unchanged general structures as compared with the native insulin structures in their respective crystals. A detailed analysis suggests that residue A21 can affect insulin receptor binding by interaction with the B25 side chain and the B chain C-terminal segment to assist the B25 side chain rearranging into the 'active' conformation.     

Distance dependent centroid to centroid force fields using high resolution decoys

Simplified force fields play an important role in protein structure prediction and de novo protein design by requiring less computational effort than detailed atomistic potentials. A side chain centroid based, distance dependent pairwise interaction potential has been developed. A linear programming based formulation was used in which non-native "decoy" conformers are forced to take a higher energy compared with the corresponding native structure. This model was trained on an enhanced and diverse protein set. High quality decoy structures were generated for approximately 1400 nonhomologous proteins using torsion angle dynamics along with restricted variations of the hydrophobic cores of the native structure. The resulting decoy set was used to train the model yielding two different side chain centroid based force fields that differ in the way distance dependence has been used to calculate energy parameters. These force fields were tested on an independent set of 148 test proteins with 500 decoy structures for each protein. The side chain centroid force fields were successful in correctly identifying approximately 86% native structures. The Z-scores produced by the proposed centroid-centroid distance dependent force fields improved compared with other distance dependent C(alpha)-C(alpha) or side chain based force fields.     

Folding and dimerization of the ionic peptide EAK 16-IV

Folding and dimerization of an ionic polyalanine-based peptide chain (EAK16-IV) are simulated with nonspecific interactions. It is found that there is a competition between two kinds of structural motifs under different strengths of electrostatic interactions. The dominance of hairpin-like structures would be realized with a strong electrostatic interaction both thermodynamically and kinetically, showing the importance of the electrostatic interaction on the formation of hairpin-like structures. Simulations on the dimerization with strong electrostatic interaction are also carried out. It is found that the concentration contributes essentially to the shape of the dimers. These studies demonstrate that the strong interactions and kinetic factors might be important for the ordered amyloid aggregates.     

Finite one-dimensional spin systems as models of biopolymers

New models are proposed for describing various properties of biopolymers, especially those of proteins and nucleic acids. Each model is constituted of a set of spins arranged on a chain, and each pair of spins produces an interaction. We examine the transitions of these spin systems between the ground state and the disordered state. It is found that the transitions of the present spin systems demonstrate various properties in response to values of the so-called interaction energy. If we define interaction energy parameters with no so-called frustration, the system exhibits two-state transitions, similar to the folding-unfolding transition of small proteins. The addition of frustrations to the model produces effects similar to those of mutations in proteins. On the other hand, if the interactions between two spins attenuate as a function of their separation along the chain, the transition of the system has characteristics similar to those of nucleic acids. Thus, the present spin systems can offer a unified view of the folding-unfolding transition of biopolymers in terms of differences in the pairwise interactions between spins. Based on our models, we propose a condition for two-state transition behavior for proteins.     

Self Assembly of Bilayers in a Lattice Model of Amphiphile and Solvent Systems

The results of a Metropolis Monte Carlo simulation of a three dimensional lattice model of an amphiphile and solvent mixture are presented. In the model each amphiphile molecule is represented as a connected chain of lattice sites with one site representing the head and the remaining chain sites representing the tail of the molecule. The remaining sites on the lattice represent solvent molecules. The amphiphiles interact through a nearest neighbour potential which includes head-solvent and tail-solvent interactions. No prior assumption is made about the structures which may be observed.

The cluster size distribution and cluster structure is studied as a function of temperature and head-solvent interaction. If the tail-solvent interaction is solvophobic and the head-solvent interaction is solvophilic, a micellar region is observed in the phase diagram. For sufficiently solvophilic head-solvent interactions the low temperature phase exhibits self assembly into castellated bilayer structures.     

Structures of heparin-derived tetrasaccharide bound to cobra cardiotoxins: heparin binding at a single protein site with diverse side chain interactions

Cobra cardiotoxins (CTXs) are three-fingered polypeptides with positively charged domains that have been shown to bind to anionic ligands of snake venom citrate, glycosaminoglycans, sulfoglycosphingolipid, and nucleotide triphosphate with various biochemical effects including toxin dimerization, cell surface retention, membrane pore formation, cell internalization and blocking of enzymatic activities of kinase and ATPase. The reported anionic binding sites, however, are found to be different among different CTX homologues for potentially different CTX activities. Herein, by NMR studies of the binding of inorganic phosphate, dATP (stable form of ATP), and heparin-derived tetrasaccharide to Naja atra CTX A1, a novel CTX molecule exhibiting in vivo necrotic activity on skeletal muscle, we demonstrate that diverse ligands binding to CTXs could also occur at a single protein site with flexible side chain interactions. The flexibility of such an interaction is also illustrated by the available heparin-CTX A3 complex structures with different heparin chain lengths binding at the same site. Our results provide a likely structural explanation on how the interaction between heparan sufate and proteins depends more on the overall charge cluster organization rather than on their fine structures. We also suggest that the ligand binding site of CTX homologues can be fine-tuned by nonconserved residues near the binding pocket because of their flexible side chain interaction and dimerization ability, even for the rigid CTX molecules tightened by four disulfide bonds.     

A possible prebiotic peptide formation from glycinamide and related compounds

We examined an experimental approach to the genesis of protocells in the primeval sea. Glycine polymers with an average chain length of 12 were formed from glycinamide in fluctuating systems (pH 7.2, 80 degrees C, 20 cycles). The resulting glycine polymers gave aggregated leaflet-like structures. A solution of the glycine polymers provided stacked disc-shaped structures in the presence of LiBr and gave sheet structures in the presence of dichloroacetic acid. The shapes of these organized structures were correlated with their molecular structures.     

A possible prebiotic peptide formation from glycinamide and related compounds

We examined an experimental approach to the genesis of protocells in the primeval sea. Glycine polymers with an average chain length of 12 were formed from glycinamide in fluctuating systems (pH 7.2, 80°C, 20 cycles). The resulting glycine polymers gave aggregated leaflet-like structures. A solution of the glycine polymers provided stacked disc-shaped structures in the presence of LiBr and gave sheet structures in the presence of dichloroacetic acid. The shapes of these organized structures were correlated with their molecular structures.     

Two polymorphs of a covalent complex between papain and a diazomethylketone inhibitor.

The three-dimensional structure of two polymorphs of a ZLFG-CH2-papain covalent complex has been determined by X-ray crystallography. The structures indicate that: (i) the methylene carbon atom of the inhibitor is covalently bound to the Sgamma atom of Cys25 of papain; (ii) the hydrophobic S2 pocket formed by Pro68, Val133, Val157, and Asp158 is occupied by the inhibitor's phenylalanyl P2 side chain; (iii) extensive hydrogen bonding and hydrophobic interactions are responsible for the interaction of the inhibitor with the enzyme. Comparison with similar structures suggests that in covalent complexes preservation of main chain-main chain interactions between the enzyme and the inhibitor may have higher priority than the P-S interactions.