Hydrogen bonds between short polar side chains and peptide backbone: Prevalence in proteins and effects on helix-forming propensities

Effects of sulfoxide and sulfone sidechain-backbone hydrogen bonding on local conformations in peptide models

We examine peptide model systems designed to probe short-range N-H⋯OS sidechain-backbone hydrogen bonding involving amino acid residues with sidechain sulfoxide or sulfone functional groups and its effects on local conformations. A strong 7-membered ring hydrogen bond of this type accompanies an intra-residue N-H⋯OC interaction and stabilizes an extended backbone conformation in preference to classical folded structures.

Peptide Stereochemistry Effects from pKa-Shift to Gold Nanoparticle Templating in a Supramolecular Hydrogel

The divergent supramolecular behavior of a series of tripeptide stereoisomers was elucidated through spectroscopic, microscopic, crystallographic, and computational techniques. Only two epimers were able to effectively self-organize into amphipathic structures, leading to supramolecular hydrogels or crystals, respectively. Despite the similarity between the two peptides' turn conformations, stereoconfiguration led to different abilities to engage in intramolecular hydrogen bonding. Self-assembly further shifted the pKa value of the C-terminal side chain. As a result, across the pH range 4-6, only one epimer predominated sufficiently as a zwitterion to reach the critical molar fraction, allowing gelation. By contrast, the differing pKa values and higher dipole moment of the other epimer favored crystallization. The four stereoisomers were further tested for gold nanoparticle (AuNP) formation, with the supramolecular hydrogel being the key to control and stabilize AuNPs, yielding a nanocomposite that catalyzed the photodegradation of a dye. Importantly, the AuNP formation occurred without the use of reductants other than the peptide, and the redox chemistry was investigated by LC-MS, NMR, and infrared scattering-type near field optical microscopy (IR s-SNOM). This study provides important insights for the rational design of simple peptides as minimalistic and green building blocks for functional nanocomposites.

Dual function of Zika virus NS2B-NS3 protease

Zika virus (ZIKV) serine protease, indispensable for viral polyprotein processing and replication, is composed of the membrane-anchored NS2B polypeptide and the N-terminal domain of the NS3 polypeptide (NS3pro). The C-terminal domain of the NS3 polypeptide (NS3hel) is necessary for helicase activity and contains an ATP-binding site. We discovered that ZIKV NS2B-NS3pro binds single-stranded RNA with a Kd of ~0.3 μM, suggesting a novel function. We tested various structural modifications of NS2B-NS3pro and observed that constructs stabilized in the recently discovered "super-open" conformation do not bind RNA. Likewise, stabilizing NS2B-NS3pro in the "closed" (proteolytically active) conformation using substrate inhibitors abolished RNA binding. We posit that RNA binding occurs when ZIKV NS2B-NS3pro adopts the "open" conformation, which we modeled using highly homologous dengue NS2B-NS3pro crystallized in the open conformation. We identified two positively charged fork-like structures present only in the open conformation of NS3pro. These forks are conserved across Flaviviridae family and could be aligned with the positively charged grove on NS3hel, providing a contiguous binding surface for the negative RNA strand exiting helicase. We propose a "reverse inchworm" model for a tightly intertwined NS2B-NS3 helicase-protease machinery, which suggests that NS2B-NS3pro cycles between open and super-open conformations to bind and release RNA enabling long-range NS3hel processivity. The transition to the closed conformation, likely induced by the substrate, enables the classical protease activity of NS2B-NS3pro.

Non-Conserved Amino Acid Residues Modulate the Thermodynamics of Zn(II) Binding to Classical ββα Zinc Finger Domains

Classical zinc fingers domains (ZFs) bind Zn(II) ion by a pair of cysteine and histidine residues to adopt a characteristic and stable ββα fold containing a small hydrophobic core. As a component of transcription factors, they recognize specific DNA sequences to transcript particular genes. The loss of Zn(II) disrupts the unique structure and function of the whole protein. It has been shown that the saturation of ZFs under cellular conditions is strictly related to their affinity for Zn(II). High affinity warrants their constant saturation, while medium affinity results in their transient structurization depending on cellular zinc availability. Therefore, there must be factors hidden in the sequence and structure of ZFs that impact Zn(II)-to-protein affinities to control their function. Using molecular dynamics simulations and experimental spectroscopic and calorimetric approaches, we showed that particular non-conserved residues derived from ZF sequences impact hydrogen bond formation. Our in silico and in vitro studies show that non-conserved residues can alter metal-coupled folding mechanisms and overall ZF stability. Furthermore, we show that Zn(II) binding to ZFs can also be entropically driven. This preference does not correlate either with Zn(II) binding site or with the extent of the secondary structure but is strictly related to a reservoir of interactions within the second coordination shell, which may loosen or tighten up the structure. Our findings shed new light on how the functionality of ZFs is modulated by non-coordinating residues diversity under cellular conditions. Moreover, they can be helpful for systematic backbone alteration of native ZF ββα scaffold to create artificial foldamers and proteins with improved stability.

Backbone and side-chain resonance assignments of the NISTmAb-scFv and antigen-binding study

Monoclonal antibodies (mAbs) therapeutics are the largest and fastest growing class of biologic drugs, amongst which, the vast majority are immunoglobulin G1 (IgG1). Their antigen binding abilities are used for the treatment of immunologic diseases, cancer therapy, reversal of drug effects, and targeting viruses and bacteria. The high importance of therapeutic mAbs and their derivatives has called for the generation of well-characterized standards for method development and calibration. One such standard, the NISTmAb RM 8621 based on the antibody motavizumab, has been developed by the National Institute of Standards and Technologies (NIST) in the US. Here, we present the resonance assignment of the single chain variable fragment, NISTmAb-scFv, that was engineered by linking the variable domains of the heavy and light chains of the NISTmAb. Also, addition of a peptide, corresponding to the target antigen of motavizumab, to samples of NISTmAb-scFv has induced chemical shift perturbations on residues lining the antigen binding interface thereby indicating proper folding of the NISTmAb-scFv.

SM-COLSARSPROT: Highly Immunogenic Supramutational Synthetic Peptides Covering the World's Population

