Nonpolar Residues Research Articles

Destabilisation of Alzheimer's amyloid-β protofibrils by Baicalein: mechanistic insights from all-atom molecular dynamics simulations.

Alzheimer's disease (AD) is the most common form of dementia and the fifth leading cause of death globally. Aggregation and deposition of neurotoxic Aβ fibrils in the neural tissues of the brain is a key hallmark in AD pathogenesis. Destabilisation studies of the amyloid-peptide by various natural molecules are highly relevant due to their neuroprotective and therapeutic potential for AD. We performed molecular dynamics (MD) simulation to investigate the destabilisation mechanism of amyloidogenic protofilament intermediate by Baicalein (BCL), a naturally occurring flavonoid. We found that the BCL molecule formed strong hydrophobic contacts with non-polar residues, specifically F19, A21, V24, and I32 of Chain A and B of the pentameric protofibril. Upon binding, it competed with the native hydrophobic contacts of the Aβ protein. BCL loosened the tight packing of the hydrophobic core by disrupting the hydrogen bonds and the prominent D23-K28 inter-chain salt bridges of the protofibril. The decrease in the structural stability of Aβ protofibrils was confirmed by the increased RMSD, radius of gyration, solvent accessible surface area (SASA), and reduced β-sheet content. PCA indicated that the presence of the BCL molecule intensified protofibril motions, particularly affecting residues in Chain A and B regions. Our findings propose that BCL would be a potent destabiliser of Aβ protofilament, and may be considered as a therapeutic agent in treating AD.

Statistical analysis of the unique characteristics of secondary structures in proteins

Protein folding is a complex process influenced by the primary sequence of amino acids. Early studies focused on understanding whether the specificity or the conservation of properties of amino acids was crucial for folding into secondary structures such as α-helices, β-sheets, turns, and coils. However, with the advent of artificial intelligence (AI) and machine learning (ML), the emphasis has shifted towards the precise nature and occurrence of specific amino acids. In our study, we analyzed a large set of proteins from diverse organisms to identify unique features of secondary structures, particularly in terms of the distribution of polar, non-polar, and charged amino acid residues. We found that α-helices tend to have a higher proportion of charged and non-polar groups compared to other secondary structures and that the presence of oppositely charged amino acid residues in helices stabilizes them, facilitating the formation of longer helices. These characteristics are distinct to α-helices. This study offers valuable insights for researchers in the field of protein design, enabling the de-novo creation of short helical peptides for a range of applications. We have also developed a web server for extensive analysis of proteins from different databases. The web server is housed at https://proseqanalyser.iitgn.ac.in/

Polymerized whey protein-SDS interactions at their high concentrations

Protein-surfactant interactions have been an ongoing topic of interest for many decades. Applications involving complexes of proteins and surfactants are relevant in food, pharmaceuticals, hygiene, molecular characterization, and other fields. In this study, the interactions of polymerized whey proteins (PWP) and sodium dodecyl sulfate (SDS) at high concentrations are investigated. Different characterization techniques are used, including electrical conductivity, turbidity, isothermal titration calorimetry, dynamic light scattering, dilute solution viscometry, rheology, and surface hydrophobicity to elucidate information on the modes and extent of interactions. Results indicate that PWP-SDS interactions produce highly extended, worm like micelles, with SDS decorating PWP chains and covering non-polar residues. PWP can host SDS up to quite high surfactant to protein ratios (SPR), producing solutions that are highly viscous with shear thickening properties, yet with no networking or gelation. Interestingly, dilution of high viscosity PWP-SDS solutions leads to even smaller size of PWP-SDS molecular complex as compared with PWP without SDS. The current study extends the vision of protein surfactant interactions by examining concentration range beyond that found in literature. The results reveal insights that can help expand studies on other systems and find applications in various fields including coatings, cosmetics, food ingredients, drug transport, and disease treatment.

Enzyme Tunnel Dynamics and Catalytic Mechanism of Norcoclaurine Synthase: Insights from a Combined LiGaMD and DFT Study.

