Evolutionarily evolved discriminators in the 3-TPR domain of the Toc64 family involved in protein translocation at the outer membrane of chloroplasts and mitochondria

Transport of polypeptides across membranes is a general and essential cellular process utilised by molecular machines. At least one component of these complexes contains a domain composed of three tetratricopeptide repeat (3-TPR) motifs. We have focussed on the receptor Toc64 to elucidate the evolved functional specifications of its 3-TPR domain. Toc64 is a component of the Toc core complex and functionally replaces Tom70 at the outer membrane of mitochondria in plants. Its 3-TPR domain recognises the conserved C-terminus of precursor-bound chaperones. We built homology models of the 3-TPR domain of chloroplastic Toc64 from different species and of the mitochondrial isoform from Arabidopsis. Guided by modelling, we identified residues essential for functional discrimination of the differently located isoforms to be located almost exclusively on the convex surface of the 3-TPR domain. The only exception is at568Ser/ps557Met, which is positioned in the ligand-binding groove. The functional implications of the homology models are discussed. Figure Motion contained within the 2nd eigenvector of the Calpha covariance matrix of the 3-TPR domain of atToc64-V indicated by a porcupine plot Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Folding mechanism of an ankyrin repeat protein: scaffold and active site formation of human CDK inhibitor p19(INK4d)

The p19(INK4d) protein consists of five ankyrin repeats (ANK) and controls the human cell cycle by inhibiting the cyclin D-dependent kinases (CDK) 4 and 6. We investigated the folding of p19(INK4d) by urea-induced unfolding transitions, kinetic analyses of unfolding and refolding, including double-mixing experiments and a special assay for folding intermediates. Folding is a sequential two-step reaction via a hyperfluorescent on-pathway intermediate. This intermediate is present under all conditions, during unfolding, refolding and at equilibrium. The folding mechanism was confirmed by a quantitative global fit of a consistent set of equilibrium and kinetic data revealing the thermodynamics and intrinsic folding rates of the different states. Surprisingly, the N<-->I transition is much faster compared to the I<-->U transition. The urea-dependence of the intrinsic folding rates causes population of the intermediate at equilibrium close to the transition midpoint. NMR detected hydrogen/deuterium exchange and the analysis of truncated variants showed that the C-terminal repeats ANK3-5 are already folded in the on-pathway intermediate, whereas the N-terminal repeats 1 and 2 are not folded. We suggest that during refolding, repeats ANK3-ANK5 first form the scaffold for the subsequent assembly of repeats ANK1 and ANK2. The binding function of p19(INK4d) resides in the latter repeats. We propose that the graded stability and the facile unfolding of repeats 1 and 2 is a prerequisite for the down-regulation of the inhibitory activity of p19(INK4d) during the cell-cycle.     

Sequential unfolding of ankyrin repeats in tumor suppressor p16

The ANK repeat is a ubiquitous 33-residue motif that adopts a beta hairpin helix-loop-helix fold. Multiple tandem repeats stack in a linear manner to produce an elongated structure that is stabilized predominantly by short-range interactions between residues close in sequence. The tumor suppressor p16(INK4) consists of four repeats and represents the minimal ANK folding unit. We found from Phi value analysis that p16 unfolded sequentially. The two N-terminal ANK repeats, which are distorted from the canonical ANK structure in all INK4 proteins and which are important for functional specificity, were mainly unstructured in the rate-limiting transition state for folding/unfolding, while the two C-terminal repeats were fully formed. A sequential unfolding mechanism could have implications for the cellular fate of wild-type and cancer-associated mutant p16 proteins.     

Ab initio calculations on the thermodynamic properties of azaborospiropentanes

Following our recent studies of the thermodynamic properties of azaspiropentane and borospiropentane, in consideration of their usefulness as new potential high energy materials, we follow up with ab initio calculations on the thermodynamic properties of azaborospiropentanes. Properties reported in this study include optimized structural parameters, vibrational frequencies, enthalpies of formation, specific enthalpies of combustion, proton affinities, and hydride affinities. Our results indicate that azatriborospiropentane gives off most energy when combusted, as evidenced by its specific enthalpy of combustion of about −52 kJ per gram. Figure Optimized geometry for R-azatriborospiropentane (10) Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Application of the PM6 method to modeling the solid state

The applicability of the recently developed PM6 method for modeling various properties of a wide range of organic and inorganic crystalline solids has been investigated. Although the geometries of most systems examined were reproduced with good accuracy, severe errors were found in the predicted structures of a small number of solids. The origin of these errors was investigated, and a strategy for improving the method proposed. Figure Detail of Structure of Dihydrogen Phosphate in KH₂PO₄ (upper pair) and in (CH₃)₄NH₂PO₄. (Footnote): X-ray structures on left, PM6 structure on right. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Aromaticity and electronic properties of Heterosuperbenzene (Heterohexabenzocoronene)