Fifty ~20-amino acid (aa)-long peptides were selected from functionally relevant SARS-CoV-2 S, M, and E proteins for trial B-21 and another 53 common ones, plus some new ones derived from the virus' main genetic variants for complementary trial C-21. Peptide selection was based on tremendous SARS-CoV-2 genetic variability for analysing them concerning vast human immunogenetic polymorphism for developing the first supramutational, Colombian SARS-protection (SM-COLSARSPROT), peptide mixture. Specific physicochemical rules were followed, i.e., aa predilection for polyproline type II left-handed (PPIIL) formation, replacing β-branched, aromatic aa, short-chain backbone H-bond-forming residues, π-π interactions (n→π* and π-CH), aa interaction with π systems, and molecular fragments able to interact with them, disrupting PPIIL propensity formation. All these modified structures had PPIIL formation propensity to enable target peptide interaction with human leukocyte antigen-DRβ1* (HLA-DRβ1*) molecules to mediate antigen presentation and induce an appropriate immune response. Such modified peptides were designed for human use; however, they induced high antibody titres against S, M, and E parental mutant peptides and neutralising antibodies when suitably modified and chemically synthesised for immunising 61 major histocompatibility complex class II (MHCII) DNA genotyped Aotus monkeys (matched with their corresponding HLA-DRβ1* molecules), predicted to cover 77.5% to 83.1% of the world's population. Such chemically synthesised peptide mixture represents an extremely pure, stable, reliable, and cheap vaccine for COVID-19 pandemic control, providing a new approach for a logical, rational, and soundly established methodology for other vaccine development.

Design, synthesis, and evaluation of the self-assembled antimicrobial peptides based on the ovalbumin-derived peptide TK913

The preparation, self-assembly, and antimicrobial activity of peptides based on TK913 is described. TK9Z4 incorporating a Pro-Pro motif exhibited self-assembly but no cytotoxicity. However, peptide TKZ3 (obtained by changing the amino acid sequence of TK9Z4) showed morphological changes at different concentrations, potent antimicrobial activity, low cytotoxicity, and trypsin resistance. Accordingly, TKZ3 is proposed as new AMP derived from ovalbumin-derived peptides.

Identification of an α-MoRF in the Intrinsically Disordered Region of the Escargot Transcription Factor

Molecular recognition features (MoRFs) are common in intrinsically disordered proteins (IDPs) and intrinsically disordered regions (IDRs). MoRFs are in constant order-disorder structural transitions and adopt well-defined structures once they are bound to their targets. Here, we study Escargot (Esg), a transcription factor in Drosophila melanogaster that regulates multiple cellular functions, and consists of a disordered N-terminal domain and a group of zinc fingers at its C-terminal domain. We analyzed the N-terminal domain of Esg with disorder predictors and identified a region of 45 amino acids with high probability to form ordered structures, which we named S2. Through 54 μs of molecular dynamics (MD) simulations using CHARMM36 and implicit solvent (generalized Born/surface area (GBSA)), we characterized the conformational landscape of S2 and found an α-MoRF of ∼16 amino acids stabilized by key contacts within the helix. To test the importance of these contacts in the stability of the α-MoRF, we evaluated the effect of point mutations that would impair these interactions, running 24 μs of MD for each mutation. The mutations had mild effects on the MoRF, and in some cases, led to gain of residual structure through long-range contacts of the α-MoRF and the rest of the S2 region. As this could be an effect of the force field and solvent model we used, we benchmarked our simulation protocol by carrying out 32 μs of MD for the (AAQAA)₃ peptide. The results of the benchmark indicate that the global amount of helix in shorter peptides like (AAQAA)₃ is reasonably predicted. Careful analysis of the runs of S2 and its mutants suggests that the mutation to hydrophobic residues may have nucleated long-range hydrophobic and aromatic interactions that stabilize the MoRF. Finally, we have identified a set of residues that stabilize an α-MoRF in a region still without functional annotations in Esg.

Conformational Switch to a β-Turn in a Staphylococcal Quorum Sensing Signal Peptide Causes a Dramatic Increase in Potency

We report the solution-phase structures of native signal peptides and related analogs capable of either strongly agonizing or antagonizing the AgrC quorum sensing (QS) receptor in the emerging pathogen Staphylococcus epidermidis. Chronic S. epidermidis infections are often recalcitrant to traditional therapies due to antibiotic resistance and formation of robust biofilms. The accessory gene regulator (agr) QS system plays an important role in biofilm formation in this opportunistic pathogen, and the binding of an autoinducing peptide (AIP) signal to its cognate transmembrane receptor (AgrC) is responsible for controlling agr. Small molecules or peptides capable of modulating this binding event are of significant interest as probes to investigate both the agr system and QS as a potential antivirulence target. We used NMR spectroscopy to characterize the structures of the three native S. epidermidis AIP signals and five non-native analogs with distinct activity profiles in the AgrC-I receptor from S. epidermidis. These studies revealed a suite of structural motifs critical for ligand activity. Interestingly, a unique β-turn was present in the macrocycles of the two most potent AgrC-I modulators, in both an agonist and an antagonist, which was distinct from the macrocycle conformation in the less-potent AgrC-I modulators and in the native AIP-I itself. This previously unknown β-turn provides a structural rationale for these ligands' respective biological activity profiles. Development of analogs to reinforce the β-turn resulted in our first antagonist with subnanomolar potency in AgrC-I, while analogs designed to contain a disrupted β-turn were dramatically less potent relative to their parent compounds. Collectively, these studies provide new insights into the AIP:AgrC interactions crucial for QS activation in S. epidermidis and advance the understanding of QS at the molecular level.

Lysine/RNA-interactions drive and regulate biomolecular condensation

Cells form and use biomolecular condensates to execute biochemical reactions. The molecular properties of non-membrane-bound condensates are directly connected to the amino acid content of disordered protein regions. Lysine plays an important role in cellular function, but little is known about its role in biomolecular condensation. Here we show that protein disorder is abundant in protein/RNA granules and lysine is enriched in disordered regions of proteins in P-bodies compared to the entire human disordered proteome. Lysine-rich polypeptides phase separate into lysine/RNA-coacervates that are more dynamic and differ at the molecular level from arginine/RNA-coacervates. Consistent with the ability of lysine to drive phase separation, lysine-rich variants of the Alzheimer's disease-linked protein tau undergo coacervation with RNA in vitro and bind to stress granules in cells. Acetylation of lysine reverses liquid-liquid phase separation and reduces colocalization of tau with stress granules. Our study establishes lysine as an important regulator of cellular condensation.

