GPU-accelerated molecular dynamics simulation of solid covalent crystals

Graphics processing unit (GPU) is becoming a powerful computational tool in science and engineering. In this paper, different from previous molecular dynamics (MD) simulation with pair potentials and many-body potentials, two MD simulation algorithms implemented on a single GPU are presented to describe a special category of many-body potentials – bond order potentials used frequently in solid covalent materials, such as the Tersoff potentials for silicon crystals. The simulation results reveal that the performance of GPU implementations is apparently superior to their CPU counterpart. Furthermore, the proposed algorithms are generalised, transferable and scalable, and can be extended to the simulations with general many-body interactions such as Stillinger–Weber potential and so on.     

The multiple forms of bovine seminal ribonuclease: Structure and stability of a C-terminal swapped dimer

Bovine seminal ribonuclease (BS-RNase) acquires an interesting anti-tumor activity associated with the swapping on the N-terminal. The first direct experimental evidence on the formation of a C-terminal swapped dimer (C-dimer) obtained from the monomeric derivative of BS-RNase, although under non-native conditions, is here reported. The X-ray model of this dimer reveals a quaternary structure different from that of the C-dimer of RNase A, due to the presence of three mutations in the hinge peptide 111–116. The mutations increase the hinge peptide flexibility and decrease the stability of the C-dimer against dissociation. The biological implications of the structural data are also discussed.     

Immobilization of Bacillus licheniformis L-Arabinose Isomerase for Semi-Continuous L-Ribulose Production

Bacillus licheniformis L-arabinose isomerase (BLAI) with a broad pH range, high substrate specificity, and high catalytic efficiency for L-arabinose was immobilized on various supports. Eupergit C, activated-carboxymethylcellulose, CNBr-activated agarose, chitosan, and alginate were tested as supports, and Eupergit C was selected as the most effective. After determination of the optimum enzyme concentration, the effects of pH and temperature were investigated using a response surface methodology. The immobilized BLAI enzyme retained 86.4% of the activity of the free enzyme. The optimal pH for the immobilized BLAI was 8.0, and immobilization improved the optimal temperature from 50 °C (free enzyme) to a range between 55 and 65 °C. The half life improved from 2 at 50 °C to 212 h at 55 °C following immobilization. The immobilized BLAI was used for semi-continuous production of L-ribulose. After 8 batch cycles, 95.1% of the BLAI activity was retained. This simple immobilization procedure and the high stability of the final immobilized BLAI on Eupergit C provide a promising solution for large-scale production of L-ribulose from an inexpensive L-arabinose precursor.     

Low‐Overpotential Electrocatalytic Water Splitting with Noble‐Metal‐Free Nanoparticles Supported in a sp3 N‐Rich Flexible COF

Covalent organic frameworks (COFs) are crystalline organic polymers with tunable structures. Here, a COF is prepared using building units with highly flexible tetrahedral sp³ nitrogens. This flexibility gives rise to structural changes which generate mesopores capable of confining very small (<2 nm sized) non‐noble‐metal‐based nanoparticles (NPs). This nanocomposite shows exceptional activity toward the oxygen‐evolution reaction from alkaline water with an overpotential of 258 mV at a current density of 10 mA cm^?2. The overpotential observed in the COF‐nanoparticle system is the best in class, and is close to the current record of ≈200 mV for any noble‐metal‐free electrocatalytic water splitting system—the Fe–Co–Ni metal‐oxide‐film system. Also, it possesses outstanding kinetics (Tafel slope of 38.9 mV dec^?1) for the reaction. The COF is able to stabilize such small‐sized NP in the absence of any capping agent because of the COF–Ni(OH)₂ interactions arising from the N‐rich backbone of the COF. Density‐functional‐theory modeling of the interaction between the hexagonal Ni(OH)₂ nanosheets and the COF shows that in the most favorable configuration the Ni(OH)₂ nanosheets are sandwiched between the sp³ nitrogens of the adjacent COF layers and this can be crucial to maximizing their synergistic interactions.     

Covalent vs. Non‐Covalent Inhibition: Tackling Drug Resistance in EGFR – A Thorough Dynamic Perspective

A persistent challenge in the treatment of non‐small cell lung cancer (NSCLC) with EGFR is the emergence of drug‐resistant caused by somatic mutations. The EGFR L858R/T790 M double mutant (EGFR^DM) was found to be the most alarming variant. Despite the development of a wide range of inhibitors, none of them could inhibit EGFR^DM effectively. Recently, 11h and 45a , have been found to be potent inhibitors against EGFR^DM through two distinctive mechanisms, non‐covalent and covalent binding, respectively. However, the structural and dynamic implications of the two modes of inhibitions remain unexplored. Herein, two molecular dynamics simulation protocols, coupled with free‐energy calculations, were applied to gain insight into the atomistic nature of each binding mode. The comparative analysis confirmed that there is a significant difference in the binding free energy between 11h and 45a (ΔΔG_bind=?21.17 kcal/mol). The main binding force that governs the binding of both inhibitors is vdW, with a higher contribution for 45a . Two residues ARG841 and THR854 were found to have curtailed role in the binding of 45a to EGFR^DM by stabilizing its flexible alcohol chain. The 45a binding to EGFR^DM induces structural rearrangement in the active site to allow easier accessibility of 45a to target residue CYS797. The findings of this work can substantially shed light on new strategies for developing novel classes of covalent and non‐covalent inhibitors with increased specificity and potency.     

Catechol-based inhibitors of bacterial urease

