Synthesis and Structure-Activity Relationship of [Nle10]Neurokinin A (4–10) Analogs with Constraint in the Backbone and at Position Six

Steric requirements of binding [Nle¹⁰]NKA(4–10) to NK-2 receptor were studied by introducing conformationally constrained amino acid analogs into its sequence. Two series of [Nle¹⁰]NKA(4–10) analogs were synthesized to investigate (i) the significance of a putative β-turn in the receptor-ligand interaction by insertion of either (S)- or (R)-Gly⁸{ANC-2}Leu⁹ γ-lactams to mimic a β-turn constraint, and (ii) the effect of hindered rotation in the Φ, χ¹ and χ² dihedral angle space of the crucially important Phe⁶ which was replaced systematically with d-Phe, d- and l-Tyr, as well as with their conformationally constrained analogs, Tic, HOTic and β-MePhe. Competition binding experiments with [³H]NKA were performed using cloned human NK-2 receptors expressed in CHO cells. The analog possessing only an (R)-Gly⁸{ANC-2}Leu⁹ constraint, had the same binding affinity as that of the parent peptide. The rank order of potency of the other analogs showed a cumulative effect of different structural modifications in decreasing the binding affinity, i.e., when changing the configuration of the lactam ring to S, replacing Phe⁶ with constrained analogs, Tic or β-MePhe, changing the configuration of the amino acid at position six to d, and introducing a hydroxyl group on the aromatic ring. Ferenc ?tv?s and Dmitry S. Gembitsky - Made an equal contribution. Abbreviations used for amino acids and peptides follow the recommendations of the IUPAC-IUB Commission of Biochemical Nomenclature, Eur. J. Biochem. (1984) 138, 9–37     

Antiamoebin I in methanol solution: rapid exchange between right-handed and left-handed 3(10)-helical conformations

Antiamoebin I (Aam-I) is a membrane-active peptaibol antibiotic isolated from fungal species belonging to the genera Cephalosporium, Emericellopsis, Gliocladium, and Stilbella. Antiamoebin I has the amino acid sequence: Ac-Phe(1)-Aib-Aib-Aib-Iva-Gly-Leu-Aib(8)-Aib-Hyp-Gln-Iva-Hyp-Aib-Pro-Phl(16). By using the uniformly (13)C,(15)N-labeled sample of Aam-I, the set of conformationally dependent J couplings and (3h)J(NC) couplings through H-bonds were measured. Analysis of these data along with the data on magnetic nonequivalence of the (13)C(beta) nuclei (Deltadelta((13)C(beta))) in Aib and Iva residues allowed us to draw the univocal conclusion that the N-terminal part (Phe(1)-Gly(6)) of Aam-I in MeOH solution is in fast exchange between the right-handed and left-handed 3(10)-helical conformations, with an approximately equal population of both states. An additional conformational exchange process was found at the Aib(8) residue. The (15)N-NMR-relaxation and CD-spectroscopy measurements confirmed these findings. Molecular modeling and Monte Carlo simulations revealed that both exchange processes are correlated and coupled with significant hinge-bending motions around the Aib(8) residue. Our results explain relatively low activity of Aam-I with respect to other 15-amino acid residue peptaibols (for example, zervamicin) in functional and biological tests. The high dynamic 'propensity' possibly prevents both initial binding of the antiamoebin to the membrane and subsequent formation of stable ionic channels according to the barrel-stave mechanism.     

‘Birth defects’ of photosystem II make it highly susceptible to photodamage during chloroplast biogenesis

High solar flux is known to diminish photosynthetic growth rates, reducing biomass productivity and lowering disease tolerance. Photosystem II (PSII) of plants is susceptible to photodamage (also known as photoinactivation) in strong light, resulting in severe loss of water oxidation capacity and destruction of the water‐oxidizing complex (WOC). The repair of damaged PSIIs comes at a high energy cost and requires de novo biosynthesis of damaged PSII subunits, reassembly of the WOC inorganic cofactors and membrane remodeling. Employing membrane‐inlet mass spectrometry and O₂‐polarography under flashing light conditions, we demonstrate that newly synthesized PSII complexes are far more susceptible to photodamage than are mature PSII complexes. We examined these ‘PSII birth defects’ in barley seedlings and plastids (etiochloroplasts and chloroplasts) isolated at various times during de‐etiolation as chloroplast development begins and matures in synchronization with thylakoid membrane biogenesis and grana membrane formation. We show that the degree of PSII photodamage decreases simultaneously with biogenesis of the PSII turnover efficiency measured by O₂‐polarography, and with grana membrane stacking, as determined by electron microscopy. Our data from fluorescence, Q_B‐inhibitor binding, and thermoluminescence studies indicate that the decline of the high‐light susceptibility of PSII to photodamage is coincident with appearance of electron transfer capability Q_A^? → Q_B during de‐etiolation. This rate depends in turn on the downstream clearing of electrons upon buildup of the complete linear electron transfer chain and the formation of stacked grana membranes capable of longer‐range energy transfer.     

