Interface residues of transient protein-protein complexes have extensive intra-protein interactions apart from inter-protein interactions

Interfacial residues in protein-protein complexes are in the eyes of the beholder

Interactions between proteins are vital in almost all biological processes. The characterization of protein-protein interactions helps us understand the mechanistic basis of biological processes, thereby enabling the manipulation of proteins for biotechnological and clinical purposes. The interface residues of a protein-protein complex are assumed to have the following two properties: (a) they always interact with a residue of a partner protein, which forms the basis for distance-based interface residue identification methods, and (b) they are solvent-exposed in the isolated form of the protein and become buried in the complex form, which forms the basis for Accessible Surface Area (ASA)-based methods. The study interrogates this popular assumption by recognizing interface residues in protein-protein complexes through these two methods. The results show that a few residues are identified uniquely by each method, and the extent of conservation, propensities, and their contribution to the stability of protein-protein interaction varies substantially between these residues. The case study analyses showed that interface residues, unique to distance, participate in crucial interactions that hold the proteins together, whereas the interface residues unique to the ASA method have a potential role in the recognition, dynamics, and specificity of the complex and can also be a hotspot. Overall, the study recommends applying both distance and ASA methods so that some interface residues missed by either method but crucial to the stability, recognition, dynamics, and function of protein-protein complexes are identified in a complementary manner.

Predictive Models and Impact of Interfacial Contacts and Amino Acids on Protein-Protein Binding Affinity

Protein-protein interactions (PPIs) play a central role in nearly all cellular processes. The strength of the binding in a PPI is characterized by the binding affinity (BA) and is a key factor in controlling protein-protein complex formation and defining the structure-function relationship. Despite advancements in understanding protein-protein binding, much remains unknown about the interfacial region and its association with BA. New models are needed to predict BA with improved accuracy for therapeutic design. Here, we use machine learning approaches to examine how well different types of interfacial contacts can be used to predict experimentally determined BA and to reveal the impact of the specific amino acids at the binding interface on BA. We create a series of multivariate linear regression models incorporating different contact features at both residue and atomic levels and examine how different methods of identifying and characterizing these properties impact the performance of these models. Particularly, we introduce a new and simple approach to predict BA based on the quantities of specific amino acids at the protein-protein interface. We found that the numbers of specific amino acids at the protein-protein interface were correlated with BA. We show that the interfacial numbers of amino acids can be used to produce models with consistently good performance across different data sets, indicating the importance of the identities of interfacial amino acids in underlying BA. When trained on a diverse set of complexes from two benchmark data sets, the best performing BA model was generated with an explicit linear equation involving six amino acids. Tyrosine, in particular, was identified as the key amino acid in controlling BA, as it had the strongest correlation with BA and was consistently identified as the most important amino acid in feature importance studies. Glycine and serine were identified as the next two most important amino acids in predicting BA. The results from this study further our understanding of PPIs and can be used to make improved predictions of BA, giving them implications for drug design and screening in the pharmaceutical industry.

A deep learning-based framework for predicting survival-associated groups in colon cancer by integrating multi-omics and clinical data

Precise prognostic classification of patients and identifying survival subgroups and their associated genes can be important clinical references when designing treatment strategies for cancer patients. Multi-omics and data integration techniques are powerful tools to achieve this goal. This study aimed to introduce a machine learning method to integrate three types of biological data, and investigate the performance of two other methods, in identifying the survival dependency of patients. The data included TCGA RNA-seq gene expression, DNA methylation, and clinical data from 368 patients with colon cancer also we use an independent external validation data set, containing 232 samples. Three methods including, hyper-parameter optimized autoencoders (HPOAE), normal autoencoder, and penalized principal component analysis (PPCA) were used for simultaneous data integration and estimation under a COX hazards model. The HPOAE was thought to outperform other methods. The HPOAE had the Log Rank Mantel-Cox value of 14.27 ± 2, and a Breslow-Generalized Wilcoxon value of 13.13 ± 1. Ten miRNA, 11 methylated genes, and 28 mRNA all by (importance of marginal cutoff > 0.95) were identified. The study demonstrated that hsa-miR-485-5p targets both ZMYM1 and tp53, the latter of which has been previously associated with cancer in numerous studies. Furthermore, compared to other methods, the HPOAE exhibited a greater capacity for identifying survival subgroups and the genes associated with them in patients with colon cancer. However, all of the results were obtained by computational methods, and clinical and experimental studies are needed to validate these results.

