A Finite Element Model of Cerebral Vascular Injury for Predicting Microbleeds Location.

Finite Element (FE) models of brain mechanics have improved our understanding of the brain response to rapid mechanical loads that produce traumatic brain injuries. However, these models have rarely incorporated vasculature, which limits their ability to predict the response of vessels to head impacts. To address this shortcoming, here we used high-resolution MRI scans to map the venous system anatomy at a submillimetre resolution. We then used this map to develop an FE model of veins and incorporated it in an anatomically detailed FE model of the brain. The model prediction of brain displacement at different locations was compared to controlled experiments on post-mortem human subject heads, yielding over 3,100 displacement curve comparisons, which showed fair to excellent correlation between them. We then used the model to predict the distribution of axial strains and strain rates in the veins of a rugby player who had small blood deposits in his white matter, known as microbleeds, after sustaining a head collision. We hypothesised that the distribution of axial strain and strain rate in veins can predict the pattern of microbleeds. We reconstructed the head collision using video footage and multi-body dynamics modelling and used the predicted head accelerations to load the FE model of vascular injury. The model predicted large axial strains in veins where microbleeds were detected. A region of interest analysis using white matter tracts showed that the tract group with microbleeds had 95th percentile peak axial strain and strain rate of 0.197 and 64.9 s-1 respectively, which were significantly larger than those of the group of tracts without microbleeds (0.163 and 57.0 s-1). This study does not derive a threshold for the onset of microbleeds as it investigated a single case, but it provides evidence for a link between strain and strain rate applied to veins during head impacts and structural damage and allows for future work to generate threshold values. Moreover, our results suggest that the FE model has the potential to be used to predict intracranial vascular injuries after TBI, providing a more objective tool for TBI assessment and improving protection against it.

Smoothed particle hydrodynamic modelling of the cerebrospinal fluid for brain biomechanics: Accuracy and stability.

The Cerebrospinal Fluid (CSF) can undergo shear deformations under head motions. Finite Element (FE) models, which are commonly used to simulate biomechanics of the brain, including traumatic brain injury, employ solid elements to represent the CSF. However, the limited number of elements paired with shear deformations in CSF can decrease the accuracy of their predictions. Large deformation problems can be accurately modelled using the mesh-free Smoothed Particle Hydrodynamics (SPH) method, but there is limited previous work on using this method for modelling the CSF. Here we explored the stability and accuracy of key modelling parameters of an SPH model of the CSF when predicting relative brain/skull displacements in a simulation of an in vivo mild head impact in human. The Moving Least Squares (MLS) SPH formulation and Ogden rubber material model were found to be the most accurate and stable. The strain and strain rate in the brain differed across the SPH and FE models of CSF. The FE mesh anchored the gyri, preventing them from experiencing the level of strains seen in the in vivo brain experiments and predicted by the SPH model. Additionally, SPH showed higher levels of strains in the sulci compared to the FE model. However, tensile instability was found to be a key challenge of the SPH method, which needs to be addressed in future. Our study provides a detailed investigation of the use of SPH and shows its potential for improving the accuracy of computational models of brain biomechanics.

The effects of topology and relative density of lattice liners on traumatic brain injury mitigation.

This paper evaluates the effects of topology and relative density of helmet lattice liners on mitigating Traumatic Brain Injury (TBI). Finite Element (FE) models of new lattice liners with prismatic and tetrahedral topologies were developed. A typical frontal head impact in motorcycle accidents was simulated, and linear and rotational accelerations of the head were recorded. A high-fidelity FE model of TBI was loaded with the accelerations to predict the brain response during the accident. The results show that prismatic lattices have better performance in preventing TBI than tetrahedral lattices and EPS that is typically used in helmets. Moreover, varying the cell size through the thickness of the liner improves its performance, but this effect was marginal. The relative density also has a significant effect, with lattices with lower relative densities providing better protection. Across different lattices studied here, the prismatic lattice with a relative density of 6% had the best performance and reduced the peak linear and rotational accelerations, Head Injury Criterion (HIC), brain strain and strain rate by 48%, 37%, 49%, 32% and 65% respectively, compared to the EPS liner. These results can be used to guide the design of lattice helmet liners for better mitigation of TBI.

Enzyme-substrate complexes of allosteric citrate synthase: evidence for a novel intermediate in substrate binding.

