Machine learning for reparameterization of four-site water models: TIP4P-BG and TIP4P-BGT.

Parameterizing an effective water model is a challenging issue because of the difficulty in maintaining a comprehensive balance among the diverse physical properties of water with a limited number of parameters. The advancement in machine learning provides a promising path to search for a reliable set of parameters. Based on the TIP4P water model, hence, about 6000 molecular dynamics (MD) simulations for pure water at 1 atm and in the range of 273-373 K are conducted here as the training data. The back-propagation (BP) neural network is then utilized to construct an efficient mapping between the model parameters and four crucial physical properties of water, including the density, vaporization enthalpy, self-diffusion coefficient and viscosity. Without additional time-consuming MD simulations, this mapping operation could result in sufficient and accurate data for high-population genetic algorithm (GA) to optimize the model parameters as much as possible. Based on the proposed parameterizing strategy, TIP4P-BG (a conventional four-site water model) and TIP4P-BGT (an advanced model with temperature-dependent parameters) are established. Both the water models exhibit excellent performance with a reasonable balance among the four crucial physical properties. The relevant mean absolute percentage errors are 3.53% and 3.08%, respectively. Further calculations on the temperature of maximum density, isothermal compressibility, thermal expansion coefficient, radial distribution function and surface tension are also performed and the resulting values are in good agreement with the experimental values. Through this water modeling example, the potential of the proposed data-driven machine learning procedure has been demonstrated for parameterizing a MD-based material model.

Water diffusion inside carbon nanotubes: mutual effects of surface and confinement.

The mutual effects of two crucial features of carbon nanotubes (CNTs) (surface and confinement) on the temperature-dependent water diffusion are studied through molecular dynamics simulations. A two-stage diffusion mechanism is detected in the CNTs of diameter smaller than 12.2 Å, which becomes obscure as the temperature increases. This peculiar phenomenon can be ascribed to the cooperation of the small confinement and the periodic surface. The diffusion coefficient of the confined water exhibits a nonmonotonic dependence on the confinement size and an unexpected increase inside the large CNTs (compared to that of bulk water). These anomalous behaviors can be attributed to the competition of the smooth surface and the small confinement. Considering the mutual effects, an empirical formula is proposed on the basis of two groups of numerical examples, whose results indicate that the confinement effect will dominate over the surface effect until the CNT diameter increases up to ~16 Å, whereas thereafter the surface effect becomes dominant and finally both of them vanish gradually.

A simulation study on nanoscale holes generated by gold nanoparticles on negative lipid bilayers.

Understanding the interactions of gold nanoparticles (AuNPs) with cellular compartments, especially cell membranes, is of fundamental importance in obtaining their control in biomedical applications. An effort is made in this paper to investigate the interactions of 2.2 nm core AuNPs with negative model bilayer membranes by coarse-grained (CG) molecular dynamics (MD) simulation. The CG model of lipid bilayer was taken from Marrink et al. ( J. Phys. Chem. B 2004, 108, 750-760 ), whereas the CG AuNPs model was developed on the basis of both atomistic MD simulations and experimental data. It was found that AuNPs functionalized with cationic ligands penetrated into the negative bilayer membranes and generated significant disruptions on bilayers. The lipids surrounding the nanoparticle were highly disordered and the bulk surface of the bilayer exhibits some defective areas. Most importantly, it is observed that a nanoscale hole can be formed and expanded spontaneously on the peripheral regions of the 20 × 20 nm bilayer. The expansion of the hole is on the time scale of hundreds of nanosceonds. The fully expanded hole had a radius of ∼5.5 nm and could transport water molecules at a rate of up to ∼1100 molecule/ns. However holes could not be formed on a larger bilayer (28 × 28 nm). The factors that can eliminate hole formation on the bilayer also include the decrease of cationic lignads on the AuNP, the reduction of negative lipids in the bilayer, the release of bilayer surface tension, the lowering of temperature, and the addition of a high concentration of salt. The results suggest that a hole can only be formed on living cell membranes under extreme conditions.

