Coarse-grained molecular dynamics-guided immunoinformatics to explain the binder and non-binder classification of Cytotoxic T-cell epitope for SARS-CoV-2 peptide-based vaccine discovery.

Epitope-based peptide vaccine can elicit T-cell immunity against SARS-CoV-2 to clear the infection. However, finding the best epitope from the whole antigen is challenging. A peptide screening using immunoinformatics usually starts from MHC-binding peptide, immunogenicity, cross-reactivity with the human proteome, to toxicity analysis. This pipeline classified the peptides into three categories, i.e., strong-, weak-, and non-binder, without incorporating the structural aspect. For this reason, the molecular detail that discriminates the binders from non-binder is interesting to be investigated. In this study, five CTL epitopes against HLA-A*02:01 were identified from the coarse-grained molecular dynamics-guided immunoinformatics screening. The strong binder showed distinctive activities from the non-binder in terms of structural and energetic properties. Furthermore, the second residue from the nonameric peptide was most important in the interaction with HLA-A*02:01. By understanding the nature of MHC-peptide interaction, we hoped to improve the chance of finding the best epitope for a peptide vaccine candidate.

Unveiling the Anti-Cancer Potential of Onoceranoid Triterpenes from Lansium domesticum Corr. cv. kokosan: An In Silico Study against Estrogen Receptor Alpha.

Breast cancer is a significant global concern, with tamoxifen, the standard treatment, raising long-term safety issues due to side effects. In this study, we evaluated the potential of five onoceranoid triterpenes from Lansium domesticum Corr. cv. kokosan against estrogen receptor alpha (ERα) using in silico techniques. Utilizing molecular docking, Lipinski's rule of five, in silico ADMET, and molecular dynamics simulations, we assessed the potency of five onoceranoid triterpenes against ERα. Molecular docking indicated competitive binding energies for these triterpenes relative to the active form of tamoxifen (4OHT) and estradiol, an ERα native ligand. Three triterpenes met drug-likeness criteria with favorable ADMET profiles. Notably, 2 demonstrated superior binding affinity in molecular dynamics simulations, outperforming estradiol, closely followed by 3 and 4. Hierarchical clustering on principal components (HCPC) and the spatial distribution of contact surface area (CSA) analyses suggest that these triterpenes, especially 2, may act as antagonist ligands akin to 4OHT. These findings highlight the potential of onoceranoid triterpenes in treating ERα-related breast cancer.

Multilayer Model of Gold Nanoparticles (AuNPs) and Its Application in the Classical Molecular Dynamics Simulation of Citrate-Capped AuNPs.

Studies on the interaction between gold nanoparticles (AuNPs) and functional proteins have been useful in developing diagnostic and therapeutic agents. Such studies require a realistic computational model of AuNPs for successful molecular design works. This study offers a new multilayer model of AuNPs to address the inconsistency between its molecular mechanics' interpretation and AuNP's plasmonic nature. We performed partial charge quantum calculation of AuNPs using Au13 and Au55 models. The result showed that it has partial negative charges on the surface and partial positive charges on the inner part, indicating that the AuNP model should be composed of multiatom types. We tested the partial charge parameters of these gold (Au) atoms in classical molecular dynamics simulation (CMD) of AuNPs. The result showed that our parameters performed better in simulating the adsorption of Na+ and dicarboxy acetone in terms of consistency with surface charge density than the zero charges Au in the interface force field (IFF). We proposed that the multiple-charged AuNP model can be developed further into a simpler four-atom type of Au in a larger AuNP size.

The nanomolar affinity of C-phycocyanin from virtual screening of microalgal bioactive as potential ACE2 inhibitor for COVID-19 therapy.

