A proteome-scale analysis of vertebrate protein amino acid occurrence: Thermoadaptation shows a correlation with protein solvation but less so with dynamics.

Despite differences in behaviors and living conditions, vertebrate organisms share the great majority of proteins, often with subtle differences in amino acid sequence. Here, we present a simple way to analyze the difference in amino acid occurrence by comparing highly homologous proteins on a subproteome level between several vertebrate model organisms. Specifically, we use this method to identify a pattern of amino acid conservation as well as a shift in amino acid occurrence between homeotherms (warm-blooded species) and poikilotherms (cold-blooded species). Importantly, this general analysis and a specific example further establish a broad correlation, if not likely connection between the thermal adaptation of protein sequences and two of their physical features: on average a change in their protein dynamics and, even more strongly, in their solvation. For poikilotherms, such as frog and fish, the lower body temperature is expected to increase the protein-protein interaction due to a decrease in protein internal dynamics. In order to counteract the tendency for enhanced binding caused by low temperatures, poikilotherms enhance the solvation of their proteins by favoring polar amino acids. This feature appears to dominate over possible changes in dynamics for some proteins. The results suggest that a general trend for amino acid choice is part of the mechanism for thermoadaptation of vertebrate organisms at the molecular level.

Computational analysis of interior mutations of SARS-CoV-2 Spike protein suggest a balance of protein stability and S2: S1 separation propensity.

SARS-CoV-2 variants often include surface mutations in the Spike protein that are important for viruses to recognize host receptors and evade antibody neutralization. The Spike protein also has mutations in the interior of the protein likely to affect the Spike protein S1 - S2 subunit's separation propensity, the most important of which is the D614G mutation. Remarkably, the Omicron variant contains a large number of internal mutations at the S2: S1 interface, which have not been investigated yet. In this study, we examined the effects of such interfacial mutations on the S2: S1 and subunit domain interactions and on the subunit's dissociation process. We found that the interaction with S2 is mainly contributed by the three encapsulation domains, named INT, ED1 and ED2 of S1, which are sandwiched between the S1 RBD and N-terminal NTD domain. We found that D614 is the strongest contributor for the S2: S1 interaction which is greatly weakened by the D614G mutation. Surprisingly, we found that, mutations T547K, H655Y, N764K, N856K, N969K, L981F in the Omicron variant largely enhance the S2: ED1 interaction, partially compensating the loss of S2: ED2 interaction due to the D614G mutation. Lastly, these results, together with biological considerations, allow us to suggest that in addition to the binding strength of between the RBD and ACE2, the stability of the Spike protein and the propensity of Spike protein S2: S1 separation are critical factors which likely exist in a balance for a particular infectivity and pathogenicity of the virus.

Computational studies of the principle of dynamic-change-driven protein interactions.

Dynamic allostery emphasizes a role of entropy change manifested as a sole change in protein fluctuations without structural changes. This kind of entropy-driven effect remains largely understudied. The most significant examples involve protein-ligand interactions, leaving protein-protein interactions, which are critical in signaling and other cellular events, largely unexplored. Here we study an example of how protein-protein interaction (binding of Ras to the Ras binding domain [RBD] of the effector protein Raf) affects a subsequent protein association process (Ras dimerization) by quenching Ras internal motions through dynamic allostery. We also investigate the influence of point mutations or ambient temperature, respectively, on the protein dynamics and interaction of two other systems: in adenylate kinase (ADK) and in the EphA2 SAM:Ship2 SAM complex. Based on these examples, we postulate that there are different ways in which dynamic-change-driven protein interactions are manifested and that it is likely a general biological phenomenon.

Neuropilin-1 assists SARS-CoV-2 infection by stimulating the separation of Spike protein S1 and S2.

