The interrelationships of side-chain and main-chain conformations in proteins

Rotavirus A Infection Prevalence and Spatio-Temporal Genotype Shift among Under-Five Children in Amhara National Regional State, Ethiopia: A Multi-Center Cross-Sectional Study

Background: Globally, rotavirus (RV) A (RVA) is the most common cause of severe and sometimes fatal diarrhea in young children. It is also the major cause of acute gastroenteritis among children in Ethiopia. Currently, the WHO has prequalified four RVA vaccines for universal childhood immunization. Ethiopia introduced the monovalent Rotarix vaccine into its national immunization program in 2013. Since then, only a few studies on the burden and genotype distribution of RVA infection post-vaccine introduction have been conducted (mostly at sentinel surveillance sites). Therefore, this study aimed to assess RVA prevalence and genotype distribution among children under five years in Ethiopia (February 2021-December 2022). Methods: This multi-center hospital-based cross-sectional study involved 537 diarrheic children under-five years old. Rotavirus A detection was conducted using a one-step reverse-transcriptase polymerase chain reaction (RT-PCR). Genotyping was conducted by Sanger sequencing of the VP7 (complete) and VP4 (partial) genes. Descriptive analysis and Pearson's chi-squared test were carried out using SPSS version 29. Phylogenetic analysis with 1000 bootstrap replicates was performed using MEGA version 11 software. Statistical significance was set at p < 0.05 for all analyses. Results: The prevalence of RVA infection among diarrheic children was 17.5%. The most prevalent G-types identified were G3 (37%), the previously uncommon G12 (28%), and G1 (20%). The predominant P-types were P[8] (51%), P[6] (29%), and P[4] (14%). The three major G/P combinations observed were G3P[8] (32.8%), G12P[6] (28.4%), and G1P[8] (19.4%). Phylogenetic analysis revealed clustering of Ethiopian strains with the globally reported strains. Many strains exhibited amino acid differences in the VP4 (VP8* domain) and VP7 proteins compared to vaccine strains, potentially affecting virus neutralization. Conclusions: Despite the high RVA vaccination rate, the prevalence of RVA infection remains significant among diarrheic children in Ethiopia. There is an observable shift in circulating RVA genotypes from G1 to G3, alongside the emergence of unusual G/P genotype combinations such as G9P[4]. Many of these circulating RVA strains have shown amino acid substitutions that may allow for neutralization escape. Therefore, further studies are warranted to comprehend the emergence of these unusual RVA strains and the diverse factors influencing the vaccine's diminished effectiveness in developing countries.

Structural study for substrate recognition of human N-terminal glutamine amidohydrolase 1 in the arginine N-degron pathway

The N-degron pathway determines the half-life of proteins by selectively destabilizing the proteins bearing N-degrons. N-terminal glutamine amidohydrolase 1 (NTAQ1) plays an essential role in the arginine N-degron (Arg/N-degron) pathway as an initializing enzyme via the deamidation of the N-terminal (Nt) glutamine (Gln). However, the Nt-serine-bound conformation of hNTAQ1 according to the previously identified crystal structure suggests the possibility of other factors influencing the recognition of Nt residues by hNTAQ1. Hence, in the current study, we aimed to further elucidate the substrate recognition of hNTAQ1; specifically, we explored 12 different substrate-binding conformations of hNTAQ1 depending on the subsequent residue of Nt-Gln. Results revealed that hNTAQ1 primarily interacts with the protein Nt backbone, instead of the side chain, for substrate recognition. Here, we report that the Nt backbone of proteins appears to be a key component of hNTAQ1 function and is the main determinant of substrate recognition. Moreover, not all second residues from Nt-Gln, but rather distinctive and charged residues, appeared to aid in detecting substrate recognition. These new findings define the substrate-recognition process of hNTAQ1 and emphasize the importance of the subsequent Gln residue in the Nt-Gln degradation system. Our extensive structural and biochemical analyses provide insights into the substrate specificity of the N-degron pathway and shed light on the mechanism underlying hNTAQ1 substrate recognition. An improved understanding of the protein degradation machinery could aid in developing therapies to promote overall health through enhanced protein regulation, such as targeted protein therapies.

Protein feature engineering framework for AMPylation site prediction

AMPylation is a biologically significant yet understudied post-translational modification where an adenosine monophosphate (AMP) group is added to Tyrosine and Threonine residues primarily. While recent work has illuminated the prevalence and functional impacts of AMPylation, experimental identification of AMPylation sites remains challenging. Computational prediction techniques provide a faster alternative approach. The predictive performance of machine learning models is highly dependent on the features used to represent the raw amino acid sequences. In this work, we introduce a novel feature extraction pipeline to encode the key properties relevant to AMPylation site prediction. We utilize a recently published dataset of curated AMPylation sites to develop our feature generation framework. We demonstrate the utility of our extracted features by training various machine learning classifiers, on various numerical representations of the raw sequences extracted with the help of our framework. Tenfold cross-validation is used to evaluate the model's capability to distinguish between AMPylated and non-AMPylated sites. The top-performing set of features extracted achieved MCC score of 0.58, Accuracy of 0.8, AUC-ROC of 0.85 and F1 score of 0.73. Further, we elucidate the behaviour of the model on the set of features consisting of monogram and bigram counts for various representations using SHapley Additive exPlanations.

Data-Efficient Generation of Protein Conformational Ensembles with Backbone-to-Side-Chain Transformers

