Measures of residue density in protein structures

Proteomics: An In-Depth Review on Recent Technical Advances and Their Applications in Biomedicine

Proteins hold pivotal importance since many diseases manifest changes in protein activity. Proteomics techniques provide a comprehensive exploration of protein structure, abundance, and function in biological samples, enabling the holistic characterization of overall changes in organisms. Nowadays, the breadth of emerging methodologies in proteomics is unprecedentedly vast, with constant optimization of technologies in sample processing, data collection, data analysis, and its scope of application is steadily transitioning from the bench to the clinic. Here, we offer an insightful review of the technical developments in proteomics and its applications in biomedicine over the past 5 years. We focus on its profound contributions in profiling disease spectra, discovering new biomarkers, identifying promising drug targets, deciphering alterations in protein conformation, and unearthing protein-protein interactions. Moreover, we summarize the cutting-edge technologies and potential breakthroughs in the proteomics pipeline and provide the principal challenges in proteomics. Based on these, we aspire to broaden the applicability of proteomics and inspire researchers to enhance our understanding of complex biological systems by utilizing such techniques.

Multiple, simultaneous, independent gradients for a versatile multidimensional liquid chromatography. Part II: Application 3 - Scouting optimization strategies for separation of monoclonal antibodies by dual simultaneous independent gradients of pH & salt on a weak cation exchange stationary phase

In previous publications we have described the pISep dual simultaneous, independent gradients (DSIGs) liquid chromatography (LC) for uncoupling gradients of non-buffering solute (NaCl, urea or acetonitrile) from externally generated pH gradients. In DSIGs the shape and slope of the [salute] gradient does not depend on the shape and slope of the pH gradient. The technique allows in a single run true simultaneous two dimensional LC separation of complex protein mixtures on various stationary phases including anion, cation exchangers (AEX, CEX), reversed phase (RP), mixed mode and mixed bed. Using a humanized IgG1 (HIgG1) monoclonal antibody (MAb) and a variety of pH & [NaCl] DSIGs, we show that most of MAb isoforms can be successfully separated from each other. These experimental observations are supported by an initial theoretical argument presented here predicting an overall improvement of all MAb isoforms separation by DSIGs of pH & [NaCl]. Theoretical calculations predict that, in general, there exists an optimal non-zero isocratic salt concentration in a pH gradient separation that will resolve isoforms close in binding energy, but a wide range of salt concentrations will be required for acceptable resolution of all isoforms. Theory also predicts better separation of weaker rather than stronger binding isoforms. Experimentally, we have found that no one set of DSIGs LC conditions could optimally baseline resolve all identifiable MAb isoforms in a single run of reasonable duration. The versatility and simplicity of the pH & [NaCl] pISep DSIGs LC allows fast, automated scouting of protein separations over any range of pH from 2.4 to 10.8 and [NaCl] from 0 to 1 M without changing the chemistry of the buffering system. Due to the universal applicability of the pISep buffering system in IEX LC, the researcher is given a powerful tool to easily develop pH & [NaCl] DSIGs protocols that vary mobile phase compositions to achieve high resolution separations of targeted proteins.

Quantitative Structural Proteomics Unveils the Conformational Changes of Proteins under the Endoplasmic Reticulum Stress

Protein structures are decisive for their activities and interactions with other molecules. Global analysis of protein structures and conformational changes cannot be achieved by commonly used abundance-based proteomics. Here, we integrated cysteine covalent labeling, selective enrichment, and quantitative proteomics to study protein structures and structural changes on a large scale. This method was applied to globally investigate protein structures in HEK293T cells and protein structural changes in the cells with the tunicamycin (Tm)-induced endoplasmic reticulum (ER) stress. We quantified several thousand cysteine residues, which contain unprecedented and valuable information of protein structures. Combining this method with pulsed stable isotope labeling by amino acids in cell culture, we further analyzed the folding state differences between pre-existing and newly synthesized proteins in cells under the Tm treatment. Besides newly synthesized proteins, unexpectedly, many pre-existing proteins were found to become unfolded upon ER stress, especially those related to gene transcription and protein translation. Furthermore, the current results reveal that N-glycosylation plays a more important role in the folding process of the tertiary and quaternary structures than the secondary structures for newly synthesized proteins. Considering the importance of cysteine in protein structures, this method can be extensively applied in the biological and biomedical research fields.

Training data composition affects performance of protein structure analysis algorithms

The three-dimensional structures of proteins are crucial for understanding their molecular mechanisms and interactions. Machine learning algorithms that are able to learn accurate representations of protein structures are therefore poised to play a key role in protein engineering and drug development. The accuracy of such models in deployment is directly influenced by training data quality. The use of different experimental methods for protein structure determination may introduce bias into the training data. In this work, we evaluate the magnitude of this effect across three distinct tasks: estimation of model accuracy, protein sequence design, and catalytic residue prediction. Most protein structures are derived from X-ray crystallography, nuclear magnetic resonance (NMR), or cryo-electron microscopy (cryo-EM); we trained each model on datasets consisting of either all three structure types or of only X-ray data. We Find that across these tasks, models consistently perform worse on test sets derived from NMR and cryo-EM than they do on test sets of structures derived from X-ray crystallography, but that the difference can be mitigated when NMR and cryo-EM structures are included in the training set. Importantly, we show that including all three types of structures in the training set does not degrade test performance on X-ray structures, and in some cases even increases it. Finally, we examine the relationship between model performance and the biophysical properties of each method, and recommend that the biochemistry of the task of interest should be considered when composing training sets.

Correlation between the secondary structure and surface activity of β-sheet forming cationic amphiphilic peptides and their anticancer activity

Cancer is one of the main causes of death worldwide. The current cancer treatment strategies often lack selectivity for cancer cells resulting in dose-limiting adverse effects and reduced quality of life. Recently, anticancer peptides (ACPs) have emerged as an alternative treatment with higher selectivity, less adverse effects, and lower propensity for drug resistance. However, most of the current studies on the ACPs are focused on α-helical ACPs and there is lack of systematic studies on β-sheet forming ACPs. Herein we report the development of a new series of rationally designed short cationic amphiphilic β-sheet forming ACPs and their structure activity relationship. The peptides had the general formula (XY1XY2)3, with X representing hydrophobic amino acids (isoleucine (I) or leucine (L)), Y1 and Y2 representing cationic amino acids (arginine (R) or lysine (K)). The cytotoxicity of the designed ACPs in HCT 116 colorectal cancer, HeLa cervical cancer and human dermal fibroblast (HDF) cells was assessed by MTT test. The physicochemical properties of the peptides were characterized by various techniques including RP-HPLC, LC-MS, and Circular Dichroism (CD) spectroscopy. The surface activity of the peptides at the air-water interface and their interaction with the lipid monolayers as models for cell membranes were studied by Langmuir trough. The peptides consisting of I with R and K had selective anticancer activity while the combination of L and R diminished the anticancer activity of the peptides but rendered them more toxic to HDFs. The anticancer activity of the peptides was directed by their surface activity (amphiphilicity) and their secondary structure in hydrophobic surfaces including cancer cell membranes. The selectivity of the peptides for cancer cells was a result of their higher penetration into cancer cell membranes compared to normal cell membranes. The peptides exerted their anticancer activity by disrupting the mitochondrial membranes and eventually apoptosis. The results presented in this study provide an insight into the structure-activity relationship of this class of ACPs which can be employed as guidance to design new ACPs with improved anticancer activity and lower toxicity against normal cells.