The cell surface receptor Neuropilin-1 (Nrp1) was recently identified as a host factor for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) entry. The Spike protein of SARS-CoV-2 is cleaved into two segments, the S1 (residues (res.) 1-685) and the S2 (res. 686-1273) domains by furin protease. Nrp1 predominantly binds to the C-terminal RRAR amino acid motif (res. 682-685) of the S1 domain. In this study, we firstly modeled the association of an Nrp1 protein (consisting of domains a2-b1-b2) with the Spike protein. Next, we studied the separation of S2 from the S1 domain, with and without Nrp1 bound, by utilizing molecular dynamics pulling simulations. During the separation, Nrp1 stabilizes the S1 C-terminal region (res. 640-685) and thereby assists the detachment of S2 N-terminal region (res. 686-700). Without Nrp1 bound, S1 tends to become stretched, whereas the bound Nrp1 stimulates an earlier separation of S2 from the S1 domain. The liberated S2 domain is known to mediate the fusion of virus and host membranes; thus, Nrp1 likely increases virus infectivity by facilitating the S1 and S2 separation. We further analyzed the possible topological structure of the SARS-CoV-2 Spike protein when bound with Nrp1 and angiotensin-converting enzyme 2 (ACE2). Understanding of such an Nrp1-assisted viral infection opens the gate for the generation of protein-protein inhibitors, such as antibodies, which could attenuate the infection mechanism and protect certain cells in a future Nrp1-ACE2 targeted combination therapy.

Beyond history and "on a roll": The list of the most well-studied human protein structures and overall trends in the protein data bank.

Of the roughly 20,000 canonical human protein sequences, as of January 20, 2021, 7,077 proteins have had their full or partial, medium- to high-resolution structures determined by x-ray crystallography or other methods. Which of these proteins dominate the protein data bank (the PDB) and why? In this paper, we list the 273 top human protein structures based on the number of their PDB entries. This set of proteins accounts for more than 40% of all available human PDB entries and represent past trends as well as current status for protein structural biology. We briefly discuss the relationship which some of the prominent protein structures have with protein research as a whole and mention their relevance to human diseases. The top-10 soluble and membrane proteins are all well-known (most of their first structures being deposited more than 30 years ago). Overall, there is no dramatic change in recent trends in the PDB. Remarkably, the number of structure depositions has grown nearly exponentially over the last 10 or more years (with a doubling time of 7 years for proteins, obtained from any organism). Growth in human protein structures is slightly faster (at 5.9 years). The information in this paper may be informative to senior scientists but also inspire researchers who are new to protein science, providing the year 2021 snap-shot for the state of protein structural biology.

Neuropilin-1 Assists SARS-CoV-2 Infection by Stimulating the Separation of Spike Protein Domains S1 and S2.

The cell surface receptor Neuropilin-1 (Nrp1) was recently identified as a host factor for SARS-CoV-2 entry. As the Spike protein of SARS-CoV-2 is cleaved into the S1 and the S2 domain by furin protease, Nrp1 binds to the newly created C-terminal RRAR amino acid sequence of the S1 domain. In this study, we model the association of a Nrp1 (a2-b1-b2) protein with the Spike protein computationally and analyze the topological constraints in the SARS-CoV-2 Spike protein for binding with Nrp1 and ACE2. Importantly, we study the exit mechanism of S2 from the S1 domain with the assistance of ACE2 as well as Nrp1 by molecular dynamics pulling simulations. In the presence of Nrp1, by binding the S1 more strongly to the host membrane, there is a high probability of S2 being pulled out, rather than S1 being stretched. Thus, Nrp1 binding could stimulate the exit of S2 from the S1 domain, which will likely increase virus infectivity as the liberated S2 domain mediates the fusion of virus and host membranes. Understanding of such a Nrp1-assisted viral infection opens the gate for the generation of protein-protein inhibitors, such as antibodies, which could attenuate the infection mechanism and protect certain cells in a future combination therapy.

Plexin-Bs enhance their GAP activity with a novel activation switch loop generating a cooperative enzyme.

Plexins receive guidance cues from semaphorin ligands and transmit their signal through the plasma membrane. This family of proteins is unique amongst single-pass transmembrane receptors as their intracellular regions interact directly with several small GTPases, which regulate cytoskeletal dynamics and cell adhesion. Here, we characterize the GTPase Activating Protein (GAP) function of Plexin-B1 and find that a cooperative GAP activity towards the substrate GTPase, Rap1b, is associated with the N-terminal Juxtamembrane region of Plexin-B1. Importantly, we unveil an activation mechanism of Plexin-B1 by identifying a novel functional loop which partially blocks Rap1b entry into the plexin GAP domain. Consistent with the concept of allokairy developed for other systems, Plexin-B activity is increased by an apparent substrate-mediated cooperative effect. Simulations and mutagenesis suggest the repositioned JM conformation is stabilized by the new activation switch loop when the active site is occupied, giving rise to faster enzymatic turnover and cooperative behavior. The biological implications, essentially those of a threshold behavior for cell migration, are discussed.

