Difference contact maps: From what to why in the analysis of the conformational flexibility of proteins

ProteinReDiff: Complex-based ligand-binding proteins redesign by equivariant diffusion-based generative models

Proteins, serving as the fundamental architects of biological processes, interact with ligands to perform a myriad of functions essential for life. Designing functional ligand-binding proteins is pivotal for advancing drug development and enhancing therapeutic efficacy. In this study, we introduce ProteinReDiff, an diffusion framework targeting the redesign of ligand-binding proteins. Using equivariant diffusion-based generative models, ProteinReDiff enables the creation of high-affinity ligand-binding proteins without the need for detailed structural information, leveraging instead the potential of initial protein sequences and ligand SMILES strings. Our evaluations across sequence diversity, structural preservation, and ligand binding affinity underscore ProteinReDiff's potential to advance computational drug discovery and protein engineering.

Structural and functional effects of the L84S mutant in the SARS-COV-2 ORF8 dimer based on microsecond molecular dynamics study

The L84S mutation has been observed frequently in the ORF8 protein of SARS-CoV-2, which is an accessory protein involved in various important functions such as virus propagation, pathogenesis, and evading the immune response. However, the specific effects of this mutation on the dimeric structure of ORF8 and its impacts on interactions with host components and immune responses are not well understood. In this study, we performed one microsecond molecular dynamics (MD) simulation and analyzed the dimeric behavior of the L84S and L84A mutants in comparison to the native protein. The MD simulations revealed that both mutations caused changes in the conformation of the ORF8 dimer, influenced protein folding mechanisms, and affected the overall structural stability. In particular, the 73YIDI76 motif has found to be significantly affected by the L84S mutation, leading to structural flexibility in the region connecting the C-terminal β4 and β5 strands. This flexibility might be responsible for virus immune modulation.  The free energy landscape (FEL) and principle component analysis (PCA) have also supported our investigation. Overall, the L84S and L84A mutations affect the ORF8 dimeric interfaces by reducing the frequency of protein-protein interacting residues (Arg52, Lys53, Arg98, Ile104, Arg115, Val117, Asp119, Phe120, and Ile121) in the ORF8 dimer.  Our findings provide detail insights for further research in designing structure-based therapeutics against the SARS-CoV-2.Communicated by Ramaswamy H. Sarma.

In vitro and in silico investigation of effects of antimicrobial peptides from Solanaceae plants against rice sheath blight pathogen Rhizoctinia solani

Rhizoctonia solani, the causative agent of sheath blight disease in rice, poses a significant threat to agricultural productivity. Traditional management approaches involving chemical fungicides have been effective but come with detrimental consequences for the ecosystem. This study aimed to investigate sustainable alternatives in the form of antifungal peptides derived from Solanaceous plant species as potential agents against R. solani. Peptide extracts were obtained using an optimized antimicrobial peptide (AMP) extraction method and desalted using the solid-phase extraction technique. The antifungal potential of peptide-rich extracts from Solanum tuberosum and Capsicum annum was assessed through in vitro tests employing the agar well diffusion method. Furthermore, peptide-protein docking analysis was performed on HPEPDOCK and HDOCK server; and molecular dynamics simulations (MDS) of 100 ns period were performed using the Gromacs 2020.4. The results demonstrated significant inhibition zones for both extracts at concentrations of 100 mg/mL. Additionally, the extracts of Solanum tuberosum and Capsicum annum had minimum inhibitory concentrations of 50 mg/mL and 25 mg/mL, respectively with minimum fungicidal concentrations of 25 mg/mL. Insights into the potential mechanisms of key peptides inhibiting R. solani targets were gleaned from in-silico studies. Notably, certain AMPs exhibited favorable free energy of binding against pathogenicity-related targets, including histone demethylase, sortin nexin, and squalene synthase, in protein-peptide docking simulations. Extended molecular dynamics simulations lasting 100 ns and MM-PBSA calculations were performed on select protein-peptide complexes. AMP10 displayed the most favorable binding free energy against all target proteins, with AMP3, AMP12b, AMP6, and AMP15 also exhibiting promising results against specific targets of R. solani. These findings underscore the potential of peptide extracts from S. tuberosum and C. annum as effective antifungal agents against rice sheath blight caused by R. solani.