Excitement at the prospect of using data-driven generative models to sample configurational ensembles of biomolecular systems stems from the extraordinary success of these models on a diverse set of high-dimensional sampling tasks. Unlike image generation or even the closely related problem of protein structure prediction, there are currently no data sources with sufficient breadth to parametrize generative models for conformational ensembles. To enable discovery, a fundamentally different approach to building generative models is required: models should be able to propose rare, albeit physical, conformations that may not arise in even the largest data sets. Here we introduce a modular strategy to generate conformations based on "backmapping" from a fixed protein backbone that (1) maintains conformational diversity of the side chains and (2) couples the side-chain fluctuations using global information about the protein conformation. Our model combines simple statistical models of side-chain conformations based on rotamer libraries with the now ubiquitous transformer architecture to sample with atomistic accuracy. Together, these ingredients provide a strategy for rapid data acquisition and hence a crucial ingredient for scalable physical simulation with generative neural networks.

Design rules for reciprocal coupling in chemically fueled assembly

Biology regulates the function and assembly of proteins through non-equilibrium reaction cycles. Reciprocally, the assembly of proteins can influence the reaction rates of these cycles. Such reciprocal coupling between assembly and reaction cycle is a prerequisite for behavior like dynamic instabilities, treadmilling, pattern formation, and oscillations between morphologies. While assemblies regulated by chemical reaction cycles gained traction, the concept of reciprocal coupling is under-explored. In this work, we provide two molecular design strategies to tweak the degree of reciprocal coupling between the assembly and reaction cycle. The strategies involve spacing the chemically active site away from the assembly or burying it into the assembly. We envision that design strategies facilitate the creation of reciprocally coupled and, by extension, dynamic supramolecular materials in the future.

Development of Orally Bioavailable Peptides Targeting an Intracellular Protein: From a Hit to a Clinical KRAS Inhibitor

Cyclic peptides as a therapeutic modality are attracting a lot of attention due to their potential for oral absorption and accessibility to intracellular tough targets. Here, starting with a drug-like hit discovered using an mRNA display library, we describe a chemical optimization that led to the orally available clinical compound known as LUNA18, an 11-mer cyclic peptide inhibitor for the intracellular tough target RAS. The key findings are as follows: (i) two peptide side chains were identified that each increase RAS affinity over 10-fold; (ii) physico-chemical properties (PCP) including Clog P can be adjusted by side-chain modification to increase membrane permeability; (iii) restriction of cyclic peptide conformation works effectively to adjust PCP and improve bio-activity; (iv) cellular efficacy was observed in peptides with a permeability of around 0.4 × 10^-6 cm/s or more in a Caco-2 permeability assay; and (v) while keeping the cyclic peptide's main-chain conformation, we found one example where the RAS protein structure was changed dramatically through induced-fit to our peptide side chain. This study demonstrates how the chemical optimization of bio-active peptides can be achieved without scaffold hopping, much like the processes for small molecule drug discovery that are guided by Lipinski's rule of five. Our approach provides a versatile new strategy for generating peptide drugs starting from drug-like hits.

The relevance of short peptides for an understanding of unfolded and intrinsically disordered proteins

Over the last thirty years the unfolded state of proteins has attracted considerable interest owing to the discovery of intrinsically disordered proteins which perform a plethora of functions despite resembling unfolded proteins to a significant extent. Research on both, unfolded and disordered proteins has revealed that their conformational properties can deviate locally from random coil behavior. In this context results from work on short oligopeptides suggest that individual amino acid residues sample the sterically allowed fraction of the Ramachandran plot to a different extent. Alanine has been found to exhibit a peculiarity in that it has a very high propensity for adopting polyproline II like conformations. This Perspectives article reviews work on short peptides aimed at exploring the Ramachandran distributions of amino acid residues in different contexts with experimental and computational means. Based on the thus provided overview the article discussed to what extent short peptides can serve as tools for exploring unfolded and disordered proteins and as benchmarks for the development of a molecular dynamics force field.

Improved NMR-data-compliant protein structure modeling captures context-dependent variations and expands the scope of functional inference

Nuclear magnetic resonance (NMR) spectroscopy can reveal conformational states of a protein in physiological conditions. However, sparsely available NMR data for a protein with large degrees of freedom can introduce structural artifacts in the built models. Currently used state-of-the-art methods deriving protein structure and conformation from NMR deploy molecular dynamics (MD) coupled with simulated annealing for building models. We provide an alternate graph-based modeling approach, where we first build substructures from NMR-derived distance-geometry constraints combined in one shot to form the core structure. The remaining molecule with inadequate data is modeled using a hybrid approach respecting the observed distance-geometry constraints. One-shot structure building is rarely undertaken for large and sparse data systems, but our data-driven bottom-up approach makes this uniquely feasible by suitable partitioning of the problem. A detailed comparison of select models with state-of-art methods reveals differences in the secondary structure regions wherein the correctness of our models is confirmed by NMR data. Benchmarking of 106 protein-folds covering 38-282 length structures shows minimal experimental-constraint violations while conforming to other structure quality parameters such as the proper folding, steric clash, and torsion angle violation based on Ramachandran plot criteria. Comparative MD studies using select protein models from a state-of-art method and ours under identical experimental parameters reveal distinct conformational dynamics that could be attributed to protein structure-function. Our work is thus useful in building enhanced NMR-evidence-based models that encapsulate the contextual secondary and tertiary structure variations present during the experimentation and expand the scope of functional inference.

Engagement with the TCR induces plasticity in antigenic ligands bound to MHC class I and CD1 molecules

