Although S-nitrosylation of cysteines is a common protein posttranslational modification, little is known about its three-dimensional structural features. This paper describes a systematic survey of the data available in the Protein Data Bank. Several interesting observations could be made. (1) As a result of radiation damage, S-nitrosylated cysteines (Snc) are frequently reduced, at least partially. (2) S-nitrosylation may be a protection against irreversible thiol oxidation; because the NO group of Snc is relatively accessible to the solvent, it may act as a cork to protect the sulfur atoms of cysteines from oxidation by molecular oxygen to sulfenic, sulfinic, and sulfonic acid; moreover, Snc are frequently found at the start or end of helices and strands and this might shield secondary structural elements from unfolding.

Solvent paramagnetic relaxation enhancement (sPRE) is a versatile nuclear magnetic resonance (NMR)-based method that allows characterization of the structure and dynamics of biomolecular systems through providing quantitative experimental information on solvent accessibility of NMR-active nuclei. Addition of soluble paramagnetic probes to the solution of a biomolecule leads to paramagnetic relaxation enhancement in a concentration-dependent manner. Here we review recent progress in the sPRE-based characterization of structural and dynamic properties of biomolecules and their complexes, and aim to deliver a comprehensive illustration of a growing number of applications of the method to various biological systems. We discuss the physical principles of sPRE measurements and provide an overview of available co-solute paramagnetic probes. We then explore how sPRE, in combination with complementary biophysical techniques, can further advance biomolecular structure determination, identification of interaction surfaces within protein complexes, and probing of conformational changes and low-population transient states, as well as deliver insights into weak, nonspecific, and transient interactions between proteins and co-solutes. In addition, we present examples of how the incorporation of solvent paramagnetic probes can improve the sensitivity of NMR experiments and discuss the prospects of applying sPRE to NMR metabolomics, drug discovery, and the study of intrinsically disordered proteins.

OBJECTIVE: Hydration forces between surfactant bilayers can be assessed using water sorption isotherms of surfactants. For a quantitative description, a water sorption model that relates water activity to water content in surfactant-based systems should be proposed.
UNASSIGNED: A water sorption model for nonionic surfactant systems based on the idea on partial solvent accessibility is proposed. The model contains only two parameters: one describes the strength of interactions, the other describes the fraction of surfactant available for water. For comparison, molecular dynamics simulations of bilayers of n-octyl β-d-glucoside with different water contents are presented.
RESULTS: The model provides an excellent fit of experimental data on water sorption isotherms of two sugar surfactants. The results of the fitting are compared with molecular dynamics simulations and show a good correlation between simulations and the theory proposed. Analysis of interaction energies shows weakly endothermic hydration both in the simulations and in the sorption model, which agrees with calorimetric data on hydration. The model also shows a non-exponential decay of hydration forces with respect to the distance between bilayers; an expression for the decay length is derived.

Development of donor-specific human leukocyte antigen (HLA) antibodies (DSA) remains a major risk factor for graft loss following organ transplantation, where DSA are directed towards patches on the three-dimensional structure of the respective organ donor\'s HLA proteins. Matching donors and recipients based on HLA epitopes appears beneficial for the avoidance of DSA. Defining surface epitopes however remains challenging and the concepts underlying their characterization are not fully understood. Based on our recently implemented computational deep learning pipeline to define HLA Class I protein-specific surface residues, we hypothesized a correlation between the number of HLA protein-specific solvent-accessible interlocus amino acid mismatches (arbitrarily called Snowflake) and the incidence of DSA. To validate our hypothesis, we considered two cohorts simultaneously. The kidney transplant cohort (KTC) considers 305 kidney-transplanted patients without DSA prior to transplantation. During the follow-up, HLA antibody screening was performed regularly to identify DSA. The pregnancy cohort (PC) considers 231 women without major sensitization events prior to pregnancy who gave live birth. Post-delivery serum was screened for HLA antibodies directed against the child\'s inherited paternal haplotype (CSA). Based on the involved individuals\' HLA typings, the numbers of interlocus-mismatched antibody-verified eplets (AbvEPS), the T cell epitope PIRCHE-II model and Snowflake were calculated locus-specific (HLA-A, -B and -C), normalized and pooled. In both cohorts, Snowflake numbers were significantly elevated in recipients/mothers that developed DSA/CSA. Univariable regression revealed significant positive correlation between DSA/CSA and AbvEPS, PIRCHE-II and Snowflake. Snowflake numbers showed stronger correlation with numbers of AbvEPS compared to Snowflake numbers with PIRCHE-II. Our data shows correlation between Snowflake scores and the incidence of DSA after allo-immunization. Given both AbvEPS and Snowflake are B cell epitope models, their stronger correlation compared to PIRCHE-II and Snowflake appears plausible. Our data confirms that exploring solvent accessibility is a valuable approach for refining B cell epitope definitions.