Macromolecular crowding effects on electrostatic binding affinity: Fundamental insights from theoretical, idealized models

The crowded cellular environment can affect biomolecular binding energetics, with specific effects depending on the properties of the binding partners and the local environment. Often, crowding effects on binding are studied on particular complexes, which provide system-specific insights but may not provide comprehensive trends or a generalized framework to better understand how crowding affects energetics involved in molecular recognition. Here, we use theoretical, idealized molecules whose physical properties can be systematically varied along with samplings of crowder placements to understand how electrostatic binding energetics are altered through crowding and how these effects depend on the charge distribution, shape, and size of the binding partners or crowders. We focus on electrostatic binding energetics using a continuum electrostatic framework to understand effects due to depletion of a polar, aqueous solvent in a crowded environment. We find that crowding effects can depend predictably on a system's charge distribution, with coupling between the crowder size and the geometry of the partners' binding interface in determining crowder effects. We also explore the effect of crowder charge on binding interactions as a function of the monopoles of the system components. Finally, we find that modeling crowding via a lowered solvent dielectric constant cannot account for certain electrostatic crowding effects due to the finite size, shape, or placement of system components. This study, which comprehensively examines solvent depletion effects due to crowding, complements work focusing on other crowding aspects to help build a holistic understanding of environmental impacts on molecular recognition.

Fine-Tuning of Alkaline Residues on the Hydrophilic Face Provides a Non-toxic Cationic α-Helical Antimicrobial Peptide Against Antibiotic-Resistant ESKAPE Pathogens

Antibiotic-resistant ESKAPE pathogens (Enterococcus faecium, Staphylococcus aureus, Klebsiella pneumoniae, Acinetobacter baumannii, Pseudomonas aeruginosa, and Enterobacter species) has become a serious threat to public health worldwide. Cationic α-helical antimicrobial peptides (CαAMPs) have attracted much attention as promising solutions in post-antibiotic era. However, strong hemolytic activity and in vivo inefficacy have hindered their pharmaceutical development. Here, we attempt to address these obstacles by investigating BmKn2 and BmKn2-7, two scorpion-derived CαAMPs with the same hydrophobic face and a distinct hydrophilic face. Through structural comparison, mutant design and functional analyses, we found that while keeping the hydrophobic face unchanged, increasing the number of alkaline residues (i.e., Lys + Arg residues) on the hydrophilic face of BmKn2 reduces the hemolytic activity and broadens the antimicrobial spectrum. Strikingly, when keeping the total number of alkaline residues constant, increasing the number of Lys residues on the hydrophilic face of BmKn2-7 significantly reduces the hemolytic activity but does not influence the antimicrobial activity. BmKn2-7K, a mutant of BmKn2-7 in which all of the Arg residues on the hydrophilic face were replaced with Lys, showed the lowest hemolytic activity and potent antimicrobial activity against antibiotic-resistant ESKAPE pathogens. Moreover, in vivo experiments indicate that BmKn2-7K displays potent antimicrobial efficacy against both the penicillin-resistant S. aureus and the carbapenem- and multidrug-resistant A. baumannii, and is non-toxic at the antimicrobial dosages. Taken together, our work highlights the significant functional disparity of Lys vs Arg in the scorpion-derived antimicrobial peptide BmKn2-7, and provides a promising lead molecule for drug development against ESKAPE pathogens.

Unevolved De Novo Proteins Have Innate Tendencies to Bind Transition Metals

Life as we know it would not exist without the ability of protein sequences to bind metal ions. Transition metals, in particular, play essential roles in a wide range of structural and catalytic functions. The ubiquitous occurrence of metalloproteins in all organisms leads one to ask whether metal binding is an evolved trait that occurred only rarely in ancestral sequences, or alternatively, whether it is an innate property of amino acid sequences, occurring frequently in unevolved sequence space. To address this question, we studied 52 proteins from a combinatorial library of novel sequences designed to fold into 4-helix bundles. Although these sequences were neither designed nor evolved to bind metals, the majority of them have innate tendencies to bind the transition metals copper, cobalt, and zinc with high nanomolar to low-micromolar affinity.

Long circulating genetically encoded intrinsically disordered zwitterionic polypeptides for drug delivery

The clinical utility of many peptide and protein drugs is limited by their short in-vivo half-life. To address this limitation, we report a new class of polypeptide-based materials that have a long plasma circulation time. The design of these polypeptides is motivated by the hypothesis that incorporating a zwitterionic sequence, within an intrinsically disordered polypeptide motif, would impart "stealth" behavior to the polypeptide and increase its plasma residence time, a behavior akin to that of synthetic stealth polymers. We designed these zwitterionic polypeptides (ZIPPs) with a repetitive (VPX₁X₂G)_n motif, where X₁ and X₂ are cationic and anionic amino acids, respectively, and n is the number of repeats. To test this hypothesis, we synthesized a set of ZIPPs with different pairs of cationic and anionic residues with varied chain length. We show that a combination of lysine and glutamic acid in the ZIPP confer superior pharmacokinetics, for both intravenous and subcutaneous administration, compared to uncharged control polypeptides. Finally, to demonstrate their clinical utility, we fused the best performing ZIPP sequence to glucagon-like peptide-1 (GLP1), a peptide drug used for treatment of type-2 diabetes and show that the ZIPP-GLP1 fusion outperforms an uncharged polypeptide of the same molecular weight in a mouse model of type-2 diabetes.

Errors in protein synthesis increase the level of saturated fatty acids and affect the overall lipid profiles of yeast