This study conducts a systematic investigation into the catalytic mechanism of norcoclaurine synthase (NCS), a key enzyme in the biosynthesis of tetrahydroisoquinolines (THIQs) with therapeutic applications. By integration of LiGaMD and DFT calculations, the reaction pathway of NCS is mapped, providing detailed insights into its catalytic activity and selectivity. Our findings underscore the critical role of E103 in substrate capture and reveal the hitherto unappreciated influence of nonpolar residues M183 and L76 on tunnel dynamics. A prominent discovery is the identification of a high-energy barrier (44.2 kcal/mol) associated with the aromatic electrophilic attack, which pinpoints the rate-limiting step. Moreover, we disclose the existence of dual transition states leading to different products with the energetically favored six-membered ring formation consistent with experimental evidence. These mechanistic revelations not only refine our understanding of NCS but also advocate for a renewed emphasis on enzyme tunnel engineering for optimizing THIQs biosynthesis. The research sets the stage for translating these findings into practical enzyme modifications. Our results highlight the potential of NCS as a biocatalyst to overcome the limitations of current synthetic methodologies, such as low yields and environmental impacts, and provide a theoretical contribution to the efficient, eco-friendly production of THIQs-based pharmaceuticals.

The effect of phosphorylation efficiency on the oncogenic properties of the protein E7 from high-risk HPV

The Human papillomavirus (HPV) causes tumors in part by hijacking the host cell cycle and forcing uncontrolled cellular division. While there are >200 genotypes of HPV, 15 are classified as high-risk and have been shown to transform infected cells and contribute to tumor formation. The remaining low-risk genotypes are not considered oncogenic and result in benign skin lesions. In high-risk HPV, the oncoprotein E7 contributes to the dysregulation of cell cycle regulatory mechanisms. High-risk E7 is phosphorylated in cells at two conserved serine residues by Casein Kinase 2 (CK2) and this phosphorylation event increases binding affinity for cellular proteins such as the tumor suppressor retinoblastoma (pRb). While low-risk E7 possesses similar serine residues, it is phosphorylated to a lesser degree in cells and has decreased binding capabilities. When E7 binding affinity is decreased, it is less able to facilitate complex interactions between proteins and therefore has less capability to dysregulate the cell cycle. By comparing E7 protein sequences from both low- and high-risk HPV variants and using site-directed mutagenesis combined with NMR spectroscopy and cell-based assays, we demonstrate that the presence of two key nonpolar valine residues within the CK2 recognition sequence, present in low-risk E7, reduces serine phosphorylation efficiency relative to high-risk E7. This results in significant loss of the ability of E7 to degrade the retinoblastoma tumor suppressor protein, thus also reducing the ability of E7 to increase cellular proliferation and reduce senescence. This provides additional insight into the differential E7-mediated outcomes when cells are infected with high-risk verses low-risk HPV. Understanding these oncogenic differences may be important to developing targeted treatment options for HPV-induced cancers.

Single Point Mutation Abolishes Water Capture in Germacradien-4-ol Synthase.

The high-fidelity sesquiterpene cyclase (-)-germacradien-4-ol synthase (GdolS) converts farnesyl diphosphate into the macrocyclic alcohol (-)-germacradien-4-ol. Site-directed mutagenesis was used to decipher the role of key residues in the water control mechanism. Replacement of Ala176, located in the G1/2 helix, with non-polar aliphatic residues of increasing size (valine, leucine, isoleucine and methionine) resulted in the accumulation of the non-hydroxylated products germacrene A and germacrene D. In contrast, hydroxylation was maintained when the polar residues threonine, glutamine or aspartate replaced Ala176. Additionally, although a contribution of His150 to the nucleophilic water addition could be ruled out, the imidazole ring of His150 appears to assist carbocation stabilisation. The results presented here shed light on how hydroxylating sesquiterpene synthases can be engineered to design modified sesquiterpene synthases to reduce the need for further steps in the biocatalytic production of oxygenated sesquiterpenoids.

A molecular synergy of non-covalent interactions facilitated via theobromine modifying the self-assembly process of type I collagen