Conserved Residues Lys57 and Lys401 of Protein Disulfide Isomerase Maintain an Active Site Conformation for Optimal Activity: Implications for Post-Translational Regulation

Despite its study since the 1960's, very little is known about the post-translational regulation of the multiple catalytic activities performed by protein disulfide isomerase (PDI), the primary protein folding catalyst of the cell. This work identifies a functional role for the highly conserved CxxC-flanking residues Lys⁵⁷ and Lys⁴⁰¹ of human PDI in vitro. Mutagenesis studies have revealed these residues as modulating the oxidoreductase activity of PDI in a pH-dependent manner. Non-conservative amino acid substitutions resulted in enzyme variants upwards of 7-fold less efficient. This attenuated activity was found to translate into a 2-fold reduction of the rate of electron shuttling between PDI and the intraluminal endoplasmic reticulum oxidase, ERO1α, suggesting a functional significance to oxidative protein folding. In light of this, the possibility of lysine acetylation at residues Lys⁵⁷ and Lys⁴⁰¹ was assessed by in vitro treatment using acetylsalicylic acid (aspirin). A total of 28 acetyllysine residues were identified, including acLys⁵⁷ and acLys⁴⁰¹. The kinetic behavior of the acetylated protein form nearly mimicked that obtained with a K57/401Q double substitution variant providing an indication that acetylation of the active site-flanking lysine residues can act to reversibly modulate PDI activity.

Bioengineered porous composite curcumin/silk scaffolds for cartilage regeneration

Articular cartilage repair is a challenge due to its limited self-repair capacity. Cartilage tissue engineering supports to overcome following injuries or degenerative diseases. Herein, we fabricated the scaffold composed of curcumin and silk fibroin as an appropriate clinical replacement for defected cartilage. The scaffolds were designed to have adequate pore size and mechanical strength for cartilage repair. Cell proliferation, sulfated glycosaminoglycan (sGAG) content and mRNA expression analysis indicated that chondrocytes remained viable and showed its growth ability in the curcumin/silk scaffolds. Especially, in 1mg/ml curcumin/silk scaffold showed higher cell viability rate and extracellular matrix formation than other experimental groups. Furthermore, curcumin/silk scaffold showed its biocompatibility and favorable environment for cartilage repair after transplantation in vivo, as indicated in histological examination results. Overall, the functional composite curcumin/silk scaffold can be applied in cartilage tissue engineering and promising substrate for cartilage repair.

Determinants of the pKa values of ionizable residues in an intrinsically disordered protein

Intrinsically disordered proteins (IDPs) are prevalent in eukaryotes; in humans, they are often associated with diseases. The protein NUPR1 is a multifunctional IDP involved in the development and progression of pancreatic cancer; therefore, it constitutes a target for drug design. In an effort to contribute to the understanding of the conformational features of NUPR1 and to provide clues on amino acid interactions in disordered states of proteins, we measured the pKa values of all its acidic groups (aspartic and glutamic residues, and backbone C terminus) by using NMR spectroscopy at low (100 mM) and high (500 mM) NaCl concentration. At low ionic strength, the pKa values were similar to those reported for random-coil models, except for Glu18 and Asp19, suggesting electrostatic interactions around these residues. Molecular modelling and simulation indicate an additional, significant role of nearby proline residues in determining the polypeptide conformational features and water accessibility in the region around Glu18, modulating the titration properties of these amino acids. In the other acidic residues of NUPR1, the small deviations of pKa values (compared to those expected for a random-coil) are likely due to electrostatic interactions with charged adjacent residues, which should be reduced at high NaCl concentrations. In fact, at high ionic strength, the pKa values of the aspartic residues were similar to those in a random coil, but there were still small differences for those of glutamic acids, probably due to hydrogen-bond formation. The overall findings suggest that local interactions and hydrophobic effects play a major role in determining the electrostatic features of NUPR1, whereas long-range charge contributions appear to be of lesser importance.

Site-Specific Internal Motions in GB1 Protein Microcrystals Revealed by 3D ²H-¹³C-¹³C Solid-State NMR Spectroscopy

(2)H quadrupolar line shapes deliver rich information about protein dynamics. A newly designed 3D (2)H-(13)C-(13)C solid-state NMR magic angle spinning (MAS) experiment is presented and demonstrated on the microcrystalline β1 immunoglobulin binding domain of protein G (GB1). The implementation of (2)H-(13)C adiabatic rotor-echo-short-pulse-irradiation cross-polarization (RESPIRATION CP) ensures the accuracy of the extracted line shapes and provides enhanced sensitivity relative to conventional CP methods. The 3D (2)H-(13)C-(13)C spectrum reveals (2)H line shapes for 140 resolved aliphatic deuterium sites. Motional-averaged (2)H quadrupolar parameters obtained from the line-shape fitting identify side-chain motions. Restricted side-chain dynamics are observed for a number of polar residues including K13, D22, E27, K31, D36, N37, D46, D47, K50, and E56, which we attribute to the effects of salt bridges and hydrogen bonds. In contrast, we observe significantly enhanced side-chain flexibility for Q2, K4, K10, E15, E19, N35, N40, and E42, due to solvent exposure and low packing density. T11, T16, and T17 side chains exhibit motions with larger amplitudes than other Thr residues due to solvent interactions. The side chains of L5, V54, and V29 are highly rigid because they are packed in the core of the protein. High correlations were demonstrated between GB1 side-chain dynamics and its biological function. Large-amplitude side-chain motions are observed for regions contacting and interacting with immunoglobulin G (IgG). In contrast, rigid side chains are primarily found for residues in the structural core of the protein that are absent from protein binding and interactions.

Structure of CrgA, a cell division structural and regulatory protein from Mycobacterium tuberculosis, in lipid bilayers

