T cell epitope based vaccine design while targeting outer capsid proteins of rotavirus strains infecting neonates: an immunoinformatics approach.

Gastrointestinal diarrhea is majorly caused by the rotavirus (RV) in the children who generally are under the age group of 5 years. WHO estimates that ∼95% of the children contract RV infection, by this age. The disease is highly contagious; notably in many cases, it is proven fatal with high mortality rates especially in the developing countries. In India alone, an estimated 145,000 yearly deaths occurs due to RV related gastrointestinal diarrhea. WHO pre-qualified vaccines that are available for RV are all live attenuated vaccines with modest efficacy range between 40 and 60%. Further, the risk of intussusceptions has been reported in some children on RV vaccination. Thus, in a quest to develop alternative candidate to overcome challenges associated with these oral vaccines, we chose immunoinformatics approach to design a multi-epitope vaccine (MEV) while targeting the outer capsid viral proteinsVP4 and VP7 of the neonatal strains of rotavirus. Interestingly, ten epitopes, that is, six CD8+T-cells and four CD4+T-cell epitopes were identified which were predicted to be antigenic, non-allergic, non-toxic and stable. These epitopes were then linked to adjuvants, linkers, and PADRE sequences to create a multi-epitope vaccine for RV. The in silico designed RV-MEV and human TLR5 complex displayed stable interactions during molecular dynamics simulations. Further, the immune simulation studies of RV-MEV corroborated that the vaccine candidate emerges as a promising immunogen. Future investigations while performing in vitro and in vivo analyses with designed RV-MEV construct are highly desirable to warrant the potential of this vaccine candidate in protective immunity against different strains of RVs infecting neonates.Communicated by Ramaswamy H. Sarma.

Effects of Conformational Constraint on Peptide Solubility Limits.

Liquid-liquid phase separation of proteins preferentially involves intrinsically disordered proteins or disordered regions. Understanding the solution chemistry of these phase separations is key to learning how to quantify and manipulate systems that involve such processes. Here, we investigate the effect of cyclization on the liquid-liquid phase separation of short polyglycine peptides. We simulated separate aqueous systems of supersaturated cyclic and linear GGGGG and observed spontaneous liquid-liquid phase separation in each of the solutions. The cyclic GGGGG phase separates less robustly than linear GGGGG and has a higher aqueous solubility, even though linear GGGGG has a more favorable single molecule solvation free energy. The versatile and abundant interpeptide contacts formed by the linear GGGGG stabilize the condensed droplet phase, driving the phase separation in this system. In particular, we find that van der Waals close contact interactions are enriched in the droplet phase as opposed to electrostatic interactions. An analysis of the change in backbone conformational entropy that accompanies the phase transition revealed that cyclic peptides lose significantly less entropy in this process as expected. However, we find that the enhanced interaction enthalpy of linear GGGGG in the droplet phase is enough to compensate for a larger decrease in conformational entropy.

Mechanistic Insights into SARS-CoV-2 Main Protease Inhibition Reveals Hotspot Residues.

The main protease (Mpro) is a key enzyme responsible for severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) replication that causes the spread of the global pandemic novel coronavirus (nCOVID-19) infection. In the present study, multiple computational approaches such as docking, long-range molecular dynamics (MD) simulations, and binding free-energy (BFE) estimation techniques were employed to investigate the mechanistic basis of the high-affinity inhibitors─GC-376, Calpain XII, and Calpain II (hereafter Calpain as Cal) from the literature─binding to Mpro. Redocking GC-376 and docking Cal XII and Cal II inhibitors to Mpro were able to reproduce all crucial interactions like the X-ray conformation. Subsequently, the apo (ligand-free) and three holo (ligand-bound) complexes were subjected to extensive MD simulations, which revealed that the ligand binding did not alter the overall Mpro structural features, whereas the heatmap analysis showed that the residues located in subsites S1 and S2, the catalytic dyad, and the 45TSEDMLN51 loop in Mpro exhibit a conformational deviation. Moreover, the BFE estimation method was used to elucidate the crucial thermodynamic properties, which revealed that Coulomb, solvation surface accessibility (Solv_SA), and lipophilic components contributed significant energies for complex formation. The decomposition of the total BFE to per-residue showed that H41, H163, M165, Q166, and Q189 residues contributed maximum energies. The overall results from the current investigation might be valuable for designing novel anti-Mpro inhibitors.

Structural Insights in Mammalian Sialyltransferases and Fucosyltransferases: We Have Come a Long Way, but It Is Still a Long Way Down.

Mammalian cell surfaces are modified with complex arrays of glycans that play major roles in health and disease. Abnormal glycosylation is a hallmark of cancer; terminal sialic acid and fucose in particular have high levels in tumor cells, with positive implications for malignancy. Increased sialylation and fucosylation are due to the upregulation of a set of sialyltransferases (STs) and fucosyltransferases (FUTs), which are potential drug targets in cancer. In the past, several advances in glycostructural biology have been made with the determination of crystal structures of several important STs and FUTs in mammals. Additionally, how the independent evolution of STs and FUTs occurred with a limited set of global folds and the diverse modular ability of catalytic domains toward substrates has been elucidated. This review highlights advances in the understanding of the structural architecture, substrate binding interactions, and catalysis of STs and FUTs in mammals. While this general understanding is emerging, use of this information to design inhibitors of STs and FUTs will be helpful in providing further insights into their role in the manifestation of cancer and developing targeted therapeutics in cancer.