Prediction of transient and permanent protein interactions using AI methods

Protein-protein interactions (PPIs) can be classified as permanent or transient interactions based on their stability or lifetime. Understanding the precise details of such protein interactions will pave the way for the discovery of inhibitors and for understanding the nature and function of PPIs. In the present work, 43 relevant physicochemical, geometrical and structural features were calculated for a curated dataset from the literature, comprising of 402 protein-protein complexes of permanent and transient categories, and 5 different Supervised Machine Learning models were developed with Scikit-learn to predict transient and permanent PPI. Additionally, deep learning method with Artificial Neural Network was also performed using Tensor Flow and Keras. Predicted models achieved accuracy ranging from 76.54% to 82.71% and k-NN has achieved the highest accuracy. Detailed analysis of these methods revealed that Interface areas such as Percent interface accessible area, Interface accessible area and Total interface area and the parameters defining the shape of the PPI interface such as Planarity, Eccentricity and Circularity are the most discriminating factors between these two categories. The present method could serve as an effective tool to understand the mechanism of protein association and to predict the transient and permanent interactions, which could supplement the costly and time-consuming experimental techniques.

Protein-protein interactions in plant antioxidant defense

The regulation of reactive oxygen species (ROS) levels in plants is ensured by mechanisms preventing their over accumulation, and by diverse antioxidants, including enzymes and nonenzymatic compounds. These are affected by redox conditions, posttranslational modifications, transcriptional and posttranscriptional modifications, Ca²⁺, nitric oxide (NO) and mitogen-activated protein kinase signaling pathways. Recent knowledge about protein-protein interactions (PPIs) of antioxidant enzymes advanced during last decade. The best-known examples are interactions mediated by redox buffering proteins such as thioredoxins and glutaredoxins. This review summarizes interactions of major antioxidant enzymes with regulatory and signaling proteins and their diverse functions. Such interactions are important for stability, degradation and activation of interacting partners. Moreover, PPIs of antioxidant enzymes may connect diverse metabolic processes with ROS scavenging. Proteins like receptor for activated C kinase 1 may ensure coordination of antioxidant enzymes to ensure efficient ROS regulation. Nevertheless, PPIs in antioxidant defense are understudied, and intensive research is required to define their role in complex regulation of ROS scavenging.

Identification of potential candidate vaccines against Mycobacterium ulcerans based on the major facilitator superfamily transporter protein

Buruli ulcer is a neglected tropical disease that is characterized by non-fatal lesion development. The causative agent is Mycobacterium ulcerans (M. ulcerans). There are no known vectors or transmission methods, preventing the development of control methods. There are effective diagnostic techniques and treatment routines; however, several socioeconomic factors may limit patients' abilities to receive these treatments. The Bacillus Calmette-Guérin vaccine developed against tuberculosis has shown limited efficacy, and no conventionally designed vaccines have passed clinical trials. This study aimed to generate a multi-epitope vaccine against M. ulcerans from the major facilitator superfamily transporter protein using an immunoinformatics approach. Twelve M. ulcerans genome assemblies were analyzed, resulting in the identification of 11 CD8+ and 7 CD4+ T-cell epitopes and 2 B-cell epitopes. These conserved epitopes were computationally predicted to be antigenic, immunogenic, non-allergenic, and non-toxic. The CD4+ T-cell epitopes were capable of inducing interferon-gamma and interleukin-4. They successfully bound to their respective human leukocyte antigens alleles in in silico docking studies. The expected global population coverage of the T-cell epitopes and their restricted human leukocyte antigens alleles was 99.90%. The population coverage of endemic regions ranged from 99.99% (Papua New Guinea) to 21.81% (Liberia). Two vaccine constructs were generated using the Toll-like receptors 2 and 4 agonists, LprG and RpfE, respectively. Both constructs were antigenic, non-allergenic, non-toxic, thermostable, basic, and hydrophilic. The DNA sequences of the vaccine constructs underwent optimization and were successfully in-silico cloned with the pET-28a(+) plasmid. The vaccine constructs were successfully docked to their respective toll-like receptors. Molecular dynamics simulations were carried out to analyze the binding interactions within the complex. The generated binding energies indicate the stability of both complexes. The constructs generated in this study display severable favorable properties, with construct one displaying a greater range of favorable properties. However, further analysis and laboratory validation are required.

