Blind, one-eyed, or eagle-eyed? pKa calculations during blind predictions with staphylococcal nuclease

In the current contribution, the performance of Poisson-Boltzmann-based pK(a) calculations of SNase mutants as part of a blind prediction exercise facilitated by the pK(a) cooperative ("pK(a) _coop") is described. A one parameter setting ("quick&dirty" approach) is used to provide an industry perspective where strong time constraints are frequently encountered. On the one hand, results are analyzed in terms of root mean square deviation performance. Furthermore, the pK(a) calculations are assessed for their ability to properly assign protonation state. For this purpose, a new measure called BIPS (binary protonation state at physiological pH) is introduced. Significant differences were found with both comparison measures based on the class of residues examined. In addition, the performance of PROPKA3 as well as the NULL model is examined on the same data set. Finally, pK(a) calculations on SNase mutants with available structural information have been performed and provide support for our calculation methods. The performance on this subset is better than on the pK(a) cooperative mutation data. In the pK(a) _coop data, no structural information on the generated mutants is available. This suggests the occurrence of a substantial structural rearrangement on the insertion of additional charged groups into SNase, which leads to improved prediction quality.     

Parametrization scheme with accuracy and transferability for tight-binding electronic structure calculations with extended Hückel approximation and molecular dynamics simulations

A transferable tight-binding parametrization procedure for extended Hückel approximation is proposed, with the charge self-consistent scheme, that could be applied to the quantum molecular dynamics (MD) simulation for long-time dynamics of large-scale systems. In this procedure, either a target molecule is divided into small molecules or another realistic set of small molecules characterizing chemical bonds in the complicated target molecule is adopted. Then, the parameters for these small molecules are adjusted and compared with reference results of energy levels and wave functions by, for example, density functional theory. Upon application to the large target molecule, these parameters are then readjusted directly in the target molecule. An example is demonstrated with MD simulation applied to the ionic liquid molecule N-methyl-N-propylpiperidinium bis trifluoromethanesulfonyl imide (PP13-TFSI). The origin and stability of HOMO–LUMO gap are discussed.     

HR-pQCT-based homogenised finite element models provide quantitative predictions of experimental vertebral body stiffness and strength with the same accuracy as μFE models

This study validated two different high-resolution peripheral quantitative computer tomography (HR-pQCT)-based finite element (FE) approaches, enhanced homogenised continuum-level (hFE) and micro-finite element (μFE) models, by comparing them with compression test results of vertebral body sections. Thirty-five vertebral body sections were prepared by removing endplates and posterior elements, scanned with HR-pQCT and tested in compression up to failure. Linear hFE and μFE models were created from segmented and grey-level CT images, and apparent model stiffness values were compared with experimental stiffness as well as strength results. Experimental and numerical apparent elastic properties based on grey-level/segmented CT images (N=35) correlated well for μFE (r2=0.748/0.842) and hFE models (r2=0.741/0.864). Vertebral section stiffness values from the linear μFE/hFE models estimated experimental ultimate apparent strength very well (r2=0.920/0.927). Calibrated hFE models were able to predict quantitatively apparent stiffness with the same accuracy as μFE models. However, hFE models needed no back-calculation of a tissue modulus or any kind of fitting and were computationally much cheaper.     

Comparing protein structures and inferring functions with a novel three-dimensional Yau–Hausdorff method

Abstract

Structures and functions of proteins play various essential roles in biological processes. The functions of newly discovered proteins can be predicted by comparing their structures with that of known-functional proteins. Many approaches have been proposed for measuring the protein structure similarity, such as the template-modeling (TM)-score method, GRaphlet (GR)-Align method as well as the commonly used root-mean-square deviation (RMSD) measures. However, the alignment comparisons between the similarity of protein structure cost much time on large dataset, and the accuracy still have room to improve. In this study, we introduce a new three-dimensional (3D) Yau–Hausdorff distance between any two 3D objects. The (3D) Yau–Hausdorff distance can be used in particular to measure the similarity/dissimilarity of two proteins of any size and does not need aligning and superimposing two structures. We apply structural similarity to study function similarity and perform phylogenetic analysis on several datasets. The results show that (3D) Yau–Hausdorff distance could serve as a more precise and effective method to discover biological relationships between proteins than other methods on structure comparison.

