Divergent evolution of a bifunctional de novo protein

Primordial proteins, the evolutionary ancestors of modern sequences, are presumed to have been minimally active and nonspecific. Following eons of selective pressure, these early progenitors evolved into highly active and specific proteins. While evolutionary trajectories from poorly active and multifunctional generalists toward highly active specialists likely occurred many times in evolutionary history, such pathways are difficult to reconstruct in natural systems, where primordial sequences are lost to time. To test the hypothesis that selection for enhanced activity leads to a loss of promiscuity, we evolved a de novo designed bifunctional protein. The parental protein, denoted Syn‐IF, was chosen from a library of binary patterned 4‐helix bundles. Syn‐IF was shown previously to rescue two different auxotrophic strains of E. coli: ΔilvA and Δfes. These two strains contain deletions for proteins with very different biochemical functions; IlvA is involved in isoleucine biosynthesis, while Fes is involved in iron assimilation. In two separate experiments, Syn‐IF, was evolved for faster rescue of either ΔilvA or Δfes. Following multiple rounds of mutagenesis, two new proteins were selected, each capable of rescuing the selected function significantly faster than the parental protein. In each case, the evolved protein also lost the ability to rescue the unselected function. In both evolutionary trajectories, the original bifunctional generalist was evolved into a monofunctional specialist with enhanced activity.     

Membrane‐induced structure of novel human tachykinin hemokinin‐1 (hHK1)

PPT‐C encoded hemokinin‐1(hHK‐1) of Homo sapiens (TGKASQFFGLM) is a structurally distinct neuropeptide among the tachykinin family that participate in the NK‐1 receptor downstream signaling processes. Subsequently, signal transduction leads to execution of various effector functions which includes aging, immunological, and central nervous system (CNS) regulatory actions. However the conformational pattern of ligand receptor binding is unclear. The three‐dimensional structure of the hemokinin‐1 in aqueous and micellar environment has been studied by one and two‐dimensional proton nuclear magnetic resonance (2D 1H‐NMR spectroscopy) and distance geometry calculations. Data shows that hemokinin‐1 was unstructured in aqueous environment; anionic detergent SDS induces α‐helix formation. Proton NMR assignments have been carried out with the aid of correlation spectroscopy (gradient‐COSY and TOCSY) and nuclear Overhauser effect spectroscopy (NOESY and ROESY) experiments. The inter proton distances and dihedral angle constraints obtained from the NMR data have been used in torsion angle dynamics algorithm for NMR applications (CYANA) to generate a family of structures, which have been refined using restrained energy minimization and dynamics. Helical conformation is observed from residue K3‐M11. The conformational range of the peptide revealed by NMR studies has been analyzed in terms of characteristic secondary features. Observed conformational features have been compared to that of Substance P potent NK1 agonist. Thus the report provides a structural insight to study hHK‐1‐NK1 interaction that is essential for hHK1 based signaling events. © 2015 Wiley Periodicals, Inc. Biopolymers 103: 702–710, 2015.     

Conformational regulation of the EGFR kinase core by the juxtamembrane and C-terminal tail: a molecular dynamics study

The catalytic domain of epidermal growth factor receptor (EGFR) is activated by dimerization, which requires allosteric coupling between distal dimerization and catalytic sites. Although crystal structures of EGFR kinases, solved in various conformational states, have provided important insights into EGFR activation by dimerization, the atomic details of how dimerization signals are dynamically coupled to catalytic regions of the kinase core are not fully understood. In this study, we have performed unrestrained and targeted molecular dynamics simulations on the active and inactive states of EGFR, followed by principal component analysis on the simulated trajectories, to identify correlated motions in the EGFR kinase domain upon dimerization. Our analysis reveals that the conformational changes associated with the catalytic functions of the kinase core are highly correlated with motions in the juxtamembrane (JM) and C-terminal tail, two flexible structural elements that play an active role in EGFR kinase activation and dimerization. In particular, the opening and closing of the ATP binding lobe relative to the substrate binding lobe is highly correlated with motions in the JM and C-terminal tail, suggesting that ATP and substrate binding can be coordinated with dimerization through conformational changes in the JM and C-terminal tail. Our study pinpoints key residues involved in this conformational coupling, and provides new insights into the role of the JM and C-terminal tail segments in EGFR kinase functions.     