A global electrophilicity parameter and the aromaticity of some heterocyclic polyaromatic hydrocarbons were evaluated on the basis of DFT calculations. The substitution of carbon atoms by nitrogen atoms dramatically changes the global electrophilicity of the molecules, with the fully substituted molecule being the most electrophilic with a reactivity very close to that of fullerene. Figure Fully substituted heterohexabenzocoronene Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Three-dimensional quantitative structure-activity relationship (3D-QSAR) studies of various benzodiazepine analogues of γ-secretase inhibitors

A 3D QSAR analysis has been performed on a series of 67 benzodiazepine analogues reported as γ-secretase inhibitors using molecular field analysis (MFA), with G/PLS to predict steric and electrostatic molecular field interaction for the activity. The MFA study was carried out using a training set of 54 compounds. The predictive ability of model developed was assessed using a test set of 13 compounds ( as high as 0.729). The analyzed MFA model has demonstrated a good fit, having r² value of 0.858 and cross validated coefficient, value as 0.790. The analysis of the best MFA model provided insight into possible modification of the molecules for better activity.   Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structural basis for recognition of 2',5'-linked oligoadenylates by human ribonuclease L

An interferon-induced endoribonuclease, ribonuclease L (RNase L), is implicated in both the molecular mechanism of action of interferon and the fundamental control of RNA stability in mammalian cells. RNase L is catalytically active only after binding to an unusual activator molecule containing a 5'-phosphorylated 2',5'-linked oligoadenylate (2-5A), in the N-terminal half. Here, we report the crystal structure of the N-terminal ankyrin repeat domain (ANK) of human RNase L complexed with the activator 2-5A. This is the first structural view of an ankyrin repeat structure directly interacting with a nucleic acid, rather than with a protein. The ANK domain folds into eight ankyrin repeat elements and forms an extended curved structure with a concave surface. The 2-5A molecule is accommodated at a concave site and directly interacts with ankyrin repeats 2-4. Interestingly, two structurally equivalent 2-5A binding motifs are found at repeats 2 and 4. The structural basis for 2-5A recognition by ANK is essential for designing stable 2-5As with a high likelihood of activating RNase L.     

Effects of tryptophan residue fluorination on streptavidin stability and biotin–streptavidin interactions via molecular dynamics simulations

Due to its highly specific and very strong binding, the (strept)avidin–biotin system forms the basis for numerous applications in the life sciences: immunoassays, DNA detection systems, affinity chromatography, etc. Fine-tuning of the ligand binding abilities of this system might provide new technologies with relevance to nanoscale research. Here, we report our computational investigations on wild type (WT) and modified streptavidin (SAV), assessing the impact of fluorination of tryptophan residues on biotin binding ability. Complexes of biotin with four SAV protein variants (WT-SAV, 4fW-SAV, 5fW-SAV and 6fW-SAV) were studied. We found that protein stability and folding are predicted to be weakly affected by fluorination. The host protein binding pocket decreases its ability to form numerous hydrogen bonds to biotin in the case of the 4fW-SAV variant. Conversely, the 5fW-SAV mutant is predicted to have an even more stable ligand–host hydrogen bonding network than WT-SAV. Thermodynamic perturbation investigations predict a decrease in biotin binding free energy from 3.0 to 6.5 kcal/mol per tetrameric host, with the 5fW-SAV mutant being least affected. Overall, the computational findings indicate that 6fW-SAV and, especially, 5fW-SAV to be promising variants of streptavidin for potential modifiable picomolar binding of the biotin ligand family. Figure Hydrogen bonding framework of the biotin–streptavidin system Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A minimum folding unit in the ankyrin repeat protein p16(INK4)