Oral microbiome and pancreatic cancer.

BACKGROUND: Microbiota profiles differ between patients with pancreatic cancer and healthy people, and understanding these differences may help in early detection of pancreatic cancer. Saliva sampling is an easy and cost-effective way to determine microbiota profiles compared to fecal and tissue sample collection. AIM: To investigate the saliva microbiome distribution in patients with pancreatic adenocarcinoma (PDAC) and the role of oral microbiota profiles in detection and risk prediction of pancreatic cancer. METHODS: We conducted a prospective study of patients with pancreatic cancer (n = 41) and healthy individuals (n = 69). Bacterial taxa were identified by 16S ribosomal ribonucleic acid gene sequencing, and a linear discriminant analysis effect size algorithm was used to identify differences in taxa. Operational taxonomic unit values of all selected taxa were converted into a normalized Z-score, and logistic regressions were used to calculate risk prediction of pancreatic cancer. RESULTS: Compared with the healthy control group, carriage of Streptococcus and Leptotrichina (z-score) was associated with a higher risk of PDAC [odds ratio (OR) = 5.344, 95% confidence interval (CI): 1.282-22.282, P = 0.021 and OR = 6.886, 95%CI: 1.423-33.337, P = 0.016, respectively]. Veillonella and Neisseria (z-score) were considered a protective microbe that decreased the risk of PDAC (OR = 0.187, 95%CI: 0.055-0.631, P = 0.007 and OR = 0.309, 95%CI: 0.100-0.952, P = 0.041, respectively). Among the patients with PDAC, patients reporting bloating have a higher abundance of Porphyromonas (P = 0.039), Fusobacterium (P = 0.024), and Alloprevotella (P = 0.041); while patients reporting jaundice had a higher amount of Prevotella (P = 0.008); patients reporting dark brown urine had a higher amount of Veillonella (P = 0.035). Patients reporting diarrhea had a lower amount of Neisseria and Campylobacter (P = 0.024 and P = 0.034), and patients reporting vomiting had decreased Alloprevotella (P = 0.036). CONCLUSION: Saliva microbiome was able to distinguish patients with pancreatic cancer and healthy individuals. Leptotrichia may be specific for patients living in Sichuan Province, southwest China. Symptomatic patients had different bacteria profiles than asymptomatic patients. Combined symptom and microbiome evaluation may help in the early detection of pancreatic cancer.

Modified Potential Functions Result in Enhanced Predictions of a Protein Complex by All-Atom Molecular Dynamics Simulations, Confirming a Stepwise Association Process for Native Protein-Protein Interactions.

The relative prevalence of native protein-protein interactions (PPIs) are the cornerstone for understanding the structure, dynamics and mechanisms of function of protein complexes. In this study, we develop a scheme for scaling the protein-water interaction in the CHARMM36 force field, in order to better fit the solvation free energy of amino acids side-chain analogues. We find that the molecular dynamics simulation with the scaled force field, CHARMM36s, as well as a recently released version, CHARMM36m, effectively improve on the overly sticky association of proteins, such as ubiquitin. We investigate the formation of a heterodimer protein complex between the SAM domains of the EphA2 receptor and the SHIP2 enzyme by performing a combined total of 48 µs simulations with the different potential functions. While the native SAM heterodimer is only predicted at a low rate of 6.7% with the original CHARMM36 force field, the yield is increased to 16.7% with CHARMM36s, and to 18.3% with CHARMM36m. By analyzing the 25 native SAM complexes formed in the simulations, we find that their formation involves a preorientation guided by Coulomb interactions, consistent with an electrostatic steering mechanism. In 12 cases, the complex could directly transform to the native protein interaction surfaces with only small adjustments in domain orientation. In the other 13 cases, orientational and/or translational adjustments are needed to reach the native complex. Although the tendency for non-native complexes to dissociate has nearly doubled with the modified potential functions, a dissociation followed by a reassociation to the correct complex structure is still rare. Instead, the remaining non-native complexes undergo configurational changes/surface searching, which, however, rarely leads to native structures on a time scale of 250 ns. These observations provide a rich picture of the mechanisms of protein-protein complex formation and suggest that computational predictions of native complex PPIs could be improved further.