Type I Collagen is a multifunctional fibrous protein present significantly in the connective tissues. Generally, fibril formation process of an interstitial collagen begins with a triple helix alignment to form a pentafibril, which sequentially self-assembles into a microfibril. Typically, this sequential formation can be achieved under an in vitro condition as well by maintaining certain specific physiological parameters to mimic the interstitial microenvironment of collagen. The necessity to study such a kinetically driven process is due to its implications in template fabrication particularly as a scaffold to enhance the growth of specific stem cells or accelerating nanoparticle synthesis. Drawing inspiration from this, theobromine, a methylxanthine that contributes towards the hydrogen bond network via dimerization and stacking interaction was used as a small molecule in modifying the molecular level interactions of a self-assembled collagen. Initially, rheological studies revealed certain variations in the viscous properties of collagen on addition of theobromine. Turbidity assay showed an increase in the rate of fibril formation and its topographical changes were observed through the high-resolution scanning electron microscope. The stability and probable physico-chemical interactions between collagen and theobromine were comprehended from the vibrational frequency changes recorded using (FT-IR) (ATR) and conformational changes using Circular Dichroism (CD) studies. Our results show possible multibinding sites on theobromine, which would compete to bind with the polar and non-polar residues exposed over the monomeric surfaces of triple helix causing a local conformational fluctuation modifying protein-protein interaction.

Effect of the Lys62Ala Mutation on the Thermal Stability of BstHPr Protein by Molecular Dynamics.

We analyzed the thermal stability of the BstHPr protein through the site-directed point mutation Lys62 replaced by Ala residue using molecular dynamics simulations at five different temperatures: 298, 333, 362, 400, and 450 K, for periods of 1 μs and in triplicate. The results from the mutant thermophilic BstHPrm protein were compared with those of the wild-type thermophilic BstHPr protein and the mesophilic BsHPr protein. Structural and molecular interaction analyses show that proteins lose stability as temperature increases. Mutant and wild-type proteins behave similarly up to 362 K. However, at 400 K the mutant protein shows greater structural instability, losing more buried hydrogen bonds and exposing more of its non-polar residues to the solvent. Therefore, in this study, we confirmed that the salt bridge network of the Glu3-Lys62-Glu36 triad, made up of the Glu3-Lys62 and Glu36-Lys62 ion pairs, provides thermal stability to the thermophilic BstHPr protein.

Structural Determinants of Peptide Nanopore Formation.

We have evolved the nanopore-forming macrolittin peptides from the bee venom peptide melittin using successive generations of synthetic molecular evolution. Despite their sequence similarity to the broadly membrane permeabilizing cytolytic melittin, the macrolittins have potent membrane selectivity. They form nanopores in synthetic bilayers made from 1-palmitoyl, 2-oleoyl-phosphatidylcholine (POPC) at extremely low peptide concentrations and yet have essentially no cytolytic activity against any cell membrane, even at high concentration. Here, we explore the structural determinants of macrolittin nanopore stability in POPC bilayers using atomistic molecular dynamics simulations and experiments on macrolittins and single-site variants. Simulations of macrolittin nanopores in POPC bilayers show that they are stabilized by an extensive, cooperative hydrogen bond network comprised of the many charged and polar side chains interacting with each other via bridges of water molecules and lipid headgroups. Lipid molecules with unusual conformations participate in the H-bond network and are an integral part of the nanopore structure. To explore the role of this H-bond network on membrane selectivity, we swapped three critical polar residues with the nonpolar residues found in melittin. All variants have potency, membrane selectivity, and cytotoxicity that were intermediate between a cytotoxic melittin variant called MelP5 and the macrolittins. Simulations showed that the variants had less organized H-bond networks of waters and lipids with unusual structures. The membrane-spanning, cooperative H-bond network is a critical determinant of macrolittin nanopore stability and membrane selectivity. The results described here will help guide the future design and optimization of peptide nanopore-based applications.

Comprehensive analysis of CXXX sequence space reveals that Saccharomyces cerevisiae GGTase-I mainly relies on a2X substrate determinants.

Many proteins undergo a post-translational lipid attachment, which increases their hydrophobicity, thus strengthening their membrane association properties or aiding in protein interactions. Geranylgeranyltransferase-I (GGTase-I) is an enzyme involved in a 3-step post-translational modification (PTM) pathway that attaches a 20-carbon lipid group called geranylgeranyl at the carboxy-terminal cysteine of proteins ending in a canonical CaaL motif (C-cysteine, a-aliphatic, L-often leucine, but can be phenylalanine, isoleucine, methionine, or valine). Genetic approaches involving 2 distinct reporters were employed in this study to assess Saccharomyces cerevisiae GGTase-I specificity, for which limited data exist, toward all 8,000 CXXX combinations. Orthogonal biochemical analyses and structure-based alignments were also performed to better understand the features required for optimal target interaction. These approaches indicate that yeast GGTase-I best modifies the Cxa[L/F/I/M/V] sequence that resembles but is not an exact match for the canonical CaaL motif. We also observed that minor modification of noncanonical sequences is possible. A consistent feature associated with well-modified sequences was the presence of a nonpolar a2 residue and a hydrophobic terminal residue, which are features recognized by mammalian GGTase-I. These results thus support that mammalian and yeast GGTase-I exhibit considerable shared specificity.