Unidirectional Self-Driving Liquid Droplet Transport on a Monolayer Graphene-Covered Textured Substrate.

Controllable directional transport of liquid droplets on a functionalized surface has been a challenge in the field of microfluidics because it does not require energy supply, and the physical mechanism of such self-driving transport exhibits extraordinary contribution to fundamental understanding of some biological processes and the design of microfluidic apparatus. In this paper, we report a novel design of a surface microstructure that can realize unidirectional self-driving liquid mercury (Hg) droplet transport on a graphene-covered copper (Cu) substrate with a three-dimensional surface microstructure. We have demonstrated that a liquid Hg droplet spontaneously propagates on a grooved Cu substrate covered by a monolayer graphene without any external force fields. Classical molecular dynamics results provide a profound insight on the self-driving process of Hg droplets. It shows that the Hg droplet undergoes acceleration, deceleration, and return stages successively from the narrow to wide ends of the gradient groove. Intriguingly, Hg droplets can move continuously and unidirectionally on the three-dimensional graphene-covered surface microstructure when they accumulate enough kinetic energy from the gradient groove to break the energy barrier at the step junctions between the two neighboring unit cells. The design of the zigzag textured surface covered by a monolayer graphene artfully uses the facts; (1) the monolayer graphene can effectively reduce the droplet pinning on the textured surface, (2) the hydrophobic graphene layer reduces the friction between Hg droplets and the substrate, and (3) the textured surface can permeably interact with the droplets through the monolayer graphene to achieve a continuous self-driving process. The findings reported here open a door to explore the graphene-covered functional surface to directional transport of liquid droplets and provide an in-depth understanding of the self-driving mechanism for liquid droplets on graphene-covered textured substrates.

Two distinct domains of the beta subunit of Aquifex aeolicus leucyl-tRNA synthetase are involved in tRNA binding as revealed by a three-hybrid selection.

The Aquifex aeolicus alphabeta-LeuRS is the only known heterodimeric class Ia aminoacyl-tRNA synthetase. In this study, we investigated the function of the beta subunit which is believed to bind tRNA(Leu). A yeast three-hybrid system was constructed on the basis of the interaction of the beta subunit with its cognate tRNA(Leu). Then, seven mutated beta subunits exhibiting impaired tRNA binding capacities were selected out from a randomly mutated library. Two mutations were identified in the class Ia-helix-bundle-domain, which might interact with the D-hairpin of the tRNA analogous to other class Ia tRNA:synthetases complexes. The five other mutations were found in the LeuRS-specific C-terminal domain of which the folding is still unknown. tRNA affinity measurements and kinetic analyses performed on the isolated beta subunits and on the co-expressed alphabeta-heterodimers showed for all the mutants an effect in tRNA affinity in the ground state. In addition, an effect on the transition state of the aminoacylation reaction was observed for a 21-residues deletion mutant of the C-terminal end. These results show that the genetic approach of the three hybrid system is widely applicable and is a powerful tool for the investigation of tRNA:synthetase interactions.

[Species-associated differences in foliage-root coupling soil-reinforcement and anti-erosion].

This paper took four kinds of common soil and water conservation plants of the study area, Caragana microphylla, Salix psammophila, Artemisia sphaerocephala and Hippophae rhamnides at ages of 4 as the research object. Thirteen indicators, i.e., single shrub to reduce wind velocity ration, shelterbelt reducing wind velocity ration, community reducing wind velocity ration, taproot tensile strength, representative root constitutive properties, representative root elasticity modulus, lateral root branch tensile strength, accumulative surface area, root-soil interface sheer strength, interface friction coefficient, accumulative root length, root-soil composite cohesive, root-soil composite equivalent friction angle, reflecting the characteristics of windbreak and roots, were chose to evaluate the differences of foliage-root coupling soil-reinforcement and anti-erosion among four kinds of plants by analytic hierarchy process (AHP) under the condition of spring gale and summer rainstorm, respectively. The results showed the anti-erosion index of foliage-root coupling was in the sequence of S. psammophila (0.841) > C. microphylla (0.454) > A. sphaerocephala (-0.466) > H. rhamnides (-0.829) in spring gale, and C. microphylla (0.841) > S. psammophila (0. 474) > A. sphaerocephala (-0.470) > H. rhamnides (-0.844) in summer rainstorm. S. psammophila could be regarded as one of the most important windbreak and anti-erosion species, while C. microphylla could be the most valuable soil and water conservation plant for the study area.