A "Tug of War" Maintains a Dynamic Protein-Membrane Complex: Molecular Dynamics Simulations of C-Raf RBD-CRD Bound to K-Ras4B at an Anionic Membrane.

Association of Raf kinase with activated Ras triggers downstream signaling cascades toward regulating transcription in the cells' nucleus. Dysregulation of Ras-Raf signaling stimulates cancers. We investigate the C-Raf RBD and CRD regions when bound to oncogenic K-Ras4B at the membrane. All-atom molecular dynamics simulations suggest that the membrane plays an integral role in regulating the configurational ensemble of the complex. Remarkably, the complex samples a few states dynamically, reflecting a competition between C-Raf CRD- and K-Ras4B- membrane interactions. This competition arises because the interaction between the RBD and K-Ras is strong while the linker between the RBD and CRD is short. Such a mechanism maintains a modest binding for the overall complex at the membrane and is expected to facilitate fast signaling processes. Competition of protein-membrane contacts is likely a common mechanism for other multiprotein complexes, if not multidomain proteins at membranes.

The RNA-Binding Site of Poliovirus 3C Protein Doubles as a Phosphoinositide-Binding Domain.

Some viruses use phosphatidylinositol phosphate (PIP) to mark membranes used for genome replication or virion assembly. PIP-binding motifs of cellular proteins do not exist in viral proteins. Molecular-docking simulations revealed a putative site of PIP binding to poliovirus (PV) 3C protein that was validated using nuclear magnetic resonance spectroscopy. The PIP-binding site was located on a highly dynamic α helix, which also functions in RNA binding. Broad PIP-binding activity was observed in solution using a fluorescence polarization assay or in the context of a lipid bilayer using an on-chip, fluorescence assay. All-atom molecular dynamics simulations of the 3C protein-membrane interface revealed PIP clustering and perhaps PIP-dependent conformations. PIP clustering was mediated by interaction with residues that interact with the RNA phosphodiester backbone. We conclude that 3C binding to membranes will be determined by PIP abundance. We suggest that the duality of function observed for 3C may extend to RNA-binding proteins of other viruses.

Computational Modeling Reveals that Signaling Lipids Modulate the Orientation of K-Ras4A at the Membrane Reflecting Protein Topology.

The structural, dynamical, and functional characterization of the small GTPase K-Ras has become a research area of intense focus due to its high occurrence in human cancers. Ras proteins are only fully functional when they interact with the plasma membrane. Here we present all-atom molecular dynamics simulations (totaling 5.8 µs) to investigate the K-Ras4A protein at membranes that contain anionic lipids (phosphatidyl serine or phosphatidylinositol bisphosphate). We find that similarly to the homologous and highly studied K-Ras4B, K-Ras4A prefers a few distinct orientations at the membrane. Remarkably, the protein surface charge and certain lipids can strongly modulate the orientation preference. In a novel analysis, we reveal that the electrostatic interaction (attraction but also repulsion) between the protein's charged residues and anionic lipids determines the K-Ras4A orientation, but that this is also influenced by the topology of the protein, reflecting the geometry of its surfaces.

Effect of vasoactive intestinal peptide on IL-1β expression after corneal alkali burn in rats / 中国免疫学杂志

Objective:To investigate the effect of vasoactive intestinal peptide on IL-1β expression after corneal alkali burn.Methods:27 healthy adult female Wistar rats were randomly divided into group A(normal blank control group),group B(alkali burn treated with saline group),and group C (alkali burn treated with vasoactive intestinal peptide group).Group A were not processed.Group B,C were made right eyes′moderate corneal alkali burn models.After inducing burn,group B and group C were followed respectively by 0.9% saline and 0.25 g/L vasoactive intestinal peptide eye drops,4 times per day.The corneal epithelial healing rates were calculated by measuring the corneal staining area by methylene blue staining at 3 d,7 d and 14 d after inducing burn.At 3 d,7 d and 21 d,the corneal tissues were removed and divided into two equal parts.The corneal tissues were observed by staining,and the expression of IL-1β were detected by immunohistochemistry.Results:The healing rates of corneal epithelium in group C were greater than those in group B at 7 d and 14 d(P<0.05).The mean optical density of IL-1β was measured by immunohistochem-istry.The mean optical density in group B and group C at 3 d were(1.62 ± 0.96) ×10-2,(0.98 ± 0.45) ×10-2,respectively,which at 7 d were (1.42 ± 0.62) ×10-2,(0.71 ± 0.22) ×10-2,respectively,which at 21 d were(0.37 ± 0.24) ×10-2,(0.13 ± 0.16) × 10-2,respectively.The mean optical density of corneal IL-1β in group C were lower than those in group B(P<0.05).Conclusion:Va-soactive intestinal peptide can decrease the expression of IL-1β after corneal alkali burn,down-regulate the level of immune stress, reduce the inflammatory reaction and promote the recovery of cornea after alkali burn.