Biography of a scientist with strength, substance, sincerity and service: Late N. Srinivasan (1962-2021)

Late N. Srinivasan belongs to the GN Ramachandran lineage of protein structural analysts. His role in the advancement of the structure based understanding of signal transduction, protein kinase analyses and host-pathogen interactions both developing and using Bioinformatics tools for protein-protein interactions, protein dynamics, remote homology detection and polypeptide stereochemistry is well documented in the literature. Thus, his contribution to the understanding of protein function through structural analysis, using computational models and tools, is exceptional.

Cryo-EM structure of a type IV secretion system

Bacterial conjugation is the fundamental process of unidirectional transfer of DNAs, often plasmid DNAs, from a donor cell to a recipient cell¹. It is the primary means by which antibiotic resistance genes spread among bacterial populations^2,3. In Gram-negative bacteria, conjugation is mediated by a large transport apparatus—the conjugative type IV secretion system (T4SS)—produced by the donor cell and embedded in both its outer and inner membranes. The T4SS also elaborates a long extracellular filament—the conjugative pilus—that is essential for DNA transfer^4,5. Here we present a high-resolution cryo-electron microscopy (cryo-EM) structure of a 2.8 megadalton T4SS complex composed of 92 polypeptides representing 8 of the 10 essential T4SS components involved in pilus biogenesis. We added the two remaining components to the structural model using co-evolution analysis of protein interfaces, to enable the reconstitution of the entire system including the pilus. This structure describes the exceptionally large protein–protein interaction network required to assemble the many components that constitute a T4SS and provides insights on the unique mechanism by which they elaborate pili.

Cryo-electron microscopy structures of a 2.8 megadalton bacterial type IV secretion system encoded by the plasmid R388 and comprising 92 polypeptides provide insights into the stepwise mechanism of pilus assembly.

Fever as an evolutionary agent to select immune complexes interfaces

We herein analyzed all available protein–protein interfaces of the immune complexes from the Protein Data Bank whose antigens belong to pathogens or cancers that are modulated by fever in mammalian hosts. We also included, for comparison, protein interfaces from immune complexes that are not significantly modulated by the fever response. We highlight the distribution of amino acids at these viral, bacterial, protozoan and cancer epitopes, and at their corresponding paratopes that belong strictly to monoclonal antibodies. We identify the “hotspots”, i.e. residues that are highly connected at such interfaces, and assess the structural, kinetic and thermodynamic parameters responsible for complex formation. We argue for an evolutionary pressure for the types of residues at these protein interfaces that may explain the role of fever as a selective force for optimizing antibody binding to antigens.

Deciphering Isoniazid Drug Resistance Mechanisms on Dimeric Mycobacterium tuberculosis KatG via Post-molecular Dynamics Analyses Including Combined Dynamic Residue Network Metrics

Resistance mutations in Mycobacterium tuberculosis (Mtb) catalase peroxidase protein (KatG), an essential enzyme in isoniazid (INH) activation, reduce the sensitivity of Mtb to first-line drugs, hence presenting challenges in tuberculosis (TB) management. Thus, understanding the mutational imposed resistance mechanisms remains of utmost importance in the quest to reduce the TB burden. Herein, effects of 11 high confidence mutations in the KatG structure and residue network communication patterns were determined using extensive computational approaches. Combined traditional post-molecular dynamics analysis and comparative essential dynamics revealed that the mutant proteins have significant loop flexibility around the heme binding pocket and enhanced asymmetric protomer behavior with respect to wild-type (WT) protein. Heme contact analysis between WT and mutant proteins identified a reduction to no contact between heme and residue His270, a covalent bond vital for the heme-enabled KatG catalytic activity. Betweenness centrality calculations showed large hub ensembles with new hubs especially around the binding cavity and expanded to the dimerization domain via interface in the mutant systems, providing possible compensatory allosteric communication paths for the active site as a result of the mutations which may destabilize the heme binding pocket and the loops in its vicinity. Additionally, an interesting observation came from Eigencentrality hubs, most of which are located in the C-terminal domain, indicating relevance of the domain in the protease functionality. Overall, our results provide insight toward the mechanisms involved in KatG-INH resistance in addition to identifying key regions in the enzyme functionality, which can be used for future drug design.

