Multi-level scoring approach to discover multi-targeting potency of medicinal plant phytochemicals against protein targets in SARS-CoV-2 and human ACE-2 receptor

SARS-CoV-2 pandemic has become a major threat to human healthcare and world economy. Due to the rapid spreading and deadly nature of infection, we are in a situation to develop quick therapeutics to combat SARS-CoV-2. In this study, we have adopted a multi-level scoring approach to identify multi-targeting potency of bioactive compounds in selected medicinal plants and compared its efficacy with two reference drugs, Nafamostat and Acalabrutinib which are under clinical trials to treat SARS-CoV-2. In particular, we employ molecular docking and implicit solvent free energy calculations (as implemented in the Molecular Mechanics -Generalized Born Surface Area approach) and QM fragmentation approach for validating the potency of bioactive compounds from the selected medicinal plants against four different viral targets and one human receptor (Angiotensin-converting enzyme 2 -ACE-2) which facilitates the SARS-CoV-2entry into the cell. The protein targets considered for the study are viral 3CL main protease (3CLpro), papain-like protease (PLpro), RNA dependent RNA polymerase (RdRp), and viral spike protein-human hACE-2 complex (Spike:hACE2)including human protein target (hACE-2). Herein, thereliable multi-level scoring approach was used to validate the mechanism behind the multi-targeting potency of selected phytochemicals from medicinal plants. The present study evidenced that the phytochemicals Chebulagic acid, Stigmosterol, Repandusinic acid and Geranin exhibited efficient inhibitory activity against PLpro while Chebulagic acid was highly active against 3CLpro. Chebulagic acid andGeranin also showed excellent target specific activity against RdRp.Luteolin, Quercetin, Chrysoeriol and Repandusinic acid inhibited the interaction of viral spike protein with human ACE-2 receptor. Moreover Piperlonguminine and Piperine displayed significant inhibitory activity against human ACE-2 receptor. Therefore, the identified compounds namely Chebulagic acid, Geranin and Repandusinic acid can serve as potent multi-targeting phytomedicine for treating COVID-1

Computer-Aided vaccine design for selected positive-sense single stranded RNA viruses

Spontaneous mutations and lack of replication fidelity in positive-sense single stranded RNA viruses (+ssRNA virus) result in emergence of genetic variants with diverse viral morphogenesis and surface proteins that affect its antigenicity. This high mutability in +ssRNA viruses has induced antiviral drug resistance and ability to overcome vaccines that subsequently resulted in rapid viral evolution and high mortality rate in human and livestock. Computer aided vaccine design and immunoinformatics play a crucial role in expediting the vaccine production protocols, antibody production and identifying suitable immunogenic regions or epitopes from the genome sequences of the pathogens. T cell and B cell epitopes can be identified in pathogens by immunoinformatics algorithms and methods that enhance the analysis of protective immunity, vaccine safety, immunity modelling and vaccine efficacy. This rapid and cost-effective computational vaccine design promotes development of potential vaccine that could induce immune response in host against rapidly mutating pathogens like +ssRNA viruses. Epitope-based vaccine is a striking concept that has been widely employed in recent years to construct vaccines targeting rapidly mutating +ssRNA viruses. Therefore, the present review provides an overview about the current progress and methodology in computer-aided vaccine design for the most notable +ssRNA viruses namely Hepatitis C virus, Dengue virus, Chikungunya virus and Coronaviruses. This review also highlights the applications of various immunoinformatics tools for vaccine design and for modelling immune response against +ssRNA viruses

Temperature-induced dynamical conformational disorder in 4-vinyl benzoic acid molecular crystals: A molecular simulation study