Determination of chlorothalonil in vegetables/fruits using liquid chromatography-tandem mass spectrometry with a discharge adapter interface.

Chlorothalonil (CHT), a broad-spectrum fungicide, has been employed widely to control foliar diseases, whereas with a major metabolite of polar 4-hydroxychlorothalonil (CHT-4-OH), only an acceptable nonpolar CHT residue is allowed by most countries. This study involves the method development for CHT residue in vegetables/fruits using liquid chromatography-tandem mass spectrometry (LC-MS/MS) with a novel modified discharge-adaptor (DA) interface. CHT residue was analyzed using LC-MS/MS with DA interface (LC-DA-MS/MS), developed in our previous works. A DA was placed on the electrospray tip to switch the ionization modes. A modified quick, easy, cheap, effective, rugged, and safe (QuEChERS) method was applied to extract CHT residue of vegetables/fruits efficiently with less sample preparation time and analysis cost. CHT and CHT-4-OH spiked in four different vegetables/fruits were extracted using the modified QuEChERS method. After LC with isocratic elution, CHT and CHT-4-OH were separated within 3min. Using LC-DA-MS/MS, the ion signals of CHT were improved two to three times, and the limit of quantification of 5ng/g and linearity (r2 > 0.99) in the range of 5-200 ng/g were achieved using 10g of vegetables/fruits. The precision and accuracy were within 15% each. The modified QuEChERS and LC-DA-MS/MS were applied to examine eight field-grown vegetables/fruits; 9.5 and 2588.9ng/g of CHT were detected in two vegetables/fruits. LC-DA-MS/MS combined with modified QuEChERS was successfully applied to determine CHT residue <10ng/g in vegetables/fruits and with satisfied validation results. The developed method could reduce both analysis cost and time, attributing to simplifications in modified QuEChERS, isocratic elution, and DA interface in LC-DA-MS/MS.

Microfibril-associated glycoprotein 4 forms octamers that mediate interactions with elastogenic proteins and cells

Microfibril-associated glycoprotein 4 (MFAP4) is a 36-kDa extracellular matrix glycoprotein with critical roles in organ fibrosis, chronic obstructive pulmonary disease, and cardiovascular disorders, including aortic aneurysms. MFAP4 multimerises and interacts with elastogenic proteins, including fibrillin-1 and tropoelastin, and with cells via integrins. Structural details of MFAP4 and its potential interfaces for these interactions are unknown. Here, we present a cryo-electron microscopy structure of human MFAP4. In the presence of calcium, MFAP4 assembles as an octamer, where two sets of homodimers constitute the top and bottom halves of each octamer. Each homodimer is linked together by an intermolecular disulphide bond. A C34S missense mutation prevents disulphide-bond formation between monomers but does not prevent octamer assembly. The atomic model, built into the 3.55 Å cryo-EM map, suggests that salt-bridge interactions mediate homodimer assembly, while non-polar residues form the interface between octamer halves. In the absence of calcium, an MFAP4 octamer dissociates into two tetramers. Binding studies with fibrillin-1, tropoelastin, LTBP4, and small fibulins show that MFAP4 has multiple surfaces for protein-protein interactions, most of which depend upon MFAP4 octamer assembly. The C34S mutation does not affect these protein interactions or cell interactions. MFAP4 assemblies with fibrillin-1 abrogate MFAP4 interactions with cells.

Molecular basis for differential recognition of an allosteric inhibitor by receptor tyrosine kinases.