Emerging roles of the HECT-type E3 ubiquitin ligases in hematological malignancies

Ubiquitination-mediated proteolysis or regulation of proteins, ultimately executed by E3 ubiquitin ligases, control a wide array of cellular processes, including transcription, cell cycle, autophagy and apoptotic cell death. HECT-type E3 ubiquitin ligases can be distinguished from other subfamilies of E3 ubiquitin ligases because they have a C-terminal HECT domain that directly catalyzes the covalent attachment of ubiquitin to their substrate proteins. Deregulation of HECT-type E3-mediated ubiquitination plays a prominent role in cancer development and chemoresistance. Several members of this subfamily are indeed frequently deregulated in human cancers as a result of genetic mutations and altered expression or activity. HECT-type E3s contribute to tumorigenesis by regulating the ubiquitination rate of substrates that function as either tumour suppressors or oncogenes. While the pathological roles of the HECT family members in solid tumors are quite well established, their contribution to the pathogenesis of hematological malignancies has only recently emerged. This review aims to provide a comprehensive overview of the involvement of the HECT-type E3s in leukemogenesis.

Stable and Temporary Enzyme Complexes and Metabolons Involved in Energy and Redox Metabolism

Significance: Alongside well-characterized permanent multimeric enzymes and multienzyme complexes, relatively unstable transient enzyme-enzyme assemblies, including metabolons, provide an important mechanism for the regulation of energy and redox metabolism. Critical Issues: Despite the fact that enzyme-enzyme assemblies have been proposed for many decades and experimentally analyzed for at least 40 years, there are very few pathways for which unequivocal evidence for the presence of metabolite channeling, the most frequently evoked reason for their formation, has been provided. Further, in contrast to the stronger, permanent interactions for which a deep understanding of the subunit interface exists, the mechanism(s) underlying transient enzyme-enzyme interactions remain poorly studied. Recent Advances: The widespread adoption of proteomic and cell biological approaches to characterize protein-protein interaction is defining an ever-increasing number of enzyme-enzyme assemblies as well as enzyme-protein interactions that likely identify factors which stabilize such complexes. Moreover, the use of microfluidic technologies provided compelling support of a role for substrate-specific chemotaxis in complex assemblies. Future Directions: Embracing current and developing technologies should render the delineation of metabolons from other enzyme-enzyme complexes more facile. In parallel, attempts to confirm that the findings reported in microfluidic systems are, indeed, representative of the cellular situation will be critical to understanding the physiological circumstances requiring and evoking dynamic changes in the levels of the various transient enzyme-enzyme assemblies of the cell. Antioxid. Redox Signal. 35, 788-807.

Fuzzy binding model of molecular interactions between tropoelastin and integrin alphaVbeta3

Tropoelastin is the highly flexible monomer subunit of elastin, required for the resilience of the extracellular matrix in elastic tissues. To elicit biological signaling, multiple sites on tropoelastin bind to cell surface integrins in a poorly understood multifactorial process. We constructed a full atomistic molecular model of the interactions between tropoelastin and integrin α_vβ₃ using ensemble-based computational methodologies. Conformational changes of integrin α_vβ₃ associated with outside-in signaling were more frequently facilitated in an ensemble in which tropoelastin bound the integrin's α1 helix rather than the upstream canonical binding site. Our findings support a model of fuzzy binding, whereby many tropoelastin conformations and defined sites cooperatively interact with multiple α_vβ₃ regions. This model explains prior experimental binding to distinct tropoelastin regions, domains 17 and 36, and points to the cooperative participation of domain 20. Our study highlights the utility of ensemble-based approaches in helping to understand the interactive mechanisms of functionally significant flexible proteins.

The Odd Faces of Oligomers: The Case of TRAF2-C, A Trimeric C-Terminal Domain of TNF Receptor-Associated Factor