Reduction of mitochondrial tRNALeu(UUR) aminoacylation by some MELAS-associated mutations.

The mitochondrial myopathy, encephalopathy, lactic acidosis and stroke-like episodes syndrome (MELAS) is a rare congenital disorder of mitochondrial DNA. Five single nucleotide substitutions within the human mitochondrial tRNALeu(UUR) gene have been reported to be associated with MELAS. Here, we provide in vitro evidence that the aminoacylation capacities of these five hmtRNALeu(UUR) transcripts are reduced to different extents relative to the wild-type hmtRNALeu(UUR) transcript. A thermal denaturation experiment showed that the A3243G and T3291C mutants, which were the least charged by LeuRS, have fragile structures. In addition, the T3291C mutant can inhibit aminoacylation of the wild-type hmtRNALeu(UUR), indicating that it may act as an inhibitor in the mitochondrial heteroplasmic environment.

Expression of gene encoding GL-7ACA acylase in Escherichia coli.

Glutaryl 7-amino cephalosporanic acid acylase (GL-7ACA acylase) from Pseudomonas sp.130 catalyzes hydrolysis of glutaryl 7-amino cephalosporanic acid to produce 7-amino cephalosporanic acid (7-ACA). 7-ACA is the starting material for the industrial production of most cephalosparonic derivatives. Six plasmids for expression of GL-7ACA acylase were constructed and these recombin ant plasmids presented different expression characteristics in Escherichia coli. The acylase gene from plasmid pKKCA1 was inserted into plasmid pMFT7-5 and the resulting plasmid pMFT7CA1 has higher expression in E.coli. The specific activity of the crude extract of the transformant JM109(DE3)/pMFT7CA1 was near 5 u/g, so the over produced enzyme was easily purified by a single-step anion exchange column chromatography. The enzyme could be purified by immobilized ion affinity chromatography after fused by 6xHis in the N-terminal of its alpha-subunit. Because plasmid pSMLCA1 brings tc(R) and p15A origin, it is special useful plasmid in fermentation. Two secretory expression plasmids, pSUCA1S and pETCA1pelB, could secrete the acylase to periplasmic space of bacteria. The whole cells containing the secretory expression plasmid may be used for production of 7-ACA directly.

Unique residues crucial for optimal editing in yeast cytoplasmic Leucyl-tRNA synthetase are revealed by using a novel knockout yeast strain.

Leucyl-tRNA synthetase (LeuRS) contains an editing domain that discriminates leucine from noncognate amino acids to ensure translational fidelity. In this study, a knock-out strain for Saccharomyces cerevisiae LeuRS was constructed to analyze in vivo the tRNA aminoacylation properties of S. cerevisiae and human cytoplasmic LeuRSs. The activities of several editing-defective mutants of ycLeuRS were determined in vitro and compared with those obtained in vivo in a complementation assay performed in the knock-out strain. The editing activities of these mutants were analyzed in the presence of either norvaline, a leucine analogue, or AN2690, a specific inhibitor that targets the editing active site. In general, the in vivo data are consistent with those obtained in vitro. Our results show that ycLeuRS post-transfer editing plays a crucial role in the establishment of the aminoacylation fidelity. When impaired, the viability of cells bearing editing-defective mutants is drastically decreased in the presence of noncognate amino acid. This study also emphasizes the crucial function of some semi-conserved residues around the editing site in modulating the editing efficiency. The assay system can be used to test the effect of compounds that potentially target the aminoacylation or editing active site of fungal LeuRS.

Two forms of human cytoplasmic arginyl-tRNA synthetase produced from two translation initiations by a single mRNA.

