RASPA: molecular simulation software for adsorption and diffusion in flexible nanoporous materials

A new software package, RASPA, for simulating adsorption and diffusion of molecules in flexible nanoporous materials is presented. The code implements the latest state-of-the-art algorithms for molecular dynamics and Monte Carlo (MC) in various ensembles including symplectic/measure-preserving integrators, Ewald summation, configurational-bias MC, continuous fractional component MC, reactive MC and Baker's minimisation. We show example applications of RASPA in computing coexistence properties, adsorption isotherms for single and multiple components, self- and collective diffusivities, reaction systems and visualisation. The software is released under the GNU General Public License.     

A biomechanical study of the relationship between running velocity and three-dimensional lumbosacral kinetics

Faster running is not performed with proportional increase in all joint torque/work exertions. Although previous studies have investigated lumbopelvic kinetics for a single velocity, it is unclear whether each lumbopelvic torque should increase for faster running. We examined the relationship between running velocity and lumbopelvic kinetics. We calculated the three-dimensional lumbosacral kinetics of 10 male sprinters during steady-state running on a temporary indoor running track at five target velocities: 3.0 (3.20 ± 0.16), 4.5 (4.38 ± 0.18), 6.0 (5.69 ± 0.47), 7.5 (7.30 ± 0.41), and maximal sprinting (9.27 ± 0.36 m/s). The lumbosacral axial rotation torque increased more markedly (from 0.37 ± 0.06 to 1.99 ± 0.46 Nm/kg) than the extension and lateral flexion torques. The increase in the axial rotation torque was larger above 7.30 m/s. Conversely, the extension and lateral flexion torques plateaued when running velocity increased above 7.30 m/s. Similar results were observed for mechanical work. The results indicate that faster running required larger lumbosacral axial rotation torque. Conversely, the extension and lateral flexion torques were relatively invariant to running velocity above 7 m/s, implying that faster running below 7 m/s might increase the biomechanical loads causing excessive pelvic posterior tilt and excessive pelvic drop which has the potential to cause pain/injury related to lumbopelvic extensors and lateral flexors, whereas these biomechanical loads might not relate with running velocity above 7 m/s.     

Efficient tuning of potential parameters for liquid–solid interactions

Spherical and cylindrical water droplets on silicon surface are studied to tune the silicon–oxygen interaction. We use molecular dynamics simulations to estimate the contact angle of two different shaped droplets. We found that the cylindrical droplets are independent of the line tension as their three phases curvature is equal zero. Additionally, we compare an analytical model, taking into account or not the Tolman length and we show that for spherical small size droplets, this length is important to be included, in contrast to cylindrical droplets in which the influence of the Tolman length is negligible. We demonstrate that the usual convenient way to exclude linear tension in the general case can give wrong results. Here, we consider cylindrical droplets, since their contact angle does not depend on the droplet size in the range of few to 10ths of nanometres. The droplets are stabilised due to the periodic boundary conditions. This allows us to propose a new parameterisation for nanoscale droplets, which is independent the size of the droplets or its shape, minimising at the same time the calculation procedure. With the proposed methodology, we can extract the epsilon parameter of the interaction potential between a liquid and a solid from the nanoscaled molecular simulation with only as input the macrosized experimental wetting angle for a given temperature.     

Substrate docking and molecular dynamic simulation for prediction of fungal enzymes from Trichoderma species-assisted extraction of nanocellulose from oil palm leaves

Abstract

Fungi of the Trichoderma species are valued industrial enzymes in support of the ‘zero-waste’ technology to convert agro-industrial biomass into valuable products, i.e. nanocellulose (NC). In this study, an in silico approach using substrate docking and molecular dynamic (MD) simulation was used to predict the order of which the multilayers of cellulosic polymers, i.e. lignin, hemicellulose and cellulose in oil palm leaves (OPL) are degraded by fungal enzymes, endocellulase and exocellulase. The study aimed to establish the catalytic tendencies of the enzymes to optimally degrade the cellulosic components of OPL for high yield production of NC. Energy minimized endocellulase and exocellulase models revealed satisfactory scores of PROCHECK (90.0% and 91.2%), Verify3D (97.23% and 98.85%) and ERRAT (95.24% and 91.00%) assessments. Active site prediction by blind docking, COACH meta-server and multiple sequence alignment indicated the catalytic triads for endocellulase and exocellulase were Ser116–His205–Glu249 and Ser382–Arg124–Asp385, respectively. Binding energy of endocellulase docked with hemicellulose (?6.0 ? kcal mol^?1) was the most favourable followed by lignin (?5.6 ? kcal mol^?1) and cellulose (?4.4 ? kcal mol^?1). Exocellulase, contrarily, bonded favorably with lignin (?8.7 ? kcal mol^?1), closely followed by cellulose (?8.5 ? kcal mol^?1) and hemicellulose (?8.4 ? kcal mol^?1). MDs simulations showed that interactions of complexes, endocellulase–hemicellulose and the exocellulase–cellulose being the most stable. Thus, the findings of the study successfully identified the specific actions of sugar-acting enzymes for NC production.

Communicated by Ramaswamy H. Sarma     

Endocytosis and Exocytosis Events Regulate Vesicle Traffic in Endothelial Cells