The occurrence of protein synthesis errors (mistranslation) above the typical mean mistranslation level of 10⁻⁴ is mostly deleterious to yeast, zebrafish and mammal cells. Previous yeast studies have shown that mistranslation affects fitness and deregulates genes related to lipid metabolism, but there is no experimental proof that such errors alter yeast lipid profiles. We engineered yeast strains to misincorporate serine at alanine and glycine sites on a global scale and evaluated the putative effects on the lipidome. Lipids from whole cells were extracted and analysed by thin layer chromatography (TLC), liquid chromatography-mass spectrometry(LC-MS) and gas chromatography (GC). Oxidative damage, fatty acid desaturation and membrane fluidity changes were screened to identify putative alterations in lipid profiles in both logarithmic (fermentative) and post-diauxic shift (respiratory) phases. There were alterations in several lipid classes, namely lyso-phosphatidylcholine, phosphatidic acid, phosphatidylethanolamine, phosphatidylinositol, phosphatidylserine, and triglyceride, and in the fatty acid profiles, namely C16:1, C16:0, C18:1 and C18:0. Overall, the relative content of lipid species with saturated FA increased in detriment of those with unsaturated fatty acids. The expression of the OLE1 mRNA was deregulated, but phospholipid fluidity changes were not observed. These data expand current knowledge of mistranslation biology and highlight its putative roles in human diseases.

Availability of NHS-biotin labeling to identify free protein lysine revealed by experiment and MD simulation

It is known that the crystallizability of protein molecules may be improved by replacing their surface lysine residues with other residue types. Here an experimental method to identify surface lysine residues by NHS-biotin chemical modification combined with MALDI-TOF MS was proposed and was evaluated using PH1033 protein from Pyrococcus horikoshii. Interestingly, the biotinylation experiment with a protein-reagent molar ratio of 1:1 revealed that only seven of twenty-two lysine residues in the protein comprising 144 residues were labeled. To investigate the result, we analyzed structures from a molecular-dynamics simulation mimicking the experiment. A logistic regression analysis revealed that the biotinylation was significantly correlated with four factors relevant to the local environment of lysine residues: the solvent accessibility, the electrostatic energy, the number of hydrogen bonds, and the estimated pKa value. This result is overall in agreement with that from the same analysis on the crystal structure. However, reflecting the flexibility of the protein molecule in solution state, the factors except for the electrostatic energy were highly variable in the MD structures depending upon the protonation state of Tyr87. The present procedure of biotin-labeling can avoid lysine residues with extensive intramolecular interactions that are incompatible with the rational design of protein crystals.

In proteins, the structural responses of a position to mutation rely on the Goldilocks principle: not too many links, not too few

A disease has distinct genetic and molecular hallmarks such as sequence variants that are likely to produce the alternative protein structures accountable for individual responses to drugs and disease development.

Designing better diffracting crystals of biotin carboxyl carrier protein from Pyrococcus horikoshii by a mutation based on the crystal-packing propensity of amino acids

An alternative rational approach to improve protein crystals by using single-site mutation of surface residues is proposed based on the results of a statistical analysis using a compiled data set of 918 independent crystal structures, thereby reflecting not only the entropic effect but also other effects upon protein crystallization. This analysis reveals a clear difference in the crystal-packing propensity of amino acids depending on the secondary-structural class. To verify this result, a systematic crystallization experiment was performed with the biotin carboxyl carrier protein from Pyrococcus horikoshii OT3 (PhBCCP). Six single-site mutations were examined: Ala138 on the surface of a β-sheet was mutated to Ile, Tyr, Arg, Gln, Val and Lys. In agreement with prediction, it was observed that the two mutants (A138I and A138Y) harbouring the residues with the highest crystal-packing propensities for β-sheet at position 138 provided better crystallization scores relative to those of other constructs, including the wild type, and that the crystal-packing propensity for β-sheet provided the best correlation with the ratio of obtaining crystals. Two new crystal forms of these mutants were obtained that diffracted to high resolution, generating novel packing interfaces with the mutated residues (Ile/Tyr). The mutations introduced did not affect the overall structures, indicating that a β-sheet can accommodate a successful mutation if it is carefully selected so as to avoid intramolecular steric hindrance. A significant negative correlation between the ratio of obtaining amorphous precipitate and the crystal-packing propensity was also found.

Protein stability: a crystallographer's perspective

Protein stability is a topic of major interest for the biotechnology, pharmaceutical and food industries, in addition to being a daily consideration for academic researchers studying proteins. An understanding of protein stability is essential for optimizing the expression, purification, formulation, storage and structural studies of proteins. In this review, discussion will focus on factors affecting protein stability, on a somewhat practical level, particularly from the view of a protein crystallographer. The differences between protein conformational stability and protein compositional stability will be discussed, along with a brief introduction to key methods useful for analyzing protein stability. Finally, tactics for addressing protein-stability issues during protein expression, purification and crystallization will be discussed.

An in silico approach to investigate the source of the controversial interpretations about the phenotypic results of the human AhR-gene G1661A polymorphism

Aryl hydrocarbon receptor (AhR) acts as an enhancer binding ligand-activated intracellular receptor. Chromatin remodeling components and general transcription factors such as TATA-binding protein (TBP) are evoked on AhR-target genes by interaction with its flexible transactivation domain (TAD). AhR-G1661A single nucleotide polymorphism (SNP: rs2066853) causes an arginine to lysine substitution in the acidic sub-domain of TAD at position 554 (R554K). Although, numerous studies associate the SNP with some abnormalities such as cancer, other reliable investigations refuse the associations. Consequently, the interpretation of the phenotypic results of G1661A-transition has been controversial. In this study, an in silico analysis were performed to investigate the possible effects of the transition on AhR-mRNA, protein structure, interaction properties and modifications. The analysis revealed that the R554K substitution affects secondary structure and solvent accessibility of adjacent residues. Also, it causes to decreasing of the AhR stability; altering the hydropathy features of the local sequence and changing the pattern of the residues at the binding site of the TAD-acidic sub-domain. Generating of new sites for ubiquitination and acetylation for AhR-K554 variant respectively at positions 544 and 560 was predicted. Our findings intensify the idea that the AhR-G1661A transition may affects AhR-TAD interactions, especially with the TBP, which influence AhR-target genes expression. However, the previously reported flexibility of the modular TAD could act as an intervening factor, moderate the SNP effects and causes distinct outcomes in different individuals and tissues.

MamA as a Model Protein for Structure-Based Insight into the Evolutionary Origins of Magnetotactic Bacteria

MamA is a highly conserved protein found in magnetotactic bacteria (MTB), a diverse group of prokaryotes capable of navigating according to magnetic fields – an ability known as magnetotaxis. Questions surround the acquisition of this magnetic navigation ability; namely, whether it arose through horizontal or vertical gene transfer. Though its exact function is unknown, MamA surrounds the magnetosome, the magnetic organelle embedding a biomineralised nanoparticle and responsible for magnetotaxis. Several structures for MamA from a variety of species have been determined and show a high degree of structural similarity. By determining the structure of MamA from Desulfovibrio magneticus RS-1 using X-ray crystallography, we have opened up the structure-sequence landscape. As such, this allows us to perform structural- and phylogenetic-based analyses using a variety of previously determined MamA from a diverse range of MTB species across various phylogenetic groups. We found that MamA has remained remarkably constant throughout evolution with minimal change between different taxa despite sequence variations. These findings, coupled with the generation of phylogenetic trees using both amino acid sequences and 16S rRNA, indicate that magnetotaxis likely did not spread via horizontal gene transfer and instead has a significantly earlier, primordial origin.