The 93-residue transmembrane protein CrgA in Mycobacterium tuberculosis is a central component of the divisome, a large macromolecular machine responsible for cell division. Through interactions with multiple other components including FtsZ, FtsQ, FtsI (PBPB), PBPA, and CwsA, CrgA facilitates the recruitment of the proteins essential for peptidoglycan synthesis to the divisome and stabilizes the divisome. CrgA is predicted to have two transmembrane helices. Here, the structure of CrgA was determined in a liquid-crystalline lipid bilayer environment by solid-state NMR spectroscopy. Oriented-sample data yielded orientational restraints, whereas magic-angle spinning data yielded interhelical distance restraints. These data define a complete structure for the transmembrane domain and provide rich information on the conformational ensembles of the partially disordered N-terminal region and interhelical loop. The structure of the transmembrane domain was refined using restrained molecular dynamics simulations in an all-atom representation of the same lipid bilayer environment as in the NMR samples. The two transmembrane helices form a left-handed packing arrangement with a crossing angle of 24° at the conserved Gly39 residue. This helix pair exposes other conserved glycine and alanine residues to the fatty acyl environment, which are potential sites for binding CrgA's partners such as CwsA and FtsQ. This approach combining oriented-sample and magic-angle spinning NMR spectroscopy in native-like lipid bilayers with restrained molecular dynamics simulations represents a powerful tool for structural characterization of not only isolated membrane proteins, but their complexes, such as those that form macromolecular machines.

ω-Turn: a novel β-turn mimic in globular proteins stabilized by main-chain to side-chain C−H···O interaction

Mimicry of structural motifs is a common feature in proteins. The 10-membered hydrogen-bonded ring involving the main-chain C − O in a β-turn can be formed using a side-chain carbonyl group leading to Asx-turn. We show that the N − H component of hydrogen bond can be replaced by a C(γ) -H group in the side chain, culminating in a nonconventional C − H···O interaction. Because of its shape this β-turn mimic is designated as ω-turn, which is found to occur ∼ three times per 100 residues. Three residues (i to i + 2) constitute the turn with the C − H···O interaction occurring between the terminal residues, constraining the torsion angles ϕi + 1, ψi + 1, ϕi + 2 and χ'1(i + 2) (using the interacting C(γ) atom). Based on these angles there are two types of ω-turns, each of which can be further divided into two groups. C(β) -branched side-chains, and Met and Gln have high propensities to occur at i + 2; for the last two residues the carbonyl oxygen may participate in an additional interaction involving the S and amino group, respectively. With Cys occupying the i + 1 position, such turns are found in the metal-binding sites. N-linked glycosylation occurs at the consensus pattern Asn-Xaa-Ser/Thr; with Thr at i + 2, the sequence can adopt the secondary structure of a ω-turn, which may be the recognition site for protein modification. Location between two β-strands is the most common occurrence in protein tertiary structure, and being generally exposed ω-turn may constitute the antigenic determinant site. It is a stable scaffold and may be used in protein engineering and peptide design.

OGlcNAcylation and phosphorylation have similar structural effects in α-helices: post-translational modifications as inducible start and stop signals in α-helices, with greater structural effects on threonine modification

OGlcNAcylation and phosphorylation are the major competing intracellular post-translational modifications of serine and threonine residues. The structural effects of both post-translational modifications on serine and threonine were examined within Baldwin model α-helical peptides (Ac-AKAAAAKAAAAKAAGY-NH₂ or Ac-YGAKAAAAKAAAAKAA-NH₂). At the N-terminus of an α-helix, both phosphorylation and OGlcNAcylation stabilized the α-helix relative to the free hydroxyls, with a larger induced structure for phosphorylation than for OGlcNAcylation, for the dianionic phosphate than for the monoanionic phosphate, and for modifications on threonine than for modifications on serine. Both phosphoserine and phosphothreonine resulted in peptides more α-helical than alanine at the N-terminus, with dianionic phosphothreonine the most α-helix-stabilizing residue here. In contrast, in the interior of the α-helix, both post-translational modifications were destabilizing with respect to the α-helix, with the greatest destabilization seen for threonine OGlcNAcylation at residue 5 and threonine phosphorylation at residue 10, with peptides containing either post-translational modification existing as random coils. At the C-terminus, both OGlcNAcylation and phosphorylation were destabilizing with respect to the α-helix, though the induced structural changes were less than in the interior of the α-helix. In general, the structural effects of modifications on threonine were greater than the effects on serine, because of both the lower α-helical propensity of Thr and the more defined induced structures upon modification of threonine than serine, suggesting threonine residues are particularly important loci for structural effects of post-translational modifications. The effects of serine and threonine post-translational modifications are analogous to the effects of proline on α-helices, with the effects of phosphothreonine being greater than those of proline throughout the α-helix. These results provide a basis for understanding the context-dependent structural effects of these competing protein post-translational modifications.

OGlcNAcylation and phosphorylation have opposing structural effects in tau: phosphothreonine induces particular conformational order

Phosphorylation and OGlcNAcylation are dynamic intracellular protein post-translational modifications that frequently are alternatively observed on the same serine and threonine residues. Phosphorylation and OGlcNAcylation commonly occur in natively disordered regions of proteins, and often have opposing functional effects. In the microtubule-associated protein tau, hyperphosphorylation is associated with protein misfolding and aggregation as the neurofibrillary tangles of Alzheimer's disease, whereas OGlcNAcylation stabilizes the soluble form of tau. A series of peptides derived from the proline-rich domain (residues 174-251) of tau was synthesized, with free Ser/Thr hydroxyls, phosphorylated Ser/Thr (pSer/pThr), OGlcNAcylated Ser/Thr, and diethylphosphorylated Ser/Thr. Phosphorylation and OGlcNAcylation were found by CD and NMR to have opposing structural effects on polyproline helix (PPII) formation, with phosphorylation favoring PPII, OGlcNAcylation opposing PPII, and the free hydroxyls intermediate in structure, and with phosphorylation structural effects greater than OGlcNAcylation. For tau196-209, phosphorylation and OGlcNAcylation had similar structural effects, opposing a nascent α-helix. Phosphomimic Glu exhibited PPII-favoring structural effects. Structural changes due to Thr phosphorylation were greater than those of Ser phosphorylation or Glu, with particular conformational restriction as the dianion, with mean (3)JαN = 3.5 Hz (pThr) versus 5.4 Hz (pSer), compared to 7.2, 6.8, and 6.2 Hz for Thr, Ser, and Glu, respectively, values that correlate with the backbone torsion angle ϕ. Dianionic phosphothreonine induced strong phosphothreonine amide protection and downfield amide chemical shifts (δmean = 9.63 ppm), consistent with formation of a stable phosphate-amide hydrogen bond. These data suggest potentially greater structural importance of threonine phosphorylation than serine phosphorylation due to larger induced structural effects.

