GROMACS 2020在ROCm平台上的移植与优化

GROMACS是应用广泛的开源分子动力学模拟软件,当前主要通过CUDA使用NVIDIA GPU进行加速计算。ROCm是一个开源的高性能异构计算平台。基于ROCm平台的HIP编程语言,首次实现了GROMACS 2020系列在ROCm平台上的完整移植。在MI50 GPU上,以一个复杂离子液体模拟算例为目标,使用GPU性能分析工具rocprof对移植代码进行了性能分析。针对MI50硬件特性,先后对成键力核函数、静电力的PME核函数和短程非成键力核函数进行了优化,优化后运行目标算例的性能相比初始版本整体上获得了约2.8倍的加速比,在 MI50上的性能高于GROMACS原版OpenCL代码60.5%,相对纯CPU版本有约2.7倍的加速比。在另外2个具有代表性算例的单结点测试以及离子液体算例的多结点扩展性测试中,优化后的代码也达到了较好的性能提升,这表明所采用的优化操作具有一定的通用性。     

A MARTINI extension for Pseudomonas aeruginosa PAO1 lipopolysaccharide

We report a course-grained, large scale simulation of the outer membrane from Pseudomonas aeruginosa. Using the MARTINI force field approach of 4-to-1 atom mapping, we simulate an asymmetrically constructed bilayer with over 1100 rough lipopolysaccharide (LPS) and 3100 16:0-18:1-phosphatidylethanolamine. We achieve 90-fold improvement in computational efficiency on a system much larger than reasonable for all-atom simulation. We also compare a coarse-grained LPS/LPS bilayer simulation with known parameters determined from neutron diffraction.     

The Flexible Rare Event Sampling Harness System (FRESHS)

