Hierarchical Visual Analysis and Steering Framework for Astrophysical Simulations

A framework for accelerating modern long-running astrophysical simulations is presented, which is based on a hierarchical architecture where computational steering in the high-resolution run is performed under the guide of knowledge obtained in the gradually refined ensemble analyses. Several visualization schemes for facilitating ensemble management, error analysis, parameter grouping and tuning are also integrated owing to the pluggable modular design. The proposed approach is prototyped based on the Flash code, and it can be extended by introducing userdefined visualization for specific requirements. Two real-world simulations, i.e., stellar wind and supernova remnant, are carried out to verify the proposed approach.

Free energy calculation of single molecular interaction using Jarzynski's identity method:the case of HIV-1 protease inhibitor system

Jarzynski＇ identity （JI） method was suggested a promising tool for reconstructing free energy landscape of biomolecular interactions in numerical simulations and ex- periments. However, JI method has not yet been well tested in complex systems such as ligand-receptor molecular pairs. In this paper, we applied a huge number of steered molec- ular dynamics （SMD） simulations to dissociate the protease of human immunodeficiency type I virus （HIV-1 protease） and its inhibitors. We showed that because of intrinsic com- plexity of the ligand-receptor system, the energy barrier pre- dicted by JI method at high pulling rates is much higher than experimental results. However, with a slower pulling rate and fewer switch times of simulations, the predictions of JI method can approach to the experiments. These results sug- gested that the JI method is more appropriate for reconstruct- ing free energy landscape using the data taken from experi- ments, since the pulling rates used in experiments are often much slower than those in SMD simulations. Furthermore, we showed that a higher loading stiffness can produce higher precision of calculation of energy landscape because it yields a lower mean value and narrower bandwidth of work distri- bution in SMD simulations.

Inhibiting the oligomerization of mycobacterial DNA-directed RNA polymerase(RNAP)using natural compound via in-silico techniques

Mycobacterium tuberculosis(Mtb)is responsible for the spread of tuberculosis(TB).The current study employed virtual screening of 2569 natural compounds against the DNA-directed RNA polymerase(RNAP)of Mtb to identify the possible binders that can inhibit its function.The in-silico methodology included molecular docking to the compounds,further,the stability and flexibility of the best complexes were studied using molecular dynamics simulation,the MM/GBSA binding free energy technique with energy decomposition,PCA,FEL,steered MD simulation,and umbrella sampling.Individual virtual screenings were conducted for the five RNAP subunits(chains A,B,C,D,and E)to identify a compound capable of inhibiting RNAP oligomerization.A promising compound,isoestradiol 3-benzoate,exhibited a low binding score(−7.28kcal/mol to−8.21kcal/mol)and showed binding ability with all subunits of the protein.Thus,the five complexes with isoestradiol 3-benzoate were selected for molecular dynamics simulation analysis.Furthermore,RMSD showed that isoestradiol 3-benzoate bound with chain E showed the lowest RMSD of 0.49nm,while with chains A and B it had the most stable and consistent conformations with RMSD of 1.75nm and 1.2nm,respectively.The H-bond between isoestradiol 3-benzoate and chains C and E showed the highest occupancy(58.27%,45.33%,and 50.80%,42.25%,11.75%).Moreover,the MMPBSA technique showed that isoestradiol 3-benzoate had a strong association with chains B and C(ΔGbind=−126.25±2.03 and−129.27±2.25).Additionally,free energy decomposition,PCA,FEL-steered MD simulation,and umbrella sampling were also performed to validate the association of the ligand with the protein.Isoestradiol 3-benzoate binds strongly to chains B and E;therefore,it should be considered as viable candidate for inhibiting the formation of RNAP protein complex,concluded in this study.