Channelrhodopsins: a bioinformatics perspective

Channelrhodopsins are microbial-type rhodopsins that function as light-gated cation channels. Understanding how the detailed architecture of the protein governs its dynamics and specificity for ions is important, because it has the potential to assist in designing site-directed channelrhodopsin mutants for specific neurobiology applications. Here we use bioinformatics methods to derive accurate alignments of channelrhodopsin sequences, assess the sequence conservation patterns and find conserved motifs in channelrhodopsins, and use homology modeling to construct three-dimensional structural models of channelrhodopsins. The analyses reveal that helices C and D of channelrhodopsins contain Cys, Ser, and Thr groups that can engage in both intra- and inter-helical hydrogen bonds. We propose that these polar groups participate in inter-helical hydrogen-bonding clusters important for the protein conformational dynamics and for the local water interactions. This article is part of a Special Issue entitled: Retinal Proteins - You can teach an old dog new tricks.

The intrinsic conformational features of amino acids from a protein coil library and their applications in force field development

The local conformational (φ, ψ, χ) preferences of amino acid residues remain an active research area, which are important for the development of protein force fields. In this perspective article, we first summarize spectroscopic studies of alanine-based short peptides in aqueous solution. While most studies indicate a preference for the P(II) conformation in the unfolded state over α and β conformations, significant variations are also observed. A statistical analysis from various coil libraries of high-resolution protein structures is then summarized, which gives a more coherent view of the local conformational features. The φ, ψ, χ distributions of the 20 amino acids have been obtained from a protein coil library, considering both backbone and side-chain conformational preferences. The intrinsic side-chain χ(1) rotamer preference and χ(1)-dependent Ramachandran plot can be generally understood by combining the interaction of the side-chain Cγ/Oγ atom with two neighboring backbone peptide groups. Current all-atom force fields such as AMBER ff99sb-ILDN, ff03 and OPLS-AA/L do not reproduce these distributions well. A method has been developed by combining the φ, ψ plot of alanine with the influence of side-chain χ(1) rotamers to derive the local conformational features of various amino acids. It has been further applied to improve the OPLS-AA force field. The modified force field (OPLS-AA/C) reproduces experimental (3)J coupling constants for various short peptides quite well. It also better reproduces the temperature-dependence of the helix-coil transition for alanine-based peptides. The new force field can fold a series of peptides and proteins with various secondary structures to their experimental structures. MD simulations of several globular proteins using the improved force field give significantly less deviation (RMSD) to experimental structures. The results indicate that the local conformational features from coil libraries are valuable for the development of balanced protein force fields.

Ser/Thr motifs in transmembrane proteins: conservation patterns and effects on local protein structure and dynamics

We combined systematic bioinformatics analyses and molecular dynamics simulations to assess the conservation patterns of Ser and Thr motifs in membrane proteins, and the effect of such motifs on the structure and dynamics of α-helical transmembrane (TM) segments. We find that Ser/Thr motifs are often present in β-barrel TM proteins. At least one Ser/Thr motif is present in almost half of the sequences of α-helical proteins analyzed here. The extensive bioinformatics analyses and inspection of protein structures led to the identification of molecular transporters with noticeable numbers of Ser/Thr motifs within the TM region. Given the energetic penalty for burying multiple Ser/Thr groups in the membrane hydrophobic core, the observation of transporters with multiple membrane-embedded Ser/Thr is intriguing and raises the question of how the presence of multiple Ser/Thr affects protein local structure and dynamics. Molecular dynamics simulations of four different Ser-containing model TM peptides indicate that backbone hydrogen bonding of membrane-buried Ser/Thr hydroxyl groups can significantly change the local structure and dynamics of the helix. Ser groups located close to the membrane interface can hydrogen bond to solvent water instead of protein backbone, leading to an enhanced local solvation of the peptide.

A propensity scale for type II polyproline helices (PPII): aromatic amino acids in proline-rich sequences strongly disfavor PPII due to proline-aromatic interactions

Type II polyproline helices (PPII) are a fundamental secondary structure of proteins, common in globular and nonglobular regions and important in cellular signaling. We developed a propensity scale for PPII using a host-guest system with sequence Ac-GPPXPPGY-NH(2), where X represents any amino acid. We found that proline has the highest PPII propensity, but most other amino acids display significant PPII propensities. The PPII propensity of leucine was the highest of all propensities of non-proline residues. Alanine and residues with linear side chains displayed the next highest PPII propensities. Three classes of residues displayed lower PPII propensities: β-branched amino acids (Thr, Val, and Ile), short amino acids with polar side chains (Asn, protonated Asp, Ser, Thr, and Cys), and aromatic amino acids (Phe, Tyr, and Trp). tert-Leucine particularly disfavored PPII. The basis of the low PPII propensities of aromatic amino acids in this context was significant cis-trans isomerism, with proline-rich peptides containing aromatic residues exhibiting 45-60% cis amide bonds, due to Pro-cis-Pro-aromatic and aromatic-cis-Pro amide bonds.

Asparagine and glutamine differ in their propensities to form specific side chain-backbone hydrogen bonded motifs in proteins