The ankyrin repeat is an abundant, 33 residue sequence motif that forms a consecutive beta-hairpin-helix-loop-helix (beta(2)alpha(2)) fold. Most ankyrin repeat proteins consist of four or more complete repeats, which provide stabilizing interactions between adjacent modules. The cyclin-dependent kinase inhibitor and tumor suppressor p16(INK4) (p16) is one of the smallest ankyrin repeat proteins with a known structure. It consists of four complete repeats plus short N and C-terminal flanking regions that are unstructured in solution. On the basis of preliminary proteolysis studies and predictions using a computer algorithm for identifying autonomous folding units, we have identified a fragment consisting of the third and fourth ankyrin repeats of p16, called p16C, that can fold independently, without the rest of the protein. Far-UV circular dichroism studies showed that p16C has a significant level of alpha-helical secondary structure, and two proline substitutions that disrupt the alpha-helical secondary structure in wild-type p16 disrupt the secondary structure in p16C. The thermal denaturation of p16C is cooperative and reversible, with a midpoint of transition at 30. 5(+/-1) degrees C. From urea-induced denaturation studies, the free energy of unfolding for p16C was estimated to be 1.7(+/-0.3) kcal/mol at 20 degrees C. (1)H-(15)N 2D NMR studies suggest that the ankyrin repeats in p16C are likely to fold into a structure similar to that of full-length p16. In order to define the minimum autonomous folding unit in p16, we have further dissected p16C into two complementary peptides, each containing a single ankyrin repeat. These peptides are unstructured in solution. Thus, p16C is the smallest ankyrin repeat module that is known to fold independently and, in general, we believe that the two-ankyrin repeat fold could be the minimum structural unit for all ankyrin repeat proteins. We further discuss the significance of p16C in protein folding and engineering.     

Solution structure of the human oncogenic protein gankyrin containing seven ankyrin repeats and analysis of its structure--function relationship

