Long range molecular dynamics study of regulation of eukaryotic glucosamine-6-phosphate synthase activity by UDP-GlcNAc

Glucosamine-6-phosphate (GlcN-6-P) synthase catalyses the first and practically irreversible step in hexosamine metabolism. The final product of this pathway, uridine 5’ diphospho N-acetyl-D-glucosamine (UDP-GlcNAc), is an essential substrate for assembly of bacterial and fungal cell walls. Moreover, the enzyme is involved in phenomenon of hexosamine induced insulin resistance in type II diabetes, which makes it a potential target for antifungal, antibacterial and antidiabetic therapy. The crystal structure of the isomerase domain of GlcN-6-P synthase from human pathogenic fungus Candida albicans, in complex with UDP-GlcNAc has been solved recently but it has not revealed the molecular mechanism of inhibition taking place under UDP-GlcNAc influence, the unique feature of the eukaryotic enzyme. UDP-GlcNAc is a physiological inhibitor of GlcN-6-P synthase, binding about 1 nm away from the active site of the enzyme. In the present work, comparative molecular dynamics simulations of the free and UDP-GlcNAc-bounded structures of GlcN-6-P synthase have been performed. The aim was to complete static X-ray structural data and detect possible changes in the dynamics of the two structures. Results of the simulation studies demonstrated higher mobility of the free structure when compared to the liganded one. Several amino acid residues were identified, flexibility of which is strongly affected upon UDP-GlcNAc binding. Importantly, the most fixed residues are those related to the inhibitor binding process and to the catalytic reaction. The obtained results constitute an important step toward understanding of mechanism of GlcN-6-P synthase inhibition by UDP-GlcNAc molecule

Is TEA an inhibitor for human Aquaporin-1?

Excessive water uptake through aquaporins can be life threatening, and disregulation of water permeability causes many diseases. Therefore, reversible aquaporin inhibitors are highly desired. In this paper, we identified the binding site for tetraethylammonium (TEA) of the membrane water channel aquaporin-1 by a combined molecular docking and molecular dynamics simulation approach. The binding site identified from docking studies was independently confirmed with an unbiased molecular dynamics simulation of an aquaporin tetramer embedded in a lipid membrane, surrounded by a 100-mM tetraethylammonium solution in water. A third independent assessment of the binding site was obtained by umbrella sampling simulations. These simulations, in addition, revealed a binding affinity of more than 17 kJ/mol, corresponding to an IC50 value of << 3 mM. Finally, we observed in our simulations a 50% reduction of the water flux in the presence of TEA, in agreement with water permeability measurements on aquaporin expressed in oocytes. These results confirm TEA as a putative lead for an aquaporin-1 inhibitor

A-Site Residues Move Independently from P-Site Residues in all-Atom Molecular Dynamics Simulations of the 70S Bacterial Ribosome

The ribosome is a large macromolecular machine, and correlated motion between residues is necessary for coordinating function across multiple protein and RNA chains. We ran two all-atom, explicit solvent molecular dynamics simulations of the bacterial ribosome and calculated correlated motion between residue pairs by using mutual information. Because of the short timescales of our simulation (ns), we expect that dynamics are largely local fluctuations around the crystal structure. We hypothesize that residues that show coupled dynamics are functionally related, even on longer timescales. We validate our model by showing that crystallographic B-factors correlate well with the entropy calculated as part of our mutual information calculations. We reveal that A-site residues move relatively independently from P-site residues, effectively insulating A-site functions from P-site functions during translation

Detection of maturity and ligament injury using magic angle directional imaging

Purpose: To investigate whether magnetic field–related anisotropies of collagen may be correlated with postmortem findings in animal models. Methods: Optimized scan planning and new MRI data‐processing methods were proposed and analyzed using Monte Carlo simulations. Six caprine and 10 canine knees were scanned at various orientations to the main magnetic field. Image intensities in segmented voxels were used to compute the orientation vectors of the collagen fibers. Vector field and tractography plots were computed. The Alignment Index was defined as a measure of orientation distribution. The knees were subsequently assessed by a specialist orthopedic veterinarian, who gave a pathological diagnosis after having dissected and photographed the joints. Results: Using 50% less scans than reported previously can lead to robust calculation of fiber orientations in the presence of noise, with much higher accuracy. The 6 caprine knees were found to range from very immature ( 3 years). Mature specimens exhibited significantly more aligned collagen fibers in their patella tendons compared with the immature ones. In 2 of the 10 canine knees scanned, partial cranial caudal ligament tears were identified from MRI and subsequently confirmed with encouragingly high consistency of tractography, Alignment Index, and dissection results. Conclusion: This method can be used to detect injury such as partial ligament tears, and to visualize maturity‐related changes in the collagen structure of tendons. It can provide the basis for new, noninvasive diagnostic tools in combination with new scanner configurations that allow less‐restricted field orientations

Fusing simulation and experiment: The effect of mutations on the structure and activity of the influenza fusion peptide

During the infection process, the influenza fusion peptide (FP) inserts into the host membrane, playing a crucial role in the fusion process between the viral and host membranes. In this work we used a combination of simulation and experimental techniques to analyse the molecular details of this process, which are largely unknown. Although the FP structure has been obtained by NMR in detergent micelles, there is no atomic structure information in membranes. To answer this question, we performed bias-exchange metadynamics (BE-META) simulations, which showed that the lowest energy states of the membrane-inserted FP correspond to helical-hairpin conformations similar to that observed in micelles. BE-META simulations of the G1V, W14A, G12A/G13A and G4A/G8A/G16A/G20A mutants revealed that all the mutations affect the peptide's free energy landscape. A FRET-based analysis showed that all the mutants had a reduced fusogenic activity relative to the WT, in particular the mutants G12A/G13A and G4A/G8A/G16A/G20A. According to our results, one of the major causes of the lower activity of these mutants is their lower membrane affinity, which results in a lower concentration of peptide in the bilayer. These findings contribute to a better understanding of the influenza fusion process and open new routes for future studies