Human p38α mitogen-activated protein kinase in the Asp168-Phe169-Gly170-in (DFG-in) state can bind allosteric inhibitor Doramapimod

Doramapimod (BIRB-796) is widely recognized as one of the most potent and selective type II inhibitors of human p38α mitogen-activated protein kinase (MAPK); however, the understanding of its binding mechanism remains incomplete. Previous studies indicated high affinity of the ligand to a so-called allosteric pocket revealed only in the ‘out’ state of the DFG motif (i.e. Asp168-Phe169-Gly170) when Phe169 becomes fully exposed to the solvent. The possibility of alternative binding in the DFG-in state was hypothesized, but the molecular mechanism was not known. Methods of bioinformatics, docking and long-time scale classical and accelerated molecular dynamics have been applied to study the interaction of Doramapimod with the human p38α MAPK. It was shown that Doramapimod can bind to the protein even when the Phe169 is fully buried inside the allosteric pocket and the kinase activation loop is in the DFG-in state. Orientation of the inhibitor in such a complex is significantly different from that in the known crystallographic complex formed by the kinase in the DFG-out state; however, the Doramapimod’s binding is followed by the ligand-induced conformational changes, which finally improve accommodation of the inhibitor. Molecular modelling has confirmed that Doramapimod combines the features of type I and II inhibitors of p38α MAPK, i.e. can directly and indirectly compete with the ATP binding. It can be concluded that optimization of the initial binding in the DFG-in state and the final accommodation in the DFG-out state should be both considered at designing novel efficient type II inhibitors of MAPK and homologous proteins.

Communicated by Ramaswamy H. Sarma     

A phylogenetic study of krill (Crustacea: Euphausiacea) reveals new taxa and co-evolution of morphological characters

The first comprehensive phylogenetic study of Euphausiacea (all 86 valid species) is presented. It is based on four molecular markers and 168 morphological characters (including 58 characters of the petasma). Phylogenetic analyses support the monophyly and robustness of the families Bentheuphausidae and Euphausiidae and reveal three major clades for which we erect three new subfamilies: Thysanopodinae, Euphausiinae and Nematoscelinae. All genus-level clades are statistically supported (except Thysanopoda in molecular analyses), deeply nested within the subfamily-level clades, and encompass 14 new species groups. Copulatory structures have a major impact on tree topology in the morphological analysis, the removal of which resulted in only half the number of supported clades and genera. We revealed three groups of morphological characters, which are probably coupled with the same biological role and thus interlinked evolutionarily: (i) antennular peduncle and petasma (copulation); (ii) eyes and anterior thoracopods (feeding); and (iii) shape of carapace and pleon (defence). We analysed the evolutionary pathways of the clades into main oceanic biotopes and compared them with morphological adaptations most likely to be coupled with this process.     

Rapid directed molecular evolution of fluorescent proteins in mammalian cells

In vivo imaging of model organisms is heavily reliant on fluorescent proteins with high intracellular brightness. Here we describe a practical method for rapid optimization of fluorescent proteins via directed molecular evolution in cultured mammalian cells. Using this method, we were able to perform screening of large gene libraries containing up to 2 × 10⁷ independent random genes of fluorescent proteins expressed in HEK cells, completing one iteration of directed evolution in a course of 8 days. We employed this approach to develop a set of green and near‐infrared fluorescent proteins with enhanced intracellular brightness. The developed near‐infrared fluorescent proteins demonstrated high performance for fluorescent labeling of neurons in culture and in vivo in model organisms such as Caenorhabditis elegans, Drosophila, zebrafish, and mice. Spectral properties of the optimized near‐infrared fluorescent proteins enabled crosstalk‐free multicolor imaging in combination with common green and red fluorescent proteins, as well as dual‐color near‐infrared fluorescence imaging. The described method has a great potential to be adopted by protein engineers due to its simplicity and practicality. We also believe that the new enhanced fluorescent proteins will find wide application for in vivo multicolor imaging of small model organisms.     