The role of flexibility and molecular shape in the crystallization of proteins by surface mutagenesis

Proteins are dynamic systems and interact with their environment. The analysis of crystal contacts in the most accurately determined protein structures (d < 1.5 Å) reveals that in contrast to current views, static disorder and high side-chain entropy are common in the crystal contact area. These observations challenge the validity of the theory that presumes that the occurrence of well ordered patches of side chains at the surface is an essential prerequisite for a successful crystallization event. The present paper provides evidence in support of the approach for understanding protein crystallization as a process dependent on multiple factors, each with its relative contribution, rather than a phenomenon driven by a few dominant physicochemical characteristics. The role of the molecular shape as a factor in the crystallization of proteins by surface mutagenesis is discussed.

Crystal structure and functional analysis of MiD49, a receptor for the mitochondrial fission protein Drp1

Mitochondrial fission requires recruitment of dynamin-related protein 1 (Drp1) to the mitochondrial surface, where assembly leads to activation of its GTP-dependent scission function. MiD49 and MiD51 are two receptors on the mitochondrial outer membrane that can recruit Drp1 to facilitate mitochondrial fission. Structural studies indicated that MiD51 has a variant nucleotidyl transferase fold that binds an ADP co-factor essential for activation of Drp1 function. MiD49 shares sequence homology with MiD51 and regulates Drp1 function. However, it is unknown if MiD49 binds an analogous co-factor. Because MiD49 does not readily crystallize, we used structural predictions and biochemical screening to identify a surface entropy reduction mutant that facilitated crystallization. Using molecular replacement, we determined the atomic structure of MiD49 to 2.4 Å. Like MiD51, MiD49 contains a nucleotidyl transferase domain; however, the electron density provides no evidence for a small-molecule ligand. Structural changes in the putative nucleotide-binding pocket make MiD49 incompatible with an extended ligand like ADP, and critical nucleotide-binding residues found in MiD51 are not conserved. MiD49 contains a surface loop that physically interacts with Drp1 and is necessary for Drp1 recruitment to the mitochondrial surface. Our results suggest a structural basis for the differential regulation of MiD51- versus MiD49-mediated fission.

Energetics-based methods for protein folding and stability measurements

Over the past 15 years, a series of energetics-based techniques have been developed for the thermodynamic analysis of protein folding and stability. These techniques include Stability of Unpurified Proteins from Rates of amide H/D Exchange (SUPREX), pulse proteolysis, Stability of Proteins from Rates of Oxidation (SPROX), slow histidine H/D exchange, lysine amidination, and quantitative cysteine reactivity (QCR). The above techniques, which are the subject of this review, all utilize chemical or enzymatic modification reactions to probe the chemical denaturant- or temperature-induced equilibrium unfolding properties of proteins and protein-ligand complexes. They employ various mass spectrometry-, sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE)-, and optical spectroscopy-based readouts that are particularly advantageous for high-throughput and in some cases multiplexed analyses. This has created the opportunity to use protein folding and stability measurements in new applications such as in high-throughput screening projects to identify novel protein ligands and in mode-of-action studies to identify protein targets of a particular ligand.

StableIsotope Labeling with Amino Acids in Cell Culture (SILAC)-based strategy for proteome-wide thermodynamic analysis of protein-ligand binding interactions

Described here is a quantitative mass spectrometry-based proteomics method for the large-scale thermodynamic analysis of protein-ligand binding interactions. The methodology utilizes a chemical modification strategy termed, Stability of Proteins from Rates of Oxidation (SPROX), in combination with a Stable Isotope Labeling with Amino Acids in Cell Culture (SILAC) approach to compare the equilibrium folding/unfolding properties of proteins in the absence and presence of target ligands. The method, which is general with respect to ligand, measures the ligand-induced changes in protein stability associated with protein-ligand binding. The methodology is demonstrated in a proof-of-principle study in which the well-characterized protein-drug interaction between cyclosporine A (CsA) and cyclophilin A was successfully analyzed in the context of a yeast cell lysate. A control experiment was also performed to assess the method's false positive rate of ligand discovery, which was found to be on the order of 0.4 - 3.5%. The new method was utilized to characterize the adenosine triphosphate (ATP)-interactome in Saccharomyces cerevisiae using the nonhydrolyzable ATP analog, adenylyl imidodiphosphate (AMP-PNP), and the proteins in a yeast cell lysate. The new methodology enabled the interrogation of 526 yeast proteins for interactions with ATP using 2035 peptide probes. Ultimately, 325 peptide hits from 139 different proteins were identified. Approximately 70% of the hit proteins identified in this work were not previously annotated as ATP binding proteins. However, nearly two-thirds of the newly discovered ATP interacting proteins have known interactions with other nucleotides and co-factors (e.g. NAD and GTP), DNA, and RNA based on GO-term analyses. The current work is the first proteome-wide profile of the yeast ATP-interactome, and it is the largest proteome-wide profile of any ATP-interactome generated, to date, using an energetics-based method. The data is available via ProteomeXchange with identifiers PXD000858, DOI 10.6019/PXD000858, and PXD000860.

Controlled apoptosis by a thermally toggled nanoscale amplifier of cellular uptake

Internalization into cancer cells is a significant challenge in the delivery of many anticancer therapeutics. Drug carriers can address this challenge by facilitating cellular uptake of cytotoxic cargo in the tumor, while preventing cellular uptake in healthy tissues. Here we describe an extrinsically controlled drug carrier, a nanopeptifier, that amplifies cellular uptake by modulating the activity of cell-penetrating peptides with thermally toggled self-assembly of a genetically encoded polypeptide nanoparticle. When appended with a proapoptotic peptide, the nanopeptifier creates a cytotoxic switch, inducing apoptosis only in its self-assembled state. The nanopeptifier provides a new approach to tune the cellular uptake and activity of anticancer therapeutics by an extrinsic thermal trigger.

A penalized-likelihood method to estimate the distribution of selection coefficients from phylogenetic data