Complementarity-determining regions (CDRs) of αβ T-cell receptors (TCRs) sense peptide-bound MHC (pMHC) complexes via chemical interactions, thereby mediating antigen specificity and MHC restriction. Flexible finger-like movement of CDR loops contributes to the establishment of optimal interactions with pMHCs. In contrast, peptide ligands captured in MHC molecules are considered more static because of the rigid hydrogen-bond network that stabilizes peptide ligands in the antigen-binding groove of MHC molecules. An array of crystal structures delineating pMHC complexes in TCR-docked and TCR-undocked forms is now available, which enables us to assess TCR engagement-induced conformational changes in peptide ligands. In this short review, we overview conformational changes in MHC class I-bound peptide ligands upon TCR docking, followed by those for CD1-bound glycolipid ligands. Finally, we analyze the co-crystal structure of the TCR:lipopeptide-bound MHC class I complex that we recently reported. We argue that TCR engagement-induced conformational changes markedly occur in lipopeptide ligands, which are essential for exposure of a primary T-cell epitope to TCRs. These conformational changes are affected by amino acid residues, such as glycine, that do not interact directly with TCRs. Thus, ligand recognition by specific TCRs involves not only T-cell epitopes but also non-epitopic amino acid residues. In light of their critical function, we propose to refer to these residues as non-epitopic residues affecting ligand plasticity and antigenicity (NR-PA).

Proline-Proline Dyad in the Fusion Peptide of the Murine β-Coronavirus Spike Protein's S2 Domain Modulates Its Neuroglial Tropism

The β-Coronavirus mouse hepatitis virus (MHV-A59)-RSA59 has a patent stretch of fusion peptide (FP) containing two consecutive central prolines (PP) in the S2 domain of the Spike protein. Our previous studies compared the PP-containing fusogenic-demyelinating strain RSA59(PP) to its one proline-deleted mutant strain RSA59(P) and one proline-containing non-fusogenic non-demyelinating parental strain RSMHV2(P) to its one proline inserted mutant strain RSMHV2(PP). These studies highlighted the crucial role of PP in fusogenicity, hepato-neuropathogenesis, and demyelination. Computational studies combined with biophysical data indicate that PP at the center of the FP provides local rigidity while imparting global fluctuation to the Spike protein that enhances the fusogenic properties of RSA59(PP) and RSMHV2(PP). To elaborate on the understanding of the role of PP in the FP of MHV, the differential neuroglial tropism of the PP and P mutant strains was investigated. Comparative studies demonstrated that PP significantly enhances the viral tropism for neurons, microglia, and oligodendrocytes. PP, however, is not essential for viral tropism for either astroglial or oligodendroglial precursors or the infection of meningeal fibroblasts in the blood-brain and blood-CSF barriers. PP in the fusion domain is critical for promoting gliopathy, making it a potential region for designing antivirals for neuro-COVID therapy.

Conformational entropy limits the transition from nucleation to elongation in amyloid aggregation

The formation of β-sheet rich amyloid fibrils in Alzheimer's disease and other neurodegenerative disorders is limited by a slow nucleation event. To understand the initial formation of β-sheets from disordered peptides, we used all-atom simulations to parameterize a lattice model that treats each amino acid as a binary variable with β and non-β states. We show that translational and conformational entropy give the nascent β-sheet an anisotropic surface tension which can be used to describe the nucleus with two-dimensional Classical Nucleation Theory. Since translational entropy depends on concentration, the aspect ratio of the critical β-sheet changes with protein concentration. Our model explains the transition from the nucleation phase to elongation as the point where the β-sheet core becomes large enough to overcome the conformational entropy cost to straighten the terminal molecule. At this point the β-strands in the nucleus spontaneously elongate, which results in a larger binding surface to capture new molecules. These results suggest that nucleation is relatively insensitive to sequence differences in co-aggregation experiments because the nucleus only involves a small portion of the peptide.

Arginine-containing dipeptides decrease affinity of gut trypsins and compromise soybean pest development

The design and production of molecules capable of mimicking the binding or/and functional sites of proteins inhibitors represent a promising strategy for the exploration and modulation of gut trypsin function in insect pests, specifically Lepidoptera. Here, for the first time, we characterized the trypsin activity present in the gut, performance and development of Anticarsia gemmatalis (Lepidoptera: Noctuidae) larvae when exposed to arginine-containing dipeptides. In silico assessment showed that arginine-containing dipeptides have a greater affinity for the active site of A. gemmatalis trypsins than lysine-containing peptides due to the presence of the double-charged guanidinium group that enhances the interaction at the S1 subsite of trypsins. Furthermore, the inhibitory and anti-insect potential of the peptides was demonstrated through kinetic and larval life cycle parameters, respectively. These dipeptides showed structural stability, binding to the active site, corroborated in vitro (competitive inhibition), and significant reduction of trypsin enzyme activity in the gut, survival, and weight of the A. gemmatalis larvae. Our findings reinforce the idea that small peptides are promising candidates for lepidopteran pest management. The optimization of DI2 and DI1 peptides, enhancing uptake and affinity to trypsins, may turn the use of these molecules feasible in agriculture.

Two Consecutive Prolines in the Fusion Peptide of Murine β-Coronavirus Spike Protein Predominantly Determine Fusogenicity and May Be Essential but Not Sufficient to Cause Demyelination

Combined in silico, in vitro, and in vivo comparative studies between isogenic-recombinant Mouse-Hepatitis-Virus-RSA59 and its proline deletion mutant, revealed a remarkable contribution of centrally located two consecutive prolines (PP) from Spike protein fusion peptide (FP) in enhancing virus fusogenic and hepato-neuropathogenic potential. To deepen our understanding of the underlying factors, we extend our studies to a non-fusogenic parental virus strain RSMHV2 (P) with a single proline in the FP and its proline inserted mutant, RSMHV2 (PP). Comparative in vitro and in vivo studies between virus strains RSA59(PP), RSMHV2 (P), and RSMHV2 (PP) in the FP demonstrate that the insertion of one proline significantly resulted in enhancing the virus fusogenicity, spread, and consecutive neuropathogenesis. Computational studies suggest that the central PP in Spike FP induces a locally ordered, compact, and rigid structure of the Spike protein in RSMHV2 (PP) compared to RSMHV2 (P), but globally the Spike S2-domain is akin to the parental strain RSA59(PP), the latter being the most flexible showing two potential wells in the energy landscape as observed from the molecular dynamics studies. The critical location of two central prolines of the FP is essential for fusogenicity and pathogenesis making it a potential site for designing antiviral.