Pairs of Mutually Compensatory Frameshifting Mutations Contribute to Protein Evolution

Insertions and deletions of lengths not divisible by 3 in protein-coding sequences cause frameshifts that usually induce premature stop codons and may carry a high fitness cost. However, this cost can be partially offset by a second compensatory indel restoring the reading frame. The role of such pairs of compensatory frameshifting mutations (pCFMs) in evolution has not been studied systematically. Here, we use whole-genome alignments of protein-coding genes of 100 vertebrate species, and of 122 insect species, studying the prevalence of pCFMs in their divergence. We detect a total of 624 candidate pCFM genes; six of them pass stringent quality filtering, including three human genes: RAB36, ARHGAP6, and NCR3LG1. In some instances, amino acid substitutions closely predating or following pCFMs restored the biochemical similarity of the frameshifted segment to the ancestral amino acid sequence, possibly reducing or negating the fitness cost of the pCFM. Typically, however, the biochemical similarity of the frameshifted sequence to the ancestral one was not higher than the similarity of a random sequence of a protein-coding gene to its frameshifted version, indicating that pCFMs can uncover radically novel regions of protein space. In total, pCFMs represent an appreciable and previously overlooked source of novel variation in amino acid sequences.     

Thermodynamic, kinetic and structural basis for recognition and repair of abasic sites in DNA by apurinic/apyrimidinic endonuclease from human placenta

X-ray analysis of enzyme–DNA interactions is very informative in revealing molecular contacts, but provides neither quantitative estimates of the relative importance of these contacts nor information on the relative contributions of specific and nonspecific interactions to the total affinity of enzymes for specific DNA. A stepwise increase in the ligand complexity approach is used to estimate the relative contributions of virtually every nucleotide unit of synthetic DNA containing abasic sites to its affinity for apurinic/apyrimidinic endonuclease (APE1) from human placenta. It was found that APE1 interacts with 9–10 nt units or base pairs of single-stranded and double-stranded ribooligonucleotides and deoxyribooligonucleotides of different lengths and sequences, mainly through weak additive contacts with internucleotide phosphate groups. Such nonspecific interactions of APE1 with nearly every nucleotide within its DNA-binding cleft provides up to seven orders of magnitude (ΔG° ~ −8.7 to −9.0 kcal/mol) of the enzyme affinity for any DNA substrate. In contrast, interactions with the abasic site together with other specific APE1–DNA interactions provide only one order of magnitude (ΔG° ~ −1.1 to −1.5 kcal/mol) of the total affinity of APE1 for specific DNA. We conclude that the enzyme's specificity for abasic sites in DNA is mostly due to a great increase (six to seven orders of magnitude) in the reaction rate with specific DNA, with formation of the Michaelis complex contributing to the substrate preference only marginally.     

A gradient of silent substitution rate in the human pseudoautosomal region

It has been demonstrated that recombination in the human p-arm pseudoautosomal region (p-PAR) is at least twenty times more frequent than the genomic average of approximately 1 cM/Mb, which may affect substitution patterns and rates in this region. Here I report the analysis of substitution patterns and rates in 10 human, chimpanzee, gorilla, and orangutan genes across the p-PAR. Between species silent divergence in the p-PAR forms a gradient, increasing toward the telomere. The correlation of silent divergence with distance from the p-PAR boundary is highly significant (rho = 0.911, P < 0.001). After exclusion of the CpG dinucleotides this correlation is still significant (rho = 0.89, P < 0.01), thus the substitution rate gradient cannot be explained solely by the differences in the extent of methylation across the p-PAR. Frequent recombination in the PAR may result in a relatively strong effect of biased gene conversion (BGC), which, because of the increased probability of fixation of the G or C nucleotides at (A or T)/(G or C) segregating sites, may affect substitution rates. BGC, however, does not seem to be the factor creating the substitution rate gradient in the p-PAR, because the only gradient is still detactable if only A<-->T and G<-->C substitutions are taken into account (rho = 0.82, P < 0.01). I hypothesize that the substitution rate gradient in the p-PAR is due to the mutagenic effect of recombination, which is very frequent in the distal human p-PAR and might be lower near the p-PAR boundary.     

On the correct use of ‘data’

Acknowledgement  Thanks are due to Drs. Paul Johnston (P.Johnston@ exeter.ac.uk) and E. William Beese (ebeese@t-online.de) for kindly providing dictionary definitions and for sharing their thoughts about the correct use of ‘data’.