Understanding kinase-inhibitor selectivity continues to be a major objective in kinase drug discovery. We probe the molecular basis of selectivity of an allosteric inhibitor (MSC1609119A-1) of the insulin-like growth factor-I receptor kinase (IGF1RK), which has been shown to be ineffective for the homologous insulin receptor kinase (IRK). Specifically, we investigated the structural and energetic basis of the allosteric binding of this inhibitor to each kinase by combining molecular modeling, molecular dynamics (MD) simulations, and thermodynamic calculations. We predict the inhibitor conformation in the binding pocket of IRK and highlight that the charged residues in the histidine-arginine-aspartic acid (HRD) and aspartic acid-phenylalanine-glycine (DFG) motifs and the nonpolar residues in the binding pocket govern inhibitor interactions in the allosteric pocket of each kinase. We suggest that the conformational changes in the IGF1RK residues M1054 and M1079, movement of the ⍺C-helix, and the conformational stabilization of the DFG motif favor the selectivity of the inhibitor toward IGF1RK. Our thermodynamic calculations reveal that the observed selectivity can be rationalized through differences observed in the electrostatic interaction energy of the inhibitor in each inhibitor/kinase complex and the hydrogen bonding interactions of the inhibitor with the residue V1063 in IGF1RK that are not attained with the corresponding residue V1060 in IRK. Overall, our study provides a rationale for the molecular basis of recognition of this allosteric inhibitor by IGF1RK and IRK, which is potentially useful in developing novel inhibitors with improved affinity and selectivity.

Binding characteristics of pheromone-binding protein 1 in Glyphodes pyloalis to organophosphorus insecticides: Insights from computational and experimental approaches.

Glyphodes pyloalis (Lepidoptera: Pyralidae) is one of the major pests in mulberry production in China, which has developed resistance to various insecticides. Chemoreception is one of the most crucial physiological tactics in insects, playing a pivotal role in recognizing chemical stimuli in the environment, including noxious stimuli such as insecticides. Herein, we obtained recombinant pheromone-binding protein 1 (GpylPBP1) that exhibited antennae-biased expression in G. pyloalis. Ligand-binding assays indicated that GpylPBP1 had the binding affinities to two organophosphorus insecticides, with a higher binding affinity to chlorpyrifos than to phoxim. Computational simulations showed that a mass of nonpolar amino acid residues formed the binding pocket of GpylPBP1 and contributed to the hydrophobic interactions in the bindings of GpylPBP1 to both insecticides. Furthermore, the binding affinities of three GpylPBP1 mutants (F12A, I52A, and F118A) to both insecticides were all significantly reduced compared to those of the GpylPBP1-wild type, suggesting that Phe12, Ile52, and Phe118 residues were crucial binding sites and played crucial roles in the bindings of GpylPBP1 to both insecticides. Our findings can be instrumental in elucidating the effects of insecticides on olfactory recognition in moths and facilitating the development of novel pest management strategies using PBPs as targets based on insect olfaction.

GRAS gene family in rye (Secale cereale L.): genome-wide identification, phylogeny, evolutionary expansion and expression analyses

BackgroundThe GRAS transcription factor family plays a crucial role in various biological processes in different plants, such as tissue development, fruit maturation, and environmental stress. However, the GRAS family in rye has not been systematically analyzed yet.ResultsIn this study, 67 GRAS genes in S. cereale were identified and named based on the chromosomal location. The gene structures, conserved motifs, cis-acting elements, gene replications, and expression patterns were further analyzed. These 67 ScGRAS members are divided into 13 subfamilies. All members include the LHR I, VHIID, LHR II, PFYRE, and SAW domains, and some nonpolar hydrophobic amino acid residues may undergo cross-substitution in the VHIID region. Interested, tandem duplications may have a more important contribution, which distinguishes them from other monocotyledonous plants. To further investigate the evolutionary relationship of the GRAS family, we constructed six comparative genomic maps of homologous genes between rye and different representative monocotyledonous and dicotyledonous plants. The response characteristics of 19 ScGRAS members from different subfamilies to different tissues, grains at filling stages, and different abiotic stresses of rye were systematically analyzed. Paclobutrazol, a triazole-based plant growth regulator, controls plant tissue and grain development by inhibiting gibberellic acid (GA) biosynthesis through the regulation of DELLA proteins. Exogenous spraying of paclobutrazol significantly reduced the plant height but was beneficial for increasing the weight of 1000 grains of rye. Treatment with paclobutrazol, significantly reduced gibberellin levels in grain in the filling period, caused significant alteration in the expression of the DELLA subfamily gene members. Furthermore, our findings with respect to genes, ScGRAS46 and ScGRAS60, suggest that these two family members could be further used for functional characterization studies in basic research and in breeding programmes for crop improvement.ConclusionsWe identified 67 ScGRAS genes in rye and further analysed the evolution and expression patterns of the encoded proteins. This study will be helpful for further analysing the functional characteristics of ScGRAS genes.

