Nonadiabatic dynamics simulations have rapidly become an indispensable tool for understanding ultrafast photochemical processes in complex systems. Here, we present our recently developed on-the-fly nonadiabatic dynamics package, JADE, which allows researchers to perform nonadiabatic excited-state dynamics simulations of polyatomic systems at an all-atomic level. The nonadiabatic dynamics is based on Tully's surface-hopping approach. Currently, several electronic structure methods (CIS, TDHF, TDDFT(RPA/TDA), and ADC(2)) are supported, especially TDDFT, aiming at performing nonadiabatic dynamics on medium- to large-sized molecules. The JADE package has been interfaced with several quantum chemistry codes, including Turbomole, Gaussian, and Gamess (US). To consider environmental effects, the Langevin dynamics was introduced as an easy-to-use scheme into the standard surface-hopping dynamics. The JADE package is mainly written in Fortran for greater numerical performance and Python for flexible interface construction, with the intent of providing open-source, easy-to-use, well-modularized, and intuitive software in the field of simulations of photochemical and photophysical processes. To illustrate the possible applications of the JADE package, we present a few applications of excited-state dynamics for various polyatomic systems, such as the methaniminium cation, fullerene (C20), p-dimethylaminobenzonitrile (DMABN) and its primary amino derivative aminobenzonitrile (ABN), and 10-hydroxybenzo[h]quinoline (10-HBQ).

The anisotropic effects and short-range quantum effects are essential characters in the formation of halogen bonds. Since there are an array of applications of halogen bonds and much difficulty in modeling them in classical force fields, the current research reports solely the polarizable ellipsoidal force field (PEff) for halogen bonds. The anisotropic charge distribution was represented with the combination of a negative charged sphere and a positively charged ellipsoid. The polarization energy was incorporated by the induced dipole model. The resulting force field is "physically motivated," which includes separate, explicit terms to account for the electrostatic, repulsion/dispersion, and polarization interaction. Furthermore, it is largely compatible with existing, standard simulation packages. The fitted parameters are transferable and compatible with the general AMBER force field. This PEff model could correctly reproduces the potential energy surface of halogen bonds at MP2 level. Finally, the prediction of the halogen bond properties of human Cathepsin L (hcatL) has been found to be in excellent qualitative agreement with the cocrystal structures.

Cellulase Cel7A from Trichoderma reesei is one of the most abundant and effective cellulases. Structural studies have established that Cel7A is a retaining glycosidase and it can processively hydrolyze cellobiose units from the reducing end of a cellulose chain. Here, to elucidate the mechanism of enzymatic catalysis of cellulase Cel7A, we carried out a multisized level theoretical study by performing MD, QM, and QM/MM calculations. At the accurate level of theory, we showed the mechanism details of the catalytic cycle, which involves the configuration inversion of the anomeric center twice: the first results in the glycosidic bond cleavage and the formation of covalent glycosyl-enzyme intermediate, and the second restores the anomeric carbon to its original configuration. Calculated results have provided detailed structural and energetic information about these two processes, both of which proceed according to a S(N)2-type-like mechanism via loose transition state structures. It is clearly indicated that the glycosidic bond hydrolysis involves the formation of a covalent glycosyl-enzyme intermediate, which has been identified as the minimum on the potential energy surface. At the catalytic active region, hydrogen bond interactions exist throughout the whole process of the catalytic cycle, which are of special importance for stabilizing the glycosyl-enzyme intermediate. The present results provide a clear paradigm of the mechanisms of general glycosidases, which hydrolyze the glycosidic bonds with net retention of the anomeric configuration.

The current COVID-19 pandemic has already had a devastating impact across the world. SARS-CoV-2 (the virus causing COVID-19) is known to use its surface spike (S) protein's receptor binding domain (RBD) to interact with the angiotensin-converting enzyme 2 (ACE2) receptor expressed on many human cell types. The RBD–ACE2 interaction is a crucial step to mediate the host cell entry of SARS-CoV-2. Recent studies indicate that the ACE2 interaction with the SARS-CoV-2 S protein has higher affinity than its binding with the structurally identical S protein of SARS-CoV-1, the virus causing the 2002-2004 SARS epidemic. However, the biophysical mechanism behind such binding affinity difference is unclear. This study utilizes a combined single-molecule force spectroscopy and steered molecular dynamics (SMD) simulation approach to quantify the specific interactions between CoV-2 or CoV-1 RBD and ACE2. Depending on the loading rates, the unbinding forces between CoV-2 RBD and ACE2 range from 70 to 110 pN, and are 30-50% higher than those of CoV-1 RBD and ACE2 under similar loading rates. SMD results indicate that CoV-2 RBD interacts with the N-linked glycan on Asn90 of ACE2. This interaction is mostly absent in the CoV-1 RBD–ACE2 complex. During the SMD simulations, the extra RBD-N-glycan interaction contributes to a greater force and prolonged interaction lifetime. The observation is confirmed by our experimental force spectroscopy study. After the removal of N-linked glycans on ACE2, its mechanical binding strength with CoV-2 RBD decreases to a similar level of the CoV-1 RBD–ACE2 interaction. Together, the study uncovers the mechanism behind the difference in ACE2 binding between SARS-CoV-2 and SARS-CoV-1, and could aid in the development of new strategies to block SARS-CoV-2 entry.

By employing ab initio quantum mechanical/molecular mechanical (QM/MM) and molecular dynamics (MD) simulations, we have provided further evidence against the previously proposed hydroperoxylation or hydroxylation mechanism of hydroxyethylphosphonate dioxygenase (HEPD). HEPD employs an interesting catalytic cycle based on concatenated bifurcations. The first bifurcation is based on the abstraction of hydrogen atoms from the substrate, which leads to a distal or proximal hydroperoxo species (Fe-OOH or Fe-(OH)O). The second and the third bifurcations refer to the carbon-carbon bond cleavage reaction. And this is achieved through a tridentate intermediate, or employing a proton-shuttle assisted mechanism, in which the residue Glu(176) or the Fe(IV)=O group serves as a general base. The reaction directions seem to be tunable and show significant environment dependence. This mechanism can provide a comprehensive interpretation for the seemingly contradicting experimental evidences and provide insight into the development of biochemistry and material sciences.