We develop a maximum penalized-likelihood (MPL) method to estimate the fitnesses of amino acids and the distribution of selection coefficients (S = 2Ns) in protein-coding genes from phylogenetic data. This improves on a previous maximum-likelihood method. Various penalty functions are used to penalize extreme estimates of the fitnesses, thus correcting overfitting by the previous method. Using a combination of computer simulation and real data analysis, we evaluate the effect of the various penalties on the estimation of the fitnesses and the distribution of S. We show the new method regularizes the estimates of the fitnesses for small, relatively uninformative data sets, but it can still recover the large proportion of deleterious mutations when present in simulated data. Computer simulations indicate that as the number of taxa in the phylogeny or the level of sequence divergence increases, the distribution of S can be more accurately estimated. Furthermore, the strength of the penalty can be varied to study how informative a particular data set is about the distribution of S. We analyze three protein-coding genes (the chloroplast rubisco protein, mammal mitochondrial proteins, and an influenza virus polymerase) and show the new method recovers a large proportion of deleterious mutations in these data, even under strong penalties, confirming the distribution of S is bimodal in these real data. We recommend the use of the new MPL approach for the estimation of the distribution of S in species phylogenies of protein-coding genes.

Protein-protein interactions: switch from classical methods to proteomics and bioinformatics-based approaches

Following the sequencing of the human genome and many other organisms, research on protein-coding genes and their functions (functional genomics) has intensified. Subsequently, with the observation that proteins are indeed the molecular effectors of most cellular processes, the discipline of proteomics was born. Clearly, proteins do not function as single entities but rather as a dynamic network of team players that have to communicate. Though genetic (yeast two-hybrid Y2H) and biochemical methods (co-immunoprecipitation Co-IP, affinity purification AP) were the methods of choice at the beginning of the study of protein-protein interactions (PPI), in more recent years there has been a shift towards proteomics-based methods and bioinformatics-based approaches. In this review, we first describe in depth PPIs and we make a strong case as to why unraveling the interactome is the next challenge in the field of proteomics. Furthermore, classical methods of investigation of PPIs and structure-based bioinformatics approaches are presented. The greatest emphasis is placed on proteomic methods, especially native techniques that were recently developed and that have been shown to be reliable. Finally, we point out the limitations of these methods and the need to set up a standard for the validation of PPI experiments.

Recovery of a chemically synthesized Japanese encephalitis virus reveals two critical adaptive mutations in NS2B and NS4A

A full-length genome infectious clone is a powerful tool for functional assays in virology. In this study, using a chemical synthesized complete genome of Japanese encephalitis virus (JEV) strain SA14 (GenBank accession no. U14163), we constructed a full-length genomic cDNA clone of JEV. The recovered virus from the cDNA clone replicated poorly in baby hamster kidney (BHK-21) cells and in suckling mice brain. Following serial passage in BHK-21 cells, adaptive mutations within the NS2B and NS4A proteins were recovered in the passaged viruses leading to viruses with a large-plaque phenotype. Mutagenesis analysis, using a genome-length RNA and a replicon of JEV, demonstrated that the adaptive mutations restored replication to different degrees, and the restoration efficiencies were in the order: NS2B-T102M<NS4A-R79K<NS2B-T102M+NS4A-R79K. An in vivo virulence assay in mice showed that the recombinant virus containing double mutations showed similar virulence to the WT SA14 (GenBank accession no. M55506). This study reports the first chemically synthesized JEV. A reverse genetics assay demonstrated that substitutions of NS2B-T102M and NS4A-R79K altered JEV replication.

Structural determinants for protein adsorption/non-adsorption to silica surface

The understanding of the mechanisms involved in the interaction of proteins with inorganic surfaces is of major interest in both fundamental research and applications such as nanotechnology. However, despite intense research, the mechanisms and the structural determinants of protein/surface interactions are still unclear. We developed a strategy consisting in identifying, in a mixture of hundreds of soluble proteins, those proteins that are adsorbed on the surface and those that are not. If the two protein subsets are large enough, their statistical comparative analysis must reveal the physicochemical determinants relevant for adsorption versus non-adsorption. This methodology was tested with silica nanoparticles. We found that the adsorbed proteins contain a higher number of charged amino acids, particularly arginine, which is consistent with involvement of this basic amino acid in electrostatic interactions with silica. The analysis also identified a marked bias toward low aromatic amino acid content (phenylalanine, tryptophan, tyrosine and histidine) in adsorbed proteins. Structural analyses and molecular dynamics simulations of proteins from the two groups indicate that non-adsorbed proteins have twice as many π-π interactions and higher structural rigidity. The data are consistent with the notion that adsorption is correlated with the flexibility of the protein and with its ability to spread on the surface. Our findings led us to propose a refined model of protein adsorption.

Why do more divergent sequences produce smaller nonsynonymous/synonymous rate ratios in pairwise sequence comparisons?

Several studies have reported a negative correlation between estimates of the nonsynonymous to synonymous rate ratio (ω = dN/dS) and the sequence distance d in pairwise comparisons of the same gene from different species. That is, more divergent sequences produce smaller estimates of ω. Explanations for this negative correlation have included segregating nonsynonymous polymorphisms in closely related species and nonlinear dynamics of the ratio of two random variables. Here we study the statistical properties of the maximum-likelihood estimates of ω and d in pairwise alignments and explore the possibility that the negative correlation can be entirely explained by those properties. We show that the ω estimate is positively biased for small d and that the bias decreases with the increase of d. We also show that the estimates of ω and d are negatively correlated when ω < 1 and positively correlated when ω > 1. However, the bias in estimates of ω and the correlation between estimates of ω and d are not enough to explain the much stronger correlation observed in real data sets. We then explore the behavior of the estimates when the model is misspecified and suggest that the observed correlation may be due to protein-level selection that causes very different amino acids to be favored in different domains of the protein. Widely used models fail to account for such among-site heterogeneity and cause underestimation of the nonsynonymous rate and ω, with the bias being much stronger for distant sequences. We point out that tests of positive selection based on the ω ratio are invariant to the parameterization of the model and thus unaffected by bias in the ω estimates or the correlation between estimates of ω and d.

The proteomic future: where mass spectrometry should be taking us

A newcomer to the -omics era, proteomics, is a broad instrument-intensive research area that has advanced rapidly since its inception less than 20 years ago. Although the 'wet-bench' aspects of proteomics have undergone a renaissance with the improvement in protein and peptide separation techniques, including various improvements in two-dimensional gel electrophoresis and gel-free or off-gel protein focusing, it has been the seminal advances in MS that have led to the ascension of this field. Recent improvements in sensitivity, mass accuracy and fragmentation have led to achievements previously only dreamed of, including whole-proteome identification, and quantification and extensive mapping of specific PTMs (post-translational modifications). With such capabilities at present, one might conclude that proteomics has already reached its zenith; however, 'capability' indicates that the envisioned goals have not yet been achieved. In the present review we focus on what we perceive as the areas requiring more attention to achieve the improvements in workflow and instrumentation that will bridge the gap between capability and achievement for at least most proteomes and PTMs. Additionally, it is essential that we extend our ability to understand protein structures, interactions and localizations. Towards these ends, we briefly focus on selected methods and research areas where we anticipate the next wave of proteomic advances.