Conformational control of antimicrobial peptide amphiphilicity: consequences for boosting membrane interactions and antimicrobial effects of photocatalytic TiO2 nanoparticles.

This study reports on the effects of conformationally controlled amphiphilicity of antimicrobial peptides (AMPs) on their ability to coat TiO2 nanoparticles (NPs) and boost the photocatalytic antimicrobial effects of such NPs. For this, TiO2 NPs were combined with AMP EFK17 (EFKRIVQRIKDFLRNLV), displaying a disordered conformation in aqueous solution but helix formation on interaction with bacterial membranes. The membrane-bound helix is amphiphilic, with all polar and charged amino acid residues located at one side and all non-polar and hydrophobic residues on the other. In contrast, the d-enantiomer variant EFK17-d (E(dF)KR(dI)VQR(dI)KD(dF)LRNLV) is unable to form the amphiphilic helix on bacterial membrane interaction, whereas the W-residues in EFK17-W (EWKRWVQRWKDFLRNLV) boost hydrophobic interactions of the amphiphilic helix. Circular dichroism results showed the effects displayed for the free peptide, to also be present for peptide-coated TiO2 NPs, causing peptide binding to decrease in the order EFK17-W > EFK17 > EFK17-d. Notably, the formation of reactive oxygen species (ROS) by the TiO2 NPs was essentially unaffected by the presence of peptide coating, for all the peptides investigated, and the coatings stabilized over hours of UV exposure. Photocatalytic membrane degradation from TiO2 NPs coated with EFK17-W and EFK17 was promoted for bacteria-like model bilayers containing anionic phosphatidylglycerol but suppressed in mammalian-like bilayers formed by zwitterionic phosphatidylcholine and cholesterol. Structural aspects of these effects were further investigated by neutron reflectometry with clear variations observed between the bacteria- and mammalian-like model bilayers for the three peptides. Mirroring these results in bacteria-like model membranes, combining TiO2 NPs with EFK17-W and EFK17, but not with non-adsorbing EFK17-d, resulted in boosted antimicrobial effects of the resulting cationic composite NPs already in darkness, effects enhanced further on UV illumination.

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

Water as a solvent plays an important role in dictating the structure and dynamics of the large biomolecules like proteins. It has been proposed that the protein dynamics and energetics are closely coupled to the water molecules in the hydration shell. In particular, it has been shown earlier that the fluctuations in intra-protein (protein-protein) interaction energy is strongly anti-correlated with the protein-water interaction energy (solvation energy). Moreover, the dynamical signatures of water gets reflected into the protein dynamics in terms of a dominant power law behavior. In this work, we have carried out systematic molecular dynamics (MD) simulations of three different protein systems to understand the key physical factors that contribute to such anti-correlations. First, we compare the dynamics in globular folded proteins like Lysozyme with intrinsically disordered peptides (IDPs) like Amyloid beta to show that enhanced conformational fluctuations lead to significantly stronger anti-correlation. Trp-cage miniprotein shows a similar behavior distributed over two sub-ensembles of folded and unfolded states. We also investigate what type of interactions and molecular processes are responsible for the observed anti-correlation. In particular, we observe that the extent of correlation decreases significantly in the order of charged, polar and non-polar amino acid residues. Charged residues that form salt bridge interactions show a dominant anti-correlation. Surprisingly, a few non-polar (hydrophobic and buried) residues demonstrate significant anti-correlation as well. We attribute this to the structural fluctuations that lead to an equilibrium between buried and solvent exposed states of these residues, where the protein-protein and protein-water interactions compensate each other. Our studies indicate that the molecular origin of the observed anti-correlation can be vary depending on the type of the interaction. The overall coupling/correlation is dominated by the electrostatic interactions and charged residues due to a general dielectric screening effect, whereas a few buried hydrophobic residues may also show such coupling due to structural dynamics between buried and exposed states.

Comparative Analysis of Protein Surface Hydrophobicity Maps Determined by Sparse Sampling INDUS and Spatial Aggregation Propensity.