Communicated by Ramaswamy H. Sarma     

Biomechanical differences of Coflex-F and pedicle screw fixation combined with TLIF or ALIF – a finite element study

Lumbar interbody fusion is a common procedure for treating lower back pain related to degenerative disc diseases. The Coflex-F is a recently developed interspinous spacer, the makers of which claim that it can provide stabilisation similar to pedicle screw fixation. Therefore, this study compares the biomechanical behaviour of the Coflex-F device and pedicle screw fixation with transforaminal lumbar interbody fusion (TLIF) or anterior lumbar interbody fusion (ALIF) surgeries by using finite element analysis. The results show that the Coflex-F device combined with ALIF surgery can provide stability similar to the pedicle screw fixation combined with TLIF or ALIF surgery. Also, the posterior instrumentations (Coflex-F and pedicle screw fixation) combined with TLIF surgery had lower stability than when combined with ALIF surgery.     

Modelling and simulation of porcine liver tissue indentation using finite element method and uniaxial stress–strain data

We hypothesize that both compression and elongation stress–strain data should be considered for modeling and simulation of soft tissue indentation. Uniaxial stress–strain data were obtained from in vitro loading experiments of porcine liver tissue. An axisymmetric finite element model was used to simulate liver tissue indentation with tissue material represented by hyperelastic models. The material parameters were derived from uniaxial stress–strain data of compressions, elongations, and combined compression and elongation of porcine liver samples. in vitro indentation tests were used to validate the finite element simulation. Stress–strain data from the simulation with material parameters derived from the combined compression and elongation data match the experimental data best. This is due to its better ability in modeling 3D deformation since the behavior of biological soft tissue under indentation is affected by both its compressive and tensile characteristics. The combined logarithmic and polynomial model is somewhat better than the 5-constant Mooney–Rivlin model as the constitutive model for this indentation simulation.     

Modeling nutrients,oxygen and critical phosphorus loading in a shallow reservoir in China with a coupled water quality – Macrophytes model

Many macrophyte-dominated clear lakes switch to a phytoplankton-dominated turbid state when the lake becomes eutrophic. An existing Yuqiao Reservoir Water Quality Model (YRWQM) and the macrophyte submodel were coupled to simulate the effect of submerged macrophytes on nutrients and dissolve oxygen cycles in a shallow reservoir in China. The level of phosphorus loading in a transition from a clear to turbid state was addressed using the integrated model. The model runs from seedling establishment until dying out, from March 1 to July 18 in 2009. The simulations were performed for a contingent range of P loadings, starting from three different initial conditions. The results indicated that the integrated model improves accuracy of predictions compared to YRWQM. The concentrations of nutrients declined slightly during the macrophyte growth period in the reservoir and dissolved oxygen increased slightly. Although nutrient concentrations increased by submerged macrophyte release during the extinction period, the effect on the nutrients was less than that of transfer with nutrient-rich water. More released nutrients may enhance increases in substantial abundance. The critical phosphorus loading level during a switch from the clear to turbid state was estimated by these scenarios. The threshold for the switch is ∼6.1 mgP m⁻² d⁻¹ with an initial total phosphorus concentration of 160 μg l⁻¹. Moreover, the results demonstrated that the switch was also dependent on the initial total phosphorus concentration. These results suggest that the reservoir in a clear water state is at risk of a switch as nutrient levels are close to the critical levels.     

Coupled molecular dynamics and continuum electrostatic method to compute the ionization pKa’s of proteins as a function of pH. Test on a large set of proteins