Human cytoplasmic arginyl-tRNA synthetase (ArgRS) is a component of a macromolecular complex consisting of at least nine tRNA synthetases and three auxiliary proteins. In mammalian cells, ArgRS is present as a free protein as well as a component of the complex. Via an alignment of ArgRSs from different vertebrates, the genes encoding full-length human cytoplasmic ArgRS and an N-terminal 72-amino acid deletion mutant (hcArgRS and DeltaNhcArgRS, respectively) were subcloned and expressed in Escherichia coli. The two ArgRS products were expressed as a soluble protein in E. coli. The level of production of DeltaNhcArgRS in E. coli and its specific activity were higher than those for hcArgRS. By Western blot analysis, using an antibody against the purified DeltaNhcArgRS, the two forms of ArgRS were detected in three human cell types. The 5'-end cDNA sequence, as confirmed by 5'RACE (5'-rapid amplification of cDNA ends), contained three start codons. Through mutation of the three codons, the two human cytoplasmic ArgRSs were found to be produced in different amounts, indicating that they resulted from two different translation initiation events. Here we show evidence that two forms of human cytoplasmic ArgRS were produced from two translational initiations by a single mRNA.

The C-terminal appended domain of human cytosolic leucyl-tRNA synthetase is indispensable in its interaction with arginyl-tRNA synthetase in the multi-tRNA synthetase complex.

Human cytosolic leucyl-tRNA synthetase is one component of a macromolecular aminoacyl-tRNA synthetase complex. This is unlike prokaryotic and lower eukaryotic LeuRSs that exist as free soluble enzymes. There is little known about it, since the purified enzyme has been unavailable. Herein, human cytosolic leucyl-tRNA synthetase was heterologously expressed in a baculovirus system and purified to homogeneity. The molecular mass (135 kDa) of the enzyme is close to the theoretical value derived from its cDNA. The kinetic constants of the enzyme for ATP, leucine, and tRNA(Leu) in the ATP-PP(i) exchange and tRNA leucylation reactions were determined, and the results showed that it is quite active as a free enzyme. Human cytosolic leucyl-tRNA synthetase expressed in human 293 T cells localizes predominantly to the cytosol. Additionally, it is found to have a long C-terminal extension that is absent from bacterial and yeast LeuRSs. A C-terminal 89-amino acid truncated human cytosolic leucyl-tRNA synthetase was constructed and purified, and the catalytic activities, thermal stability, and subcellular location were found to be almost identical to native enzyme. In vivo and in vitro experiments, however, show that the C-terminal extension of human cytosolic leucyl-tRNA synthetase is indispensable for its interaction with the N-terminal of human cytosolic arginyl-tRNA synthetase in the macromolecular complex. Our results also indicate that the two molecules interact with each other only through their appended domains.

High-level expression and single-step purification of leucyl-tRNA synthetase from Aquifex aeolicus.

Aquifex aeolicus leucyl-tRNA synthetase is the only known heterodimeric LeuRS, consisting of two subunits with molecular masses of 74.0 and 33.5 kDa, and named alphabeta-LeuRS. The gene encoding alpha subunit was cloned into pSBET-b vector. Synthetic oligonucleotide encoding six histidine residues was also inserted in front of alpha subunit. PSBET-b vector contains argU gene, which encodes a rare Escherichia coli tRNA(Arg)(AGA/AGG). The argU gene helps A. aeolicus LeuRS, which contains AGA/AGG codons in exceptionally high frequency, express well in E. coli. The gene encoding beta subunit was inserted into pET-15b vector. E. coli BL21-CodonPlus (DE3) cells were transformed with the two recombinant plasmids to produce alphabeta-LeuRS with a His6 tag at the N-terminus of alpha subunit. The enzyme was purified by affinity chromatography on Ni-NTA Superflow. About 7 mg purified alphabeta-LeuRS was obtained from 250 ml culture. The His6-tag at the N-terminus did not affect the aminoacylation activity of the enzyme.