Protein surface hydrophobicity plays a central role in various biological processes such as protein folding and aggregation, as well as in the design and manufacturing of biotherapeutics. While the hydrophobicity of protein surface patches has been linked to their constituent residue hydropathies, recent research has shown that protein surface hydrophobicity is more complex and characterized by the response of water to these surfaces. In this work, we employ water density perturbations to map the surface hydrophobicity of a set of model proteins using sparse indirect umbrella sampling simulations (SSI). This technique is used to identify hydrophobic surface patches for the set of model proteins, and the results are compared to those obtained from the widely adopted spatial aggregation propensity (SAP) technique. While SAP-based calculations show agreement with SSI in some cases, there are several examples of disagreement. We identify four general classes of difference in behavior and study factors that contribute to these differences. We find that the SAP method can sometimes mask the effect of weakly nonpolar or isolated nonpolar residues that can lead to strong hydrophobic patches on the protein surface. In addition, hydrophobic patches identified by SAP can exhibit shifts in both position and strength on the SSI map. Our results demonstrate that the combination of topography and chemical context controls the hydrophobicity of a given patch above and beyond the intrinsic polarity of the residues present on the patch surface. The availability of more accurate protein hydrophobicity maps in concert with new classes of hydrophobic molecular descriptors may create significant opportunities for in silico prediction of protein behavior for a range of applications, such as protein design, biomanufacturability, and downstream bioprocessing.

Structural analysis of open reading frames of bovine immunodeficiency virus proteins

The goal of this study was determining the structural organization peculiarities of the ORF2 and ORF3 proteins of the bovine immunodeficiency virus (BIV). Five ORFs were determined for two BIV isolates with complete genome using the ATGpr software, which permits effective prediction of translation initiation codons with nucleotide accuracy. Phyre2 software was used to predict, analyze the secondary structure and function of proteins. PONDR-FIT software was used to search for protein fragments in a disordered or natively unfolded state. Analysis of the amino acid composition of ORF2 and ORF3 proteins of BIVisolates regarding the presence of nonpolar, polar, aromatic, and hydrophobic amino acid residues was carried out using PSIPRED software. Models of the 3D-structure of proteins were obtained by I-TASSER server. 14% of α helices, 17% of β strands and 43% of disordered structure are predicted for the ORF3 protein. 37% of α helices, 0% of β strands, and 41% of disordered structure were determined for Gag polyprotein, which is translated from ORF2. The distribution of charged amino acid residues characterizes the surface properties of proteins. Their number reaches 23.9% for ORF2 protein. The amount of Arg is 5.2%, Lys — 8.0%, Glu — 7.3%, Asp — 3.4%. The total number of charged amino acid residues of ORF3 is 23.3%. The number of Arg is 12.6%, Lys — 4.9%, Glu — 1.9%, Asp — 3.9%. Only two ORFs of five ones coincide in nucleotide length (and, therefore, in length of corresponding proteins) for the two BIV isolates. The ORF3 protein belongs to the intrinsically disordered proteins that cannot be stably folded into a unique three-dimensional structure under physiological conditions, and the Gag polyprotein, which is translated from ORF2, belongs to the class of fully structured proteins. The secondary structure of both proteins shows the presence of α-helices

Catalysis by KDM6 Histone Demethylases - A Synergy between the Non-Heme Iron(II) Center, Second Coordination Sphere, and Long-Range Interactions.

KDM6A (UTX) and KDM6B (JMJD3) are human non-heme Fe(II) and 2-oxoglutarate (2OG) dependent JmjC oxygenases that catalyze the demethylation of trimethylated lysine 27 in the N-terminal tail of histone H3, a post-translational modification that regulates transcription. A Combined Quantum Mechanics/ Molecular Mechanics (QM/MM) and Molecular Dynamics (MD) study on the catalytic mechanism of KDM6A/B reveals that the transition state for the rate-limiting hydrogen atom transfer (HAT) reaction in KDM6A catalysis is stabilized by polar (Asn217) and aromatic (Trp369)/non-polar (Pro274) residues in contrast to KDM4, KDM6B and KDM7 demethylases where charged residues (Glu, Arg, Asp) are involved. KDM6A employs both σ- and π-electron transfer pathways for HAT, whereas KDM6B employs the σ-electron pathway. Differences in hydrogen bonding of the Fe-chelating Glu252(KDM6B) contribute to the lower energy barriers in KDM6B vs. KDM6A. The study reveals a dependence of the activation barrier of the rebound hydroxylation on the Fe-O-C angle in the transition state of KDM6A. Anti-correlation of the Zn-binding domain with the active site residues is a key factor distinguishing KDM6A/B from KDM7/4s. The results reveal the importance of communication between the Fe center, second coordination sphere, and long-range interactions in catalysis by KDMs and, by implication, other 2OG oxygenases.