A computational method, to predict the pKa values of the ionizable residues Asp, Glu, His, Tyr, and Lys of proteins, is presented here. Calculation of the electrostatic free-energy of the proteins is based on an efficient version of a continuum dielectric electrostatic model. The conformational flexibility of the protein is taken into account by carrying out molecular dynamics simulations of 10 ns in implicit water. The accuracy of the proposed method of calculation of pKa values is estimated from a test set of experimental pKa data for 297 ionizable residues from 34 proteins. The pKa-prediction test shows that, on average, 57, 86, and 95% of all predictions have an error lower than 0.5, 1.0, and 1.5 pKa units, respectively. This work contributes to our general understanding of the importance of protein flexibility for an accurate computation of pKa, providing critical insight about the significance of the multiple neutral states of acid and histidine residues for pKa-prediction, and may spur significant progress in our effort to develop a fast and accurate electrostatic-based method for pKa-predictions of proteins as a function of pH.     

An α/β hydrolase and associated Per-ARNT-Sim domain comprise a bipartite sensing module coupled with diverse output domains

The RsbQ α/β hydrolase and RsbP serine phosphatase form a signaling pair required to activate the general stress factor σ(B) of Bacillus subtilis in response to energy limitation. RsbP has a predicted N-terminal Per-ARNT-Sim (PAS) domain, a central coiled-coil, and a C-terminal protein phosphatase M (PPM) domain. Previous studies support a model in which RsbQ provides an activity needed for PAS to regulate the phosphatase domain via the coiled-coil. RsbQ and the PAS domain (RsbP-PAS) therefore appear to form a sensory module. Here we test this hypothesis using bioinformatic and genetic analysis. We found 45 RsbQ and RsbP-PAS homologues encoded by adjacent genes in diverse bacteria, with PAS and a predicted coiled-coil fused to one of three output domains: PPM phosphatase (Gram positive bacteria), histidine protein kinase (Gram negative bacteria), and diguanylate cyclase (both lineages). Multiple alignment of the RsbP-PAS homologues suggested nine residues that distinguish the class. Alanine substitutions at four of these conferred a null phenotype in B. subtilis, indicating their functional importance. The F55A null substitution lay in the Fα helix of an RsbP-PAS model. F55A inhibited interaction of RsbP with RsbQ in yeast two-hybrid and pull-down assays but did not significantly affect interaction of RsbP with itself. We propose that RsbQ directly contacts the PAS domains of an RsbP oligomer to provide the activating signal, which is propagated to the phosphatase domains via the coiled-coil. A similar mechanism would allow the RsbQ-PAS module to convey a common input signal to structurally diverse output domains, controlling a variety of physiological responses.     

Red blood cell simulation using a coupled shell–fluid analysis purely based on the SPH method

Biomechanics and Modeling in Mechanobiology - In this paper, a novel 3D numerical method has been developed to simulate red blood cells (RBCs) based on the interaction between a shell-like solid...     

The transition state of coupled folding and binding for a flexible β-finger

Flexible and fully disordered protein regions that fold upon binding mediate numerous protein-protein interactions. However, little is known about their mechanism of interaction. One such coupled folding and binding occurs when a flexible region of neuronal nitric oxide synthase adopts a β-finger structure upon binding to its protein ligand, a PDZ [PSD-95 (postsynaptic density protein-95)/Discs large/ZO-1] domain from PSD-95. We have analyzed this binding reaction by protein engineering combined with kinetic experiments. Mutational destabilization of the β-finger changed mainly the dissociation rate constant of the proteins and, to a lesser extent, the association rate constant. Thus, mutation affected late events in the coupled folding and binding reaction. Our results therefore suggest that the native binding interactions of the β-finger are not present in the rate-limiting transition state for binding but form on the downhill side in a cooperative manner. However, by mutation, we could destabilize the β-finger further and change the rate-limiting step such that an initial conformational change becomes rate limiting. This switch in rate-limiting step shows that multistep binding mechanisms are likely to be found among flexible and intrinsically disordered regions of proteins.     

Synthesis of aromatic n-alkyl-glucoside esters in a coupled β-glucosidase and lipase reaction

Regioselective synthesis of e.g. 6-O-phenylbutyryl-1-n-butyl--D-glucopyranose was achieved in 21 % yield using almond--glucosidase and Candida antarctica lipase B. The -glucosidase reaction was performed in a biphasic (buffer/n-alcohol) system using free and Eupergit CTM-immobilized glucosidase. Immobilized enzyme allowed product formation even at a water content of 1 %. © Rapid Science Ltd. 1998     