Short range side chain-backbone hydrogen bonded motifs involving Asn and Gln residues have been identified from a data set of 1370 protein crystal structures (resolution ≤ 1.5 Å). Hydrogen bonds involving residues i - 5 to i + 5 have been considered. Out of 12,901 Asn residues, 3403 residues (26.4%) participate in such interactions, while out of 10,934 Gln residues, 1780 Gln residues (16.3%) are involved in these motifs. Hydrogen bonded ring sizes (C(n), where n is the number of atoms involved), directionality and internal torsion angles are used to classify motifs. The occurrence of the various motifs in the contexts of protein structure is illustrated. Distinct differences are established between the nature of motifs formed by Asn and Gln residues. For Asn, the most highly populated motifs are the C(10)(CO(δ)(i) …NH(i + 2)), C(13)(CO(δ)(i) …NH(i + 3)) and C(17)(N(δ)H(i) …CO(i - 4)) structures. In contrast, Gln predominantly forms C(16)(CO(ε)(i) …NH(i - 3)), C(12)(N(ε)H(i) …CO(i - 2)), C(15)(N(ε)H(i) …CO(i - 3)) and C(18)(N(ε)H(i) …CO(i - 4)) motifs, with only the C(18) motif being analogous to the Asn C(17) structure. Specific conformational types are established for the Asn containing motifs, which mimic backbone β-turns and α-turns. Histidine residues are shown to serve as a mimic for Asn residues in side chain-backbone hydrogen bonded ring motifs. Illustrative examples from protein structures are considered.

Why are polar residues within the membrane core evolutionary conserved?

Here, we present a study of polar residues within the membrane core of alpha-helical membrane proteins. As expected, polar residues are less frequent in the membrane than expected. Further, most of these residues are buried within the interior of the protein and are only rarely exposed to lipids. However, the polar groups often border internal water filled cavities, even if the rest of the sidechain is buried. A survey of their functional roles in known structures showed that the polar residues are often directly involved in binding of small compounds, especially in channels and transporters, but other functions including proton transfer, catalysis, and selectivity have also been attributed to these proteins. Among the polar residues histidines often interact with prosthetic groups in photosynthetic- and oxidoreductase-related proteins, whereas prolines often are required for conformational changes of the proteins. Indeed, the polar residues in the membrane core are more conserved than other residues in the core, as well as more conserved than polar residues outside the membrane. The reason is twofold; they are often (i) buried in the interior of the protein and (ii) directly involved in the function of the proteins. Finally, a method to identify which polar residues are present within the membrane core directly from protein sequences was developed. Applying the method to the set of all human membrane proteins the prediction indicates that polar residues were most frequent among active transporter proteins and GPCRs, whereas infrequent in families with few transmembrane regions, such as non-GPCR receptors.

Influence of side chain conformations on local conformational features of amino acids and implication for force field development

Statistical analysis of coil regions in protein structures has been used to obtain the local backbone phi, psi preferences of amino acids, which agree well with the NMR experiments of unfolded peptides and proteins. We analyzed the conformational features of amino acid residues in a restricted coil library of 4220 high-resolution protein crystal structures. In addition to Gly, Ala, and Pro, the phi, psi distribution (Ramachandran plot) of each amino acid is analyzed with respect to three side chain conformers: g+ (chi(1) approximately -60 degrees), g- (chi(1) approximately 60 degrees), and t (chi(1) approximately 180 degrees). The statistical study indicates that the effect of side chain conformations on phi, psi distributions is even greater than the effect of amino acid types. On the basis of the chi(1), phi, psi conformational preferences, the amino acids in addition to Gly, Pro, and Ala can be divided into five types: (1) ordinary amino acids, (2) Ser, (3) Asp and Asn, (4) Val and Ile, and (5) Thr, each with distinguished chi(1) rotamers. The alpha-helix, beta-sheet, and type-I beta-turn preferences of the different rotamers of various amino acid types can be captured by their intrinsic phi, psi preferences from our coil library. Molecular dynamics simulations of dipeptide Ac-X-NHMe and tetrapeptide Ac-A-X-A-NHMe models give nearly the same side chain rotamer distributions. However, for many amino acids, both OPLS-AA/L and AMBER-FF03 force fields give very different chi(1) rotamer distributions from the coil library. This may partially explain why dipeptide models sometimes cannot reproduce those of protein structures well. The current coil library analysis may be valuable in improving the force field for protein simulations.

Backbone statistical potential from local sequence-structure interactions in protein loops

Native proteins have been optimized by evolution simultaneously for structure and sequence. Structural databases reflect this interdependency. In this paper, we present a new statistical potential for a reduced backbone representation that has both structure and sequence characteristics as variables. We use information from structural data available in the Protein Coil Library, selected on the basis of resolution and refinement factor. In these structures, the nonlocal interactions are randomly distributed and, thus, average out in statistics, so structural propensities due to local backbone-based interactions can be studied separately. We collect data in the form of local sequence-specific phi-psi backbone dihedral pairs. From these data, we construct dihedral probability density functions (DPDFs) that quantify any adjacent phi-psi pair distribution in the context of all possible combinations of local residue types. We use a probabilistic analysis to deduce how the correlations encoded in the various DPDFs as well as in residue frequencies propagate along the sequence and can be cumulated in a statistical potential capable of efficiently scoring a loop by its backbone conformation and sequence only. Our potential is able to identify with high accuracy the native structure of a loop with a given sequence among possible alternative conformations from sets of well-constructed decoys. Conversely, the potential can also be used for sequence prediction problems and is shown to score the native sequence of a given loop structure among the most fit of the possible sequence combinations. Applications for both structure prediction and sequence design are discussed.

The importance of serine 161 in the sodium channel beta3 subunit for modulation of Na(V)1.2 gating

Voltage-gated sodium (Na) channels contribute to the regulation of cellular excitability due to their role in the generation and propagation of action potentials. They are composed of a pore-forming alpha subunit and are modulated by at least two of four distinct beta subunits (beta1-4). Recent studies have implicated a role for the intracellular domain of beta subunits in modulating Na channel gating and trafficking. In beta3, the intracellular domain contains a serine residue at position 161 that is replaced by an alanine in beta1. In this study, we have probed the functional importance of beta3S161 for modulating Na channel gating. Wild-type beta3 and point mutations beta3S161A or beta3S161E were individually co-expressed in HEK 293 cells stably expressing human Na(v)1.2. WTbeta3 expression increased Na current density, shifted steady-state inactivation in a depolarized direction, and accelerated the kinetics of recovery from inactivation of the Na current. Analogous effects were observed with beta3S161E co-expression. In contrast, beta3S161A abolished the shifts in steady-state inactivation and recovery from inactivation of the Na current, but did increase Na current density. Immunocytochemistry and Western blot experiments demonstrate membrane expression of WTbeta3, beta3S161E, and beta3S161A, suggesting that the differences in Na channel gating were not due to disruptions in beta subunit trafficking. These studies suggest that modification of beta3S161 may be important in modulating Na-channel gating.