TNF Receptor Associated Factor 2 (TRAF2) is a trimeric protein that belongs to the TNF receptor associated factor family (TRAFs). The TRAF2 oligomeric state is crucial for receptor binding and for its interaction with other proteins involved in the TNFR signaling. The monomer-trimer equilibrium of a C- terminal domain truncated form of TRAF2 (TRAF2-C), plays also a relevant role in binding the membrane, causing inward vesiculation. In this study, we have investigated the conformational dynamics of TRAF2-C through circular dichroism, fluorescence, and dynamic light scattering, performing temperature-dependent measurements. The data indicate that the protein retains its oligomeric state and most of its secondary structure, while displaying a significative increase in the heterogeneity of the tyrosines signal, increasing the temperature from ≈15 to ≈35 °C. The peculiar crowding of tyrosine residues (12 out of 18) at the three subunit interfaces and the strong dependence on the trimer concentration indicate that such conformational changes mainly involve the contact areas between each pair of monomers, affecting the oligomeric state. Molecular dynamic simulations in this temperature range suggest that the interfaces heterogeneity is an intrinsic property of the trimer that arises from the continuous, asymmetric approaching and distancing of its subunits. Such dynamics affect the results of molecular docking on the external protein surface using receptor peptides, indicating that the TRAF2-receptor interaction in the solution might not involve three subunits at the same time, as suggested by the static analysis obtainable from the crystal structure. These findings shed new light on the role that the TRAF2 oligomeric state might have in regulating the protein binding activity in vivo.

A molecular modelling approach for identifying antiviral selenium-containing heterocyclic compounds that inhibit the main protease of SARS-CoV-2: an in silico investigation

Coronavirus disease 2019 (COVID-19), an infectious disease caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has been declared a global pandemic by the World Health Organization, and the situation worsens daily, associated with acute increases in case fatality rates. The main protease (Mpro) enzyme produced by SARS-CoV-2 was recently demonstrated to be responsible for not only viral reproduction but also impeding host immune responses. The element selenium (Se) plays a vital role in immune functions, both directly and indirectly. Thus, we hypothesised that Se-containing heterocyclic compounds might curb the activity of SARS-CoV-2 Mpro. We performed a molecular docking analysis and found that several of the selected selenocompounds showed potential binding affinities for SARS-CoV-2 Mpro, especially ethaselen (49), which exhibited a docking score of -6.7 kcal/mol compared with the -6.5 kcal/mol score for GC376 (positive control). Drug-likeness calculations suggested that these compounds are biologically active and possess the characteristics of ideal drug candidates. Based on the binding affinity and drug-likeness results, we selected the 16 most effective selenocompounds as potential anti-COVID-19 drug candidates. We also validated the structural integrity and stability of the drug candidate through molecular dynamics simulation. Using further in vitro and in vivo experiments, we believe that the targeted compound identified in this study (ethaselen) could pave the way for the development of prospective drugs to combat SARS-CoV-2 infections and trigger specific host immune responses.

Mechanistic differences in the effects of sucrose and sucralose on the phase stability of lysozyme solutions

The effect of two disaccharide analogues, sucrose and sucralose, on the phase stability of aqueous lysozyme solutions has been addressed from a mechanistic viewpoint by a combination of experiment and molecular dynamics (MD) simulations. The influence of the added low molecular weight salts (NaBr, NaI and NaNO₃) was considered as well. The cloud-point temperature measurements revealed a larger stabilizing effect of sucralose. Upon increasing sugar concentration, the protein solutions became more stable and differences in the effect of sucralose and sucrose amplified. It was confirmed that the addition of either of the two sugars imposed no secondary structure changes of the lysozyme. Enthalpies of lysozyme-sugar mixing were exothermic and a larger effect was recorded for sucralose. MD simulations indicated that acidic, basic and polar amino acid residues play predominant roles in the sugar-protein interactions, mainly through hydrogen bonding. Such sugar mediated protein-protein interactions are thought to be responsible for the biopreserative nature of sugars. Our observations hint at mechanistic differences in sugar-lysozyme interactions: while sucrose does not interact directly with the protein's surface for the most part (in line with the preferential hydration hypothesis), sucralose forms hydrogen bonds with acidic, basic and polar amino acid residues at the lysozyme's surface (in line with the water replacement hypothesis).