Targeted covalent inhibitors of urease were developed on the basis of the catechol structure. Forty amide and ester derivatives of 3,4-dihydroxyphenylacetic acid, caffeic acid, ferulic acid and gallic acid were obtained and screened against Sporosarcinia pasteurii urease. The most active compound, namely propargyl ester of 3,4-dihydroxyphenylacetic acid exhibited IC₅₀?=?518?nM and$k_{\mathit{inact}}/K_i$?=?1379?M^?1?s^?1. Inhibitory activity of this compound was better and toxicity lower than those obtained for the starting compound – catechol. The molecular modelling studies revealed a mode of binding consistent with structure-activity relationships.     

A novel hydrophilic support, CoFoam, for enzyme immobilization

Lipases from different sources (porcine pancreas, Mucor miehei and Candida antarctica B) were covalently immobilized on a hydrophilic polyurethane composite (CoFoam). Their hydrolytic activities assayed with tributyrin were 0.55, 2.1 and 447 U g(-1), respectively. The activity of the C. antarctica B composite in the synthesis of methyl oleate in hexane was 8.8 U g(-1) compared to 60.6 U g(-1) for commercial Novozyme 435. The advantage of the CoFoam composite lies in the low pressure drop in a packed-bed reactor at fairly large flow rates. For example, at flow rates of 10-12 l min(-1), the pressure drop over 15 cm is typically 3 kPa.     

Dynamics and DNA substrate recognition by the catalytic domain of lambda integrase

Bacteriophage lambda integrase (lambda-Int) is the prototypical member of a large family of enzymes that catalyze site-specific DNA recombination via the formation of a Holliday junction intermediate. DNA strand cleavage by lambda-Int is mediated by nucleophilic attack on the scissile phosphate by a conserved tyrosine residue, forming an intermediate with the enzyme covalently attached to the 3'-end of the cleaved strand via a phosphotyrosine linkage. The crystal structure of the catalytic domain of lambda-Int (C170) obtained in the absence of DNA revealed the tyrosine nucleophile at the protein's C terminus to be located on a beta-hairpin far from the other conserved catalytic residues and adjacent to a disordered loop. This observation suggested that a conformational change in the C terminus of the protein was required to generate the active site in cis, or alternatively, that the active site could be completed in trans by donation of the tyrosine nucleophile from a neighboring molecule in the recombining synapse. We used NMR spectroscopy together with limited proteolysis to examine the dynamics of the lambda-Int catalytic domain in the presence and absence of DNA half-site substrates with the goal of characterizing the expected conformational change. Although the C terminus is indeed flexible in the absence of DNA, we find that conformational changes in the tyrosine-containing beta-hairpin are not coupled to DNA binding. To gain structural insights into C170/DNA complexes, we took advantage of mechanistic conservation with Cre and Flp recombinases to model C170 in half-site and tetrameric Holliday junction complexes. Although the models do not reveal the nature of the conformational change required for cis cleavage, they are consistent with much of the available experimental data and provide new insights into the how trans complementation could be accommodated.     

Optimal covalent immobilization of glucose oxidase-containing liposomes for highly stable biocatalyst in bioreactor

The glucose oxidase-containing liposomes (GOL) were prepared by entrapping glucose oxidase (GO) in the liposomes composed of phosphatidylcholine (PC), dimyristoyl L-alpha-phosphatidylethanolamine (DMPE), and cholesterol (Chol) and then covalently immobilized in the glutaraldehyde-activated chitosan gel beads. The immobilized GOL gel beads (IGOL) were characterized to obtain a highly stable biocatalyst applicable to bioreactor. At first, the glutaraldehyde concentration used in the gel beads activation as well as the immobilizing temperature and time were optimized to enhance the immobilization yield of the GOL to the highest extent. The liposome membrane composition and liposome size were then optimized to obtain the greatest possible immobilization yield of the GOL, the highest possible activity efficiency of the IGOL, and the lowest possible leakage of the entrapped GO during the GOL immobilization. As a result, the optimal immobilization conditions were found to be as follows: the liposome composition, PC/DMPE/Chol = 65/5/30 (molar percentage); the liposome size, 100 nm; the glutaraldehyde concentration, 2% (w/v); the immobilizing temperature, 4 degrees C; and the immobilizing time, 10 h. Furthermore, the optimal IGOL prepared were characterized by its rapidly increasing effective GO activity to the externally added substrate (glucose) with increasing temperature from 20 to 40 degrees C, and also by its high stability at 40 degrees C against not only the thermal denaturation in a long-term (7 days) incubation but also the bubbling stress in a bubble column. Finally, compared to the conventionally immobilized glucose oxidase (IGO), the higher operational stability of the optimal IGOL was verified by using it either repeatedly (4 times) or for a long time (7 days) to catalyze the glucose oxidation in a small-scale airlift bioreactor.     

Extraction of membrane proteins from a living cell surface using the atomic force microscope and covalent crosslinkers

The force curve mode of the atomic force microscope (AFM) was applied to extract intrinsic membrane proteins from the surface of live cells using AFM tips modified by amino reactive bifunctional covalent crosslinkers. The modified AFM tips were individually brought into brief contact with the living cell surface to form covalent bonds with cell surface molecules. The force curves recorded during the detachment process from the cell surface were often characterized by an extension of a few hundred nanometers followed mostly by a single step jump to the zero force level. Collection and analysis of the final rupture force revealed that the most frequent force values (of the force) were in the range of 0.4–0.6 nN. The observed rupture force most likely represented extraction events of intrinsic membrane proteins from the cell membrane because the rupture force of a covalent crosslinking system was expected to be significantly larger than 1.0 nN, and the separation force of noncovalent ligand-receptor pairs to be less than 0.2 nN, under similar experimental conditions. The transfer of cell surface proteins to the AFM tip was verified by recording characteristic force curves of protein stretching between the AFM tips used on the cell surface and a silicon surface modified with amino reactive bifunctional crosslinkers. This method will be a useful addition to bionanotechnological research for the application of AFM.