We present the software package FRESHS (http://www.freshs.org) for parallel simulation of rare events using sampling techniques from the ‘splitting’ family of methods. Initially, Forward Flux Sampling (FFS) and Stochastic Process Rare Event Sampling (SPRES) have been implemented. These two methods together make rare event sampling available for both quasi-static and full non-equilibrium regimes. Our framework provides a plugin system for software implementing the underlying physics of the system of interest. At present, example plugins exist for our framework to steer the popular MD packages GROMACS, LAMMPS and ESPResSo, but due to the simple interface of our plugin system, it is also easy to attach other simulation software or self-written code. Use of our framework does not require recompilation of the simulation program. The modular structure allows the flexible implementation of further sampling methods or physics engines and creates a basis for objective comparison of different sampling algorithms.     

A PM7 dynamic residue-ligand interactions energy landscape of the BACE1 inhibitory pathway by hydroxyethylamine compounds. Part I: The flap closure process

BACE1 is an enzyme of scientific interest because it participates in the progression of Alzheimer’s disease. Hydroxyethylamines (HEAs) are a family of compounds which exhibit inhibitory activity toward BACE1 at a nanomolar level, favorable pharmacokinetic properties and oral bioavailability. The first step in the inhibition of BACE1 by HEAs consists of their entrance into the protease active site and the resultant conformational change in the protein, from Apo to closed form. These two conformations differ in the position of an antiparallel loop (called the flap) which covers the entrance to the catalytic site. For BACE1, closure of this flap is vital to its catalytic activity and to inhibition of the enzyme due to the new interactions thereby formed with the ligand. In the present study a dynamic energy landscape of residue-ligand interaction energies (ReLIE) measured for 112 amino acids in the BACE1 active site and its immediate vicinity during the closure of the flap induced by 8 HEAs of different inhibitory power is presented. A total of 6.272 million ReLIE calculations, based on the PM7 semiempirical method, provided a deep and quantitative view of the first step in the inhibition of the aspartyl protease. The information suggests that residues Asp93, Asp289, Thr292, Thr293, Asn294 and Arg296 are anchor points for the ligand, accounting for approximately 45% of the total protein-ligand interaction. Additionally, flap closure improved the BACE1-HEA interaction by around 25%. Furthermore, the inhibitory activity of HEAs could be related to the capacity of these ligands to form said anchor point interactions and maintain them over time: the lack of some of these anchor interactions delayed flap closure or impeded it completely, or even caused the flap to reopen. The methodology employed here could be used as a tool to evaluate future structural modifications which lead to improvements in the favorability and stability of BACE1-HEA ReLIEs, aiding in the design of better inhibitors.     

Conformational Solvation Studies of LIGNOLs with Molecular Dynamics and Conductor-Like Screening Model

Molecular dynamics (MD) simulations were performed on sterically hindered α-conidendrin-based chiral 1,4-diols (LIGNOLs) from the naturally occurring lignan hydroxymatairesinol (HMR) using the GROMACS software. The aim of this study was to explore the conformational behaviour of the LIGNOLs in aqueous solution adopting the TIP4P model. The topologies of the LIGNOLs were constructed manually and they were modeled with the OPLS-AA force field implemented in GROMACS. The four most relevant torsional angles in the LIGNOLs were properly analyzed during the simulations. The determining property for the conformation preferred in aqueous solution was found to be the lowest energy in gas phase. The solvation effects on the LIGNOLs were also studied by quantum chemical calculations applying the COnductor-like Screening MOdel (COSMO). The hydration studies of the MD simulations showed that several of these LIGNOLs, produced from a renewable source, have a great potential of acting as chiral catalysts.     

Small molecule interaction with lipid bilayers: A molecular dynamics study of chlorhexidine

Chlorhexidine (CHX) is an effective anti-bacterial agent whose mode of action is thought to be the disruption of the cell membrane. We tested the capability of the Slipids all atom force fields using data from neutron scattering and NMR experiments on the drug chlorhexidine in a 1,2-dimyrisoyl-3-sn-phosphatidylcholine (DMPC) membrane. Since it is not known what the charge of the CHX molecule is inside an apolar environment, a neutral, as well as a +1 and +2 charge model for the molecule were created and tested at several concentrations. This study shows that the location of CHX is minorly dependent on concentration, and dominantly reliant on the charge. The effect of adding CHX to DMPC is a thinning of the membrane, thus increasing the area per lipid.     

Structural insights into type I and type II of nsp4 porcine reproductive and respiratory syndrome virus (nsp4 PRRSV) by molecular dynamics simulations

Porcine reproductive and respiratory virus (PRRSV) causes major economic concerns for the swine industry worldwide. We have performed molecular dynamics simulations (MD) and principle component analysis (PCA) to investigate the role of the catalytic triad and conformational dynamics of type I and type II of nsp4 PRRSV. The results showed that the RMSF of residues 136–142 near the active site of all models was highly flexible. Moreover, we identified the effect of single structural mutations of the catalytic triad. The percentage of residue with a 0.1 nm RMSF value and PCA results revealed that the mutations affected domain I and II suggesting the wild types were more stable than the mutants. At the catalytic triad, the distances between H39 and S118 were very flexible while the distances between H39 and D64 were very stable. H39, D64 and S118 showed high occupancy percentage of the hydrogen bond interaction with many residues that are conserved in PRRSVAS, PRRSVES, LDVC, LDVP and EAV. Moreover, S118 of wild-type protein formed H-bonds with T134 and G135 but these interactions were lost in PRRSVAV (S118A) and PRRSVES (S117A) indicating that the substitution of important H-bond interaction in PRRSVAS (S118A) and PRRSVES (S117A) affected the flexibility around the catalytic triad, conformation and proteolytic activity. Overall, our study may provide the structural basic of the catalytic triad and be useful for testing the protein activity in future experiments.     

Long-range interactions and parallel scalability in molecular simulations

Typical biomolecular systems such as cellular membranes, DNA, and protein complexes are highly charged. Thus, efficient and accurate treatment of electrostatic interactions is of great importance in computational modeling of such systems. We have employed the GROMACS simulation package to perform extensive benchmarking of different commonly used electrostatic schemes on a range of computer architectures (Pentium-4, IBM Power 4, and Apple/IBM G5) for single processor and parallel performance up to 8 nodes—we have also tested the scalability on four different networks, namely Infiniband, GigaBit Ethernet, Fast Ethernet, and nearly uniform memory architecture, i.e. communication between CPUs is possible by directly reading from or writing to other CPUs' local memory. It turns out that the particle-mesh Ewald method (PME) performs surprisingly well and offers competitive performance unless parallel runs on PC hardware with older network infrastructure are needed. Lipid bilayers of sizes 128, 512 and 2048 lipid molecules were used as the test systems representing typical cases encountered in biomolecular simulations. Our results enable an accurate prediction of computational speed on most current computing systems, both for serial and parallel runs. These results should be helpful in, for example, choosing the most suitable configuration for a small departmental computer cluster.     

GROMACS软件并行计算性能分析

分子动力学模拟是对微观分子原子体系在时间与空间上的运动模拟,是从微观本质上认识体系宏观性质的有力方法.针对如何提升分子动力学并行模拟性能的问题,本文以著名软件GROMACS为例,分析其在分子动力学模拟并行计算方面的实现策略,结合分子动力学模拟关键原理与测试实例,提出MPI+OpenMP并行环境下计算性能的优化策略,为并行计算环境下实现分子动力学模拟的最优化计算性能提供理论和实践参考.对GPU异构并行环境下如何进行MPI、OpenMP、GPU搭配选择以达到性能最优,本文亦给出了一定的理论和实例参考.     

A molecular dynamics study of the BACE1 conformational change from Apo to closed form induced by hydroxyethylamine derived compounds

BACE1 is an aspartyl protease which is a therapeutic target for Alzheimer’s disease (AD) because of its participation in the rate-limiting step in the production of Aβ-peptide, the accumulation of which produces senile plaques and, in turn, the neurodegenerative effects associated with AD. The active site of this protease is composed in part by two aspartic residues (Asp93 and Asp289). Additionally, the catalytic site has been found to be covered by an antiparallel hairpin loop called the flap. The dynamics of this flap are fundamental to the catalytic function of the enzyme. When BACE1 is inactive (Apo), the flap adopts an open conformation, allowing a substrate or inhibitor to access the active site. Subsequent interaction with the ligand induces flap closure and the stabilization of the macromolecular complex. Further, the protonation state of the aspartic dyad is affected by the chemical nature of the species entering the active site, so that appropriate selection of protonation states for the ligand and the catalytic residues will permit the elucidation of the inhibitory pathway for BACE1. In the present study, comparative analysis of different combinations of protonation states for the BACE1-hydroxyethylamine (HEA) system is reported. HEAs are potent inhibitors of BACE1 with favorable pharmacological and kinetic properties, as well as oral bioavailability. The results of Molecular Dynamics (MD) simulations and population density calculations using 8 different parameters demonstrate that the LnAsp289 configuration (HEA with a neutral amine and the Asp289 residue protonated) is the only one which permits the expected conformational change in BACE1, from apo to closed form, after flap closure. Additionally, differences in their capacities to establish and maintain interactions with residues such as Asp93, Gly95, Thr133, Asp289, Gly291, and Asn294 during this step allow differentiation among the inhibitory activities of the HEAs. The results and methodology here reported will serve to elucidate the inhibitory pathway of other families of compounds that act as BACE1 inhibitors, as well as the design of better leader compounds for the treatment of AD.