On the induced‐fit mechanism of substrate‐enzyme binding structures of nylon‐oligomer hydrolase

We present a detailed computational investigation of the induced‐fit motion in a nylon‐oligomer hydrolase (NylB) upon substrate binding. To this aim, we resort on the recently introduced parallel cascade selection molecular dynamics approach, allowing for an accelerated access to the set of conformational changes from an open‐ to a closed‐state structure to form the enzyme‐substrate complex in a specific induce‐fit mechanism. The structural investigation is quantitatively complemented by free energy analyses within the umbrella sampling algorithm accompanied by weighted histogram analysis. We find that the stabilization free energy is about 1.4 kcal/mol, whereas the highest free energy barrier to be overcome is about 2.3 kcal/mol. Conversely, the energetic contribution for the substrate binding is about 20 kcal/mol, as estimated from Generalized Born/Surface Area. This means that the open‐close induced‐fit motion could occur frequently once the substrate binds to the open state of NylB. © 2014 Wiley Periodicals, Inc.     

Simulations of remote mutants of dihydrofolate reductase reveal the nature of a network of residues coupled to hydride transfer

Recent experimental and theoretical studies have proposed that enzymes involve networks of coupled residues throughout the protein that participate in motions accompanying chemical barrier crossing. Here, we have examined portions of a proposed network in dihydrofolate reductase (DHFR) using quantum mechanics/molecular mechanics simulations. The simulations use a hybrid quantum mechanics‐molecular mechanics approach with a recently developed semiempirical AM1‐SRP Hamiltonian that provides accurate results for this reaction. The simulations reproduce experimentally determined catalytic rates for the wild type and distant mutants of E. coli DHFR, underscoring the accuracy of the simulation protocol. Additionally, the simulations provide detailed insight into how residues remote from the active site affect the catalyzed chemistry, through changes in the thermally averaged properties along the reaction coordinate. The mutations do not greatly affect the structure of the transition state near the bond activation, but we observe differences somewhat removed from the point of C? H cleavage that affect the rate. The mutations have global effects on the thermally averaged structure that propagate throughout the enzyme and the current simulations highlight several interactions that appear to be particularly important. © 2014 Wiley Periodicals, Inc.     

19-取代雄烯二酮芳构酶抑制剂的定量构效关系研究

用分子动力学和量子化学AM1经验方法对系列19-位取代的雄烯二酮芳构酶抑制 剂的定量构效关系进行了研究，结果表明该位置的取代基必须符合P450芳构酶活性 位点的立体效应和疏水作用的要求。此类化合物作为较好的电子供体与芳构酶结合 ，给电子的主要基团为21-位含S取代基。这些研究结果为设计活性更高的芳构酶抑 制剂提供了有益的参考信息。     

葡萄糖水溶液氢键结构和动力分析

利用分子动力学模拟方法研究了不同浓度下葡萄糖水溶液的氢键结构和氢键生存周期. 分析了参与i个氢键(分子内、分子间、所有类型)的葡萄糖分子和水分子的百分比分布. 研究发现存在一个特征数N, 参与N个氢键的分子的比例最高, 当iN时, 参与i个氢键的分子的比例随着浓度的增加而减小. 还分析了不同类型氢键(葡萄糖分子内、葡萄糖分子间、水分子间、葡萄糖分子与水分子间、所有类型)的连续和截断自相关函数, 并计算了对应的氢键生存周期.     

Identification,characterization and evolution of non-local quasi-Lagrangian structures in turbulence

The recent progress on non-local Lagrangian and quasi-Lagrangian structures in turbulence is reviewed. The quasi-Lagrangian structures, e.g., vortex surfaces in vis-cous flow, gas-liquid interfaces in multi-phase flow, and flame fronts in premixed combustion, can show essential Lagrangian following properties, but they are able to have topological changes in the temporal evolution. In addition, they can represent or influence the turbulent flow field. The challenges for the investigation of the non-local structures include their identification, characterization, and evolution. The improving understanding of the quasi-Lagrangian struc-tures is expected to be helpful to elucidate crucial dynamics and develop structure-based predictive models in turbulence.     

Immersed boundary method for unsteady kinetic model equations