Serine 477 plays a crucial role in the interaction of the SARS-CoV-2 spike protein with the human receptor ACE2

Since the worldwide outbreak of the infectious disease COVID-19, several studies have been published to understand the structural mechanism of the novel coronavirus SARS-CoV-2. During the infection process, the SARS-CoV-2 spike (S) protein plays a crucial role in the receptor recognition and cell membrane fusion process by interacting with the human angiotensin-converting enzyme 2 (hACE2) receptor. However, new variants of these spike proteins emerge as the virus passes through the disease reservoir. This poses a major challenge for designing a potent antigen for an effective immune response against the spike protein. Through a normal mode analysis (NMA) we identified the highly flexible region in the receptor binding domain (RBD) of SARS-CoV-2, starting from residue 475 up to residue 485. Structurally, the position S477 shows the highest flexibility among them. At the same time, S477 is hitherto the most frequently exchanged amino acid residue in the RBDs of SARS-CoV-2 mutants. Therefore, using MD simulations, we have investigated the role of S477 and its two frequent mutations (S477G and S477N) at the RBD during the binding to hACE2. We found that the amino acid exchanges S477G and S477N strengthen the binding of the SARS-COV-2 spike with the hACE2 receptor.

Molecular basis for protein-protein interactions

This minireview provides an overview on the current knowledge of protein-protein interactions, common characterisation methods to characterise them, and their role in protein complex formation with some examples. A deep understanding of protein-protein interactions and their molecular interactions is important for a number of applications, including drug design. Protein-protein interactions and their discovery are thus an interesting avenue for understanding how protein complexes, which make up the majority of proteins, work.

Recifin A, Initial Example of the Tyr-Lock Peptide Structural Family, Is a Selective Allosteric Inhibitor of Tyrosyl-DNA Phosphodiesterase I

Tyrosyl-DNA phosphodiesterase 1 (TDP1) is a molecular target for the sensitization of cancer cells to the FDA-approved topoisomerase inhibitors topotecan and irinotecan. High-throughput screening of natural product extract and fraction libraries for inhibitors of TDP1 activity resulted in the discovery of a new class of knotted cyclic peptides from the marine sponge Axinella sp. Bioassay-guided fractionation of the source extract resulted in the isolation of the active component which was determined to be an unprecedented 42-residue cysteine-rich peptide named recifin A. The native NMR structure revealed a novel fold comprising a four strand antiparallel β-sheet and two helical turns stabilized by a complex disulfide bond network that creates an embedded ring around one of the strands. The resulting structure, which we have termed the Tyr-lock peptide family, is stabilized by a tyrosine residue locked into three-dimensional space. Recifin A inhibited the cleavage of phosphodiester bonds by TDP1 in a FRET assay with an IC₅₀ of 190 nM. Enzyme kinetics studies revealed that recifin A can specifically modulate the enzymatic activity of full-length TDP1 while not affecting the activity of a truncated catalytic domain of TDP1 lacking the N-terminal regulatory domain (Δ1-147), suggesting an allosteric binding site for recifin A on the regulatory domain of TDP1. Recifin A represents both the first of a unique structural class of knotted disulfide-rich peptides and defines a previously unseen mechanism of TDP1 inhibition that could be productively exploited for potential anticancer applications.

Can COVID-19 pandemic boost the epidemic of neurodegenerative diseases?

The pandemic of Coronavirus Disease 2019 (COVID-19) presents the world with the medical challenge associated with multifactorial nature of this pathology. Indeed COVID-19 affects several organs and systems and presents diversified clinical picture. COVID-19 affects the brain in many ways including direct infection of neural cells with SARS-CoV-2, severe systemic inflammation which floods the brain with pro-inflammatory agents thus damaging nervous cells, global brain ischaemia linked to a respiratory failure, thromboembolic strokes related to increased intravascular clotting and severe psychological stress. Often the COVID-19 is manifested by neurological and neuropsychiatric symptoms that include dizziness, disturbed sleep, cognitive deficits, delirium, hallucinations and depression. All these indicate the damage to the nervous tissue which may substantially increase the incidence of neurodegenerative diseases and promote dementia.

Interaction of Isocitrate Lyase with Proteins Involved in the Energetic Metabolism in Paracoccidioides lutzii