Replacing thymine with a strongly pairing fifth Base: A combined quantum mechanics and molecular dynamics study.

The non-natural ethynylmethylpyridone C-nucleoside (W), a thymidine (T) analogue that can be incorporated in oligonucleotides by automated synthesis, has recently been reported to form a high fidelity base pair with adenosine (A) and to be well accommodated in B-DNA duplexes. The enhanced binding affinity for A of W, as compared to T, makes it an ideal modification for biotechnological applications, such as efficient probe hybridization for the parallel detection of multiple DNA strands. In order to complement the experimental study and rationalize the impact of the non-natural W nucleoside on the structure, stability and dynamics of DNA structures, we performed quantum mechanics (QM) calculations along with molecular dynamics (MD) simulations. Consistently with the experimental study, our QM calculations show that the A:W base pair has an increased stability as compared to the natural A:T pair, due to an additional CH-π interaction. Furthermore, we show that mispairing between W and guanine (G) causes a distortion in the planarity of the base pair, thus explaining the destabilization of DNA duplexes featuring a G:W pair. MD simulations show that incorporation of single or multiple consecutive A:W pairs in DNA duplexes causes minor changes to the intra- and inter-base geometrical parameters, while a moderate widening/shrinking of the major/minor groove of the duplexes is observed. QM calculations applied to selected stacks from the MD simulations also show an increased stacking energy for W, over T, with the neighboring bases.

Occurrence and stability of lone pair-π and OH-π interactions between water and nucleobases in functional RNAs.

We identified over 1000 instances of water-nucleobase stacking contacts in a variety of RNA molecules from a non-redundant set of crystal structures with resolution ≤3.0 Å. Such contacts may be of either the lone pair-π (lp-π) or the OH-π type, in nature. The distribution of the distances of the water oxygen from the nucleobase plane peaks at 3.5 Å for A, G and C, and approximately at 3.1-3.2 Å for U. Quantum mechanics (QM) calculations confirm, as expected, that the optimal energy is reached at a shorter distance for the lp-π interaction as compared to the OH-π one (3.0 versus 3.5 Å). The preference of each nucleobase for either type of interaction closely correlates with its electrostatic potential map. Furthermore, QM calculations show that for all the nucleobases a favorable interaction, of either the lp-π or the OH-π type, can be established at virtually any position of the water molecule above the nucleobase skeleton, which is consistent with the uniform projection of the OW atoms over the nucleobases ring we observed in the experimental occurrences. Finally, molecular dynamics simulations of a model system for the characterization of water-nucleobase stacking contacts confirm the stability of these interactions also under dynamic conditions.

Nanoscale Dynamics and Energetics of Proteins and Protein-Nucleic Acid Complexes in Classical Molecular Dynamics Simulations.

The present article describes techniques for classical simulations of proteins and protein-nucleic acid complexes, revealing their dynamics and protein-substrate binding energies. The approach is based on classical atomistic molecular dynamics (MD) simulations of the experimentally determined structures of the complexes. MD simulations can provide dynamics of complexes in realistic solvents on microsecond timescales, and the free energy methods are able to provide Gibbs free energies of binding of substrates, such as nucleic acids, to proteins. The chapter describes methodologies for the preparation of computer models of biomolecular complexes and free energy perturbation methodology for evaluating Gibbs free energies of binding. The applications are illustrated with examples of snapshots of proteins and their complexes with nucleic acids, as well as the precise Gibbs free energies of binding.

Computational Model and Dynamics of Monomeric Full-Length APOBEC3G.

APOBEC3G (A3G) is a restriction factor that provides innate immunity against HIV-1 in the absence of viral infectivity factor (Vif) protein. However, structural information about A3G, which can aid in unraveling the mechanisms that govern its interactions and define its antiviral activity, remains unknown. Here, we built a computer model of a full-length A3G using docking approaches and molecular dynamics simulations, based on the available X-ray and NMR structural data for the two protein domains. The model revealed a large-scale dynamics of the A3G monomer, as the two A3G domains can assume compact forms or extended dumbbell type forms with domains visibly separated from each other. To validate the A3G model, we performed time-lapse high-speed atomic force microscopy (HS-AFM) experiments enabling us to get images of a fully hydrated A3G and to directly visualize its dynamics. HS-AFM confirmed that A3G exists in two forms, a globular form (∼84% of the time) and a dumbbell form (∼16% of the time), and can dynamically switch from one form to the other. The obtained HS-AFM results are in line with the computer modeling, which demonstrates a similar distribution between two forms. Furthermore, our simulations capture the complete process of A3G switching from the DNA-bound state to the closed state. The revealed dynamic nature of monomeric A3G could aid in target recognition including scanning for cytosine locations along the DNA strand and in interactions with viral RNA during packaging into HIV-1 particles.