Intraprotein Interactions Research Articles

Protein solvation: Site-specific hydrophilicity, hydrophobicity, counter ions, and interaction entropy.

Solvation free energy is a driving force that plays an important role in the stability of biomolecular conformations. Currently, the implicit solvent model is widely used to calculate solvation energies of biomolecules such as proteins. However, for proteins, the implicit solvent calculation does not provide much detailed information since a protein is highly inhomogeneous on its surface. In this study, we develop an explicit solvent approach to protein solvation, which allows us to investigate detailed site-specific hydrophilicity and hydrophobicity, including the role of counter ions and intra-protein interactions. This approach facilitates the analysis of specific residue interactions with solvent molecules, extending the understanding of protein solubility to the energetic impacts of site-specific residue-solvent interactions. Our study showed that specific residue-solvent interactions are strongly influenced by the electrostatic environment created by its nearby residues, especially charged residues. In particular, charged residues on the protein surface are mainly responsible for the heterogeneity of the electrostatic environment of the protein surface, and they significantly affect the local distribution of water. In addition, counter ions change the local electrostatic environment and alter specific residue-water interactions. Neutral residues also interact with water, with polar residues being more prominent than nonpolar ones but contributing less to solvation energy than charged residues. This study illustrates an explicit solvent approach to protein solvation, which gives residue-specific contributions to protein solvation and provides detailed information on site-specific hydrophilicity and hydrophobicity.

Transient Non-local Interactions Dominate the Dynamics of Measles Virus NTAIL.

The RNA genome of measles virus is encapsidated by the nucleoprotein within a helical nucleocapsid that serves as template for both transcription and replication. The intrinsically disordered domain of the nucleoprotein (NTAIL), partly protruding outward from the nucleocapsid, is essential for binding the polymerase complex responsible for viral transcription and replication. As for many IDPs, binding of NTAIL occurs through a short molecular recognition element (MoRE) that folds upon binding, with the majority of NTAIL remaining disordered. Though NTAIL regions far from the MoRE influence the binding affinity, interactions between them and the MoRE have not been investigated in depth. Using an integrated approach, relying on photo-induced electron transfer (PET) experiments between tryptophan and cysteine pairs placed at different positions in the protein under varying salt and pH conditions, combined with simulations and analytical models, we identified transient interactions between two disordered regions distant in sequence, which dominate NTAIL dynamics, and regulate the conformational preferences of both the MoRE and the entire NTAIL domain. Co-evolutionary analysis corroborates our findings, and suggests an important functional role for the same intramolecular interactions. We propose mechanisms by which these non-local interactions may regulate binding to the phosphoprotein, polymerase recruitment, and ultimately viral transcription and replication. Our findings may be extended to other IDPs, where non-local intra-protein interactions affect the conformational preferences of intermolecular binding sites.

Mutations in the Receptor Binding Domain of Severe Acute Respiratory Coronavirus-2 Omicron Variant Spike Protein Significantly Stabilizes Its Conformation.

The Omicron variant and its sub-lineages are the only current circulating SARS-CoV-2 viruses worldwide. In this study, the conformational stability of the isolated Receptor Binding Domain (RBD) of Omicron's spike protein is examined in detail. The parent Omicron lineage has over ten mutations in the ACE2 binding region of the RBD that are specifically associated with its β hairpin loop domain. It is demonstrated through biophysical molecular computations that the mutations in the β hairpin loop domain significantly increase the intra-protein interaction energies of intra-loop and loop-RBD interactions. The interaction energy increases include the formation of new hydrogen bonds in the β hairpin loop domain that help stabilize this critical ACE2 binding region. Our results also agree with recent experiments on the stability of Omicron's core β barrel domain, outside of its loop domain, and help demonstrate the overall conformational stability of the Omicron RBD. It is further shown here through dynamic simulations that the unbound state of the Omicron RBD remains closely aligned with the bound state configuration, which was not observed for the wild-type RBD. Overall, these studies demonstrate the significantly increased conformational stability of Omicron over its wild-type configuration and raise a number of questions on whether conformational stability could be a positive selection feature of SARS-CoV-2 viral mutational changes.

Quantitative entropy-enthalpy compensation in intraprotein interactions from model compound data.