Thermodynamics of globular protein native structure

The temperature dependence of the partial heat capacity of the native protein structure in an aqueous solution has been analyzed. It is shown that the strictly linear temperature dependence is due to the contributions of the vibrational and conformational components, which indicates volume consistensy and the absence of conformational transitions up to the main two-state transition. The two-level structural and functional organization of the protein three-dimensional structure are considered in relation to the energy and conformational entropy properties in accordance with the principles of the organization of the protein macromolecule.Communicated by Ramaswamy H. Sarma.

Study on the Potential Distribution of Leptinotarsa decemlineata and Its Natural Enemy Picromerus bidens Under Climate Change

The Colorado potato beetle (CPB), scientifically known as Leptinotarsa decemlineata, is a destructive quarantine pest that has invaded more than 40 countries and regions worldwide. It causes a 20–100% reduction in plant production, leading to severe economic losses. Picromerus bidens L. is a predatory insect that preys on CPB. This study used the MaxEnt model to predict the current and future potential distribution areas of CPB and P. bidens under different climatic scenarios to determine the possibility of using P. bidens as a natural enemy to control CPB. The possible introduction routes of CPB and P. bidens were subsequently predicted by combining their potential distribution with the current distribution of airports and ports. Notably, the potential distribution area of P. bidens was similar to that of CPB, suggesting that P. bidens could be used as a natural enemy to control CPB. Future changes in the suitable growth areas of CPB under different climate scenarios increased and decreased but were insignificant, while those of P. bidens decreased. Consequently, a reduction of the suitable habitats of P. bidens may cause a decrease in its population density, leading to a lack of adequate and timely prevention and control of invasive pests. Active measures should thus be enacted to minimize global warming and protect biodiversity. This study provides a theoretical basis and data support for early warning, monitoring, and control of the CPB spread.

Remarkable diastereomeric effect on thermoresponsive behavior of polyurethane based on lysine and tartrate ester derivatives

This study describes the long-distance diastereomeric effect on thermoresponsive properties in water-soluble diastereomeric polyurethanes (PUs) composed of an l-lysine ethyl ester diisocyanate and a trimethylene glycol l-/d-tartrate ester, which have differences in spatial arrangements of the ethyl esters in the mirror image. The PUs based on l-lysine and l-/d-tartrate ester, named l-PU and d-PU, were synthesized with various number average molecular weights from 4700 to 13 100. In turbidimetry, l-PU showed a steep phase transition from 100%T to 0%T within about 10 °C at 4 g L⁻¹, whereas d-PU did not change completely to 0%T transmittance even at 80 °C at 4 g L⁻¹. In addition, the thermoresponsive properties of l-PU were less affected by concentration than those of d-PU. This long-distance diastereomeric effect on thermoresponsive behavior between l-PU and d-PU appeared in common among 6 samples with 4700 to 13 100 number average molecular weight. In the dynamic light scattering experiments at each transmittance, the hydrodynamic diameter (D_h) of l-PU increased up to 1000 nm, while the D_h of d-PU remained almost at 200–300 nm. The C Created by potrace 1.16, written by Peter Selinger 2001-2019 ]]> O stretching vibration of FT-IR spectra showed that d-PU has more hydrogen-bonded ester groups than L-PU. Thus, we speculated that the difference in the retention of polymer chains in the micelle to promote intermicellar bridging generates the long-distance diastereomeric effect.

The long-distance diastereomeric effect on thermoresponsive properties in a polyurethane system consisting of chiral monomers was reported.

Evolutionary selectivity of amino acid is inspired from the enhanced structural stability and flexibility of the folded protein

AIM

The present study attempts to decipher the site-specific amino acid alterations at certain positions experiencing preferential selectivity and their effect on proteins' stability and flexibility. The study examines the selection preferences by considering pair-wise non-bonded interaction energies of adjacent and interacting amino acids present at the interacting site, along with their evolutionary history.

MATERIALS AND METHODS

For the study, variations in the interacting residues of spike protein (S-Protein) receptor-binding domain (RBD) of different coronaviruses were examined. The MD simulation trajectory analysis revealed that, though all the variants studied were structurally stable at their native and bound confirmations, the RBD of 2019-nCoV/SARS-CoV-2 was found to be more flexible and more dynamic. Furthermore, a noticeable change observed in the non-bonded interaction energies of the amino acids interacting with the receptor corroborated their selection at respective positions.

KEY FINDINGS

The conformational changes exerted by the altered amino acids could be the reason for a broader range of interacting receptors among the selected proteins.

SIGNIFICANCE

The results envisage a strong indication that the residue selection at certain positions is governed by a well-orchestrated feedback mechanism, which follows increased stability and flexibility in the folded structure compared to its evolutionary predecessor.

Spike protein fusion loop controls SARS-CoV-2 fusogenicity and infectivity