Determination of protein fold class from Raman or Raman optical activity spectra using random forests

Knowledge of the fold class of a protein is valuable because fold class gives an indication of protein function and evolution. Fold class can be accurately determined from a crystal structure or NMR structure, though these methods are expensive, time-consuming, and inapplicable to all proteins. In contrast, vibrational spectra [infra-red, Raman, or Raman optical activity (ROA)] are rapidly obtained for proteins under wide range of biological molecules under diverse experimental and physiological conditions. Here, we show that the fold class of a protein can be determined from Raman or ROA spectra by converting a spectrum into data of 10 cm⁻¹ bin widths and applying the random forest machine learning algorithm. Spectral data from 605 and 1785 cm⁻¹ were analyzed, as well as the amide I, II, and III regions in isolation and in combination. ROA amide II and III data gave the best performance, with 33 of 44 proteins assigned to one of the correct four top-level structural classification of proteins (SCOP) fold class (all α, all β, α and β, and disordered). The method also shows which spectral regions are most valuable in assigning fold class.     

Activity optimization of an undecapeptide analogue derived from a frog-skin antimicrobial peptide

While natural antimicrobial peptides are potential therapeutic agents, their physicochemical properties and bioactivity generally need to be enhanced for clinical and commercial development. We have previously developed a cationic, amphipathic α-helical, 11-residue peptide (named herein GA-W2: FLGWLFKWASK-NH₂) with potent antimicrobial and hemolytic activity, which was derived from a 24-residue, natural antimicrobial peptide isolated from frog skin. Here, we attempted to optimize peptide bioactivity by a rational approach to sequence modification. Seven analogues were generated from GA-W2, and their activities were compared with that of a 12-residue peptide, omiganan, which is being developed for clinical and commercial applications. Most of the modifications reported here improved antimicrobial activity. Among them, the GA-K4AL (FAKWAFKWLKK-NH₂) peptide displayed the most potent antimicrobial activity with negligible hemolytic activity, superior to that of omiganan. The therapeutic index of GA-K4AL was improved more than 53- and more than 31-fold against Gram-negative and Gram-positive bacteria, respectively, compared to that of the starting peptide, GA-W2. Given its relatively shorter length and simpler amino acid composition, our sequence-optimized GA-K4AL peptide may thus be a potentially useful antimicrobial peptide agent.     

Organization and assembly of the TRAPPII complex

Current models suggest that TRAPP tethering complexes exist in two forms. Whereas the seven-subunit TRAPPI complex mediates ER-to-Golgi transport, TRAPPII contains three additional subunits (Trs65, Trs120 and Trs130) and is required for distinct tethering events at Golgi membranes. It is not clear how TRAPPII assembly is regulated. Here, we show that Tca17 is a fourth TRAPPII-specific component, and that Trs65 and Tca17 interact with distinct domains of Trs130 and make different contributions to complex assembly. Whereas Tca17 promotes the stable association of TRAPPII-specific subunits with the core complex, Trs65 stabilizes TRAPPII in an oligomeric form. We show that Trs85, which was previously reported to be a subunit of both TRAPPI and TRAPPII, is not associated with the TRAPPII complex in yeast. However, we find that proteins related to Trs85, Trs65 and Tca17 are part of the same TRAPP complex in mammalian cells. These findings have implications for models of TRAPP complex formation and suggest that TRAPP complexes may be organized differently in yeast and mammals.     

Determination of alpha-helix and beta-sheet stability in the solid state: a solid-state NMR investigation of poly(L-alanine)

