Theoretical study of N<Subscript>4</Subscript>X (X = O,S, Se) systems

A series of N₄X (X = O, S, Se) compounds have been examined with ab initio and density functional theory (DFT) methods. To our knowledge, these compounds, except for the C_2v ring and the C_3v towerlike isomers of N₄O, are first reported here. The ring structures are the most energetically favored for N₄X (X = O and S) systems. For N₄Se, the cagelike structure is the most energetically favored. Several decomposition and isomerization pathways for the N₄X species have been investigated. The dissociation of C_2v ring N₄O and N₄S structures via ring breaking and the barrier height are only 1.1 and −0.2 kcal mol⁻¹ at the CCSD(T)/6-311+G*//MP2/6-311+G* level of theory. The dissociation of the cagelike N₄X species is at a cost of 12.1–16.2 kcal mol⁻¹. As for the towerlike and triangle bipyramidal isomers, their decomposition or isomerization barrier heights are all lower than 10.0 kcal mol⁻¹. Although the C_S cagelike N₄S isomer has a moderate isomerization barrier (18.3–29.1 kcal mol⁻¹), the low dissociation barrier (−1.0 kcal mol⁻¹) indicates that it will disappear when going to the higher CCSD(T) level. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Musculoskeletal model choice influences hip joint load estimations during gait

The prevalence of musculoskeletal modeling studies investigating hip contact forces and the number of models used to conduct such investigations has increased in recent years. However, the consistency between models remain unknown and differences in model predicted hip contact forces between studies are difficult to distinguish from natural inter-individual differences. The purpose of this study was therefore to evaluate differences in hip joint contact forces during gait between four OpenSim models. These models included the generic models gait2392 and the Arnold Lower Limb Model, as well as the hip specific models hip2372 and London Lower Limb Model. Data from four individuals who have had a total hip replacement with instrumented hip implants performing slow, normal, and fast walking trials were taken from the HIP98 database to evaluate the various models effectiveness at estimating hip loads. Muscle forces were estimated using static optimization and hip contact forces were calculated using the JointReaction analysis in OpenSim. Results indicated that, for gait, the hip specific London Lower Limb Model consistently predicted peak push-off hip joint contact forces with lower magnitude and timing errors compared to the other models. Likewise, root mean square error values were lowest and correlation coefficients were highest for the London Lower Limb Model. These results suggest that the London Lower Limb Model is the most appropriate model for investigations focused on hip joint loading.     

Identification and characterization of carbapenem binding sites within the RND-transporter AcrB

Understanding the molecular determinants for recognition, binding and transport of antibiotics by multidrug efflux systems is important for basic research and useful for the design of more effective antimicrobial compounds. Imipenem and meropenem are two carbapenems whose antibacterial activity is known to be poorly and strongly affected by MexAB-OprM, the major efflux pump transporter in Pseudomonas aeruginosa. However, not much is known regarding recognition and transport of these compounds by AcrAB-TolC, which is the MexAB-OprM homologue in Escherichia coli and by definition the paradigm model for structural studies on efflux pumps. Prompted by this motivation, we unveiled the molecular details of the interaction of imipenem and meropenem with the transporter AcrB by combining computer simulations with biophysical experiments. Regarding the interaction with the two main substrate binding regions of AcrB, the so-called access and deep binding pockets, molecular dynamics simulations revealed imipenem to be more mobile than meropenem in the former, while comparable mobilities were observed in the latter. This result is in line with isothermal titration calorimetry, differential scanning experiments, and binding free energy calculations, indicating a higher affinity for meropenem than imipenem at the deep binding pocket, while both sharing similar affinities at the access pocket. Our findings rationalize how different physico-chemical properties of compounds reflect on their interactions with AcrB. As such, they constitute precious information to be exploited for the rational design of antibiotics able to evade efflux pumps.     

Price of disorder in the lac repressor hinge helix

The Lac system of genes has been pivotal in understanding gene regulation. When the lac repressor protein binds to the correct DNA sequence, the hinge region of the protein goes through a disorder to order transition. The structure of this region of the protein is well understood when it is in this bound conformation, but less so when it is not. Structural studies show that this region is flexible. Our simulations show this region is extremely flexible in solution; however, a high concentration of salt can help kinetically trap the hinge helix. Thermodynamically, disorder is more favorable without the DNA present.     

Modeling,simulations, and optimization of smooth muscle cell tissue engineering for the production of vascular grafts