Background: Systemic mycosis is a cause of death of immunocompromised subjects. The treatment directed to evade fungal pathogens shows severe limitations, such as time of drug exposure and side effects. The paracoccidioidomycosis (PCM) treatment depends on the severity of the infection and may last from months to years. Methods: To analyze the main interactions of Paracoccidioides lutzii isocitrate lyase (ICL) regarding the energetic metabolism through affinity chromatography, we performed blue native PAGE and co-immunoprecipitation to identify ICL interactions. We also performed in silico analysis by homology, docking, hot-spot prediction and contact preference analysis to identify the conformation of ICL complexes. Results: ICL interacted with 18 proteins in mycelium, 19 in mycelium-to-yeast transition, and 70 in yeast cells. Thirty complexes were predicted through docking and contact preference analysis. ICL has seven main regions of interaction with protein partners. Conclusions: ICL seems to interfere with energetic metabolism of P. lutzii, regulating aerobic and anaerobic metabolism as it interacts with proteins from glycolysis, gluconeogenesis, TCA and methylcitrate cycles, mainly through seven hot-spot residues.

Solvent dynamics play a decisive role in the complex formation of biologically relevant redox proteins

Electron transfer processes between proteins are vital in many biological systems. Yet, the role of the solvent in influencing these redox reactions remains largely unknown. In this study, terahertz-time domain spectroscopy (THz-TDS) is used to probe the collective hydration dynamics of flavoenzyme ferredoxin-NADP+-reductase (FNR), electron transfer protein ferredoxin-1 (PetF), and the transient complex that results from their interaction. Results reveal changes in the sub-picosecond hydration dynamics that are dependent upon the surface electrostatic properties of the individual proteins and the transient complex. Retarded solvent dynamics of 8-9 ps are observed for FNR, PetF, and the FNR:PetF transient complex. Binding of the FNR:PetF complex to the substrate NADP+ results in bulk-like solvent dynamics of 7 ps, showing that formation of the ternary complex is entropically favored. Our THz measurements reveal that the electrostatic interaction of the protein surface with water results in charge sensitive changes in the solvent dynamics. Complex formation between the positively charged FNR:NADP+ pre-complex and the negatively charged PetF is not only entropically favored, but in addition the solvent reorganization into more bulk-like water assists the molecular recognition process. The change in hydration dynamics observed here suggests that the interaction with the solvent plays a significant role in mediating electron transfer processes between proteins.

Exploring the Mechanism of Skeletal Muscle in a Tacrolimus-Induced Posttransplantation Diabetes Mellitus Model on Gene Expression Profiles

OBJECTIVE

Posttransplantation diabetes mellitus (PTDM) is a known complication of transplantation that affects the prognosis. Tacrolimus (Tac or FK506) is a widely used immunosuppressant that has been reported to be a risk factor for PTDM and to further induce complications in heart and skeletal muscles, but the mechanism is still largely unknown. In our preliminary experiments, we found that after Tac treatment, blood glucose increased, and the weight of skeletal muscle declined. Here, we hypothesize that tacrolimus can induce PTDM and influence the atrophy of skeletal muscle.

METHODS

We designed preliminary experiments to establish a tacrolimus-induced PTDM model. Gene expression profiles in quadriceps muscle from this rat model were characterized by oligonucleotide microarrays. Then, differences in gene expression profiles in muscle from PTDM rats that received tacrolimus and control subjects were analyzed by using GeneSpring GX 11.0 software (Agilent). Functional annotation and enrichment analysis of differentially expressed genes (DEGs) helped us identify clues for the side effects of tacrolimus.

RESULTS

Our experiments found that the quadriceps in tacrolimus-induced PTDM group were smaller than those in the control group. The study identified 275 DEGs that may be responsible for insulin resistance and the progression of PTDM, including 86 upregulated genes and 199 downregulated genes. GO and KEGG functional analysis of the DEGs showed a significant correlation between PTDM and muscle development. PPI network analysis screened eight hub genes and found that they were related to troponin and tropomyosin.

CONCLUSIONS

This study explored the molecular mechanism of muscle atrophy in a tacrolimus-induced PTDM model by bioinformatics analyses. We identified 275 DEGs and identified significant biomarkers for predicting the development and progression of tacrolimus-induced PTDM.