The high SARS-CoV-2 reproductive number driving the COVID-19 pandemic has been a mystery. Our recent in vitro, and in vivo coronaviral pathogenesis studies involving Mouse Hepatitis Virus (MHV-A59) suggest a crucial role for a small host membrane-virus contact initiator region of the Spike protein, called the fusion peptide that enhances the virus fusogenicity and infectivity. Here I study the Spike from five human β-coronaviruses (HCoV) including the SARS-CoV-2, and MHV-A59 for comparison. The structural and dynamics analyses of the Spike show that its fusion loop spatially organizes three fusion peptides contiguous to each other to synergistically trigger the virus-host membrane fusion process. I propose a Contact Initiation Model based on the architecture of the Spike quaternary structure that explains the obligatory participation of the fusion loop in the initiation of the host membrane contact for the virus fusion process. Among all the HCoV Spikes in this study, SARS-CoV-2 has the most hydrophobic surface and the extent of hydrophobicity correlates with the reproductive number and infectivity of the other HCoV. Comparison between results from standard and replica exchange molecular dynamics reveal the unique physicochemical properties of the SARS-CoV-2 fusion peptides, accrued in part from the presence of consecutive prolines that impart backbone rigidity which aids the virus fusogenicity. The priming of the Spike by its cleavage and subsequent fusogenic conformational transition steered by the fusion loop may be critical for the SARS-CoV-2 spread. The importance of the fusion loop makes it an apt target for anti-virals and vaccine candidates.

Studying the helical conformations of aspereline peptides

Asperelines are short-sequence peptaibol molecules, and these peptides composed of 10 residues were isolated from the Trichoderma asperellum. In our study, a detailed structural characterization was performed on the asperelines by means of molecular dynamics methods. For the aspereline peptides, the occurrence of various secondary structural elements (i.e. β-turns and helical structures) was investigated along their entire sequences. The results derived from the simulated annealing calculations led to the observations that in the case of asperelines, the types I, III and III' β-turn structures, as well as their stabilizing i ← i+3 H-bonds appeared. However, beside the different β-turns, shorter or longer helical structures were also detected. Based on the results obtained by the molecular dynamics simulations, it was concluded that the three-dimensional structure of aspereline peptides could be characterized by helical conformations (i.e. 3₁₀ - and α-helix). Nevertheless, on the basis of individual molecular dynamics trajectories, it was observed that the asperelines could adopt not only the right-handed, but also the left-handed helical structures.

Reciprocal Coupling in Chemically Fueled Assembly: A Reaction Cycle Regulates Self-Assembly and Vice Versa

In biology, self-assembly of proteins and energy-consuming reaction cycles are intricately coupled. For example, tubulin is activated and deactivated for assembly by a guanosine triphosphate (GTP)-driven reaction cycle, and the emerging microtubules catalyze this reaction cycle by changing the microenvironment of the activated tubulin. Recently, synthetic analogs of chemically fueled assemblies have emerged, but examples in which assembly and reaction cycles are reciprocally coupled remain rare. In this work, we report a peptide that can be activated and deactivated for self-assembly. The emerging assemblies change the microenvironment of their building blocks, which consequently accelerate the rates of building block deactivation and reactivation. We quantitatively understand the mechanisms at play, and we are thus able to tune the catalysis by molecular design of the peptide precursor.

Increasing protein stability by engineering the n → π* interaction at the β-turn

Abundant n → π* interactions between adjacent backbone carbonyl groups, identified by statistical analysis of protein structures, are predicted to play an important role in dictating the structure of proteins. However, experimentally testing the prediction in proteins has been challenging due to the weak nature of this interaction. By amplifying the strength of the n → π* interaction via amino acid substitution and thioamide incorporation at a solvent exposed β-turn within the GB1 proteins and Pin 1 WW domain, we demonstrate that an n → π* interaction increases the structural stability of proteins by restricting the ϕ torsion angle. Our results also suggest that amino acid side-chain identity and its rotameric conformation play an important and decisive role in dictating the strength of an n → π* interaction.

Amino acid residues adopt a right-handed α-helical conformation with increasing strength of the n → π* interaction. We also demonstrate a direct consequence of n → π* interactions on enhancing the structural stability of proteins.

Solvation and the secondary structure of a proline-containing dipeptide: insights from VCD spectroscopy

In this study we investigate the IR and VCD spectra of the diastereomeric dipeptide Boc-Pro-Phe-(n-propyl) 1 in chloroform-d₁ (CDCl₃) and the strongly hydrogen bonding solvent dimethylsulfoxide-d₆ (DMSO-d₆). From comparison of the experimental spectra, the amide II spectral region is identified as marker signature for the stereochemistry of the dipeptide: the homochiral LL-1 features a (+/-)-pattern in the amide II region of the VCD spectrum, while the amide II signature of the diastereomer LD-1 is inverted. Computational analysis of the IR and VCD spectra of LL-1 reveals that the experimentally observed amide II signature is characteristic for a β_I-turn structure of the peptide. Likewise, the inverted pattern found for LD-1 arises from a β_II-turn structure of the dipeptide. Following a micro-solvation approach, the experimental spectra recorded in DMSO-d₆ are computationally well reproduced by considering only a single solvent molecule in a hydrogen bond with N-H groups. Considering a second solvent molecule, which would lead to a cleavage of intramolecular hydrogen bonds in 1, is found to give a significantly worse match with the experiment. Hence, the detailed computational analysis of the spectra of LL- and LD-1 recorded in DMSO-d₆ confirms that the intramolecular hydrogen bonding pattern, that stabilizes the β-turns and other conformations of LL- and LD-1 in apolar solvents, remains intact. Our findings also show that it is essential to consider solvation explicitly in the analysis of the IR and VCD spectra of dipeptides in strongly hydrogen bonding solvents. As the solute-solvent interactions affect both conformational preferences and spectral signatures, it is also demonstrated that this inclusion of solvent molecules cannot be circumvented by applying fitting procedures to non-solvated structures.

New facets of larger Nest motifs in proteins