The paper presents a transient, continuum, two-phase model of the tissue engineering in fibrous scaffolds, including transport equations for the flowing culture medium, nutrient and cell concentration with transverse and in-plane diffusion and cell migration, a novel feature of local in-plane transport across a phenomenological pore and innovative layer-by-layer cell filling approach. The model is successfully validated for the smooth muscle cell tissue engineering of a vascular graft using crosslinked, electrospun gelatin fiber scaffolds for both static and dynamic cell culture, the latter in a dynamic bioreactor with a rotating shaft on which the tubular scaffold is attached. Parametric studies evaluate the impact of the scaffold microstructure, cell dynamics, oxygen transport, and static or dynamic conditions on the rate and extent of cell proliferation and depth of oxygen accessibility. An optimized scaffold of 75% dry porosity is proposed that can be tissue engineered into a viable and still fully oxygenated graft of the tunica media of the coronary artery within 2 days in the dynamic bioreactor. Such scaffold also matches the mechanical properties of the tunica media of the human coronary artery and the suture retention strength of a saphenous vein, often used as a coronary artery graft.     

Insulin adsorption onto zinc oxide nanoparticle mediates conformational rearrangement into amyloid-prone structure with enhanced cytotoxic propensity

Background

Injection localized amyloidosis is one of the most prevalent disorders in type II diabetes mellitus (TIIDM) patients relying on insulin injections. Previous studies have reported that nanoparticles can play a role in the amyloidogenic process of proteins. Hence, the present study deals with the effect of zinc oxide nanoparticles (ZnONP) on the amyloidogenicity and cytotoxicity of insulin.

Effect of heterochiral inversions on the structure of a β-hairpin peptide

We study computationally a family of β-hairpin peptides with systematically introduced chiral inversions, in explicit water, and we investigate the extent to which the backbone structure is able to fold in the presence of heterochiral perturbations. In contrast to the recently investigated case of a helical peptide, we do not find a monotonic change in secondary structure content as a function of the number of L- to D-inversions. The effects of L- to D-inversions are instead found to be highly position-specific. Additionally, in contrast to the helical peptide, some inversions increase the stability of the folded peptide: in such cases, we compute an increase in β-sheet content in the aqueous solution equilibrium ensemble. However, the tertiary structures of the stable (folded) configurations for peptides for which inversions cause an increase in β-sheet content show differences from one another, as well as from the native fold of the nonchirally perturbed β-hairpin. Our results suggest that although some chiral perturbations can increase folding stability, chirally perturbed proteins may still underperform functionally, given the relationship between structure and function.     

Combinatorial peptide library screening for discovery of diverse α-glucosidase inhibitors using molecular dynamics simulations and binary QSAR models

Human α-glucosidase is an enzyme involved in the catalytic cleavage of the glucoside bond and involved in numerous functionalities of the organism, as well as in the insurgence of diabetes mellitus 2 and obesity. Thus, developing chemicals that inhibit this enzyme is a promising approach for the treatment of several pathologies. Small peptides such as di- and tri-peptides may be in natural organism as well as in the GI tract in high concentration, coming from the digestive process of meat, wheat and milk proteins. In this work, we reported the first tentative hierarchical structure-based virtual screening of peptides for human α-glucosidase. The goal of this work is to discover novel and diverse lead compounds that my act as inhibitors of α-glucosidase such as small peptides by performing a computer aided virtual screening and to find novel scaffolds for further development. Thus, in order to select novel candidates with original structure we performed molecular dynamics (MD) simulations among the 12 top-ranked peptides taking as comparison the MD simulations performed on crystallographic inhibitor acarbose. The compounds with the lower RMSD variability during the MD, were reserved for in vitro biological assay. The selected 4 promising structures were prepared on solid phase peptide synthesis and used for the inhibitory assay, among them compound 2 showed good inhibitory activity, which validated our method as an original strategy to discover novel peptide inhibitors. Moreover, pharmacokinetic profile predictions of these 4 peptides were also carried out with binary QSAR models using MetaCore/MetaDrug applications.     

Molecular Modelling reveals the inhibition mechanism and structure–activity relationship of curcumin and its analogues to Staphylococcal aureus Sortase A

Previous studies found that the activity of Sortase A, a bacterial surface protein from Staphylococcus aureus, was inhibited by curcumin and its analogues. To explore this inhibitory mechanism, Sortase A and its inhibitors in complex systems were studied by molecular docking, molecular modelling, binding energy decomposition calculation and steered molecular dynamics simulations. Energy decomposition analysis indicated that PRO-163, LEU-169, GLN-172, ILE-182 and ILE-199 are key residues in Sortase A-inhibitor complexes. Furthermore, interactions between the methoxyl group on the benzene ring in the conjugated molecule (curcumin, demethoxycurcumin, bisdemethoxycurcumin) and VAL-168, LEU-169 and GLN-172 induce the inhibitory activity based on the energy decomposition and distance analyses between the whole residues and inhibitors. However, because of its coiled structure, the non-conjugated molecule, tetrahydrocurcumin, with key residues in the binding sites of Sortase A, interacted weakly with SrtA, leading to the loss of inhibitory activity. Based on these results, the methoxyl group on the benzene ring in the conjugated molecule largely influenced the inhibitory activity of the Sortase A inhibitors.