Many small globular proteins exist in only two states-the physiologically relevant folded state and an inactive unfolded state. The active state is stabilized by numerous weak attractive contacts, including hydrogen bonds, other polar interactions, and the hydrophobic effect. Knowledge of these interactions is key to understanding the fundamental equilibrium thermodynamics of protein folding and stability. We focus on one such interaction, that between amide and aromatic groups. We provide a statistically convincing case for quantitative, linear entropy-enthalpy compensation in forming aromatic-amide interactions using published model compound transfer-free energy data.

Nontraditional Roles of Magnesium Ions in Modulating Sav2152: Insight from a Haloacid Dehalogenase-like Superfamily Phosphatase from Staphylococcus aureus.

Methicillin-resistant Staphylococcus aureus (MRSA) infection has rapidly spread through various routes. A genomic analysis of clinical MRSA samples revealed an unknown protein, Sav2152, predicted to be a haloacid dehalogenase (HAD)-like hydrolase, making it a potential candidate for a novel drug target. In this study, we determined the crystal structure of Sav2152, which consists of a C2-type cap domain and a core domain. The core domain contains four motifs involved in phosphatase activity that depend on the presence of Mg2+ ions. Specifically, residues D10, D12, and D233, which closely correspond to key residues in structurally homolog proteins, are responsible for binding to the metal ion and are known to play critical roles in phosphatase activity. Our findings indicate that the Mg2+ ion known to stabilize local regions surrounding it, however, paradoxically, destabilizes the local region. Through mutant screening, we identified D10 and D12 as crucial residues for metal binding and maintaining structural stability via various uncharacterized intra-protein interactions, respectively. Substituting D10 with Ala effectively prevents the interaction with Mg2+ ions. The mutation of D12 disrupts important structural associations mediated by D12, leading to a decrease in the stability of Sav2152 and an enhancement in binding affinity to Mg2+ ions. Additionally, our study revealed that D237 can replace D12 and retain phosphatase activity. In summary, our work uncovers the novel role of metal ions in HAD-like phosphatase activity.

Prediction of the antigenic regions in eight RhD variants identified by computational biology.

Changes in RHD generate variations in protein structure that lead to antigenic variants. The classical model divides them into quantitative (weak and Del) and qualitative (partial D). There are two types of protein antigens: linear and conformational. Computational biology analyses the theoretical assembly of tertiary protein structures and allows us to identify the 'topological' differences between isoforms. Our aim was to determine the theoretical antigenic differences between weak RhD variants compared with normal RhD based on structural analysis using bioinformatic techniques. We analysed the variations in secondary structures and hydrophobicity of RHD*01, RHD*01W.1, W2, W3, RHD*09.03.01, RHD*09.04, RHD*11, RHD*15 and RHD*21. We then modelled the tertiary structure and calculated their probable antigenic regions, intra-protein interactions, displacement and membrane width and compared them with Rhce. The 10 proteins are similar in their secondary structure and hydrophobicity, with the main differences observed in the exofacial coils. We identified six potential antigenic regions: one that is unique to RhD (R3), one that is common to all D (R6), three that are highly variable among RhD isoforms (R1, R2 and R4), one that they share with Rhce (R5) and two that are unique to Rhce (Ra and Rbc). The alloimmunization capacity of these subjects could be explained by the variability of the antigen pattern, which is not necessarily recognized or recognized with lower intensity by the commercially available antibodies, and not because they have a lower protein concentration in the membrane.

Design of stable circular permutants of the GroEL chaperone apical domain

Enhancing protein stability holds paramount significance in biotechnology, therapeutics, and the food industry. Circular permutations offer a distinctive avenue for manipulating protein stability while keeping intra-protein interactions intact. Amidst the creation of circular permutants, determining the optimal placement of the new N- and C-termini stands as a pivotal, albeit largely unexplored, endeavor. In this study, we employed PONDR-FIT’s predictions of disorder propensity to guide the design of circular permutants for the GroEL apical domain (residues 191–345). Our underlying hypothesis posited that a higher predicted disorder value would correspond to reduced stability in the circular permutants, owing to the increased likelihood of fluctuations in the novel N- and C-termini. To substantiate this hypothesis, we engineered six circular permutants, positioning glycines within the loops as locations for the new N- and C-termini. We demonstrated the validity of our hypothesis along the set of the designed circular permutants, as supported by measurements of melting temperatures by circular dichroism and differential scanning microcalorimetry. Consequently, we propose a novel computational methodology that rationalizes the design of circular permutants with projected stability.8ghBPrif1HpKfHNvAL1G2xVideo