The Nest is a concave-shaped structural motif in proteins formed by consecutive enantiomeric left-handed (L) and right-handed (R) helical conformation of the backbone. This important motif subsumes many turn and helix capping structures and binds electron-rich ligands. Simple Nests are either RL or LR. Larger Nests (>2 residues long) may be RLR, LRL, RLRL, and so forth, being considered as composed of overlapping simple Nests. The larger Nests remain under-explored despite their widely known contributions to protein function. In our study, we address whether the recurrence of enantiomeric geometry in the larger Nests constrains the peptide backbone such that distinct compositional and conformational preferences are seen compared to simple Nests. Our analysis reveals the critical role of the L helical torsion angle in the formation of larger Nests. This can be observed through the higher propensity of residue or secondary structure combinations in LR and LRL backbone conformation in comparison to RL or RLR, although LR/LRL is considerably lower by occurrence. We also find that the most abundant doublets and triplets in Nests have a propensity for particular secondary structures, suggesting a strong sequence-structure relationship in the larger Nest. Overall, our analysis corroborates distinct features of simple and the larger Nests. Such insights would be helpful towards in-vitro design of peptides and peptidomimetic studies.

Delineation of a new structural motif involving NHN γ-turn

Macromolecules are characterized by distinctive arrangement of hydrogen bonds. Different patterns of hydrogen bonds give rise to distinct and stable structural motifs. An analysis of 4114 non-redundant protein chains reveals the existence of a three-residue, (i - 1) to (i + 1), structural motif, having two hydrogen-bonded five-membered pseudo rings (the first, an NH···OC involving the first residue, and the second being NH∙∙∙N involving the last two residues), separated by a peptide bond. There could be an additional hydrogen bond between the side-chain at (i-1) and the main-chain NH of (i + 1). The average backbone torsion angles of -76(±21)° and - 12(±17)° at i creates a tight turn in the polypeptide chain, akin to a γ-turn. Indeed, a search of three-residue fragments with restriction on the terminal C^α ···C^α distance and the existence of the two pseudo rings on either side revealed the presence 14 846 cases of a variant, termed NHN γ-turn, distinct from the NHO γ-turn (2032 cases) that has traditionally been characterized by the presence of NHO hydrogen bond linking the terminal main-chain atoms. As in the latter, the newly identified γ-turns are also of two types-classical and inverse, occurring in the ratio of 1:6. The propensities of residues to occur in these turns and their secondary structural features have been enumerated. An understanding of these turns would be useful for structure prediction and loop modeling, and may serve as models to represent some of the unfolded state or disordered region in proteins.

Crystal structure of N-{N-[N-(tert-but-oxy-carbon-yl)-l-α-aspart-yl]-l-α-aspart-yl}-l-α-aspartic acid 14,24,34-trimethyl ester 31-2-oxo-2-phenyl-ethyl ester {Boc-[Asp(OMe)]3-OPac}

In the title homotripeptide {Boc-[Asp(OMe)]3-OPac}, C28H37N3O13, all peptide bonds adopt an s-trans conformation with respect to the N-H and C=O groups. In the crystal, N-H⋯O hydrogen bonds result in an infinite parallel β-sheet structure running along the b-axis direction. The Boc protecting group at the N-terminus of the peptide is disordered over two sites with occupancy factors of 0.504 (5) and 0.496 (5).

A proline insertion-deletion in the spike glycoprotein fusion peptide of mouse hepatitis virus strongly alters neuropathology

Fusion peptides (FPs) in spike proteins are key players mediating early events in cell-to-cell fusion, vital for intercellular viral spread. A proline residue located at the central FP region has often been suggested to have a distinctive role in this fusion event. The spike glycoprotein from strain RSA59 (PP) of mouse hepatitis virus (MHV) contains two central, consecutive prolines in the FP. Here, we report that deletion of one of these proline residues, resulting in RSA59 (P), significantly affected neural cell syncytia formation and viral titers postinfection in vitro Transcranial inoculation of C57Bl/6 mice with RSA59 (PP) or RSA59 (P) yielded similar degrees of necrotizing hepatitis and meningitis, but only RSA59 (PP) produced widespread encephalitis that extended deeply into the brain parenchyma. By day 6 postinfection, both virus variants were mostly cleared from the brain. Interestingly, inoculation with the RSA59 (P)-carrying MHV significantly reduced demyelination at the chronic stage. We also found that the presence of two consecutive prolines in FP promotes a more ordered, compact, and rigid structure in the spike protein. These effects on FP structure were due to proline's unique stereochemical properties intrinsic to its secondary amino acid structure, revealed by molecular dynamics and NMR experiments. We therefore propose that the differences in the severity of encephalitis and demyelination between RSA59 (PP) and RSA59 (P) arise from the presence or absence, respectively, of the two consecutive prolines in FP. Our studies define a structural determinant of MHV entry in the brain parenchyma important for altered neuropathogenesis.

Ramachandran maps for side chains in globular proteins

The Ramachandran plot for backbone ϕ,ψ-angles in a blocked monopeptide has played a central role in understanding protein structure. Curiously, a similar analysis for side chain χ-angles has been comparatively neglected. Instead, efforts have focused on compiling various types of side chain libraries extracted from proteins of known structure. Departing from this trend, the following analysis presents backbone-based maps of side chains in blocked monopeptides. As in the original ϕ,ψ-plot, these maps are derived solely from hard-sphere steric repulsion. Remarkably, the side chain biases exhibit marked similarities to corresponding biases seen in high-resolution protein structures. Consequently, some of the entropic cost for side chain localization in proteins is prepaid prior to the onset of folding events because conformational bias is built into the chain at the covalent level. Furthermore, side chain conformations are seen to experience fewer steric restrictions for backbone conformations in either the α or β basins, those map regions where repetitive ϕ,ψ-angles result in α-helices or strands of β-sheet, respectively. Here, these α and β basins are entropically favored for steric reasons alone; a blocked monopeptide is too short to accommodate the peptide hydrogen bonds that stabilize repetitive secondary structure. Thus, despite differing energetics, α/β-basins are favored for both monopeptides and repetitive secondary structure, underpinning an energetically unfrustrated compatibility between these two levels of protein structure.