Folding without charges

Surface charges of proteins have in several cases been found to function as "structural gatekeepers," which avoid unwanted interactions by negative design, for example, in the control of protein aggregation and binding. The question is then if side-chain charges, due to their desolvation penalties, play a corresponding role in protein folding by avoiding competing, misfolded traps? To find out, we removed all 32 side-chain charges from the 101-residue protein S6 from Thermus thermophilus. The results show that the charge-depleted S6 variant not only retains its native structure and cooperative folding transition, but folds also faster than the wild-type protein. In addition, charge removal unleashes pronounced aggregation on longer timescales. S6 provides thus an example where the bias toward native contacts of a naturally evolved protein sequence is independent of charges, and point at a fundamental difference in the codes for folding and intermolecular interaction: specificity in folding is governed primarily by hydrophobic packing and hydrogen bonding, whereas solubility and binding relies critically on the interplay of side-chain charges.

Solution structure of CyanoP from Synechocystis sp. PCC 6803: new insights on the structural basis for functional specialization amongst PsbP family proteins

The structure of the CyanoP subunit of photosystem II from the cyanobacterium Synechocystis sp. PCC 6803 has been determined in solution by Nuclear Magnetic Resonance spectroscopy. Combined with homology modeling of PsbP-like structures we have identified distinct structural differences between PsbP homologues which may account for the functional differences apparent between members of this protein family. A surface cleft containing a large number of conserved residues found only in CyanoP and PsbP-like homologues has been identified and our findings suggest that one of the potential cation binding sites found in CyanoP may be functionally significant. Evidence for the evolution and divergence of the PsbP super family is presented from a structural perspective including identification of residues which distinguish the PsbP family from unrelated proteins with a similar domain fold. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Genetic construct design and recombinant protein expression for structural biology

Obtaining diffraction quality crystals is frequently an iterative process which traditionally has involved screening large numbers of crystallization conditions. Due to advances in high-throughput gene engineering, recombinant expression, and purification, the protein of interest has now become one of the many variables routinely investigated during crystallization trials. As such, construct design is a critical step in the path toward successful crystallization. In this chapter will we address construct design strategies frequently employed to improve the solution and crystallization behavior of proteins. Topics covered include choosing a recombinant expression system and reducing disorder through truncations and surface mutagenesis. Also covered are strategies to reduce heterogeneity from posttranslational modifications, impurities, and aggregation.

Improving MM-GB/SA Scoring through the Application of the Variable Dielectric Model

A variable dielectric model based on residue types for better description of protein-ligand electrostatics in MM-GBSA scoring is reported. The variable dielectric approach provides better correlation with binding data and reduces the score dynamic range, typically observed in the standard MM-GB/SA method. The latter supports the view that exaggerated enthalpic separation between weak and potent compounds due to the lack of shielding effects in the model is greatly responsible for the wide scoring spread.

On the selectivity and efficacy of defense peptides with respect to cancer cells

Here, we review potential determinants of the anticancer efficacy of innate immune peptides (ACPs) for cancer cells. These determinants include membrane-based factors, such as receptors, phosphatidylserine, sialic acid residues, and sulfated glycans, and peptide-based factors, such as residue composition, sequence length, net charge, hydrophobic arc size, hydrophobicity, and amphiphilicity. Each of these factors may contribute to the anticancer action of ACPs, but no single factor(s) makes an overriding contribution to their overall selectivity and toxicity. Differences between the anticancer actions of ACPs seem to relate to different levels of interplay between these peptide and membrane-based factors.

Probing phenylalanine environments in oligomeric structures with pentafluorophenylalanine and cyclohexylalanine

Stabilization of protein structures and protein-protein interactions are critical in the engineering of industrially useful enzymes and in the design of pharmaceutically valuable ligands. Hydrophobic interactions involving phenylalanine residues play crucial roles in protein stability and protein-protein/peptide interactions. To establish an effective method to explore the hydrophobic environments of phenylalanine residues, we present a strategy that uses pentafluorophenylalanine (F5Phe) and cyclohexylalanine (Cha). In this study, substitution of F5Phe or Cha for three Phe residues at positions 328, 338, and 341 in the tetramerization domain of the tumor suppressor protein p53 was performed. These residues are located at the interfaces of p53-p53 interactions and are important in the stabilization of the tetrameric structure. The stability of the p53 tetrameric structure did not change significantly when F5Phe-containing peptides at positions Phe328 or Phe338 were used. In contrast, the substitution of Cha for Phe341 in the hydrophobic core enhanced the stability of the tetrameric structure with a T(m) value of 100 degrees C. Phe328 and Phe338 interact with each other through pi-interactions, whereas Phe341 is buried in the surrounding alkyl side-chains of the hydrophobic core of the p53 tetramerization domain. Furthermore, high pressure-assisted denaturation analysis indicated improvement in the occupancy of the hydrophobic core. Considerable stabilization of the p53 tetramer was achieved by filling the identified cavity in the hydrophobic core of the p53 tetramer. The results indicate the status of the Phe residues, indicating that the "pair substitution" of Cha and F5Phe is highly suitable for probing the environments of Phe residues.

An accurate feature-based method for identifying DNA-binding residues on protein surfaces

Proteins that interact with DNA play vital roles in all mechanisms of gene expression and regulation. In order to understand these activities, it is crucial to analyze and identify DNA-binding residues on DNA-binding protein surfaces. Here, we proposed two novel features B-factor and packing density in combination with several conventional features to characterize the DNA-binding residues in a well-constructed representative dataset of 119 protein-DNA complexes from the Protein Data Bank (PDB). Based on the selected features, a prediction model for DNA-binding residues was constructed using support vector machine (SVM). The predictor was evaluated using a 5-fold cross validation on above dataset of 123 DNA-binding proteins. Moreover, two independent datasets of 83 DNA-bound protein structures and their corresponding DNA-free forms were compiled. The B-factor and packing density features were statistically analyzed on these 83 pairs of holo-apo proteins structures. Finally, we developed the SVM model to accurately predict DNA-binding residues on protein surface, given the DNA-free structure of a protein. Results showed here indicate that our method represents a significant improvement of previously existing approaches such as DISPLAR. The observation suggests that our method will be useful in studying protein-DNA interactions to guide consequent works such as site-directed mutagenesis and protein-DNA docking.

Synthetic symmetrization in the crystallization and structure determination of CelA from Thermotoga maritima

