Revisiting the Ramachandran plot based on statistical analysis of static and dynamic characteristics of protein structures

Special Issue: "Molecular Dynamics Simulations and Structural Analysis of Protein Domains"

The 3D protein structure is the basis for all their biological functions [...].

Avian influenza A H5N1 hemagglutinin protein models have distinct structural patterns re-occurring across the 1959-2023 strains

Influenza A H5N1 hemagglutinin (HA) plays a crucial role in viral pathogenesis and changes in the HA receptor binding domain (RBD) have been attributed to alterations in viral pathogenesis. Mutations often occur within the HA which in-turn results in HA structural changes that consequently contribute to protein evolution. However, the possible occurrence of mutations that results to reversion of the HA protein (going back to an ancestral protein conformation) which in-turn creates distinct HA structural patterns across the 1959-2023 H5N1 viral evolution has never been investigated. Here, we generated and verified the quality of the HA models, identified similar HA structural patterns, and elucidated the possible variations in HA RBD structural dynamics. Our results show that there are 7 distinct structural patterns occurring among the 1959-2023 H5N1 HA models which suggests that reversion of the HA protein putatively occurs during viral evolution. Similarly, we found that the HA RBD structural dynamics vary among the 7 distinct structural patterns possibly affecting viral pathogenesis.

MDavocado: Analysis and Visualization of Protein Motion by Time-Dependent Angular Diagrams