Structural determinants of coiled coil mechanics

The natural abundance of coiled coil (CC) motifs in the cytoskeleton and the extracellular matrix suggests that CCs play a crucial role in the bidirectional mechanobiochemical signaling between cells and the matrix. Their functional importance and structural simplicity has allowed the development of numerous applications, such as protein-origami structures, drug delivery systems and biomaterials. With the goal of establishing CCs as nanomechanical building blocks, we investigated the importance of helix propensity and hydrophobic core packing on the mechanical stability of 4-heptad CC heterodimers. Using single-molecule force spectroscopy, we show that both parameters determine the force-induced dissociation in shear loading geometry; however, with different effects on the energy landscape. Decreasing the helix propensity lowers the transition barrier height, leading to a concomitant decrease in the distance to the transition state. In contrast, a less tightly packed hydrophobic core increases the distance to the transition state. We propose that this originates from a larger side chain dynamics, possible water intrusion at the interface as well as differences in solvation of the hydrophobic amino acids at the transition state. In conclusion, the different contributions of helix propensity and hydrophobic core packing need to be considered when tuning the mechanical properties of CCs for applications.

Can Mutational Analysis Be Used To Assist Structure Determination of Peptides?

Mutational analysis is widely used to study the relationship between sequence and structure of proteins and peptides. It is often assumed that substituting a proline with another amino acid "locks" the peptide bond in the trans conformation, allowing only a subset of the initial molecular geometries to be observed. To test this assumption, we assess the result of substituting two prolines in the bradykinin sequence with alanine using field-asymmetric ion mobility spectrometry combined with cryogenic ion spectroscopy in the gas phase. While the structure of the mutant coincides with a part of the conformational space of the original peptide, the higher flexibility of the alanine backbone compared to proline allows it to access additional structures. We conclude that proline-to-nonproline substitutions are helpful to assign structures, but they should be used in conjunction with spectroscopic techniques that allow detailed comparison of the structures of the mutant and the native peptide.

Prediction of Metal Ion Binding Sites in Proteins from Amino Acid Sequences by Using Simplified Amino Acid Alphabets and Random Forest Model

Metal binding proteins or metallo-proteins are important for the stability of the protein and also serve as co-factors in various functions like controlling metabolism, regulating signal transport, and metal homeostasis. In structural genomics, prediction of metal binding proteins help in the selection of suitable growth medium for overexpression's studies and also help in obtaining the functional protein. Computational prediction using machine learning approach has been widely used in various fields of bioinformatics based on the fact all the information contains in amino acid sequence. In this study, random forest machine learning prediction systems were deployed with simplified amino acid for prediction of individual major metal ion binding sites like copper, calcium, cobalt, iron, magnesium, manganese, nickel, and zinc.

Intrinsic Ribosome Destabilization Underlies Translation and Provides an Organism with a Strategy of Environmental Sensing

Nascent polypeptides can modulate the polypeptide elongation speed on the ribosome. Here, we show that nascent chains can even destabilize the translating Escherichia coli ribosome from within. This phenomenon, termed intrinsic ribosome destabilization (IRD), occurs in response to a special amino acid sequence of the nascent chain, without involving the release or the recycling factors. Typically, a consecutive array of acidic residues and those intermitted by alternating prolines induce IRD. The ribosomal protein bL31, which bridges the two subunits, counteracts IRD, such that only strong destabilizing sequences abort translation in living cells. We found that MgtL, the leader peptide of a Mg²⁺ transporter (MgtA), contains a translation-aborting sequence, which sensitizes the ribosome to a decline in Mg²⁺ concentration and thereby triggers the MgtA-upregulating genetic scheme. Translation proceeds at an inherent risk of ribosomal destabilization, and nascent chain-ribosome complexes can function as a Mg²⁺ sensor by harnessing IRD.

Correlated Mutation in the Evolution of Catalysis in Uracil DNA Glycosylase Superfamily

Enzymes in Uracil DNA glycosylase (UDG) superfamily are essential for the removal of uracil. Family 4 UDGa is a robust uracil DNA glycosylase that only acts on double-stranded and single-stranded uracil-containing DNA. Based on mutational, kinetic and modeling analyses, a catalytic mechanism involving leaving group stabilization by H155 in motif 2 and water coordination by N89 in motif 3 is proposed. Mutual Information analysis identifies a complexed correlated mutation network including a strong correlation in the EG doublet in motif 1 of family 4 UDGa and in the QD doublet in motif 1 of family 1 UNG. Conversion of EG doublet in family 4 Thermus thermophilus UDGa to QD doublet increases the catalytic efficiency by over one hundred-fold and seventeen-fold over the E41Q and G42D single mutation, respectively, rectifying the strong correlation in the doublet. Molecular dynamics simulations suggest that the correlated mutations in the doublet in motif 1 position the catalytic H155 in motif 2 to stabilize the leaving uracilate anion. The integrated approach has important implications in studying enzyme evolution and protein structure and function.

Making a Right or Left Choice: Chiral Self-Sorting as a Tool for the Formation of Discrete Complex Structures