Protein crystallization continues to be a major bottleneck in X-ray crystallography. Previous studies suggest that symmetric proteins, such as homodimers, might crystallize more readily than monomeric proteins or asymmetric complexes. Proteins that are naturally monomeric can be made homodimeric artificially. Our approach is to create homodimeric proteins by introducing single cysteines into the protein of interest, which are then oxidized to form a disulfide bond between the two monomers. By introducing the single cysteine at different sequence positions, one can produce a variety of synthetically dimerized versions of a protein, with each construct expected to exhibit its own crystallization behavior. In earlier work, we demonstrated the potential utility of the approach using T4 lysozyme as a model system. Here we report the successful application of the method to Thermotoga maritima CelA, a thermophilic endoglucanase enzyme with low sequence identity to proteins with structures previously reported in the Protein Data Bank. This protein had resisted crystallization in its natural monomeric form, despite a broad survey of crystallization conditions. The synthetic dimerization of the CelA mutant D188C yielded well-diffracting crystals with molecules in a packing arrangement that would not have occurred with native, monomeric CelA. A 2.4 Å crystal structure was determined by single anomalous dispersion using a seleno-methionine derivatized protein. The results support the notion that synthetic symmetrization can be a useful approach for enlarging the search space for crystallizing monomeric proteins or asymmetric complexes.

Structure of the putative dihydroorotate dehydrogenase from Streptococcus mutans

Streptococcus mutans is one of the pathogenic species involved in dental caries, especially in the initiation and development stages. Here, the crystal structure of SMU.595, a putative dihydroorotate dehydrogenase (DHOD) from S. mutans, is reported at 2.4 Å resolution. DHOD is a flavin mononucleotide-containing enzyme which catalyzes the oxidation of L-dihydroorotate to orotate, which is the fourth step and the only redox reaction in the de novo biosynthesis of pyrimidine nucleotides. The reductive lysine-methylation procedure was applied in order to improve the diffraction qualities of the crystals. Analysis of the S. mutans DHOD crystal structure shows that this enzyme is a class 1A DHOD and also suggests potential sites that could be exploited for the design of highly specific inhibitors using the structure-based chemotherapeutic design technique.

Molecular modeling of mechanosensory ion channel structural and functional features

The DEG/ENaC (Degenerin/Epithelial Sodium Channel) protein family comprises related ion channel subunits from all metazoans, including humans. Members of this protein family play roles in several important biological processes such as transduction of mechanical stimuli, sodium re-absorption and blood pressure regulation. Several blocks of amino acid sequence are conserved in DEG/ENaC proteins, but structure/function relations in this channel class are poorly understood. Given the considerable experimental limitations associated with the crystallization of integral membrane proteins, knowledge-based modeling is often the only route towards obtaining reliable structural information. To gain insight into the structural characteristics of DEG/ENaC ion channels, we derived three-dimensional models of MEC-4 and UNC-8, based on the available crystal structures of ASIC1 (Acid Sensing Ion Channel 1). MEC-4 and UNC-8 are two DEG/ENaC family members involved in mechanosensation and proprioception respectively, in the nematode Caenorhabditis elegans. We used these models to examine the structural effects of specific mutations that alter channel function in vivo. The trimeric MEC-4 model provides insight into the mechanism by which gain-of-function mutations cause structural alterations that result in increased channel permeability, which trigger cell degeneration. Our analysis provides an introductory framework to further investigate the multimeric organization of the DEG/ENaC ion channel complex.

Defining the functional domain of programmed cell death 10 through its interactions with phosphatidylinositol-3,4,5-trisphosphate

Cerebral cavernous malformations (CCM) are vascular abnormalities of the central nervous system predisposing blood vessels to leakage, leading to hemorrhagic stroke. Three genes, Krit1 (CCM1), OSM (CCM2), and PDCD10 (CCM3) are involved in CCM development. PDCD10 binds specifically to PtdIns(3,4,5)P3 and OSM. Using threading analysis and multi-template modeling, we constructed a three-dimensional model of PDCD10. PDCD10 appears to be a six-helical-bundle protein formed by two heptad-repeat-hairpin structures (α1–3 and α4–6) sharing the closest 3D homology with the bacterial phosphate transporter, PhoU. We identified a stretch of five lysines forming an amphipathic helix, a potential PtdIns(3,4,5)P3 binding site, in the α5 helix. We generated a recombinant wild-type (WT) and three PDCD10 mutants that have two (Δ2KA), three (Δ3KA), and five (Δ5KA) K to A mutations. Δ2KA and Δ3KA mutants hypothetically lack binding residues to PtdIns(3,4,5)P3 at the beginning and the end of predicted helix, while Δ5KA completely lacks all predicted binding residues. The WT, Δ2KA, and Δ3KA mutants maintain their binding to PtdIns(3,4,5)P3. Only the Δ5KA abolishes binding to PtdIns(3,4,5)P3. Both Δ5KA and WT show similar secondary and tertiary structures; however, Δ5KA does not bind to OSM. When WT and Δ5KA are co-expressed with membrane-bound constitutively-active PI3 kinase (p110-CAAX), the majority of the WT is co-localized with p110-CAAX at the plasma membrane where PtdIns(3,4,5)P3 is presumably abundant. In contrast, the Δ5KA remains in the cytoplasm and is not present in the plasma membrane. Combining computational modeling and biological data, we propose that the CCM protein complex functions in the PI3K signaling pathway through the interaction between PDCD10 and PtdIns(3,4,5)P3.

Electropositive charge in alpha-defensin bactericidal activity: functional effects of Lys-for-Arg substitutions vary with the peptide primary structure

Cationic amino acids contribute to alpha-defensin bactericidal activity. Curiously, although Arg and Lys have equivalent electropositive charges at neutral pH, alpha-defensins contain an average of nine Arg residues per Lys residue. To investigate the role of high alpha-defensin Arg content, all Arg residues in mouse Paneth cell alpha-defensin cryptdin 4 (Crp4) and rhesus myeloid alpha-defensin 4 (RMAD-4) were replaced with Lys to prepare (R/K)-Crp4 and (R/K)-RMAD-4, respectively. Lys-for-Arg replacements in Crp4 attenuated bactericidal activity and slowed the kinetics of Escherichia coli ML35 cell permeabilization, and (R/K)-Crp4 required longer exposure times to reduce E. coli cell survival. In marked contrast, Lys substitutions in RMAD-4 improved microbicidal activity against certain bacteria and permeabilized E. coli more effectively. Therefore, Arg-->Lys substitutions attenuated activity in Crp4 but not in RMAD-4, and the functional consequences of Arg-->Lys replacements in alpha-defensins are dependent on the peptide primary structure. In addition, the bactericidal effects of (R/K)-Crp4 and (R/K)-RMAD-4 were more sensitive to inhibition by NaCl than those of the native peptides, suggesting that the high Arg content of alpha-defensins may be under selection to confer superior microbicidal function under physiologic conditions.

Toward rational protein crystallization: A Web server for the design of crystallizable protein variants