The global pandemic of COVID-19 caused by SARS-CoV-2 has caused more than 400 million infections with more than 5.7 million deaths worldwide, and the number of validated therapies from natural products for treating coronavirus infections needs to be increased. Therefore, the virtual screening of bioactive compounds from natural products based on computational methods could be an interesting strategy. Among many sources of bioactive natural products, compounds from marine organisms, particularly microalgae and cyanobacteria, can be potential antiviral agents. The present study investigates bioactive antiviral compounds from microalgae and cyanobacteria as a potential inhibitor of SARS-CoV-2 by targeting Angiotensin-Converting Enzyme II (ACE2) using integrated in silico and in vitro approaches. Our in silico analysis demonstrates that C-Phycocyanin (CPC) can potentially inhibit the binding of ACE2 receptor and SARS-CoV-2 with the docking score of -9.7 kcal mol-1. This score is relatively more favorable than the native ligand on ACE2 receptor. Molecular dynamics simulation also reveals the stability interaction between both CPC and ACE2 receptor with a root mean square deviation (RMSD) value of 1.5 Å. Additionally, our in vitro analysis using the surface plasmon resonance (SPR) method shows that CPC has a high affinity for ACE2 with a binding affinity range from 5 to 125 µM, with KD 3.37 nM. This study could serve as a reference to design microalgae- or cyanobacteria-based antiviral drugs for prophylaxis in SARS-CoV-2 infections.

Residual Interactions of LL-37 with POPC and POPE:POPG Bilayer Model Studied by All-Atom Molecular Dynamics Simulation.

LL-37 is a membrane-active antimicrobial peptide (AMP) that could disrupt the integrity of bacterial membranes due to its inherent cationic and amphipathic nature. Developing a shorter derivative of a long peptide such as LL-37 is of great interest, as it can reduce production costs and cytotoxicity. However, more detailed information about the residual interaction between LL-37 and the membrane is required for further optimization. Previously, molecular dynamics simulation using mixed all-atom and united-atom force fields showed that LL-37 could penetrate the bilayer membrane. This study aimed to perform all-atom molecular dynamics simulations, highlighting the residual interaction of LL-37 with the simplest model of the bacterial membrane, POPE:POPG (2:1), and compare its interaction with the POPC, which represents the eukaryotic membrane. The result showed leucine-leucine as the leading residues of LL-37 that first contact the membrane surface. Then, the cationic peptide of LL-37 started to penetrate the membrane by developing salt bridges between positively charged amino acids, Lys-Arg, and the exposed phosphate group of POPE:POPG, which is shielded in POPC. Residues 18 to 29 are suggested as the core region of LL-37, as they actively interact with the POPE:POPG membrane, not POPC. These results could provide a basis for modifying the amino acid sequence of LL-37 and developing a more efficient design for LL-37 derivatives.

Molecular Dynamics Simulation of T10609C and C10676G Mutations of Mitochondrial ND4L Gene Associated With Proton Translocation in Type 2 Diabetes Mellitus and Cataract Patients.

The mutation rate of mitochondrial DNA (mtDNA) is 17 times higher than nuclear DNA, and these mutations can cause mitochondrial disease in 1 of 10.000 people. The T10609C mutation was identified in type 2 diabetes mellitus (T2DM) patients and the C10676G mutation in cataract patients, with both mutations occurring in the ND4L gene of mtDNA that encodes ND4L protein. ND4L protein, a subunit of complex I in the respiratory complex, has been shown to play a role in the proton translocation process. The purpose of this study was to investigate the effect of both mutations on the proton translocation mechanism. Mutation mapping showed changes in amino acids M47T (T10609C) and C69W (C10676G). The 100 ns molecular dynamics (MD) simulations performed on native and mutants of ND4L-ND6 subunits. It is revealed that the native model had a similar proton translocation pathway to that of complex I from other organisms. Interestingly, the mutant M47T and C69W showed the interruption of the translocation pathway by a hydrogen bond formation between Glu34 and Tyr157. It is observed that the mutations were restricting the passage of water molecules through the transmembrane region. These results could help to develop the computational assay for the validation of a specific genetic biomarker for T2DM and cataracts.