Given the challenges for the experimental determination of RNA tertiary structures, probing solvent accessibility has become increasingly important to gain functional insights. Among various chemical probes developed, backbone-cleaving hydroxyl radical is the only one that can provide unbiased detection of all accessible nucleotides. However, the readouts have been based on reverse transcription (RT) stop at the cleaving sites, which are prone to false positives due to PCR amplification bias, early drop-off of reverse transcriptase, and the use of random primers in RT reaction. Here, we introduced a fixed-primer method called RL-Seq by performing RtcB Ligation (RL) between a fixed 5\'-OH-end linker and unique 3\'-P-end fragments from hydroxyl radical cleavage prior to high-throughput sequencing. The application of this method to E. coli ribosomes confirmed its ability to accurately probe solvent accessibility with high sensitivity (low required sequencing depth) and accuracy (strong correlation to structure-derived values) at the single-nucleotide resolution. Moreover, a near-perfect correlation was found between the experiments with and without using unique molecular identifiers, indicating negligible PCR biases in RL-Seq. Further improvement of RL-Seq and its potential transcriptome-wide applications are discussed.

Molecular Dynamics (MD) is a method used to calculate the movement of atoms and molecules broadly applied to several aspects of science. It involves computational simulation, which makes it, at first glance, not easily accessible. The rise of several automated tools to perform molecular simulations has allowed researchers to navigate through the various steps of MD. This enables to elucidate structural properties of proteins that could not be analyzed otherwise, such as the impact of glycosylation. Glycosylation dictates the physicochemical and biological properties of a protein modulating its solubility, stability, resistance to proteolysis, interaction partners, enzymatic activity, binding and recognition. Given the high conformational and compositional diversity of the glycan chains, assessing their influence on the protein structure is challenging using conventional analytical techniques. In this manuscript, we present a step-by-step workflow to build and perform MD analysis of glycoproteins focusing on the SPIKE glycoprotein of SARS-CoV-2 to appraise the impact of glycans in structure stabilization and antibody occlusion.

Sequence-based predictors of the residue-level protein function and structure cover a broad spectrum of characteristics including intrinsic disorder, secondary structure, solvent accessibility and binding to nucleic acids. They were catalogued and evaluated in numerous surveys and assessments. However, methods focusing on a given characteristic are studied separately from predictors of other characteristics, while they are typically used on the same proteins. We fill this void by studying complementarity of a representative collection of methods that target different predictions using a large, taxonomically consistent, and low similarity dataset of human proteins. First, we bridge the gap between the communities that develop structure-trained vs. disorder-trained predictors of binding residues. Motivated by a recent study of the protein-binding residue predictions, we empirically find that combining the structure-trained and disorder-trained predictors of the DNA-binding and RNA-binding residues leads to substantial improvements in predictive quality. Second, we investigate whether diverse predictors generate results that accurately reproduce relations between secondary structure, solvent accessibility, interaction sites, and intrinsic disorder that are present in the experimental data. Our empirical analysis concludes that predictions accurately reflect all combinations of these relations. Altogether, this study provides unique insights that support combining results produced by diverse residue-level predictors of protein function and structure.

Characterizing RNA structures and functions have mostly been focused on 2D, secondary and 3D, tertiary structures. Recent advances in experimental and computational techniques for probing or predicting RNA solvent accessibility make this 1D representation of tertiary structures an increasingly attractive feature to explore. Here, we provide a survey of these recent developments, which indicate the emergence of solvent accessibility as a simple 1D property, adding to secondary and tertiary structures for investigating complex structure-function relations of RNAs.

MTHFR deficiency still deserves an investigation to associate the phenotype to protein structure variations. To this aim, considering the MTHFR wild type protein structure, with a catalytic and a regulatory domain and taking advantage of state-of-the-art computational tools, we explore the properties of 72 missense variations known to be disease associated. By computing the thermodynamic ΔΔG change according to a consensus method that we recently introduced, we find that 61% of the disease-related variations destabilize the protein, are present both in the catalytic and regulatory domain and correspond to known biochemical deficiencies. The propensity of solvent accessible residues to be involved in protein-protein interaction sites indicates that most of the interacting residues are located in the regulatory domain, and that only three of them, located at the interface of the functional protein homodimer, are both disease-related and destabilizing. Finally, we compute the protein architecture with Hidden Markov Models, one from Pfam for the catalytic domain and the second computed in house for the regulatory domain. We show that patterns of disease-associated, physicochemical variation types, both in the catalytic and regulatory domains, are unique for the MTHFR deficiency when mapped into the protein architecture.

Spike (S) glycoprotein of SARS-CoV2 exists chiefly in two conformations, open and closed. Most previous structural studies on S protein have been conducted at pH 8.0, but knowledge of the conformational propensities under both physiological and endosomal pH conditions is important to inform vaccine development. Our current study employed single-particle cryoelectron microscopy to visualize multiple states of open and closed conformations of S protein at physiological pH 7.4 and near-physiological pH 6.5 and pH 8.0. Propensities of open and closed conformations were found to differ with pH changes, whereby around 68% of S protein exists in open conformation at pH 7.4. Furthermore, we noticed a continuous movement in the N-terminal domain, receptor-binding domain (RBD), S2 domain, and stalk domain of S protein conformations at various pH values. Several key residues involving RBD-neutralizing epitopes are differentially exposed in each conformation. This study will assist in developing novel therapeutic measures against SARS-CoV2.