The citrate synthase (CS) of Escherichia coli is an allosteric hexameric enzyme specifically inhibited by NADH. The crystal structure of wild type (WT) E. coli CS, determined by us previously, has no substrates bound, and part of the active site is in a highly mobile region that is shifted from the position needed for catalysis. The CS of Acetobacter aceti has a similar structure, but has been successfully crystallized with bound substrates: both oxaloacetic acid (OAA) and an analog of acetyl coenzyme A (AcCoA). We engineered a variant of E. coli CS wherein five amino acids in the mobile region have been replaced by those in the A. aceti sequence. The purified enzyme shows unusual kinetics with a low affinity for both substrates. Although the crystal structure without ligands is very similar to that of the WT enzyme (except in the mutated region), complexes are formed with both substrates and the allosteric inhibitor NADH. The complex with OAA in the active site identifies a novel OAA-binding residue, Arg306, which has no functional counterpart in other known CS-OAA complexes. This structure may represent an intermediate in a multi-step substrate binding process where Arg306 changes roles from OAA binding to AcCoA binding. The second complex has the substrate analog, S-carboxymethyl-coenzyme A, in the allosteric NADH-binding site and the AcCoA site is not formed. Additional CS variants unable to bind adenylates at the allosteric site show that this second complex is not a factor in positive allosteric activation of AcCoA binding.

Urea as a protein destabilizing agent in electrospray ionisation.

Urea is well known as a denaturant of proteins, but there is also evidence that millimolar amounts of urea may in fact stabilize protein complexes. Advances in mass spectrometric analysis have given us the opportunity to test the effect of urea on several noncovalent complexes in buffered solutions. We expected to see lower charge states if folded proteins were more compact (and therefore more stable), and higher charge states if the proteins were denatured. We have found that mM urea interferes with some noncovalent interactions, and that the extent of interference depends on the specific protein complex. The difference seems to be related to the type of interactions, with weak ones, such as H-bonds, more sensitive to urea. Examples show that a quick check with urea may give some insights into protein stability in the mass spectrometer.

Method for estimating the isotopic distributions of metabolically labeled proteins by MALDI-TOFMS: application to NMR samples.

We have developed an efficient method of estimating metabolic incorporation of heavy isotopes into proteins, including those where a single amino acid carries the label. The protein is digested with trypsin, and the resulting peptide mixture is examined directly by MALDI-TOF mass spectrometry. Peptides are chosen for analysis if they contain one or more labeled atoms and also exhibit clearly separated mass spectra. The known atomic composition of the peptide is then used to simulate ion distributions for various proportions of heavy isotope incorporation, to obtain the best match to the observed ion distribution. We demonstrate the method by comparing simulated and observed mass spectra of tryptic peptides of Escherichia coli citrate synthase labeled with 15N in several ways and show that the method is particularly applicable when only one amino acid is isotopically labeled.

Proteomic identification of AP2 gamma as a rat placental lactogen II trophoblast cell-specific enhancer binding protein.

The factors that regulate the developmental expression of the rodent prolactin gene family in placenta remain poorly defined. We previously identified an enhancer element in the 5' flanking region of one family member, rat placental lactogen II (rPLII), which could target reporter gene expression to the placenta in transgenic mice; this enhancer functioned in the Rcho rat trophoblast cell line but not in the rat pituitary GC cell line. In further experiments to identify the factors that bind this element, we have selectively enriched for DNA binding proteins in nuclear extract from Rcho cells using magnetic beads coupled to a 43-bp enhancer oligonucleotide. Tryptic peptides of bound proteins were analyzed by HPLC coupled off-line to matrix-assisted laser desorption ionization time of flight mass spectrometry. Several peptides of AP2 gamma, a key trophoblast cell-specific transcription factor, were identified. Gel mobility shift assays using AP2 gamma-specific antiserum and mutant enhancer oligonucleotides demonstrated binding specifically to the FP2 DNase I-protected region of the element, identifying an atypical binding site for this factor. In cotransfection assays in rat pituitary GC cells, AP2 gamma transactivated the enhancer via this region. Chromatin immunoprecipitation assays confirmed AP2 gamma occupancy of the enhancer region in situ in the nuclei of Rcho giant cells. These data support a role for AP2 gamma in the placental giant cell-specific expression of the rPLII gene and provide the first direct evidence for the involvement of a placental-specific transcription factor in the regulation of a member of this gene family.