Molecular Characterization of Alpha-Amylase of Aspergillus flavus and Aspergillus oryzae: An In Silico Study

Alpha amylases (EC 3.2.1.1.) hydrolyze polysaccharides like starch and glycogen to create oligosaccharides of different sizes by random cleavage of internal 1, 4-glucosidic linkages. These enzymes are produced by many different bacteria and fungi. In this investigation, α-amylase from Aspergillus flavus and Aspergillus oryzae was taken into consideration. Amino acid sequences and protein structure were retrieved from databases. Sequences of Aspergillus oryzae possess a high number of charged residues whereas Aspergillus flavus have a higher abundance of uncharged polar and hydrophobic amino acid residues. 3 common superfamilies were observed between both species. The alignment of sequences showed their similar conservative nature. Non-covalent intra-protein interactions like the salt bridge, aromatic-aromatic interactions, aromatic-sulfur interactions and cation-pi interactions helped to increase the stability of α-amylase of Aspergillus oryzae. The tunnels and cavities also help to increase the functionality and catalytic activity of α-amylase. This study will also be helpful for protein engineering.

Elucidating the microscopic origin of observed anti-correlation between the protein-protein and protein-water interaction energies

Elucidating the microscopic origin of observed anti-correlation between the protein-protein and protein-water interaction energies

Comparative analysis of permanent and transient domain-domain interactions in multi-domain proteins.

Protein domains are structural, functional, and evolutionary units. These domains bring out the diversity of functionality by means of interactions with other co-existing domains and provide stability. Hence, it is important to study intra-protein inter-domain interactions from the perspective of types of interactions. Domains within a chain could interact over short timeframes or permanently, rather like protein-protein interactions (PPIs). However, no systematic study has been carried out between two classes, namely permanent and transient domain-domain interactions. In this work, we studied 263 two-domain proteins, belonging to either of these classes and their interfaces on the basis of several factors, such as interface area and details of interactions (number, strength, and types of interactions). We also characterized them based on residue conservation at the interface, correlation of residue motions across domains, its involvement in repeat formation, and their involvement in particular molecular processes. Finally, we could analyze the interactions arising from domains in two-domain monomeric proteins, and we observed significant differences between these two classes of domain interactions and a few similarities. This study will help to obtain a better understanding of structure-function and folding principles of multi-domain proteins.

P2Y12 receptor residues crucial for thrombosis regulation.

Thrombosis, or the formation of blood clots, can lead to serious medical conditions such as stroke, heart attack, and deep vein thrombosis. The purinoreceptor P2Y12 plays a critical role in the thrombotic pathway and is targeted for therapy to prevent clot formation. However, it is essential to balance the regulation of thrombosis to avoid adverse situations. This study focuses on the P2Y12 receptor and aims to discern the protein residue network and differentiate residues based on their intramolecular interactions. The study utilized a statistical analysis to characterize the significant residues involved in ligand interaction, which helps to identify critical residues that are essential for the function of the receptor. A parametric analysis of interactions of residues in the intraprotein interaction was conducted, which revealed significant residue-based contacts that facilitate protein interactions. By examining the interactions between residues, the mechanisms underlying protein interactions were studied and the importance of specific residues in facilitating these interactions was determined. This research provides important information on P2Y12, and the findings based on the networkbased significance of interacting residues may contribute to the development of new therapies that target the receptor to prevent clot formation while maintaining a balance in thrombosis regulation to avoid adverse outcomes. Ultimately, this study could lead to improved treatments for thrombotic disorders and better patient outcomes.

Model-Dependent Solvation of the K-18 Domain of the Intrinsically Disordered Protein Tau.