Extracting meaningful information from atomistic molecular dynamics (MD) simulations of proteins remains a challenging task due to the high-dimensionality and complexity of the data. MD simulations yield trajectories that contain the positions of thousands of atoms in millions of steps. Gaining a comprehensive understanding of local dynamical events across the entire trajectory is often difficult. Here, we present a novel approach to visualize MD trajectories in the form of time-dependent Ramachandran plots. Specialized data aggregation techniques are employed to address the challenge of plotting millions of data points on a single image, thereby ensuring that the analysis is independent of the molecule size and/or length of the MD simulation. This approach facilitates quick identification of flexible and dynamic regions, and its strength is the ability to simultaneously observe the movements of all amino acids over time. The Python program MDavocado is freely available at GitHub (https://github.com/zoranstefanic/MDavocado).

The Peptide Bond: Resonance Increases Bond Order and Complicates Fragmentation

The enhancement of the peptide bond order by a resonance in the lone pair of N and the π-bond of CO is analyzed. A decomposition of the bond order in terms of localized molecular orbitals is developed and applied to the peptide bond. A combination of two rotations of hybrid orbitals is proposed to improve the boundary treatment in the fragment molecular orbital method. The developed approach is applied to peptide bonds, and it is found crucial to retain the π orbital in the variational space of both fragments across the boundary. The interaction energies between conventional amino acid residues in Trp-cage (1L2Y) are discussed.

Structural exploration of the PfBLM Helicase-ATP Binding Domain and implications in the quest for antimalarial therapies

BACKGROUND OBJECTIVES

The battle against malaria has witnessed remarkable progress in recent years, characterized by increased funding, development of life-saving tools, and a significant reduction in disease prevalence. Yet, the formidable challenge of drug resistance persists, threatening to undo these gains.

METHODS

To tackle this issue, it is imperative to identify new effective drug candidates against the malaria parasite that exhibit minimal toxicity. This study focuses on discovering such candidates by targeting PfRecQ1, also known as PfBLM, a vital protein within the malaria parasite Plasmodium falciparum . PfRecQ1 plays a crucial role in the parasite's life cycle and DNA repair processes, making it an attractive drug development target. The study employs advanced computational techniques, including molecular modeling, structure-based virtual screening (SBVS), ADMET profiling, molecular docking, and dynamic simulations.

RESULTS

The study sources ligand molecules from the extensive MCULE database and utilizes strict filters to ensure that the compounds meet essential criteria. Through these techniques, the research identifies MCULE-3763806507-0-9 as a promising antimalarial drug candidate, surpassing the binding affinity of potential antimalarial drugs. However, it is essential to underscore that drug-like properties are primarily based on in silico experiments, and wet lab experiments are necessary to validate these candidates' therapeutic potential.

INTERPRETATION CONCLUSION

This study represents a critical step in addressing the challenge of drug resistance in the fight against malaria.

Insights into the mechanism of peptide fibril growth on gold surface

Understanding the formation of β-fibrils over the gold surface is of paramount interest in nano-bio-medicinal Chemistry. The intricate mechanism of self-assembly of neurofibrillogenic peptides and their growth over the gold surface remains elusive, as experiments are limited in unveiling the microscopic dynamic details, in particular, at the early stage of the peptide aggregation. In this work, we carried out equilibrium molecular dynamics and enhanced sampling simulations to elucidate the underlying mechanism of the growth of an amyloid-forming sequence of tau fragments over the gold surface. Our results disclose that the collective intermolecular interactions between the peptide chains and peptides with the gold surface facilitate the peptide adsorption, followed by integration, finally leading to the fibril formation.

Computer-assisted discovery and evaluation of potential ribosomal protein S6 kinase beta 2 inhibitors

S6K2 is an important protein in mTOR signaling pathway and cancer. To identify potential S6K2 inhibitors for mTOR pathway treatment, a virtual screening of 1,575,957 active molecules was performed using PLANET, AutoDock GPU, and AutoDock Vina, with their classification abilities compared. The MM/PB(GB)SA method was used to identify four compounds with the strongest binding energies. These compounds were further investigated using molecular dynamics (MD) simulations to understand the properties of the S6K2/ligand complex. Due to a lack of available 3D structures of S6K2, OmegaFold served as a reliable 3D predictive model with higher evaluation scores in SAVES v6.0 than AlphaFold, AlphaFold2, and RoseTTAFold2. The 150 ns MD simulation revealed that the S6K2 structure in aqueous solvation experienced compression during conformational relaxation and encountered potential energy traps of about 19.6 kJ mol-1. The virtual screening results indicated that Lys75 and Lys99 in S6K2 are key binding sites in the binding cavity. Additionally, MD simulations revealed that the ligands remained attached to the activation cavity of S6K2. Among the compounds, compound 1 induced restrictive dissociation of S6K2 in the presence of a flexible region, compound 8 achieved strong stability through hydrogen bonding with Lys99, compound 9 caused S6K2 tightening, and the binding of compound 16 was heavily influenced by hydrophobic interactions. This study suggests that these four potential inhibitors with different mechanisms of action could provide potential therapeutic options.

Computational screening of damaging nsSNPs in human SOD1 genes associated with amyotrophic lateral sclerosis identifies destabilising effects of G38R and G42D mutations through in silico evaluation

Amyotrophic lateral sclerosis (ALS), a complicated neurodegenerative disorder affected by hereditary and environmental variables, is a condition. In this study, the genetic makeup of ALS is investigated, with a focus on the SOD1 gene's single-nucleotide polymorphisms (SNPs) and their ability to affect disease risk. Eleven high-risk missense variations that may impair the functionality of the SOD1 protein were discovered after a thorough examination of SNPs in the SOD1 gene. These mutations were chosen using a variety of prediction approaches, highlighting their importance in the aetiology of ALS. Notably, it was discovered that the stability of the SOD1 wild-type protein structure was compromised by the G38R and G42D SOD1 variants. Additionally, Edaravone, a possible ALS medication, showed a greater affinity for binding mutant SOD1 structures, pointing to potential personalised treatment possibilities. The high-risk SNPs discovered in this investigation seem to have functional effects, especially on the stability of proteins and their interactions with other molecules. This study clarifies the complex genetics of ALS and offers insights into how these genetic variations may affect the effectiveness of therapeutic interventions, particularly in the context of edaravone. In this study advances our knowledge of the genetic mechanisms causing ALS vulnerability and prospective therapeutic strategies. Future studies are necessary to confirm these results and close the gap between individualised clinical applications and improved ALS care.

Construction of an aerolysin-based multi-epitope vaccine against Aeromonas hydrophila: an in silico machine learning and artificial intelligence-supported approach

Aeromonas hydrophila, a gram-negative coccobacillus bacterium, can cause various infections in humans, including septic arthritis, diarrhea (traveler's diarrhea), gastroenteritis, skin and wound infections, meningitis, fulminating septicemia, enterocolitis, peritonitis, and endocarditis. It frequently occurs in aquatic environments and readily contacts humans, leading to high infection rates. This bacterium has exhibited resistance to numerous commercial antibiotics, and no vaccine has yet been developed. Aiming to combat the alarmingly high infection rate, this study utilizes in silico techniques to design a multi-epitope vaccine (MEV) candidate against this bacterium based on its aerolysin toxin, which is the most toxic and highly conserved virulence factor among the Aeromonas species. After retrieval, aerolysin was processed for B-cell and T-cell epitope mapping. Once filtered for toxicity, antigenicity, allergenicity, and solubility, the chosen epitopes were combined with an adjuvant and specific linkers to create a vaccine construct. These linkers and the adjuvant enhance the MEV's ability to elicit robust immune responses. Analyses of the predicted and improved vaccine structure revealed that 75.5%, 19.8%, and 1.3% of its amino acids occupy the most favored, additional allowed, and generously allowed regions, respectively, while its ERRAT score reached nearly 70%. Docking simulations showed the MEV exhibiting the highest interaction and binding energies (-1,023.4 kcal/mol, -923.2 kcal/mol, and -988.3 kcal/mol) with TLR-4, MHC-I, and MHC-II receptors. Further molecular dynamics simulations demonstrated the docked complexes' remarkable stability and maximum interactions, i.e., uniform RMSD, fluctuated RMSF, and lowest binding net energy. In silico models also predict the vaccine will stimulate a variety of immunological pathways following administration. These analyses suggest the vaccine's efficacy in inducing robust immune responses against A. hydrophila. With high solubility and no predicted allergic responses or toxicity, it appears safe for administration in both healthy and A. hydrophila-infected individuals.

Classification of MLH1 Missense VUS Using Protein Structure-Based Deep Learning-Ramachandran Plot-Molecular Dynamics Simulations Method

Pathogenic variation in DNA mismatch repair (MMR) gene MLH1 is associated with Lynch syndrome (LS), an autosomal dominant hereditary cancer. Of the 3798 MLH1 germline variants collected in the ClinVar database, 38.7% (1469) were missense variants, of which 81.6% (1199) were classified as Variants of Uncertain Significance (VUS) due to the lack of functional evidence. Further determination of the impact of VUS on MLH1 function is important for the VUS carriers to take preventive action. We recently developed a protein structure-based method named "Deep Learning-Ramachandran Plot-Molecular Dynamics Simulation (DL-RP-MDS)" to evaluate the deleteriousness of MLH1 missense VUS. The method extracts protein structural information by using the Ramachandran plot-molecular dynamics simulation (RP-MDS) method, then combines the variation data with an unsupervised learning model composed of auto-encoder and neural network classifier to identify the variants causing significant change in protein structure. In this report, we applied the method to classify 447 MLH1 missense VUS. We predicted 126/447 (28.2%) MLH1 missense VUS were deleterious. Our study demonstrates that DL-RP-MDS is able to classify the missense VUS based solely on their impact on protein structure.

Design of a novel multi-epitopes vaccine against Escherichia fergusonii: a pan-proteome based in- silico approach

Escherichia fergusonii a gram-negative rod-shaped bacterium in the Enterobacteriaceae family, infect humans, causing serious illnesses such as urinary tract infection, cystitis, biliary tract infection, pneumonia, meningitis, hemolytic uremic syndrome, and death. Initially treatable with penicillin, antibiotic misuse led to evolving resistance, including resistance to colistin, a last-resort drug. With no licensed vaccine, the study aimed to design a multi-epitope vaccine against E. fergusonii. The study started with the retrieval of the complete proteome of all known strains and proceeded to filter the surface exposed virulent proteins. Seventeen virulent proteins (4 extracellular, 4 outer membranes, 9 periplasmic) with desirable physicochemical properties were identified from the complete proteome of known strains. Further, these proteins were processed for B-cell and T-cell epitope mapping. Obtained epitopes were evaluated for antigenicity, allergenicity, solubility, MHC-binding, and toxicity and the filtered epitopes were fused by specific linkers and an adjuvant into a vaccine construct. Structure of the vaccine candidate was predicted and refined resulting in 78.1% amino acids in allowed regions and VERIFY3D score of 81%. Vaccine construct was docked with TLR-4, MHC-I, and MHC-II, showing binding energies of -1040.8 kcal/mol, -871.4 kcal/mol, and -1154.6 kcal/mol and maximum interactions. Further, molecular dynamic simulation of the docked complexes was carried out resulting in a significant stable nature of the docked complexes (high B-factor and deformability values, lower Eigen and high variance values) in terms of intermolecular binding conformation and interactions. The vaccine was also reported to stimulate a variety of immunological pathways after administration. In short, the designed vaccine revealed promising predictions about its immune protective potential against E. fergusonii infections however experimental validation is needed to validate the results.