This review discusses chiral self-sorting-the process of choosing an interaction partner with a given chirality from a complex mixture of many possible racemic partners. Chiral self-sorting (also known as chiral self-recognition or chiral self-discrimination) is fundamental for creating functional structures in nature and in the world of chemistry because interactions between molecules of the same or the opposite chirality are characterized by different interaction energies and intrinsically different resulting structures. However, due to the similarity between recognition sites of enantiomers and common conformational lability, high fidelity homochiral or heterochiral self-sorting poses a substantial challenge. Chiral self-sorting occurs among natural and synthetic molecules that leads to the amplification of discrete species. The review covers a variety of complex self-assembled structures ranging from aggregates made of natural and racemic peptides and DNA, through artificial functional receptors, macrocyles, and cages to catalytically active metal complexes and helix mimics. The examples involve a plethora of reversible interactions: electrostatic interactions, π-π stacking, hydrogen bonds, coordination bonds, and dynamic covalent bonds. A generalized view of the examples collected from different fields allows us to suggest suitable geometric models that enable a rationalization of the observed experimental preferences and establishment of the rules that can facilitate further design.

Symmetry-Directed Self-Organization in Peptide Nanoassemblies through Aromatic π-π Interactions

Almost all biological systems are assemblies of one or more biomolecules from nano- to macrodimensions. Unlike inorganic molecules, peptide systems attune with the conceptual framework of aggregation models when forming nanoassemblies. Three significant recent theoretical models have indicated that nucleation, end-to-end association, and geometry of growth are determined primarily by the size and electrostatics of the individual basic building blocks. In this study, we tested six model systems, differentially modulating the prominence of three design variables, namely, aromatic π-π interactions, local electrostatics, and overall symmetry of the basic building unit. Our results indicate that the crucial design elements in a peptide-based nanoassembly are (a) a stable extended π-π interaction network, (b) size, and (c) overall symmetry of the basic building blocks. The six model systems represent all of the design variables in the best manner possible, considering the complexity of a biomolecule. The results provide important directives in deciding the morphology and crystallinity of peptide nanoassemblies.

Mapping the Geometric Evolution of Protein Folding Motor

Polypeptide chain has an invariant main-chain and a variant side-chain sequence. How the side-chain sequence determines fold in terms of its chemical constitution has been scrutinized extensively and verified periodically. However, a focussed investigation on the directive effect of side-chain geometry may provide important insights supplementing existing algorithms in mapping the geometrical evolution of protein chains and its structural preferences. Geometrically, folding of protein structure may be envisaged as the evolution of its geometric variables: ϕ, and ψ dihedral angles of polypeptide main-chain directed by χ1, and χ2 of side chain. In this work, protein molecule is metaphorically modelled as a machine with 4 rotors ϕ, ψ, χ1 and χ2, with its evolution to the functional fold is directed by combinations of its rotor directions. We observe that differential rotor motions lead to different secondary structure formations and the combinatorial pattern is unique and consistent for particular secondary structure type. Further, we found that combination of rotor geometries of each amino acid is unique which partly explains how different amino acid sequence combinations have unique structural evolution and functional adaptation. Quantification of these amino acid rotor preferences, resulted in the generation of 3 substitution matrices, which later on plugged in the BLAST tool, for evaluating their efficiency in aligning sequences. We have employed BLOSUM62 and PAM30 as standard for primary evaluation. Generation of substitution matrices is a logical extension of the conceptual framework we attempted to build during the development of this work. Optimization of matrices following the conventional routines and possible application with biologically relevant data sets are beyond the scope of this manuscript, though it is a part of the larger project design.

A novel secondary structure based on fused five-membered rings motif

An analysis of protein structures indicates the existence of a novel, fused five-membered rings motif, comprising of two residues (i and i + 1), stabilized by interresidue N_i+1–H∙∙∙N_i and intraresidue N_i+1–H∙∙∙O=C_i+1 hydrogen bonds. Fused-rings geometry is the common thread running through many commonly occurring motifs, such as β-turn, β-bulge, Asx-turn, Ser/Thr-turn, Schellman motif, and points to its structural robustness. A location close to the beginning of a β-strand is rather common for the motif. Devoid of side chain, Gly seems to be a key player in this motif, occurring at i, for which the backbone torsion angles cluster at ~(−90°, −10°) and (70°, 20°). The fused-rings structures, distant from each other in sequence, can hydrogen bond with each other, and the two segments aligned to each other in a parallel fashion, give rise to a novel secondary structure, topi, which is quite common in proteins, distinct from two major secondary structures, α-helix and β-sheet. Majority of the peptide segments making topi are identified as aggregation-prone and the residues tend to be conserved among homologous proteins.

Synthesis of a new β-amino acid with a 3-deoxy-L-ara furnaoside side chain: the influence of the side chain on the conformation of α/β-peptides

The important role of side chains in the stabilization of helical folds in peptidic foldamers containing C-linked carbo-β-amino acids (β-Caa), an interesting class of β-amino acids, with carbohydrate side chains has been extensively elaborated. As a pragmatic approach to alleviate the interference of substituents in the side chains on the folding propensities of the peptides, they are often modified or removed. The present study reports the synthesis of a new β-Caa with a 3-deoxy-L-ara furanoside side chain, [(R)-β-Caa(da)], from D-glucose, and its use in the synthesis of α/β-peptides in 1 : 1 alternation with D-Ala. The synthesis of peptides using (R)-β-Caa(da), was facile unlike those from (R)-β-Caa(a) having the L-ara furanoside side chain. The detailed NMR, molecular dynamics (MD) and CD studies on the new α/β-peptides showed the presence of robust left-handed 11/9-mixed helices. The study demonstrates that the new (R)-β-Caa(da), behaves differently compared to the other two related monomers, (R)-β-Caa(x) with the D-xylo furanoside side chain and (R)-β-Caa(a).