Multiple equilibria of the Escherichia coli chaperonin GroES revealed by mass spectrometry.

Nanospray time-of-flight mass spectrometry has been used to study the assembly of the heptamer of the Escherichia coli cochaperonin protein GroES, a system previously described as a monomer-heptamer equilibrium. In addition to the monomers and heptamers, we have found measurable amounts of dimers and hexamers, the presence of which suggests the following mechanism for heptamer assembly: 2 Monomers <--> Dimer; 3 Dimers <--> Hexamer; Hexamer + Monomer <--> Heptamer. Equilibrium constants for each of these steps, and an overall constant for the Monomer <--> Heptamer equilibrium, have been estimated from the data. These constants imply a standard free-energy change, DeltaG(0), of about 9 kcal/mol for each contact surface formed between GroES subunits, except for the addition of the last subunit, where DeltaG(0) = 6 kcal/mol. This lower value probably reflects the loss of entropy when the heptamer ring is formed. These experiments illustrate the advantages of electrospray mass spectrometry as a method of measuring all components of a multiple equilibrium system.

Characterization of a large subunit catalase truncated by proteolytic cleavage.

The large subunit catalase HPII from Escherichia coli can be truncated by proteolysis to a structure similar to small subunit catalases. Mass spectrometry analysis indicates that there is some heterogeneity in the precise cleavage sites, but approximately 74 N-terminal residues, 189 C-terminal residues, and a 9-11-residue internal fragment, including residues 298-308, are removed. Crystal structure refinement at 2.8 A reveals that the tertiary and quaternary structure of the native enzyme is retained with only very subtle changes despite the loss of 36% of the sequence. The truncated variant exhibits a 1.8 times faster turnover rate and enhanced sensitivity to high concentrations of H(2)O(2), consistent with easier access of the substrate to the active site. In addition, the truncated variant is more sensitive to inhibition, particularly by reagents such as aminotriazole and azide which are larger than substrate H(2)O(2). The main channel leading to the heme cavity is largely unaffected by the truncation, but the lateral channel is shortened and its entrance widened by removal of the C-terminal domain, providing an explanation for easier access to the active site. Opening of the entrance to the lateral channel also opens the putative NADPH binding site, but NADPH binding could not be demonstrated. Despite the lack of bound NADPH, the compound I species of both native and truncated HPII are reduced back to the resting state with compound II being evident in the absorbance spectrum only of the heme b-containing H392A variant.

Catalase-peroxidases (KatG) exhibit NADH oxidase activity.

Catalase-peroxidases (KatG) produced by Burkholderia pseudomallei, Escherichia coli, and Mycobacterium tuberculosis catalyze the oxidation of NADH to form NAD+ and either H2O2 or superoxide radical depending on pH. The NADH oxidase reaction requires molecular oxygen, does not require hydrogen peroxide, is not inhibited by superoxide dismutase or catalase, and has a pH optimum of 8.75, clearly differentiating it from the peroxidase and catalase reactions with pH optima of 5.5 and 6.5, respectively, and from the NADH peroxidase-oxidase reaction of horseradish peroxidase. B. pseudomallei KatG has a relatively high affinity for NADH (Km=12 microm), but the oxidase reaction is slow (kcat=0.54 min(-1)) compared with the peroxidase and catalase reactions. The catalase-peroxidases also catalyze the hydrazinolysis of isonicotinic acid hydrazide (INH) in an oxygen- and H2O2-independent reaction, and KatG-dependent radical generation from a mixture of NADH and INH is two to three times faster than the combined rates of separate reactions with NADH and INH alone. The major products from the coupled reaction, identified by high pressure liquid chromatography fractionation and mass spectrometry, are NAD+ and isonicotinoyl-NAD, the activated form of isoniazid that inhibits mycolic acid synthesis in M. tuberculosis. Isonicotinoyl-NAD synthesis from a mixture of NAD+ and INH is KatG-dependent and is activated by manganese ion. M. tuberculosis KatG catalyzes isonicotinoyl-NAD formation from NAD+ and INH more efficiently than B. pseudomallei KatG.

Defining regulatory regions in the rat prolactin gene family locus using a large P1 genomic clone.