Kinking the coiled coil--negatively charged residues at the coiled-coil interface

The coiled coil is one of the most common protein-structure motifs. It is believed to be adopted by 3-5% of all amino acids in proteins. It comprises two or more alpha-helical chains wrapped around one another. The sequences of most coiled coils are characterized by a seven-residue (heptad) repeat, denoted (abcdefg)(n). Residues at the a and d positions define the helical interface (core) and are usually hydrophobic, though about 20% are polar or charged. We show that parallel coiled-coils have a unique pattern of their negatively charged residues at the core positions: aspartic acid is excluded from these positions while glutamic acid is not. In contrast the antiparallel structures are more permissive in their amino acid usage. We show further, and for the first time, that incorporation of Asp but not Glu into the a positions of a parallel coiled coil creates a flexible hinge and that the maximal hinge angle is being directly related to the number of incorporated mutations. These new computational and experimental observations will be of use in improving protein-structure predictions, and as rules to guide rational design of novel coiled-coil motifs and coiled coil-based materials.

EPR spectroscopy and theoretical study of gamma-irradiated asparagine and aspartic acid in solid state

Aspartic acid (Asp) and asparagine (Asn) are vulnerable amino acids. One-electron addition or withdrawal reactions initiate many deleterious processes involving these amino acids. To study these redox processes we have irradiated by gamma-rays asparagine or aspartic acid in the solid state. The nature of the resulting free radicals was determined by electron paramagnetic resonance (EPR) and by calculations using DFT methods in various environments. Reactions initiated by electron transfer are different for both amino acids: Asn anion loses hydrogen atom whereas the cation undergoes decarboxylation. Conversely, Asp cation loses hydrogen atom from amine group, which triggers decarboxylation.

Statistical characterization of salt bridges in proteins

The structure and folding mechanism of a given protein are determined by many factors, including the electrostatic interactions between charged residues of protein molecules known in general as salt bridges. In this study, analyses were conducted on 10,370 salt bridges in 2017 proteins and the results compared to previous statistical surveys of 36 protein structures. Although many of the general trends remained consistent with other studies, more detailed information was illuminated by the larger dataset. In particular, it was shown that there is a strong correlation between secondary structure and salt bridge formation, and that salt bridges display preferential formation in an environment of about 30% solvent accessible surface area.

Pairwise Coupling in an Arg-Phe-Met Triplet Stabilizes α-Helical Peptide via Shared Rotamer Preferences

The hydrophobic Arg-Phe and Phe-Met side chain interactions stabilize the alpha-helix by -0.29 and -0.59 kcal/mol, respectively, when placed i, i + 4 in an alanine-based peptide. When both interactions are present simultaneously, however, they stabilize the helix by an additional -0.75 kcal/mol, nearly as much as the sum of its parts. We attribute this coupling to a shared rotamer preference, as the central Phe is t in both bonds. The energetic cost of restricting the Phe residue into a t conformation is only paid once in the triplet, rather than twice when the interactions are separate. Coupling is thus demonstrated to have large effects on protein stability.

Mimicry by asx- and ST-turns of the four main types of beta-turn in proteins

Hydrogen-bonded beta-turns in proteins occur in four categories: type I (the most common), type II, type II', and type I'. Asx-turns resemble beta-turns, in that both have an NH. . .OC hydrogen bond forming a ring of 10 atoms. Serine and threonine side chains also commonly form hydrogen-bonded turns, here called ST-turns. Asx-turns and ST-turns can be categorized into four classes, based on side chain rotamers and the conformation of the central turn residue, which are geometrically equivalent to the four types of beta-turns. We propose asx- and ST-turns be named using the type I, II, I', and II' beta-turn nomenclature. Using this, the frequency of occurrence of both asx- and ST-turns is: type II' > type I > type II > type I', whereas for beta-turns it is type I > type II > type I' > type II'. Almost all type II asx-turns occur as a recently described three residue feature named an asx-nest.

Analysis of side-chain rotamers in transmembrane proteins

We measured the frequency of side-chain rotamers in 14 alpha-helical and 16 beta-barrel membrane protein structures and found that the membrane environment considerably perturbs the rotamer frequencies compared to soluble proteins. Although there are limited experimental data, we found statistically significant changes in rotamer preferences depending on the residue environment. Rotamer distributions were influenced by whether the residues were lipid or protein facing, and whether the residues were found near the N- or C-terminus. Hydrogen-bonding interactions with the helical backbone perturbs the rotamer populations of Ser and His. Trp and Tyr favor side-chain conformations that allow their side chains to extend their polar atoms out of the membrane core, thereby aligning the side-chain polarity gradient with the polarity gradient of the membrane. Our results demonstrate how the membrane environment influences protein structures, providing information that will be useful in the structure prediction and design of transmembrane proteins.

Structure-based Design of Selective and Potent Inhibitors of Protein-tyrosine Phosphatase β

Protein-tyrosine phosphatases (PTPs) are considered important therapeutic targets because of their pivotal role as regulators of signal transduction and thus their implication in several human diseases such as diabetes, cancer, and autoimmunity. In particular, PTP1B has been the focus of many academic and industrial laboratories because it was found to be an important negative regulator of insulin and leptin signaling, and hence a potential therapeutic target in diabetes and obesity. As a result, significant progress has been achieved in the design of highly selective and potent PTP1B inhibitors. In contrast, little attention has been given to other potential drug targets within the PTP family. Guided by x-ray crystallography, molecular modeling, and enzyme kinetic analyses with wild type and mutant PTPs, we describe the development of a general, low molecular weight, non-peptide, non-phosphorus PTP inhibitor into an inhibitor that displays more than 100-fold selectivity for PTPbeta over PTP1B. Of note, our structure-based design principles, which are based on extensive bioinformatics analyses of the PTP family, are general in nature. Therefore, we anticipate that this strategy, here applied to PTPbeta, in principle can be used in the design and development of selective inhibitors of many, if not most PTPs.