Extensive molecular simulations are carried out as a function of temperature to understand and quantify the conformational disorder in molecular crystals of 4-vinyl benzoic acid. The conformational disorder is found to be dynamic and associated with a flip-flop motion of vinyl groups. The population of minor conformer is less than 3% up to 300 Kand is 13.2% at 350 K and these results are consistent with the experimental observations. At still higher temperatures, the population of minor conformer increases up to 25%. The evolution of structure at both molecular and unit-cell level of the molecular crystal as a function of temperature has been characterized by various quantities such as radial distribution functions, average cell parameters, volume,and interaction energies. The van't Hoff plot shows a nonlinear behavior at lower temperatures as it has been reported recently by O gawa and co-workers in the case of stilbene, azobenzene [J. Am. Chem.Soc. 126, 3539 (2004)], and N-(4-methylbenzylidene)-4-methylaniline[Acta Crystallogr, Sect. B. Struct. Sci. B60, 589 (2004)] molecularcrystals. A set of rigid body simulations were also carried out to quantify the effect of conformational disorder on structural quantities such as unit-cell volume and interaction energy. The anomalous shrinkage of vinyl C=C bond length as a function of temperature has been explained by combining the results of simulations and a set of constrained optimizations using ab initio electronic structure calculations for various molecular structures differing in torsional angle

What does pressure decide to cook with orientationally disordered plastic phase of cubane: An orientational glass or crystal?

The effect of pressure on the structure and reorientational motion of molecules in orientationally disordered (OD) crystalline phase of cubane has been investigated in detail using variable shape molecular simulations in constant-pressure constant-temperature ensemble. Complete orientational ordering occurs at a pressure of 1.0 GPa and the OD phase transforms to an orientationally ordered phase at this pressure. The transition is associated with a kink in the variation of structural parameters such as cell parameters, unit-cell volume, and interaction energy. This transition is also associated with an anomaly in specific heat. Above this transition pressure, the structural quantities display only smaller changes with further increase in pressure. The structure of high-pressure orientationally ordered (HPOO) phase has been characterized using radial distribution functions and orientational distribution function. From detailed analysis of the structure of HPOO phase we conclude that it is isostructural with low-temperature orientationally ordered phase. The OD phase has four times larger compressibility than the HPOO phase

Orientational melting and reorientational motion in a cubane molecular crystal: A molecular simulation study

Detailed molecular simulations are carried out to investigate the effect of temperature on orientational order in cubane molecular crystal. We report a transition from an orientationally ordered to an orientationally disordered plastic crystalline phase in the temperature range 425-450 K. This is similar to the experimentally reported transition at 395 K. The nature of this transition is first order and is associated with a 4.8% increase in unit Cell volume that is comparable to the experimentally reported unit cell volume change of 5.4% (Phys. Rev. Lett. 1997, 78, 4938). An orientational order parameter, eta(T), has been defined in terms of average angle of libration of a molecular 3-fold axis and the orientational melting has been characterized by using eta(T). The orientational melting is associated with an anomaly in specific heat at constant pressure (C-p) and compressibility (kappa). The enthalpy of transition and entropy of transition associated with this orientational melting are 20.8 J mol(-1) and 0.046 J mol(-1) K-1, respectively. The structure of crystalline as well as plastic crystalline phases is characterized by using various radial distribution functions and orientational distribution functions. The coefficient of thermal expansion of the plastic crystalline phase is more than twice that of the crystalline phase

Pressure-Induced Ordering in Adamantane: A Monte Carlo Simulation Study

Isothermal-isobaric ensemble Monte Carlo simulation of adamantane has been carried out with a variable shape simulation cell. The low-temperature crystalline phase and the room-temperature plastic crystalline phases have been studied employing the modified Williams potential. We show that at room temperature, the plastic crystalline phase transforms to the crystalline phase on increase in pressure. Further, we show that this is the same phase as the low-temperature ordered tetragonal phase of adamantane. The high-pressure ordered phase appears to be characterized by a slightly larger shift of the first peak toward a lower value of r in C-C, C-H, and H-H radial distribution functions as compared to the low-temperature tetragonal phase. The coexistence curve between the crystalline and plastic crystalline phase has been obtained approximately up to a pressure of 4 GPa