Members of the large rat prolactin gene family located on chromosome 17 are expressed in one or more placental trophoblast cell types and maternal decidua at specific times during pregnancy. Studies to identify the factors involved in these highly specific developmental expression patterns, using limited amounts of 5'-flanking DNA, have met with only partial success. Here we report the isolation and characterization of an 80-kb rat genomic clone, P1 12830, containing linked rat placental lactogen II, rat prolactin-like protein-I, and rat prolactin-like protein-B genes with substantial amounts of 5'- and 3'-flanking DNA as well as a rat placental lactogen II-related pseudogene, the first to be described in this gene family. This clone was used to create F0 transgenic mice, and the levels of expression of the three rat genes were compared with those of the endogenous mouse genes, using RT-PCR. Each rat gene was expressed differently in the same placenta, confirming the importance of sufficient flanking sequences in the expression of the individual genes. These studies emphasize the need for large genomic clones in defining the complete complement of factors that regulate the developmental expression of the rat prolactin gene locus.

Characterization of the catalase-peroxidase KatG from Burkholderia pseudomallei by mass spectrometry.

The electron density maps of the catalase-peroxidase from Burkholderia pseudomallei (BpKatG) presented two unusual covalent modifications. A covalent structure linked the active site Trp111 with Tyr238 and Tyr238 with Met264, and the heme was modified, likely by a perhydroxy group added to the vinyl group on ring I. Mass spectrometry analysis of tryptic digests of BpKatG revealed a cluster of ions at m/z 6585, consistent with the fusion of three peptides through Trp111, Tyr238, and Met264, and a cluster at m/z approximately 4525, consistent with the fusion of two peptides linked through Trp111 and Tyr238. MS/MS analysis of the major ions at m/z 4524 and 4540 confirmed the expected sequence and suggested that the multiple ions in the cluster were the result of multiple oxidation events and transfer of CH3-S to the tyrosine. Neither cluster of ions at m/z 4525 or 6585 was present in the spectrum of a tryptic digest of the W111F variant of BpKatG. The spectrum of the tryptic digest of native BpKatG also contained a major ion for a peptide in which Met264 had been converted to homoserine, consistent with the covalent bond between Tyr238 and Met264 being susceptible to hydrolysis, including the loss of the CH3-S from the methionine. Analysis of the tryptic digests of hydroperoxidase I (KatG) from Escherichia coli provided direct evidence for the covalent linkage between Trp105 and Tyr226 and indirect evidence for a covalent linkage between Tyr226 and Met252. Tryptic peptide analysis and N-terminal sequencing revealed that the N-terminal residue of BpKatG is Ser22.

Probing the roles of key residues in the unique regulatory NADH binding site of type II citrate synthase of Escherichia coli.

The citrate synthase of Escherichia coli is an example of a Type II citrate synthase, a hexamer that is subject to allosteric inhibition by NADH. In previous crystallographic work, we defined the NADH binding sites, identifying nine amino acids whose side chains were proposed to make hydrogen bonds with the NADH molecule. Here, we describe the functional properties of nine sequence variants, in which these have been replaced by nonbonding residues. All of the variants show some changes in NADH binding and inhibition and small but significant changes in kinetic parameters for catalysis. In three cases, Y145A, R163L, and K167A, NADH inhibition has become extremely weak. We have used nanospray/time-of-flight mass spectrometry, under non-denaturing conditions, to show that two of these, R163L and K167A, do not form hexamers in response to NADH binding, unlike the wild type enzyme. One variant, R109L, shows tighter NADH binding. We have crystallized this variant and determined its structure, with and without bound NADH. Unexpectedly, the greatest structural changes in the R109L variant are in two regions outside the NADH binding site, both of which, in wild type citrate synthase, have unusually high mobilities as measured by crystallographic thermal factors. In the R109L variant, both regions (residues 260 -311 and 316-342) are much less mobile and have rearranged significantly. We argue that these two regions are elements in the path of communication between the NADH binding sites and the active sites and are centrally involved in the regulatory conformational change in E. coli citrate synthase.

Insights into the evolution of allosteric properties. The NADH binding site of hexameric type II citrate synthases.