Mutational analysis of different regions in the coxsackievirus 2B protein: requirements for homo-multimerization, membrane permeabilization, subcellular localization, and virus replication

The coxsackievirus 2B protein is a small hydrophobic protein (99 amino acids) that increases host cell membrane permeability, possibly by forming homo-multimers that build membrane-integral pores. Previously, we defined the functional role of the two hydrophobic regions HR1 and HR2. Here, we investigated the importance of regions outside HR1 and HR2 for multimerization, increasing membrane permeability, subcellular localization, and virus replication through analysis of linker insertion and substitution mutants. From these studies, the following conclusions could be drawn. (i) The hydrophilic region ((58)RNHDD(62)) between HR1 and HR2 is critical for multimerization and increasing membrane permeability. Substitution analysis of Asn(61) and Asn(62) demonstrated the preference for short polar side chains (Asp, Asn), residues that are often present in turns, over long polar side chains (Glu, Gln). This finding supports the idea that the hydrophilic region is involved in pore formation by facilitating a turn between HR1 and HR2 to reverse chain direction. (ii) Studies undertaken to define the downstream boundary of HR2 demonstrated that the aromatic residues Trp(80) and Trp(82), but not the positively charged residues Arg(81), Lys(84), and Lys(86) are important for increasing membrane permeability. (iii) The N terminus is not required for multimerization but does contribute to the membrane-active character of 2B. (iv) The subcellular localization of 2B does not rely on regions outside HR1 and HR2 and does not require multimerization. (v) Virus replication requires both the membrane-active character and an additional function of 2B that is not connected to this activity.

Thermodynamic Dissection of the Binding Energetics of Proline-rich Peptides to the Abl-SH3 Domain: Implications for Rational Ligand Design

The inhibition of the interactions between SH3 domains and their targets is emerging as a promising therapeutic strategy. To date, rational design of potent ligands for these domains has been hindered by the lack of understanding of the origins of the binding energy. We present here a complete thermodynamic analysis of the binding energetics of the p41 proline-rich decapeptide (APSYSPPPPP) to the SH3 domain of the c-Abl oncogene. Isothermal titration calorimetry experiments have revealed a thermodynamic signature for this interaction (very favourable enthalpic contributions opposed by an unfavourable binding entropy) inconsistent with the highly hydrophobic nature of the p41 ligand and the Abl-SH3 binding site. Our structural and thermodynamic analyses have led us to the conclusion, having once ruled out any possible ionization events or conformational changes coupled to the association, that the establishment of a complex hydrogen-bond network mediated by water molecules buried at the binding interface is responsible for the observed thermodynamic behaviour. The origin of the binding energetics for proline-rich ligands to the Abl-SH3 domain is further investigated by a comparative calorimetric analysis of a set of p41-related ligands. The striking effects upon the enthalpic and entropic contributions provoked by conservative substitutions at solvent-exposed positions in the ligand confirm the complexity of the interaction. The implications of these results for rational ligand design are discussed.

Folding of a highly conserved diverging turn motif from the SH3 domain

Recent NMR structural characterization studies showed that a seven-residue segment (FKKGERL) from the src SH3 domain adopts the nativelike diverging type II beta-turn in aqueous solution in support of the prediction based on the I-sites library of sequence structural motifs. We study the conformational variability and folding/unfolding thermodynamics of this peptide in explicit solvent using replica-exchange molecular dynamics simulations, which greatly enhances the sampling of the conformational space. This peptide samples three main free energy basins (nativelike, intermediate, and unfolded) separated by small barriers. The nativelike basin is fractionally populated (DeltaG(300K) = 0.4 kcal/mol) with structures that satisfy a subset of the NMR-derived constraints. The intrinsic stability of the diverging turn is examined in relationship to the nature of three specific contacts: a turn-hydrogen bond, a mainchain-to-sidechain hydrogen bond, and an end-to-end hydrophobic contact. We have carried out simulations of mutants at the highly conserved GE positions in the sequence. The mutation E5D destabilizes the isolated diverging turn motif, contrary to the observation that this mutation stabilizes the fyn SH3 domain. The G4T mutation also destabilizes the isolated diverging turn; however, the extent of destabilization is smaller than that of the reverse mutation in the drk SH3.

Rotamer libraries in the 21st century

Rotamer libraries are widely used in protein structure prediction, protein design, and structure refinement. As the size of the structure data base has increased rapidly in recent years, it has become possible to derive well-refined rotamer libraries using strict criteria for data inclusion and for studying dependence of rotamer populations and dihedral angles on local structural features.

Exploring the role of alanine in the structure of the Lac repressor tetramerization domain, a ferritin-like Alacoil

We are interested in the determinants that specify the structure of antiparallel coiled coils. Antiparallel coiled coils often contain alanine as an important interfacial packing residue; structures containing alanine at certain well-defined positions in the heptad-repeating unit are referred to as Alacoils. Two types have been identified, containing alanine at either the g position of the heptad repeating unit (defined as the d position in the Richardson nomenclature), referred to as a rop-like Alacoil, or the e position (a position in the Richardson nomenclature), referred to as a ferritin-like Alacoil. The Lac repressor tetramerization domain forms an antiparallel four-chain coiled coil, which falls into the second class of Alacoils based on recent crystal structures. The role of alanine in such structures has not yet been explored experimentally. We test the importance of alanine at the e positions on the oligomeric state and stability of the isolated coiled-coil domain of Lac repressor by testing the effect of mutations at this position. We find that mutation to leucine is tolerated and its moderately stabilizing effect is most likely a consequence of plasticity of this motif. The effects on stability of the mutations to either serine or glutamine can be largely accounted for by helix propensity differences between these residues and alanine. The ability of the helices to adjust to such mutations, while maintaining the basic fold of this coiled coil, was further tested by making the same changes at the more highly exposed g position. Leucine at the g positions also causes an increase in stability, presumably by subtle rearrangement of the helices to allow partial desolvation of this side-chain.