PyRETIS 3: Conquering rare and slow events without boundaries

We present and discuss the advancements made in PyRETIS 3, the third instalment of our Python library for an efficient and user-friendly rare event simulation, focused to execute molecular simulations with replica exchange transition interface sampling (RETIS) and its variations. Apart from a general rewiring of the internal code towards a more modular structure, several recently developed sampling strategies have been implemented. These include recently developed Monte Carlo moves to increase path decorrelation and convergence rate, and new ensemble definitions to handle the challenges of long-lived metastable states and transitions with unbounded reactant and product states. Additionally, the post-analysis software PyVisa is now embedded in the main code, allowing fast use of machine-learning algorithms for clustering and visualising collective variables in the simulation data.

Highly parallelizable path sampling with minimal rejections using asynchronous replica exchange and infinite swaps

Capturing rare yet pivotal events poses a significant challenge for molecular simulations. Path sampling provides a unique approach to tackle this issue without altering the potential energy landscape or dynamics, enabling recovery of both thermodynamic and kinetic information. However, despite its exponential acceleration compared to standard molecular dynamics, generating numerous trajectories can still require a long time. By harnessing our recent algorithmic innovations-particularly subtrajectory moves with high acceptance, coupled with asynchronous replica exchange featuring infinite swaps-we establish a highly parallelizable and rapidly converging path sampling protocol, compatible with diverse high-performance computing architectures. We demonstrate our approach on the liquid-vapor phase transition in superheated water, the unfolding of the chignolin protein, and water dissociation. The latter, performed at the ab initio level, achieves comparable statistical accuracy within days, in contrast to a previous study requiring over a year.

Path sampling with memory reduction and replica exchange to reach long permeation timescales

Assessing kinetics in biological processes with molecular dynamics simulations remains a computational and conceptual challenge, given the large time and length scales involved. For kinetic transport of biochemical compounds or drug molecules, the permeability through the phospholipid membranes is a key kinetic property, but long timescales are hindering the accurate computation. Technological advances in high-performance computing therefore need to be accompanied by theoretical and methodological developments. In this contribution, the replica exchange transition interface sampling (RETIS) methodology is shown to give perspective toward observing longer permeation pathways. It is first reviewed how RETIS, a path-sampling methodology that gives in principle exact kinetics, can be used to compute membrane permeability. Next, recent and current developments in three RETIS aspects are discussed: several new Monte Carlo moves in the path-sampling algorithm, memory reduction by reducing pathlengths, and exploitation of parallel computing with CPU-imbalanced replicas. Finally, the memory reduction presenting a new replica exchange implementation, coined REPPTIS, is showcased with a permeant needing to pass a membrane with two permeation channels, either representing an entropic or energetic barrier. The REPPTIS results showed clearly that inclusion of some memory and enhancing ergodic sampling via replica exchange moves are both necessary to obtain correct permeability estimates. In an additional example, ibuprofen permeation through a dipalmitoylphosphatidylcholine membrane was modeled. REPPTIS succeeded in estimating the permeability of this amphiphilic drug molecule with metastable states along the permeation pathway. In conclusion, the presented methodological advances allow for deeper insight into membrane biophysics even if the pathways are slow, as RETIS and REPPTIS push the permeability calculations to longer timescales.

Enhanced path sampling using subtrajectory Monte Carlo moves

Path sampling allows the study of rare events, such as chemical reactions, nucleation, and protein folding, via a Monte Carlo (MC) exploration in path space. Instead of configuration points, this method samples short molecular dynamics (MD) trajectories with specific start- and end-conditions. As in configuration MC, its efficiency highly depends on the types of MC moves. Since the last two decades, the central MC move for path sampling has been the so-called shooting move in which a perturbed phase point of the old path is propagated backward and forward in time to generate a new path. Recently, we proposed the subtrajectory moves, stone-skipping (SS) and web-throwing, that are demonstrably more efficient. However, the one-step crossing requirement makes them somewhat more difficult to implement in combination with external MD programs or when the order parameter determination is expensive. In this article, we present strategies to address the issue. The most generic solution is a new member of subtrajectory moves, wire fencing (WF), that is less thrifty than the SS but more versatile. This makes it easier to link path sampling codes with external MD packages and provides a practical solution for cases where the calculation of the order parameter is expensive or not a simple function of geometry. We demonstrate the WF move in a double-well Langevin model, a thin film breaking transition based on classical force fields, and a smaller ruthenium redox reaction at the ab initio level in which the order parameter explicitly depends on the electron density.