Study of the hexameric and allosterically regulated citrate synthases (type II CS) provides a rare opportunity to gain not only an understanding of a novel allosteric mechanism but also insight into how such properties can evolve from an unregulated structural platform (the dimeric type I CS). To address both of these issues, we have determined the structure of the complex of NADH (a negative allosteric effector) with the F383A variant of type II Escherichia coli CS. This variant was chosen because its kinetics indicate it is primarily in the T or inactive allosteric conformation, the state that strongly binds to NADH. Our structural analyses show that the six NADH binding sites in the hexameric CS complex are located at the interfaces between dimer units such that most of each site is formed by one subunit, but a number of key residues are drawn from the adjacent dimer. This arrangement of interactions serves to explain why NADH allosteric regulation is a feature only of hexameric type II CS. Surprisingly, in both the wild-type enzyme and the NADH complex, the two subunits of each dimer within the hexameric conformation are similar but not identical in structure, and therefore, while the general characteristics of NADH binding interactions are similar in each subunit, the details of these are somewhat different between subunits. Detailed examination of the observed NADH binding sites indicates that both direct charged interactions and the overall cationic nature of the sites are likely responsible for the ability of these sites to discriminate between NADH and NAD(+). A particularly novel characteristic of the complex is the horseshoe conformation assumed by NADH, which is strikingly different from the extended conformation found in its complexes with most proteins. Sequence homology studies suggest that this approach to binding NADH may arise out of the evolutionary need to add an allosteric regulatory function to the base CS structure. Comparisons of the amino acid sequences of known type II CS enzymes, from different Gram-negative bacteria taxonomic groups, show that the NADH-binding residues identified in our structure are strongly conserved, while hexameric CS molecules that are insensitive to NADH have undergone key changes in the sequence of this part of the protein.

Hydroperoxidase II of Escherichia coli exhibits enhanced resistance to proteolytic cleavage compared to other catalases.

Catalase (hydroperoxidase) HPII of Escherichia coli is the largest catalase so far characterized, existing as a homotetramer of 84 kDa subunits. Each subunit has a core structure that closely resembles small subunit catalases, supplemented with an extended N-terminal sequence and compact flavodoxin-like C-terminal domain. Treatment of HPII with trypsin, chymotrypsin, or proteinase K, under conditions of limited digestion, resulted in cleavage of 72-74 residues from the N-terminus of each subunit that created a homotetramer of 76 kDa subunits with 80% of wild-type activity. Longer treatment with proteinase K removed the C-terminal domain, producing a transient 59 kDa subunit which was subsequently cleaved into two fragments, 26 and 32 kDa. The tetrameric structure was retained despite this fragmentation, with four intermediates being observed between the 336 kDa native form and the 236 kDa fully truncated form corresponding to tetramers with a decreasing complement of C-termini (4, 3, 2, and 1). The truncated tetramers retained 80% of wild-type activity. The T(m) for loss of activity during heating was decreased from 85 to 77 degrees C by removal of the N-terminal sequence and to 59 degrees C by removal of the C-terminal domain, revealing the importance of the C-terminal domain in enzyme stability. The sites of cleavage were determined by N- and C-terminal sequencing, and two were located on the surface of the tetramer with a third being exposed by removal of the C-terminal domain.

Structure of the Clade 1 catalase, CatF of Pseudomonas syringae, at 1.8 A resolution.

Catalase CatF of Pseudomonas syringae has been identified phylogenetically as a clade 1 catalase, closely related to plant catalases, a group from which no structure has been determined. The structure of CatF has been refined at 1.8 A resolution by using X-ray synchrotron data collected from a crystal flash-cooled with liquid nitrogen. The crystallographic agreement factors R and R(free) are, respectively, 18.3% and 24.0%. The asymmetric unit of the crystal contains a whole molecule that shows accurate 222-point group symmetry. The crystallized enzyme is a homotetramer of subunits with 484 residues, some 26 residues shorter than predicted from the DNA sequence. Mass spectrometry analysis confirmed the absence of 26 N-terminal residues, possibly removed by a periplasmic transport system. The core structure of the CatF subunit was closely related to seven other catalases with root-mean-square deviations (RMSDs) of 368 core Calpha atoms of 0.99-1.30 A. The heme component of CatF is heme b in the same orientation that is found in Escherichia coli hydroperoxidase II, an orientation that is flipped 180 degrees with respect the orientation of the heme in bovine liver catalase. NADPH is not found in the structure of CatF because key residues required for nucleotide binding are missing; 2129 water molecules were refined into the model. Water occupancy in the main or perpendicular channel of CatF varied among the four subunits from two to five in the region between the heme and the conserved Asp150. A comparison of the water occupancy in this region with the same region in other catalases reveals significant differences among the catalases.