Human gankyrin (226 residues, 24.4 kDa) is a liver oncoprotein that plays an important role in the development of human hepatocellular carcinomas. In this paper, its solution structure is reported, which is the largest ankyrin protein ever determined by NMR. The highly degenerate primary sequences of the seven ankyrin repeats presented a major challenge, which was overcome by combined use of TROSY experiments, perdeuterated samples, isotope-filtered NMR experiments, and residual dipolar couplings. The final structure was of high quality, with atomic rmsds for the backbone (N, C', and C(alpha)) and all heavy atoms (residues 4-224) of 0.69 +/- 0.09 and 1.04 +/- 0.09 A, respectively. Detailed analyses of NMR data suggested that the conserved TPLH motifs play important structural roles in stabilizing the repeating ankyrin scaffold. Gankyrin is conformationally more stable than the tumor suppressor p16(INK4A), possibly due to the structural roles of conserved residues evidenced by slowly exchanging backbone amides as well as hydrogen bonding networks involving labile side chain protons. Structural comparison with p16(INK4A) identified several residues of gankyrin that are potentially important for CDK4 binding, whereas observation of the thiol proton of C180 indicated a well-structured Rb-binding site in the helical region of the sixth ankyrin repeat. Interestingly, the CDK4-binding site and Rb-binding site located in N- and C-terminal regions, respectively, are separated by comparatively more stable ankyrin repeats and highly condensed positive surface charge. These results and analyses will shed light on the structural basis of the function of human gankyrin.     

Pharmacophore mapping of a series of pyrrolopyrimidines,indolopyrimidines and their congeners as multidrug-resistance-associated protein (MRP1) modulators

Pharmacophore mapping studies were undertaken for a series of molecules belonging to pyrrolopyrimidines, indolopyrimidines and their congeners as multidrug resistance-associated protein (MRP1) modulators. A five-point pharmacophore with two hydrogen bond acceptors (A), one lipophilic/hydrophobic group (H), one positive ionic feature (P) and one aromatic ring (R) as pharmacophoric features was developed. The pharmacophore hypothesis yielded a statistically significant 3D-QSAR model, with a correlation coefficient of r ² = 0.799 for training set molecules. The model generated showed excellent predictive power, with a correlation coefficient Q ² = 0.679 for an external test set of 20 molecules. The pharmacophore was further validated using four structurally diverse compounds with MRP1 modulatory activity. These compounds mapped well onto four of the five features of the pharmacophore. The pharmacophore proposed here was then utilised for the successful retrieval of active molecules with diverse chemotypes from database search. The geometry and features of pharmacophore are expected to be useful for the design of selective MRP1 inhibitors. Figure Alignment of multidrug resistance-associated protein (MRP1) inhibitors with the developed pharmacophore. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structure-based substitutions for increased solubility of a designed protein

Manipulation of protein solubility is important for many aspects of protein design and engineering. Previously, we designed a series of consensus ankyrin repeat proteins containing one, two, three and four identical repeats (1ANK, 2ANK, 3ANK and 4ANK). These proteins, particularly 4ANK, are intended for use as a universal scaffold on which specific binding sites can be constructed. Despite being well folded and extremely stable, 4ANK is soluble only under acidic conditions. Designing interactions with naturally occurring proteins requires the designed protein to be soluble at physiological pH. Substitution of six leucines with arginine on exposed hydrophobic patches on the surface of 4ANK resulted in increased solubility over a large pH range. Study of the pH dependence of stability demonstrated that 4ANK is one of the most stable ankyrin repeat proteins known. In addition, analogous leucine to arginine substitutions on the surface of 2ANK allowed the partially folded protein to assume a fully folded conformation. Our studies indicate that replacement of surface-exposed hydrophobic residues with positively charged residues can significantly improve protein solubility at physiological pH.     

Investigation of the intermolecular proton transfer in the supersystems adenine-methanol/ethanol/i-propanol: MP2 and DFT levels study

Twelve H-bonded supersystems constructed between the adenine tautomers and methanol, ethanol, and i-propanol were studied at the B3LYP and MP2 levels of theory using 6-311G(d,p) and 6-311++G(d,p) basis functions. The thermodynamic parameters of the complex formations were calculated in order to estimate the exact stability of the supersystems. It was proven that the calculated energy barriers of the alcohol-assisted proton transfers are about 60% lower than those of the intramolecular proton transfers in adenine found earlier (Gu and Leszczynski in J Phys Chem A 103:2744–2750, 1999). Figure H-bonded complex between i-propanol and adenine     

Ab initio multireference study of Hetero-Diels-Alder reaction of buta-1,3-diene with alkyl glyoxylates

Hetero-Diels-Alder (HDA) reaction of methyl glyoxylate with buta-1,3-diene has been investigated using multireference methods (complete active space SCF and multi-reference perturbation theory) and compared with several single-reference methods (including DFT) often used in calculations of catalysed [4+2] cycloadditions. Concerted and stepwise mechanisms, found in the literature, are compared. It is shown, that the stepwise mechanism may be a result of choosing unbalanced active space. Such choice leads to very close singlet and triplet states in the intermediate geometry - an artificial effect, that disappears if properly balanced active space is used (here, we use active space of 12 orbitals and 12 electrons). Conclusions concerning the mechanism and usefulness of the applied methodology are drawn, which might be important for theoretical investigation of stereoselectivity and specificity of catalysts for the HDA reaction. Figure Hetero-Diels-Alder reaction of alkyl glyoxylates with buta-1,3-diene, investigated using multi- and single-reference ab initio methods Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Ligation of Aza bases to the AgF<Subscript>2</Subscript> molecule: a theoretical study

We have analyzed the electronic structure and chemical bonding for molecular adducts of the Ag(II)F₂ molecule with various aza Lewis bases including ammonia, nitriles, secondary amines, and their derivatives exhibiting various degrees of fluorination. Density functional theory calculations indicate that a progressive shift occurs of the spin density from the Ag center towards the coordinating nitrogen atoms of aza ligands, as the ligation energy increases. Chemistry of Ag(II) might be extended with little effort beyond the known aza connections, to include nitriles, perfluorinated nitriles and perfluorinated amines. Figure Properties of a variety of novel adducts of the AgF₂ molecule with two aza bases (L), possible precursors of the AgF₂L₂ extended solids, were assessed by the DFT calculations Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. This work is dedicated to memory of Wojciech Ochmański, unforgettable person, good-hearted man, whose craftsmanship in work was second-to-none.     

Theoretical study of spin-spin coupling across the hydrogen (O-H⋯N) bond in adenosine derivatives

The study of spin-spin coupling constants across hydrogen bond provides useful information about configuration of complexes. The interesting case of such interactions was observed as a coupling across an intramolecular hydrogen bond in 8-bromo-2′,3′-O-isopropylideneadenosine between the -CH₂OH (at 5″ proton) group and the nitrogen atom of adenine. In this paper we report theoretical investigations on the ^4h J _NH coupling across the H″-C-O-H···N hydrogen bond in adenosine derivatives in various solvent models. Figure Coupling constants in 8-bromo-2′,3′-O-isopropylideneadenosine Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Structural model of an antistasin/notch-like fusion protein from the cocoon wall of the aquatic leech, Theromyzon tessulatum

The aquatic leech, Theromyzon tessulatum, secretes a proteinaceous cocoon with extraordinary physical properties (e.g., proteolytic, thermal resiliency). The deduced amino acid sequence of a major protein (Tcp—Theromyzon cocoon protein) from the T. tessulatum cocoon wall has been used to model the endogenous structure of the Tcp protein. The Tcp protein sequence comprises six internal repeats, each containing 12 ordered Cys residues. Amino acid alignments suggest that the region Cys1→6 is homologous to antistasin, a leech anticoagulant, and Cys7→12 is homologous to an epidermal growth factor-like domain found in notch-class proteins, which play critical roles in development, signaling, and adhesion throughout the Animalia. Modeling of individual domains (i.e., antistasin and notch) positions multiple hydrophobic and charged residues on the surface. When the antistasin and notch domains were fused, hydrophobic pockets appeared that may facilitate a polymerization mechanism. Collectively, the predicted features of our Tcp model are consistent with the physical properties of the leech cocoon wall.