Growing well-diffracting crystals constitutes a serious bottleneck in structural biology. A recently proposed crystallization methodology for "stubborn crystallizers" is to engineer surface sequence variants designed to form intermolecular contacts that could support a crystal lattice. This approach relies on the concept of surface entropy reduction (SER), i.e., the replacement of clusters of flexible, solvent-exposed residues with residues with lower conformational entropy. This strategy minimizes the loss of conformational entropy upon crystallization and renders crystallization thermodynamically favorable. The method has been successfully used to crystallize more than 15 novel proteins, all stubborn crystallizers. But the choice of suitable sites for mutagenesis is not trivial. Herein, we announce a Web server, the surface entropy reduction prediction server (SERp server), designed to identify mutations that may facilitate crystallization. Suggested mutations are predicted based on an algorithm incorporating a conformational entropy profile, a secondary structure prediction, and sequence conservation. Minor considerations include the nature of flanking residues and gaps between mutation candidates. While designed to be used with default values, the server has many user-controlled parameters allowing for considerable flexibility. Within, we discuss (1) the methodology of the server, (2) how to interpret the results, and (3) factors that must be considered when selecting mutations. We also attempt to benchmark the server by comparing the server's predictions with successful SER structures. In most cases, the structure yielding mutations were easily identified by the SERp server. The server can be accessed at http://www.doe-mbi.ucla.edu/Services/SER.

Crystal structures of the luciferase and green fluorescent protein from Renilla reniformis

Due to its ability to emit light, the luciferase from Renilla reniformis (RLuc) is widely employed in molecular biology as a reporter gene in cell culture experiments and small animal imaging. To accomplish this bioluminescence, the 37-kDa enzyme catalyzes the degradation of its substrate coelenterazine in the presence of molecular oxygen, resulting in the product coelenteramide, carbon dioxide, and the desired photon of light. We successfully crystallized a stabilized variant of this important protein (RLuc8) and herein present the first structures for any coelenterazine-using luciferase. These structures are based on high-resolution data measured to 1.4 A and demonstrate a classic alpha/beta-hydrolase fold. We also present data of a coelenteramide-bound luciferase and reason that this structure represents a secondary conformational form following shift of the product out of the primary active site. During the course of this work, the structure of the luciferase's accessory green fluorescent protein (RrGFP) was also determined and shown to be highly similar to that of Aequorea victoria GFP.

(NZ)CH...O contacts assist crystallization of a ParB-like nuclease

Background

The major bottleneck for determination of 3 D structures of proteins using X-rays is the production of diffraction quality crystals. Often proteins are subjected to chemical modification to improve the chances of crystallization

Results

Here, we report the successful crystallization of a nuclease employing a reductive methylation protocol. The key to crystallization was the successful introduction of 44 new cohesive (NZ) CH...O contacts (3.2 – 3.7 Å) by the addition of 2 methyl groups to the side chain amine nitrogen (NZ) of 9 lysine residues of the nuclease. The new contacts dramatically altered the crystallization properties of the protein, resulting in crystals that diffracted to 1.2 Å resolution. Analytical ultracentrifugation analysis and thermodynamics results revealed a more compact protein structure with better solvent exclusion of buried Trp residues in the folded state of the methylated protein, assisting crystallization.

Conclusion

In this study, introduction of novel cohesive (NZ)CH...O contacts by reductive methylation resulted in the crystallization of a protein that had previously resisted crystallization in spite of extensive purification and crystallization space screening. Introduction of (NZ)CH...O contacts could provide a solution to crystallization problems for a broad range of protein targets.

Toward Understanding the Cationicity of Defensins

Human defensins are a family of small antimicrobial proteins found predominantly in leukocytes and epithelial cells that play important roles in the innate and adaptive immune defense against microbial infection. The most distinct molecular feature of defensins is cationicity, manifested by abundant Arg and/or Lys residues in their sequences. Sequence analysis indicates that Arg is strongly selected over Lys in alpha-defensins but not in beta-defensins. To understand this Arg/Lys disparity in defensins, we chemically synthesized human alpha-defensin 1 (HNP1) and several HNP1 analogs where three Arg residues were replaced by each of the following six alpha-amino acids: Lys, ornithine (Orn), diaminobutyric acid (Dab), diaminopropionic acid (Dap), N,N-dimethyl-Lys ((diMe)Lys), and homo-Arg ((homo)Arg). In addition, we prepared human beta-defensin 1 (hBD1) and (Lys-->Arg)hBD1 in which all four Lys residues were substituted for Arg. Bactericidal activity assays revealed the following. 1) Arg-containing HNP1 and (Lys-->Arg)hBD1 are functionally better than Lys-HNP1 and hBD1, respectively; the difference between Arg and Lys is more evident in the alpha-defensin than in the beta-defensin and is more evident at low salt concentrations than at high salt concentrations. 2) For HNP1, the Arg/Lys disparity is much more pronounced with Staphylococcus aureus than with Escherichia coli, and the Arg-rich HNP1 kills bacteria faster than its Lys-rich analog. 3) Arg and Lys appear to have optimal chain lengths for bacterial killing as shortening Lys or lengthening Arg in HNP1 invariably becomes functionally deleterious. Our findings provide insights into the Arg/Lys disparity in defensins, and shed light on the cationicity of defensins with respect to their antimicrobial activity and specificity.

Geometry of nonbonded interactions involving planar groups in proteins

Although hydrophobic interaction is the main contributing factor to the stability of the protein fold, the specificity of the folding process depends on many directional interactions. An analysis has been carried out on the geometry of interaction between planar moieties of ten side chains (Phe, Tyr, Trp, His, Arg, Pro, Asp, Glu, Asn and Gln), the aromatic residues and the sulfide planes (of Met and cystine), and the aromatic residues and the peptide planes within the protein tertiary structures available in the Protein Data Bank. The occurrence of hydrogen bonds and other nonconventional interactions such as C-H...pi, C-H...O, electrophile-nucleophile interactions involving the planar moieties has been elucidated. The specific nature of the interactions constraints many of the residue pairs to occur with a fixed sequence difference, maintaining a sequential order, when located in secondary structural elements, such as alpha-helices and beta-turns. The importance of many of these interactions (for example, aromatic residues interacting with Pro or cystine sulfur atom) is revealed by the higher degree of conservation observed for them in protein structures and binding regions. The planar residues are well represented in the active sites, and the geometry of their interactions does not deviate from the general distribution. The geometrical relationship between interacting residues provides valuable insights into the process of protein folding and would be useful for the design of protein molecules and modulation of their binding properties.

Protein engineering

This chapter focuses on protein engineering strategies that aim to increase the chances of obtaining crystals suitable for X-ray diffraction. The chapter is divided into three main parts: one dealing with protein engineering through a bioinformatics approach, the second focusing on DNA modifications via random mutagenesis, and the third describing a nonexhaustive number of in vitro modifications based on site-directed mutagenesis.