[Interpretation for the group standards in the handling of accidental leak coronavirus disease 2019 sample]

Biosafety is an important guarantee of the new coronavirus laboratory test. The accident treatment of sample overflow and sprinkle is a necessary part of the emergency plan for testing activities. Beijing Preventive Medicine Association coordinated biosafety experts of COVID-19 laboratories from Beijing CDC, to write up "The standard for handling of accidents of corona virus disease 2019 sample (T/BPMA 0005-2020)" . The group standard was based on the guidelines of China and WHO, and combined with the practical experience of COVID-19 epidemic and the principle of "scientific, normative, applicable and feasible" . Through all kinds of risk Assessment, it included the spillover of samples caused by the packing of COVID-19 (highly pathogenic) samples, the overflow and sprinkle in the laboratory during the detection operation, and the spillage accident occurred during the transfer of samples in the same building. The standard could guide and standardize the handling methods of accidental overflow and sprinkle that may occur in the SARS-CoV-2 testing laboratories in the city.

[Study on the super-antigen genes of group A Streptococcus pyogenes strains isolated from patients with scarlet fever and pharyngeal infection, in Beijing, 2015-2017]

Objective: To analyze the characteristics of super-antigen (SAg) of group A Streptococcus pyogenes (GAS), isolated from patients with scarlet fever or pharyngeal infections in Beijing between 2015-2017. Methods: Throat swab specimens from patients with scarlet fever or pharyngeal infections were collected and tested for GAS. Eleven currently known SAg genes including SpeA, speC, speG, speH, speI, speJ, speK, speL, speM, smeZ and ssa were tested by real-time PCR while M protein genes (emm genes) were amplified and sequenced by PCR. Results: A total of 377 GAS were isolated from 6 801 throat swab specimens, with the positive rate as 5.5%. There were obvious changes noticed among speC, speG, speH and speK in three years. A total of 45 SAg genes profiles were observed, according to the SAgs inclusion. There were significant differences appeared in the frequencies among two of the highest SAg genes profiles between emm1 and emm12 strains (χ(2)=38.196, P<0.001; χ(2)=72.310, P<0.001). There also appeared significant differences in the frequencies of speA, speH, speI and speJ between emm1 and emm12 strains (χ(2)=146.154, P<0.001; χ(2)=52.31, P<0.001; χ(2)=58.43, P<0.001; χ(2)=144.70, P<0.001). Conclusions: Obvious changes were noticed among SAg genes including speC, speG, speH and speK from patients with scarlet fever or pharyngeal infections in Beijing between 2015-2017. SAg genes including speA, speH, speI and speJ appeared to be associated with the emm 1 and emm 12 strains. More kinds of SAg genes profiles were isolated form GAS but with no significant differences seen in the main SAg genes profiles, during the epidemic period.

[Open mitral valvuloplasty in 40 cases]

From Dec. 1987 to June 1991, open mitral valvuloplasty was done in 40 patients. Among them rheumatic valve disease was observed in 35 cases. The early postoperative result was excellent. We consider that the requisites for this operation are noncalcification of anterior mitral leaflet with sufficient size and flexibility, without significant fusion and shortening of subvalvar structures, and readiness of valve replacement in case of its failure.

[The double-chambered right ventricle. Report of 8 cases]

From April 1985 to 1988, of 264 patients under went repair of congenital heart defect, 8 patients with double-chambered right ventricle (DCRV) were successfully managed. All patients survived well. The catheterization and angiocardiography in the right ventricle(RV) play an important role in the preoperative diagnosis. To avoid missing and incorrect management the surgeon must keep the possibility of DCRV in mind during repair of a VSD. A dimple in the anterior wall of the RV during systole represents the underlying attachment of the aberrant muscular bundle. The right ventricular outflow tract must be explored when through a right atriotomy to repair the VSD and the tricuspid valve must be examined when through a right ventriculotomy.