Transition of care: experiences and preferences of patients across the primary/secondary interface – a qualitative study

BACKGROUND: Coordination between care providers of different disciplines is essential to improve the quality of care, in particular for patients with chronic diseases. The way in which general practitioners (GP's) and medical specialists interact has important implications for any healthcare system in which the GP plays the role of gatekeeper to specialist care. Patient experiences and preferences have proven to be increasingly important in discussing healthcare policy. The Dutch government initiated the development of a special website with information for patients on performance indicators of hospitals as well as information on illness or treatment.In the present study we focus on the transition of care at the primary - secondary interface with reference to the impact of patients' ability to make choices about their secondary care providers. The purpose of this study is to (a) explore experiences and preferences of patients regarding the transition between primary and secondary care, (b) study informational resources on illness/treatment desired by patients and (c) determine how information supplied could make it easier for the patient to choose between different options for care (hospital or specialist). METHODS: We conducted a qualitative study using semi-structured focus group interviews among 71 patients referred for various indications in the north and west of The Netherlands. RESULTS: Patients find it important that they do not have to wait, that they are taken seriously, and receive adequate and individually relevant information. A lack of continuity from secondary to primary care was experienced. The patient's desire for free choice of type of care did not arise in any of the focus groups. CONCLUSION: Hospital discharge information needs to be improved. The interval between discharge from specialist care and the report of the specialist to the GP might be a suitable performance indicator in healthcare. Patients want to receive information, tailored to their own situation. The need for information, however, is quite variable. Patients do not feel strongly about self-chosen healthcare, contrary to what administrators presently believe

Scaling behaviour for the water transport in nanoconfined geometries

The transport of water in nanoconfined geometries is different from bulk phase and has tremendous implications in nanotechnology and biotechnology. Here molecular dynamics is used to compute the self-diffusion coefficient D of water within nanopores, around nanoparticles, carbon nanotubes and proteins. For almost 60 different cases, D is found to scale linearly with the sole parameter theta as D(theta)=DB[1+(DC/DB-1)theta], with DB and DC the bulk and totally confined diffusion of water, respectively. The parameter theta is primarily influenced by geometry and represents the ratio between the confined and total water volumes. The D(theta) relationship is interpreted within the thermodynamics of supercooled water. As an example, such relationship is shown to accurately predict the relaxometric response of contrast agents for magnetic resonance imaging. The D(theta) relationship can help in interpreting the transport of water molecules under nanoconfined conditions and tailoring nanostructures with precise modulation of water mobility

Detection of Functional Modes in Protein Dynamics

Proteins frequently accomplish their biological function by collective atomic motions. Yet the identification of collective motions related to a specific protein function from, e.g., a molecular dynamics trajectory is often non-trivial. Here, we propose a novel technique termed “functional mode analysis” that aims to detect the collective motion that is directly related to a particular protein function. Based on an ensemble of structures, together with an arbitrary “functional quantity” that quantifies the functional state of the protein, the technique detects the collective motion that is maximally correlated to the functional quantity. The functional quantity could, e.g., correspond to a geometric, electrostatic, or chemical observable, or any other variable that is relevant to the function of the protein. In addition, the motion that displays the largest likelihood to induce a substantial change in the functional quantity is estimated from the given protein ensemble. Two different correlation measures are applied: first, the Pearson correlation coefficient that measures linear correlation only; and second, the mutual information that can assess any kind of interdependence. Detecting the maximally correlated motion allows one to derive a model for the functional state in terms of a single collective coordinate. The new approach is illustrated using a number of biomolecules, including a polyalanine-helix, T4 lysozyme, Trp-cage, and leucine-binding protein

Glucokinase Gene Mutations: Structural and Genotype-Phenotype Analyses in MODY Children from South Italy

BACKGROUND: Maturity onset diabetes of the young type 2 (or GCK MODY) is a genetic form of diabetes mellitus provoked by mutations in the glucokinase gene (GCK). METHODOLOGY/PRINCIPAL FINDINGS: We screened the GCK gene by direct sequencing in 30 patients from South Italy with suspected MODY. The mutation-induced structural alterations in the protein were analyzed by molecular modeling. The patients' biochemical, clinical and anamnestic data were obtained. Mutations were detected in 16/30 patients (53%); 9 of the 12 mutations identified were novel (p.Glu70Asp, p.Phe123Leu, p.Asp132Asn, p.His137Asp, p.Gly162Asp, p.Thr168Ala, p.Arg392Ser, p.Glu290X, p.Gln106_Met107delinsLeu) and are in regions involved in structural rearrangements required for catalysis. The prevalence of mutation sites was higher in the small domain (7/12: approximately 59%) than in the large (4/12: 33%) domain or in the connection (1/12: 8%) region of the protein. Mild diabetic phenotypes were detected in almost all patients [mean (SD) OGTT = 7.8 mMol/L (1.8)] and mean triglyceride levels were lower in mutated than in unmutated GCK patients (p = 0.04). CONCLUSIONS: The prevalence of GCK MODY is high in southern Italy, and the GCK small domain is a hot spot for MODY mutations. Both the severity of the GCK mutation and the genetic background seem to play a relevant role in the GCK MODY phenotype. Indeed, a partial genotype-phenotype correlation was identified in related patients (3 pairs of siblings) but not in two unrelated children bearing the same mutation. Thus, the molecular approach allows the physician to confirm the diagnosis and to predict severity of the mutation