Are indirect genetic benefits associated with polyandry? Testing predictions in a natural population of lemon sharks

Multiple mating has clear fitness benefits for males, but uncertain benefits and costs for females. We tested for indirect genetic benefits of polyandry in a natural population, by using data from a long-term genetic and demographic study of lemon sharks ( Negaprion brevirostris ) at Bimini, Bahamas. To do so, we followed the fates of individuals from six cohorts (450 age-0 and 254 age-1 fish) in relation to their individual level of genetic variation, and whether they were from polyandrous or monoandrous litters. We find that offspring from polyandrous litters did not have a greater genetic diversity or greater survival than did the offspring of monoandrous litters. We also find no evidence of positive associations between individual offspring genetic diversity metrics and our surrogate measure of fitness (i.e. survival). In fact, age-1 individuals with fewer heterozygous microsatellite loci and more genetically similar parents were more likely to survive to age-2. Thus, polyandry in female lemon sharks does not appear to be adaptive from the perspective of indirect genetic benefits to offspring. It may instead be the result of convenience polyandry, whereby females mate multiply to avoid harassment by males. Our inability to find indirect genetic benefits of polyandry despite detailed pedigree and survival information suggests the need for similar assessments in other natural populations.     

Efficient production of deoxynucleoside-5′-monophosphates using deoxynucleoside kinase coupled with a GTP-regeneration system

Deoxynucleoside-5′-monophosphates (5′-dNMPs) are the basic components of DNA and are widely used in medicine and as chemical and biochemical reagents. A large amount of effort has been expended to obtain 5′-dNMPs of high quality and at a low cost. However, these procedures are inefficient and inconvenient. In this study, deoxyadenosine-5′-monophosphate (5′-dAMP), 2,6-diaminopurine deoxynucleoside-5′-monophosphate (5′-dDAMP), and deoxycytidine-5′-monophosphate (5′-dCMP) were biosynthesized using recombinant N-deoxyribosyltransferase II (NDT-II), deoxycytidine kinase, and acetate kinase in a one-pot reaction system. The ndt-II gene from Lactobacillus delbrueckii, dck from Bacillus subtilus, and ack from Escherichia coli K12 were overexpressed in E. coli BL21 (DE3). Thymidine was used as the deoxyribose donor; GTP was used as the phosphate donor, and acetyl phosphate was used to regenerate GTP. Under optimized conditions, each 10 mM adenine, 10 mM 2,6-diaminopurine, or 10 mM cytosine were converted into 9.01 mM 5′-dAMP, 8.68 mM 5′-dDAMP, or 6.23 mM 5′-dCMP, respectively. The high yield indicated that this process of biosynthesis of 5′-dAMP, 5′-dDAMP, or 5′-dCMP was efficient and economical, and this one-pot system may also potentially be used for the preparation of other types of 5′-dNMPs.     

Dating and localizing an invasion from post-introduction data and a coupled reaction–diffusion–absorption model

Journal of Mathematical Biology - Invasion of new territories by alien organisms is of primary concern for environmental and health agencies and has been a core topic in mathematical modeling, in...     

A new colorimetric DPPH• scavenging activity method with no need for a spectrophotometer applied on synthetic and natural antioxidants and medicinal herbs

2,2-Diphenyl-1-picrylhydrazyl (DPPH^?) radical scavenging, the most commonly used antioxidant method with more than seventeen thousand articles cited, is very practical; however, as with most assays, it has the major disadvantage of dependence on a spectrophotometer. To overcome this drawback, the colorimetric determination of the antioxidant activity using a scanner and freely available Image J software was developed. In this new method, the mixtures of solutions of DPPH^? and standard antioxidants or extracts of common medicinal herbs were dropped onto TLC plates, after an incubation period. The spot images were evaluated with Image J software to determine CSC₅₀ values, the sample concentrations providing 50% colour reduction, which were very similar with the SC₅₀ values obtained with spectrophotometric method. The advantages of the new method are the use of lower amounts of reagents and solvents, no need for costly spectrophotometers, and thus significantly lowered costs, and convenient implementation in any environment and situation.