Interaction of peptides with cell membranes: insights from molecular modeling.

The investigation of the interaction of peptides with cell membranes is the focus of active research. It can enhance the understanding of basic membrane functions such as membrane transport, fusion, and signaling processes, and it may shed light on potential applications of peptides in biomedicine. In this review, we will present current advances in computational studies on the interaction of different types of peptides with the cell membrane. Depending on the properties of the peptide, membrane, and external environment, the peptide-membrane interaction shows a variety of different forms. Here, on the basis of recent computational progress, we will discuss how different peptides could initiate membrane pores, translocate across the membrane, induce membrane endocytosis, produce membrane curvature, form fibrils on the membrane surface, as well as interact with functional membrane proteins. Finally, we will present a conclusion summarizing recent progress and providing some specific insights into future developments in this field.

Protective effects of icariin-mediated SIRT1/FOXO3 signaling pathway on intestinal ischemia/reperfusion-induced acute lung injury.

Acute lung injury (ALI) is a common complication following intestinal ischemia/reperfusion (I/R) and is a major contributing factor to its high mortality rate. Sirtuin 1 (SIRT1), a NAD+-dependent deacetylase, has been reported to have an important role in apoptosis inhibition, oxidative stress resistance and cell lifespan extension through its deacetylation of forkhead box protein O3 (FOXO3). It has been demonstrated that icariin (ICA), a flavonoid extracted from Epimedium, upregulates SIRT1 expression. The aim of the present study was to examine whether ICA-mediated SIRT1/FOXO3 signaling pathway activation had a protective effect on intestinal I/R-induced ALI. The effects of ICA on intestinal I/R-induced ALI and its regulation of the SIRT1/FOXO3 signaling pathway on intestinal I/R-induced ALI were investigated in rats. The results demonstrated that ICA pretreatment markedly reduced intestinal I/R-induced ALI as indicated by histological alterations, including decreased tumor necrosis factor-α (TNF-α), interleukin 6 (IL-6), reduced oxidative stress, acetylated FOXO3 and B-cell lymphoma 2 (Bcl-2)-interacting mediator of cell death levels, and increased glutathione (GSH), GSH peroxidase, SIRT1, manganese superoxide dismutase and Bcl-2 levels in rat lung tissues. Furthermore, ICA pretreatment upregulated SIRT1 expression, which then downregulated FOXO3 acetylation. In conclusion, ICA exhibited significant protective effects in intestinal I/R-induced ALI. The protective effect of ICA may be attributed to the upregulation of SIRT1, which contributed to FOXO3 deacetylation and the modulation of downstream antioxidative and anti-apoptotic factors.

In vivo antioxidant, hypoglycemic, and anti-tumor activities of anthocyanin extracts from purple sweet potato

Anthocyanin from purple sweet potato (PSP) extracted by microwave baking (MB) and acidified electrolyzed water (AEW) exhibited antioxidant activity. After further purification by macroporous AB-8 resin, the color value of PSP anthocyanin (PSPA) reached 30.15 with a total flavonoid concentration of 932.5 mg/g. The purified extracts had more potent antioxidant activities than the crude extracts. After continuously administering the PSP extracts to 12-mo-old mice for 1 mo, the anti-aging index of the experimental group was not significantly different from that of 5-mo-old mice. To a certain degree, PSPA was also effective for controlling plasma glucose levels in male Streptozocin (STZ)-treated diabetic mice. In addition, the extracts inhibited Sarcoma S180 cell growth in ICR mice. Mice consuming the PSP extracts formed significantly fewer and smaller sarcomas than mice consuming the control diets. The highest inhibition rate was 69.03%. These results suggest that anthocyanin extracts from PSP not only exert strong antioxidant effects in vitro, but also had anti-aging, anti-hyperglycemic, and anti-tumor activities.