We used water-soluble styryl pyridinium dyes that fluoresce at the membrane-water interface to study vesicle traffic in endothelial cells. Cultured endothelial cells derived from bovine and human pulmonary microvessels were incubated in styryl probes, washed to remove dye from the plasmalemmal outer face, and observed by digital fluorescence microscopy. Vesicles that derived from plasmalemma by endocytosis were filled with the styryl dye. These vesicles were distributed throughout the cytosol as numerous particles of heterogeneous diameter and brightness. Vesicle formation was activated 2-fold following addition of extracellular albumin whereas a control protein, immunoglobulin G, had no effect. Dye uptake was abrogated by labeling at low temperatures and inhibitors of phosphoinositide-3-kinase (PI 3-kinase). Tyrosine kinase inhibitors (genistein and herbimycin A) prevented the albumin-induced vesicle formation. Cytochalasin B prevented vesicle redistribution indicating involvement of actin filaments in translocation of endosomes away from sites of vesicle formation. Styryl dye was lost from cells by exocytosis as evident by the disappearance of discrete fluorescent particles. N-ethylmaleimide and botulinum toxin types A and B caused cells to accumulate increased number of vesicles suggesting that exocytosis was regulated by NSF-dependent SNARE mechanism. The results suggest that phosphoinositide metabolism regulates endocytosis in endothelial cells and that extracellular albumin activates endocytosis by a mechanism involving tyrosine phosphorylation, whereas exocytosis is a distinct process regulated by the SNARE machinery. The results support the hypothesis that albumin regulates its internalization and release in vascular endothelial cells via activation of specific endocytic and exocytic pathways.     

Short-term variability of plant populations within a regularly disturbed habitat

The changing populations of weeds during 13 years of the Broadbalk continuous wheat experiment were analysed to investigate the extent of differences in shortterm variability of cover between species. The data were from two sections of the experiment where winter wheat was grown continuously under herbicide treatment for 13 and 6 years respectively. Logistic regressions were fitted to the data. Equisetum arvense showed significant long-term increases on both sections; long-term trends were also detected in the longer data run for Agrostis stolonifera, Cirsium arvensa, Poa trivialis, Ranunculus arvensis and Vicia sativa, and for Medicago lupulina on the shorter data run. Variation around long-term trends was low in the case of Equisetum, and, in the longer data run, for Cirsium and Tussilago farfara, and high for Poa spp. and Vicia. Cover values on the two sections were positively correlated for Alopecurus myosuroides, Equisetum, Poa annua and Tripleurospermum inodorum. There was a weak correlation between C-S-R strategy and short-term variability; the more competitive species displayed less variability than the ruderal species. Furthermore, species regenerating from persistent seed banks were more variable in the short term than those regenerating from short-lived seed or bud banks. This can be explained by differences in response to year-to-year variation in environmental conditions, those species with persistent seed banks being typically more sensitive to annual fluctuations in the environment than those without.     

Root growth dynamics linked to above-ground growth in walnut (Juglans regia)

Background and Aims Examination of plant growth below ground is relatively scant compared with that above ground, and is needed to understand whole-plant responses to the environment. This study examines whether the seasonal timing of fine root growth and the spatial distribution of this growth through the soil profile varies in response to canopy manipulation and soil temperature.Methods Plasticity in the seasonal timing and vertical distribution of root production in response to canopy and soil water manipulation was analysed in field-grown walnut (Juglans regia ‘Chandler’) using minirhizotron techniques.Key Results Root production in walnuts followed a unimodal curve, with one marked flush of root growth starting in mid-May, with a peak in mid-June. Root production declined later in the season, corresponding to increased soil temperature, as well as to the period of major carbohydrate allocation to reproduction. Canopy and soil moisture manipulation did not influence the timing of root production, but did influence the vertical distribution of roots through the soil profile. Water deficit appeared to promote root production in deeper soil layers for mining soil water. Canopy removal appeared to promote shallow root production.Conclusions The findings of this study add to growing evidence that root growth in many ecosystems follows a unimodal curve with one marked flush of root growth in coordination with the initial leaf flush of the season. Root vertical distribution appeared to have greater plasticity than timing of root production in this system, with temperature and/or carbohydrate competition constraining the timing of root growth. Effects on root distribution can have serious impacts on trees, with shallow rooting having negative impacts in years with limited soil water or positive impacts in years with wet springs, and deep rooting having positive impacts on soil water mining from deeper soil layers but negative impacts in years with wet springs.     

Analysis of non-peptidic compounds as potential malarial inhibitors against Plasmodial cysteine proteases via integrated virtual screening workflow

Falcipain-2 (FP-2) and falcipain-3 (FP-3), haemoglobin-degrading enzymes in Plasmodium falciparum, are validated drug targets for the development of effective inhibitors against malaria. However, no commercial drug-targeting falcipains has been developed despite their central role in the life cycle of the parasites. In this work, in silico approaches are used to identify key structural elements that control the binding and selectivity of a diverse set of non-peptidic compounds onto FP-2, FP-3 and homologues from other Plasmodium species as well as human cathepsins. Hotspot residues and the underlying non-covalent interactions, important for the binding of ligands, are identified by interaction fingerprint analysis between the proteases and 2-cyanopyridine derivatives (best hits). It is observed that the size and chemical type of substituent groups within 2-cyanopyridine derivatives determine the strength of protein–ligand interactions. This research presents novel results that can further be exploited in the structure-based molecular-guided design of more potent antimalarial drugs.