The relative stability of alpha-helix and beta-sheet secondary structure in the solid state was investigated using poly(L-alanine) (PLA) as a model system. Protein folding and stability has been well studied in solution, but little is known about solid-state environments, such as the core of a folded protein, where peptide packing interactions are the dominant factor in determining structural stability. (13)C cross-polarization with magic angle spinning (CPMAS) NMR spectroscopy was used to determine the backbone conformation of solid powder samples of 15-kDa and 21.4-kDa PLA before and after various sample treatments. Reprecipitation from helix-inducing solvents traps the alpha-helical conformation of PLA, although the method of reprecipitation also affects the conformational distribution. Grinding converts the secondary structure of PLA to a final steady-state mixture of 55% beta-sheet and 45% alpha-helix at room temperature regardless of the initial secondary structure. Grinding PLA at liquid nitrogen temperatures leads to a similar steady-state mixture with 60% beta-sheet and 40% alpha-helix, indicating that mechanical shear force is sufficient to induce secondary structure interconversion. Cooling the sample in liquid nitrogen or subjecting it to high pressure has no effect on secondary structure. Heating the sample without grinding results in equilibration of secondary structure to 50% alpha-helix/50% beta-sheet at 100 degrees C when starting from a mostly alpha-helical state. No change was observed upon heating a beta-sheet sample, perhaps due to kinetic effects and the different heating rate used in the experiments. These results are consistent with beta-sheet approximately 260 J/mol more stable than alpha-helix in solid-state PLA.     

The cytoplasmic domain of MT1-MMP is dispensable for migration augmentation but necessary to mediate viability of MCF-7 breast cancer cells

Membrane-type-1 Matrix Metalloproteinase (MT1-MMP) is a multifunctional protease that regulates ECM degradation, proMMP-2 activation, and varied cellular processes including migration and viability. MT1-MMP is believed to be a central mediator of tumourigenesis whose role is dictated by its functionally distinct protein domains. Both the localization and signal transduction capabilities of MT1-MMP are dependent on its cytoplasmic domain, exemplifying diverse regulatory functions. To further our understanding of the multifunctional contributions of MT1-MMP to cellular processes, we overexpressed cytoplasmic domain altered constructs in MCF-7 breast cancer cells and analyzed migration and viability in 2D culture conditions, morphology in 3D Matrigel culture, and tumorigenic ability in vivo. We found that the cytoplasmic domain was not needed for MT1-MMP mediated migration promotion, but was necessary to maintain viability during serum depravation in 2D culture. Similarly, during 3D Matrigel culture the cytoplasmic domain of MT1-MMP was not needed to initiate a protrusive phenotype, but was necessary to prevent colony blebbing when cells were serum deprived. We also tested in vivo tumorigenic potential to show that cells expressing cytoplasmic domain altered constructs demonstrated a reduced ability to vascularize tumours. These results suggest that the cytoplasmic domain regulates MT1-MMP function in a manner required for cell survival, but is dispensable for cell migration.     

Crystal structure of the RNA 2'-phosphotransferase from Aeropyrum pernix K1

In the final step of tRNA splicing, the 2'-phosphotransferase catalyzes the transfer of the extra 2'-phosphate from the precursor-ligated tRNA to NAD. We have determined the crystal structure of the 2'-phosphotransferase protein from Aeropyrum pernix K1 at 2.8 Angstroms resolution. The structure of the 2'-phosphotransferase contains two globular domains (N and C-domains), which form a cleft in the center. The N-domain has the winged helix motif, a subfamily of the helix-turn-helix family, which is shared by many DNA-binding proteins. The C-domain of the 2'-phosphotransferase superimposes well on the NAD-binding fold of bacterial (diphtheria) toxins, which catalyze the transfer of ADP ribose from NAD to target proteins, indicating that the mode of NAD binding by the 2'-phosphotransferase could be similar to that of the bacterial toxins. The conserved basic residues are assembled at the periphery of the cleft and could participate in the enzyme contact with the sugar-phosphate backbones of tRNA. The modes by which the two functional domains recognize the two different substrates are clarified by the present crystal structure of the 2'-phosphotransferase.