Development of multi epitope subunit vaccines against emerging carp viruses Cyprinid herpesvirus 1 and 3 using immunoinformatics approach

Cyprinid herpesvirus is a causative agent of a destructive disease in common and koi carp (Cyprinus carpio), which leads to substantial global financial losses in aquaculture industries. Among the strains of C. herpesvirus, C. herpesvirus 1 (CyHV-1) and C. herpesvirus 3 (CyHV-3) are known as highly pathogenic to carp fishes in Europe, Asia, and Africa. To date, no effective vaccine has been developed to combat these viruses. This study aimed to develop unique multi-epitope subunit vaccines targeting the CyHV-1 and CyHV-3 using a reverse vaccinology approach. The study began with a comprehensive literature review to identify the most critical proteins, which were then subjected to in silico analyses to predict highly antigenic epitopes. These analyses involved assessing antigenicity, transmembrane topology screening, allergenecity, toxicity, and molecular docking approaches. We constructed two multi-epitope-based vaccines incorporating a suitable adjuvant and appropriate linkers. It revealed that both the vaccines are non-toxic and immunogenic. The tertiary structures of the vaccine proteins were generated, refined, and validated to ensure their suitability. The binding affinity between the vaccine constructs and TLR3 and TLR5 receptors were assessed by molecular docking studies. Molecular dynamics simulations indicated that vaccine construct V1 exhibited greater stability with both TLR3 and TLR5 based on RMSD analysis. Hydrogen bond analysis revealed a stronger binding affinity between the vaccine constructs and TLR5 compared to TLR3. Furthermore, MM-PBSA analysis suggested that both vaccine constructs exhibited a better affinity for TLR5. Considering all aspects, the results suggest that in silico development of CyHV vaccines incorporating multiple epitopes holds promise for management of diseases caused by CyHV-1 and CyHV-3. However, further in vivo trials are highly recommended to validate the efficacies of these vaccines.

Antifungal plant flavonoids identified in silico with potential to control rice blast disease caused by Magnaporthe oryzae

Rice blast disease, caused by the fungus Magnaporthe oryzae, poses a severe threat to rice production, particularly in Asia where rice is a staple food. Concerns over fungicide resistance and environmental impact have sparked interest in exploring natural fungicides as potential alternatives. This study aimed to identify highly potent natural fungicides against M. oryzae to combat rice blast disease, using advanced molecular dynamics techniques. Four key proteins (CATALASE PEROXIDASES 2, HYBRID PKS-NRPS SYNTHETASE TAS1, MANGANESE LIPOXYGENASE, and PRE-MRNA-SPLICING FACTOR CEF1) involved in M. oryzae's infection process were identified. A list of 30 plant metabolites with documented antifungal properties was compiled for evaluation as potential fungicides. Molecular docking studies revealed that 2-Coumaroylquinic acid, Myricetin, Rosmarinic Acid, and Quercetin exhibited superior binding affinities compared to reference fungicides (Azoxystrobin and Tricyclazole). High throughput molecular dynamics simulations were performed, analyzing parameters like RMSD, RMSF, Rg, SASA, hydrogen bonds, contact analysis, Gibbs free energy, and cluster analysis. The results revealed stable interactions between the selected metabolites and the target proteins, involving important hydrogen bonds and contacts. The SwissADME server analysis indicated that the metabolites possess fungicide properties, making them effective and safe fungicides with low toxicity to the environment and living beings. Additionally, bioactivity assays confirmed their biological activity as nuclear receptor ligands and enzyme inhibitors. Overall, this study offers valuable insights into potential natural fungicides for combating rice blast disease, with 2-Coumaroylquinic acid, Myricetin, Rosmarinic Acid, and Quercetin standing out as promising and environmentally friendly alternatives to conventional fungicides. These findings have significant implications for developing crop protection strategies and enhancing global food security, particularly in rice-dependent regions.

Exploring the Alternative Conformation of a Known Protein Structure Based on Contact Map Prediction

The rapid development of deep learning-based methods has considerably advanced the field of protein structure prediction. The accuracy of predicting the 3D structures of simple proteins is comparable to that of experimentally determined structures, providing broad possibilities for structure-based biological studies. Another critical question is whether and how multistate structures can be predicted from a given protein sequence. In this study, analysis of tens of two-state proteins demonstrated that deep learning-based contact map predictions contain structural information on both states, which suggests that it is probably appropriate to change the target of deep learning-based protein structure prediction from one specific structure to multiple likely structures. Furthermore, by combining deep learning- and physics-based computational methods, we developed a protocol for exploring alternative conformations from a known structure of a given protein, by which we successfully approached the holo-state conformations of multiple representative proteins from their apo-state structures.