A known imbalance between intra-protein and protein-water interactions in many empirical force fields results in collapsed conformational ensembles of intrinsically disordered proteins in explicit solvent simulations that disagree with experiments. Multiple strategies have been introduced in the literature to modify protein-water interactions, which improve agreement between experiments and simulations. In this work, we combine simulations with standard and modified force fields with a spatially resolved analysis of solvation free energy contributions and compare the consequences of each strategy. We find that enhanced Lennard-Jones (LJ) interactions between protein atoms and water oxygens primarily improve the solvation of nonpolar functional groups of the protein. In contrast, modified electrostatics in the water model or strengthened LJ interactions between the protein and water hydrogens mainly affect the hydration of polar functional groups. Modified electrostatics further impact the average orientation of water molecules in the hydration shell. As a result, protein-water interactions with the first hydration layers are strengthened, while interactions with water molecules in higher hydration shells are weakened. Hence, distinct strategies to balance intra-protein and protein-water interactions in simulations have qualitatively different effects on protein solvation. These differences are not necessarily captured by comparisons to experiments that report on global parameters describing protein conformational ensembles, e.g., the radius of gyration, but will influence the tendency of a protein to form aggregates or phase-separated droplets.

Identification of fungus-growing termite-associated halogenated-PKS maduralactomycin a as a potential inhibitor of MurF protein of multidrug-resistant Acinetobacter baumannii.

Multidrug-resistant Acinetobacter baumannii infections have become a major public health concern globally. Inhibition of its essential MurF protein has been proposed as a potential target for broad-spectrum drugs. This study aimed to evaluate the potential of a novel ecological niche of 374 fungus-growing termite associated Natural Products (NPs). The molecular docking and computational pharmacokinetics screened four compounds, i.e., Termstrin B, Fridamycin A, Maduralactomycin A, and Natalenamide C, as potential compounds that have higher binding affinities and favourable protein-ligand interactions. The compound Maduralactomycin A induced more stability based on its lowest average RMSD value (2.31Å) and low standard deviation (0.35) supported by the consistent flexibility and β-factor during the protein's time-dependent motion. While hydrogen bond analysis indicated that Termstrin B has formed the strongest intra-protein interaction, solvent accessibility was in good agreement with Maduralactomycin A compactness. Maduralactomycin A has the strongest binding energy among all the compounds (-348.48kcal/mol) followed by Termstrin B (-321.19kcal/mol). Since these findings suggest Maduralactomycin A and Termstrin B as promising candidates for inhibition of MurF protein, the favourable binding energies of Maduralactomycin A make it a more important compound to warrant further investigation. However, experimental validation using animal models and clinical trials is recommended before reaching any final conclusions.

Comparative molecular dynamics simulations of pathogenic and non-pathogenic huntingtin protein monomers and dimers.

Polyglutamine expansion at the N-terminus of the huntingtin protein exon 1 (Htt-ex1) is closely associated with a number of neurodegenerative diseases, which result from the aggregation of the increased polyQ repeat. However, the underlying structures and aggregation mechanism are still poorly understood. We performed microsecond-long all-atom molecular dynamics simulations to study the folding and dimerization of Htt-ex1 (about 100 residues) with non-pathogenic and pathogenic polyQ lengths, and uncovered substantial differences. The non-pathogenic monomer adopts a long α-helix that includes most of the polyQ residues, which forms the interaction interface for dimerization, and a PPII-turn-PPII motif in the proline-rich region. In the pathogenic monomer, the polyQ region is disordered, leading to compact structures with many intra-protein interactions and the formation of short β-sheets. Dimerization can proceed via different modes, where those involving the N-terminal headpiece bury more hydrophobic residues and are thus more stable. Moreover, in the pathogenic Htt-ex1 dimers the proline-rich region interacts with the polyQ region, which slows the formation of β-sheets.

A novel high-prevalence antigen in the Lutheran system, LUGA (LU24), and an updated, full-length 3D BCAM model.

The basal cell adhesion molecule (BCAM) carries the antigens of the Lutheran (LU, ISBT005) system. We report a novel Lutheran antigen and propose an updated, full-length 3D model of BCAM. Red blood cell testing, antibody identification, and BCAM genomic DNA sequencing were done by standard methods. Multi-template homology modeling of BCAM used structural templates selected for coverage, highest sequence identity, and protein domain family. All variants causing the loss or gain of a Lutheran antigen were analyzed for residue accessibility and intraprotein interactions. An antibody to a high-prevalence antigen in the plasma of a pregnant woman was determined to be directed at a novel Lutheran antigen. Sequencing of BCAM found three homozygous changes: c.212G > A (p.Arg71His) and two silent, c.711C > T and c.714C > T. The model was built from the first two immunoglobulin crystallized domains of BCAM (D1, D2), three other templates (for D3, D4 and D5 with a higher sequence identity with the target than those used for the model proposed by Burton and Brady in 2008, and for the transmembrane region) and RaptorX (for the intracellular domain). All residues associated with a Lutheran antigen were found to be exposed in wild-type or variant proteins, except p.447 associated with loss of Lu13 expression. The c.212G > A change results in the loss of LUGA (LU24) antigen. Whole genome sequencing continues to reveal polymorphisms with uncertain immunogenicity. This model and demonstration that nearly all residues associated with the expression of a Lutheran antigen are exposed will help evaluate the significance of new polymorphisms.