Predicting unsteady flows and aerodynamic forces for large displacement motion of microstructures requires transient solution of Boltzmann equation with moving boundaries. For the inclusion of moving complex boundaries for these problems, three immersed boundary method flux formulations (interpolation, relaxation, and interrelaxation) are presented. These formulations are implemented in a 2‐D finite volume method solver for ellipsoidal‐statistical (ES)‐Bhatnagar‐Gross‐Krook (BGK) equations using unstructured meshes. For the verification, a transient analytical solution for free molecular 1‐D flow is derived, and results are compared with the immersed boundary (IB)‐ES‐BGK methods. In 2‐D, methods are verified with the conformal, non‐moving finite volume method, and it is shown that the interrelaxation flux formulation gives an error less than the interpolation and relaxation methods for a given mesh size. Furthermore, formulations applied to a thermally induced flow for a heated beam near a cold substrate show that interrelaxation formulation gives more accurate solution in terms of heat flux. As a 2‐D unsteady application, IB/ES‐BGK methods are used to determine flow properties and damping forces for impulsive motion of microbeam due to high inertial forces. IB/ES‐BGK methods are compared with Navier–Stokes solution at low Knudsen numbers, and it is shown that velocity slip in the transitional rarefied regime reduces the unsteady damping force. Copyright © 2015 John Wiley & Sons, Ltd.     

Strain rate effects on compressive behavior of covalently bonded CNT networks

In this study, strain rate effects on the compressive mechanical properties of randomly structured carbon nanotube (CNT) networks were examined. For this purpose, three-dimensional atomistic models of CNT networks with covalently-bonded junctions were generated. After that, molecular dynamics (MD) simulations of compressive loading were performed at five different strain rates to investigate the basic deformation characteristic mechanisms of CNT networks and determine the effect of strain rate on stress–strain curves. The simulation results showed that the strain rate of compressive loading increases, so that a higher resistance of specimens to deformation is observed. Furthermore, the local deformation characteristics of CNT segments, which are mainly driven by bending and buckling modes, and their prevalence are strongly affected by the deformation rate. It was also observed that CNT networks have superior features to metal foams such as metal matrix syntactic foams (MMSFs) and porous sintered fiber metals (PSFMs) in terms of energy absorbing capabilities.     

Effect of Vacancy Defects on the Young's Modulus and Fracture Strength of Graphene: A Molecular Dynamics Study

The Young's modulus of graphene with various rectangular and circular vacancy defects is investigated by molecular dynamics simulation. By comparing with the results calculated from an effective spring model, it is demonstrated that the Young's modulus of graphene is largely correlated to the size of vacancy defects perpendicular to the stretching direction. And a linear reduction of Young's modulus with the increasing concentration of mono‐atomic‐vacancy defects (i.e., the slope of ?0.03) is also observed. The fracture behavior of graphene, including the fracture strength, crack initiation and propagation are then studied by the molecular dynamics simulation, the effective spring model, and the quantized fracture mechanics. The blunting effect of vacancy edges is demonstrated, and the characterized crack tip radius of 4.44 Å is observed.     

Ion beam mixing of α-MoSi2

Abstract

Using molecular dynamics simulation, we study the mixing of an α-MoSi₂ crystal by 1 keV Ar ions. The observed order of magnitude of the mixing is compatible with a spike model of ion beam mixing. The influence of the target surface and of chemical effects are discussed.     

Novel Dual‐Site‐Binding Neuraminidase Inhibitor from Virtual Screening by Pharmacophore and Molecular Dynamics Methods

Neuraminidase is a significant anti‐influenza target that plays crucial role in virus replication cycle. The discovery of 150‐cavity in Group‐1 neuraminidase provides us a novel mentality of designing inhibitor which can bind with both conserved site and 150‐cavity. In order to discover novel dual‐site‐binding inhibitors, a 3D chemical‐feature‐based pharmacophore model was established to cover dual‐site in neuraminidase. The dual‐site‐binding model was consistent in predicting the binding conformation of Group‐1 neuraminidase inhibitor and applied for virtual screening of Specs database. Compound 4 (ZINC05790048) that aligned well to the model was selected after multiple filtrations for molecular dynamics simulations, indicating improved binding energy with neuraminidase. It can sever as the lead compound for a novel series of inhibitors.