An immunoinformatics and extended molecular dynamics approach for designing a polyvalent vaccine against multiple strains of Human T-lymphotropic virus (HTLV)

Human T-lymphotropic virus (HTLV), a group of retroviruses belonging to the oncovirus family, has long been associated with various inflammatory and immunosuppressive disorders. At present, there is no approved vaccine capable of effectively combating all the highly pathogenic strains of HTLV that makes this group of viruses a potential threat to human health. To combat the devastating impact of any potential future outbreak caused by this virus group, our study employed a reverse vaccinology approach to design a novel polyvalent vaccine targeting the highly virulent subtypes of HTLV. Moreover, we comprehensively analyzed the molecular interactions between the designed vaccine and corresponding Toll-like receptors (TLRs), providing valuable insights for future research on preventing and managing HTLV-related diseases and any possible outbreaks. The vaccine was designed by focusing on the envelope glycoprotein gp62, a crucial protein involved in the infectious process and immune mechanisms of HTLV inside the human body. Epitope mapping identified T cell and B cell epitopes with low binding energies, ensuring their immunogenicity and safety. Linkers and adjuvants were incorporated to enhance the vaccine's stability, antigenicity, and immunogenicity. Initially, two vaccine constructs were formulated, and among them, vaccine construct-2 exhibited superior solubility and structural stability. Molecular docking analyses also revealed strong binding affinity between the vaccine construct-2 and both targeted TLR2 and TLR4. Molecular dynamics simulations demonstrated enhanced stability, compactness, and consistent hydrogen bonding within TLR-vaccine complexes, suggesting a strong binding affinity. The stability of the complexes was further corroborated by contact, free energy, structure, and MM-PBSA analyses. Consequently, our research proposes a vaccine targeting multiple HTLV subtypes, offering valuable insights into the molecular interactions between the vaccine and TLRs. These findings should contribute to developing effective preventive and treatment approaches against HTLV-related diseases and preventing possible outbreaks. However, future research should focus on in-depth validation through experimental studies to confirm the interactions identified in silico and to evaluate the vaccine's efficacy in relevant animal models and, eventually, in clinical trials.

Molecular Dynamic Studies of Dye-Dye and Dye-DNA Interactions Governing Excitonic Coupling in Squaraine Aggregates Templated by DNA Holliday Junctions

Dye molecules, arranged in an aggregate, can display excitonic delocalization. The use of DNA scaffolding to control aggregate configurations and delocalization is of research interest. Here, we applied Molecular Dynamics (MD) to gain an insight on how dye-DNA interactions affect excitonic coupling between two squaraine (SQ) dyes covalently attached to a DNA Holliday junction (HJ). We studied two types of dimer configurations, i.e., adjacent and transverse, which differed in points of dye covalent attachments to DNA. Three structurally different SQ dyes with similar hydrophobicity were chosen to investigate the sensitivity of excitonic coupling to dye placement. Each dimer configuration was initialized in parallel and antiparallel arrangements in the DNA HJ. The MD results, validated by experimental measurements, suggested that the adjacent dimer promotes stronger excitonic coupling and less dye-DNA interaction than the transverse dimer. Additionally, we found that SQ dyes with specific functional groups (i.e., substituents) facilitate a closer degree of aggregate packing via hydrophobic effects, leading to a stronger excitonic coupling. This work advances a fundamental understanding of the impacts of dye-DNA interactions on aggregate orientation and excitonic coupling.

Immunoinformatics Approach to Design Novel Subunit Vaccine against the Epstein-Barr Virus