Fingerprinting full-length proteins through single-file translocations.

Fingerprinting full-length proteins through single-file translocations.

Role of gRNA binding in HIV-1 Gag assembly process.

Role of gRNA binding in HIV-1 Gag assembly process.

Investigating the interaction of TDP-43 with RNA using structural mass spectrometry.

Investigating the interaction of TDP-43 with RNA using structural mass spectrometry.

Two new Scianna variants causing loss of high prevalence antigens: ERMAPmodel and 3D analysis of the antigens.

Scianna (Sc) antigens, seven high and two of low prevalence, are expressed on erythrocyte membrane-associated protein (ERMAP). We investigated SC (ERMAP) in individuals who made antibodies to high prevalence Scianna antigens, and propose a 3D model for ERMAP to precisely localize the residues associated with the known antigens. Serological testing and DNA sequencing was performed by standard methods. A 3D structural model was built using a multi-template homology approach. Protein structures representing missense variants associated with the loss or gain of an antigen were generated. Residue accessibility and intraprotein interactions were compared with the wild-type protein. Two new SC alleles, one with c.349C > T (p.Arg117Cys) in a woman from South India with anti-Sc3 in her plasma, and a c.217_219delinsTGT (p.Arg73Cys) in an African-American woman with an antibody to a new high prevalence antigen, termed SCAC, were identified. Six structural templates were used to model ERMAP. 3D analysis showed that residues key for Scianna antigen expression were all exposed at the surface of the extracellular domain. The p.Arg117Cys change was predicted to abolish interactions between residues 93 and 117, with no compensating interactions. We confirm the extracellular location of Scianna residues responsible for antigen expression which predicts direct accessibility to antibodies. Loss of intraprotein interactions appear to be responsible for a Sc null and production of anti-Sc3 with p.117Cys, SC*01 N.03, and for loss of a high prevalence antigen with p.73Cys, termed SCAC for Sc Arg to Cys. Comparative modeling aids our understanding of new alleles and Scianna antigen expression.

3D analysis of CROMER (DAF) and a new antigen CRAG.

Cromer antigens are carried on decay accelerating factor (DAF, CD55), for which the crystal structure is available. We investigated two samples with an unidentified antibody to a high prevalence antigen and evaluate the location and characteristics of amino acids associated with antigens on the CD55 by 3D modelling. Antigen typing and antibody identification were by standard methods. CD55 was sequenced, and Cromer variants were generated using the protein's crystal structure (1OK3, chain A). Antigen-associated residues and intraprotein interactions were investigated in 3D (Naccess, Protein Interactions Calculator). The antibody in the sample from a woman of Kashmiri descent was identified as anti-IFC (anti-CROM7). Her RBCs were negative for high-prevalence Cromer antigens including IFC. CD55 sequencing revealed a silent c.147G>A (p.Leu49=) and c.148G>T (p.Glu50Ter) changes, designated CROM*01N.05. The antibody in the sample from a woman of Greek ancestry was only compatible with IFC- RBCs but her RBCs were positive for known high-prevalence Cromer antigens. CD55 sequencing found she was homozygous for c.173A>G (p.Asp58Gly). The high prevalence antigen was named CRAG (ISBT CROM18 or 021018) and the allele designated CROM*01.-18. By 3D analysis, all known antigen-associated residues, including the new CRAG antigen, were exposed at the protein surface. Interactions between antigen-associated residues within the same CD55 domains were identified. Identification of antibodies to high prevalence Cromer antigens can be challenging. The surface exposure of antigen-associated residues likely accounts for their immunogenicity. 3D analysis of CD55 provides insight into previous serologic observations regarding the influence of some Cromer antigens on the expression of others.