Epstein-Barr virus (EBV) is a lymphotropic virus responsible for numerous epithelial and lymphoid cell malignancies, including gastric carcinoma, Hodgkin's lymphoma, nasopharyngeal carcinoma, and Burkitt lymphoma. Hundreds of thousands of people worldwide get infected with this virus, and in most cases, this viral infection leads to cancer. Although researchers are trying to develop potential vaccines and drug therapeutics, there is still no effective vaccine to combat this virus. In this study, the immunoinformatics approach was utilized to develop a potential multiepitope subunit vaccine against the two most common subtypes of EBV, targeting three of their virulent envelope glycoproteins. Eleven cytotoxic T lymphocyte (CTL) epitopes, 11 helper T lymphocyte (HTL) epitopes, and 10 B-cell lymphocyte (BCL) epitopes were predicted to be antigenic, nonallergenic, nontoxic, and fully conserved among the two subtypes, and nonhuman homologs were used for constructing the vaccine after much analysis. Later, further validation experiments, including molecular docking with different immune receptors (e.g., Toll-like receptors [TLRs]), molecular dynamics simulation analyses (including root means square deviation [RMSD], root mean square fluctuation [RMSF], radius of gyration [Rg], principal-component analysis [PCA], dynamic cross-correlation [DCC], definition of the secondary structure of proteins [DSSP], and Molecular Mechanics Poisson-Boltzmann Surface Area [MM-PBSA]), and immune simulation analyses generated promising results, ensuring the safe and stable response of the vaccine with specific immune receptors after potential administration within the human body. The vaccine's high binding affinity with TLRs was revealed in the docking study, and a very stable interaction throughout the simulation proved the potential high efficacy of the proposed vaccine. Further, in silico cloning was also conducted to design an efficient mass production strategy for future bulk industrial vaccine production. IMPORTANCE Epstein-Barr virus (EBV) vaccines have been developing for over 30 years, but polyphyletic and therapeutic vaccines have failed to get licensed. Our vaccine surpasses the limitations of many such vaccines and remains very promising, which is crucial because the infection rate is higher than most viral infections, affecting a whopping 90% of the adult population. One of the major identifications covers a holistic analysis of populations worldwide, giving us crucial information about its effectiveness for everyone's unique immunological system. We targeted three glycoproteins that enhance the virulence of the virus to design an epitope-based polyvalent vaccine against two different strains of EBV, type 1 and 2. Our methodology in this study is nonconventional yet swift to show effective results while designing vaccines.

Extending the Horizon of Homology Detection with Coevolution-based Structure Prediction

Traditional sequence analysis algorithms fail to identify distant homologies when they lie beyond a detection horizon. In this review, we discuss how co-evolution-based contact and distance prediction methods are pushing back this homology detection horizon, thereby yielding new functional insights and experimentally testable hypotheses. Based on correlated substitutions, these methods divine three-dimensional constraints among amino acids in protein sequences that were previously devoid of all annotated domains and repeats. The new algorithms discern hidden structure in an otherwise featureless sequence landscape. Their revelatory impact promises to be as profound as the use, by archaeologists, of ground-penetrating radar to discern long-hidden, subterranean structures. As examples of this, we describe how triplicated structures reflecting longin domains in MON1A-like proteins, or UVR-like repeats in DISC1, emerge from their predicted contact and distance maps. These methods also help to resolve structures that do not conform to a "beads-on-a-string" model of protein domains. In one such example, we describe CFAP298 whose ubiquitin-like domain was previously challenging to perceive owing to a large sequence insertion within it. More generally, the new algorithms permit an easier appreciation of domain families and folds whose evolution involved structural insertion or rearrangement. As we exemplify with α1-antitrypsin, coevolution-based predicted contacts may also yield insights into protein dynamics and conformational change. This new combination of structure prediction (using innovative co-evolution based methods) and homology inference (using more traditional sequence analysis approaches) shows great promise for bringing into view a sea of evolutionary relationships that had hitherto lain far beyond the horizon of homology detection.

Dynamics-Evolution Correspondence in Protein Structures

The genotype-phenotype mapping of proteins is a fundamental question in structural biology. In this Letter, with the analysis of a large dataset of proteins from hundreds of protein families, we quantitatively demonstrate the correlations between the noise-induced protein dynamics and mutation-induced variations of native structures, indicating the dynamics-evolution correspondence of proteins. Based on the investigations of the linear responses of native proteins, the origin of such a correspondence is elucidated. It is essential that the noise- and mutation-induced deformations of the proteins are restricted on a common low-dimensional subspace, as confirmed from the data. These results suggest an evolutionary mechanism of the proteins gaining both dynamical flexibility and evolutionary structural variability.