Variational Monte Carlo study of the hydrogen electronic structure at ultra-high pressure

Die causa finalis der vorliegenden Arbeit ist das Verständnis des Phasendiagramms von Wasserstoff bei ultrahohen Drücken, welche von nichtleitendem H2 bis hin zu metallischem H reichen. Da die Voraussetzungen für ultrahohen Druck im Labor schwer zu schaffen sind, bilden Computersimulationen ein wichtiges alternatives Untersuchungsinstrument. Allerdings sind solche Berechnungen eine große Herausforderung. Eines der größten Probleme ist die genaue Auswertung des Born-Oppenheimer Potentials, welches sowohl für die nichtleitende als auch für die metallische Phase geeignet sein muss. Außerdem muss es die starken Korrelationen berücksichtigen, die durch die kovalenten H2 Bindungen und die eventuellen Phasenübergänge hervorgerufen werden. Auf dieses Problem haben unsere Anstrengungen abgezielt. Im Kontext von Variationellem Monte Carlo (VMC) ist die Shadow Wave Function (SWF) eine sehr vielversprechende Option. Aufgrund ihrer Flexibilität sowohl lokalisierte als auch delokalisierte Systeme zu beschreiben sowie ihrer Fähigkeit Korrelationen hoher Ordnung zu berücksichtigen, ist sie ein idealer Kandidat für unsere Zwecke. Unglücklicherweise bringt ihre Formulierung ein Vorzeichenproblem mit sich, was die Anwendbarkeit limitiert. Nichtsdestotrotz ist es möglich diese Schwierigkeit zu umgehen indem man die Knotenstruktur a priori festlegt. Durch diesen Formalismus waren wir in der Lage die Beschreibung der Elektronenstruktur von Wasserstoff signifikant zu verbessern, was eine sehr vielversprechende Perspektive bietet. Während dieser Forschung haben wir also die Natur des Vorzeichenproblems untersucht, das sich auf die SWF auswirkt, und dabei ein tieferes Verständnis seines Ursprungs erlangt. Die vorliegende Arbeit ist in vier Kapitel unterteilt. Das erste Kapitel führt VMC und die SWF mit besonderer Ausrichtung auf fermionische Systeme ein. Kapitel 2 skizziert die Literatur über das Phasendiagramm von Wasserstoff bei ultrahohem Druck. Das dritte Kapitel präsentiert die Implementierungen unseres VMC Programms und die erhaltenen Ergebnisse. Zum Abschluss fasst Kapitel 4 unsere Bestrebungen zur Lösung des zur SWF zugehörigen Vorzeichenproblems zusammen.

Numerical analysis of the relation between interactions and structure in a molecular fluid

Coarse graining is a popular technique used in physics to speed up the computer simulation of molecular fluids. An essential part of this technique is a method that solves the inverse problem of determining the interaction potential or its parameters from the given structural data. Due to discrepancies between model and reality, the potential is not unique, such that stability of such method and its convergence to a meaningful solution are issues.rnrnIn this work, we investigate empirically whether coarse graining can be improved by applying the theory of inverse problems from applied mathematics. In particular, we use the singular value analysis to reveal the weak interaction parameters, that have a negligible influence on the structure of the fluid and which cause non-uniqueness of the solution. Further, we apply a regularizing Levenberg-Marquardt method, which is stable against the mentioned discrepancies. Then, we compare it to the existing physical methods - the Iterative Boltzmann Inversion and the Inverse Monte Carlo method, which are fast and well adapted to the problem, but sometimes have convergence problems.rnrnFrom analysis of the Iterative Boltzmann Inversion, we elaborate a meaningful approximation of the structure and use it to derive a modification of the Levenberg-Marquardt method. We engage the latter for reconstruction of the interaction parameters from experimental data for liquid argon and nitrogen. We show that the modified method is stable, convergent and fast. Further, the singular value analysis of the structure and its approximation allows to determine the crucial interaction parameters, that is, to simplify the modeling of interactions. Therefore, our results build a rigorous bridge between the inverse problem from physics and the powerful solution tools from mathematics. rn

Factor Structure of the Neck Disability Index

BACKGROUND CONTEXT: The Neck Disability Index frequently is used to measure outcomes of the neck. The statistical rigor of the Neck Disability Index has been assessed with conflicting outcomes. To date, Confirmatory Factor Analysis of the Neck Disability Index has not been reported for a suitably large population study. Because the Neck Disability Index is not a condition-specific measure of neck function, initial Confirmatory Factor Analysis should consider problematic neck patients as a homogenous group. PURPOSE: We sought to analyze the factor structure of the Neck Disability Index through Confirmatory Factor Analysis in a symptomatic, homogeneous, neck population, with respect to pooled populations and gender subgroups. STUDY DESIGN: This was a secondary analysis of pooled data. PATIENT SAMPLE: A total of 1,278 symptomatic neck patients (67.5% female, median age 41 years), 803 nonspecific and 475 with whiplash-associated disorder. OUTCOME MEASURES: The Neck Disability Index was used to measure outcomes. METHODS: We analyzed pooled baseline data from six independent studies of patients with neck problems who completed Neck Disability Index questionnaires at baseline. The Confirmatory Factor Analysis was considered in three scenarios: the full sample and separate sexes. Models were compared empirically for best fit. RESULTS: Two-factor models have good psychometric properties across both the pooled and sex subgroups. However, according to these analyses, the one-factor solution is preferable from both a statistical perspective and parsimony. The two-factor model was close to significant for the male subgroup (p<.07) where questions separated into constructs of mental function (pain, reading headaches and concentration) and physical function (personal care, lifting, work, driving, sleep, and recreation). CONCLUSIONS: The Neck Disability Index demonstrated a one-factor structure when analyzed by Confirmatory Factor Analysis in a pooled, homogenous sample of neck problem patients. However, a two-factor model did approach significance for male subjects where questions separated into constructs of mental and physical function. Further investigations in different conditions, subgroup and sex-specific populations are warranted.

Structure determination of channel and transport proteins by high-resolution microscopy techniques

High-resolution microscopy techniques provide a plethora of information on biological structures from the cellular level down to the molecular level. In this review, we present the unique capabilities of transmission electron and atomic force microscopy to assess the structure, oligomeric state, function and dynamics of channel and transport proteins in their native environment, the lipid bilayer. Most importantly, membrane proteins can be visualized in the frozen-hydrated state and in buffer solution by cryo-transmission electron and atomic force microscopy, respectively. We also illustrate the potential of the scintillation proximity assay to study substrate binding of detergent-solubilized transporters prior to crystallization and structural characterization.

Fine structure of hippocampal mossy fiber synapses following rapid high-pressure freezing

Synapses of hippocampal neurons play important roles in learning and memory processes and are involved in aberrant hippocampal function in temporal lobe epilepsy. Major neuronal types in the hippocampus as well as their input and output synapses are well known, but it has remained an open question to what extent conventional electron microscopy (EM) has provided us with the real appearance of synaptic fine structure under in vivo conditions. There is reason to assume that conventional aldehyde fixation and dehydration lead to protein denaturation and tissue shrinkage, likely associated with the occurrence of artifacts. However, realistic fine-structural data of synapses are required for our understanding of the transmission process and for its simulation. Here, we used high-pressure freezing and cryosubstitution of hippocampal tissue that was not subjected to aldehyde fixation and dehydration in ethanol to monitor the fine structure of an identified synapse in the hippocampal CA3 region, that is, the synapse between granule cell axons, the mossy fibers, and the proximal dendrites of CA3 pyramidal neurons. Our results showed that high-pressure freezing nicely preserved ultrastructural detail of this particular synapse and allowed us to study rapid structural changes associated with synaptic plasticity.

Structural basis for the sheddase function of human meprin β metalloproteinase at the plasma membrane

Ectodomain shedding at the cell surface is a major mechanism to regulate the extracellular and circulatory concentration or the activities of signaling proteins at the plasma membrane. Human meprin β is a 145-kDa disulfide-linked homodimeric multidomain type-I membrane metallopeptidase that sheds membrane-bound cytokines and growth factors, thereby contributing to inflammatory diseases, angiogenesis, and tumor progression. In addition, it cleaves amyloid precursor protein (APP) at the β-secretase site, giving rise to amyloidogenic peptides. We have solved the X-ray crystal structure of a major fragment of the meprin β ectoprotein, the first of a multidomain oligomeric transmembrane sheddase, and of its zymogen. The meprin β dimer displays a compact shape, whose catalytic domain undergoes major rearrangement upon activation, and reveals an exosite and a sugar-rich channel, both of which possibly engage in substrate binding. A plausible structure-derived working mechanism suggests that substrates such as APP are shed close to the plasma membrane surface following an "N-like" chain trace.

Genetic Variation and Population Structure in Native Americans

We examined genetic diversity and population structure in the American landmass using 678 autosomal microsatellite markers genotyped in 422 individuals representing 24 Native American populations sampled from North, Central, and South America. These data were analyzed jointly with similar data available in 54 other indigenous populations worldwide, including an additional five Native American groups. The Native American populations have lower genetic diversity and greater differentiation than populations from other continental regions. We observe gradients both of decreasing genetic diversity as a function of geographic distance from the Bering Strait and of decreasing genetic similarity to Siberians-signals of the southward dispersal of human populations from the northwestern tip of the Americas. We also observe evidence of: (1) a higher level of diversity and lower level of population structure in western South America compared to eastern South America, (2) a relative lack of differentiation between Mesoamerican and Andean populations, (3) a scenario in which coastal routes were easier for migrating peoples to traverse in comparison with inland routes, and (4) a partial agreement on a local scale between genetic similarity and the linguistic classification of populations. These findings offer new insights into the process of population dispersal and differentiation during the peopling of the Americas.

Hydrogel-based engineered skeletal muscle grafts normalize heart function early after myocardial infarction

Tissue engineering represents an attractive approach for the treatment of congestive heart failure. The influence of the differentiation of myogenic graft for functional recovery is not defined. We engineered a biodegradable skeletal muscle graft (ESMG) tissue and investigated its functional effect after implantation on the epicardium of an infarcted heart segment. ESMGs were synthesized by mixing collagen (2 mg/mL), Matrigel (2 mg/mL), and rat skeletal muscle cells (10(6)). Qualitative and quantitative aspects of ESMGs were optimized. Two weeks following coronary ligation, the animals were randomized in three groups: ESMG glued to the epicardial surface with fibrin (ESMG, n = 7), fibrin alone (fibrin, n = 5), or sham operation (sham, n = 4). Echocardiography, histology, and immunostaining were performed 4 weeks later. A cohesive three-dimensional tissular structure formed in vitro within 1 week. Myoblasts differentiated into randomly oriented myotubes. Four weeks postimplantation, ESMGs were vascularized and invaded by granulation tissue. Mean fractional shortening (FS) was, however, significantly increased in the ESMG group as compared with preimplantation values (42 +/- 6 vs. 33 +/- 5%, P < 0.05) and reached the values of controlled noninfarcted animals (control, n = 5; 45 +/- 3%; not significant). Pre- and postimplantation FS did not change over these 4 weeks in the sham group and the fibrin-treated animals. This study showed that it is possible to improve systolic heart function following myocardial infarction through implantation of differentiated muscle fibers seeded on a gel-type scaffold despite a low rate of survival.

Vertebroplasty: experimental characterization of polymethylmethacrylate bone cement spreading as a function of viscosity, bone porosity, and flow rate

STUDY DESIGN: This is an experimental study on an artificial vertebra model and human cadaveric spine. OBJECTIVE: Characterization of polymethylmethacrylate (PMMA) bone cement distribution in the vertebral body as a function of cement viscosity, bone porosity, and injection speed. Identification of relevant parameters for improved cement flow predictability and leak prevention in vertebroplasty. SUMMARY OF BACKGROUND DATA: Vertebroplasty is an efficient procedure to treat vertebral fractures and stabilize osteoporotic bone in the spine. Severe complications result from bone cement leakage into the spinal canal or the vascular system. Cement viscosity has been identified as an important parameter for leak prevention but the influence of bone structure and injection speed remain obscure. METHODS: An artificial vertebra model based on open porous aluminum foam was used to simulate bone of known porosity. Fifty-six vertebroplasties with 4 different starting viscosity levels and 2 different injection speeds were performed on artificial vertebrae of 3 different porosities. A validation on a human cadaveric spine was executed. The experiments were radiographically monitored and the shape of the cement clouds quantitatively described with the 2 indicators circularity and mean cement spreading distance. RESULTS: An increase in circularity and a decrease in mean cement spreading distance was observed with increasing viscosity, with the most striking change occurring between 50 and 100 Pas. Larger pores resulted in significantly reduced circularity and increased mean cement spreading distance whereas the effect of injection speed on the 2 indicators was not significant. CONCLUSION: Viscosity is the key factor for reducing the risk of PMMA cement leakage and it should be adapted to the degree of osteoporosis encountered in each patient. It may be advisable to opt for a higher starting viscosity but to inject the material at a faster rate.

Investigating within-canopy variation of functional traits and cellular structure of sugar maple (Acer saccharum) leaves

Patterns of increasing leaf mass per area (LMA), area-based leaf nitrogen (Narea), and carbon isotope composition (δ13C) with increasing height in the canopy have been attributed to light gradients or hydraulic limitation in tall trees. Theoretical optimal distributions of LMA and Narea that scale with light maximize canopy photosynthesis; however, sub-optimal distributions are often observed due to hydraulic constraints on leaf development. Using observational, experimental, and modeling approaches, we investigated the response of leaf functional traits (LMA, density, thickness, and leaf nitrogen), leaf carbon isotope composition (δ13C), and cellular structure to light availability, height, and leaf water potential (Ψl) in an Acer saccharum forest to tease apart the influence of light and hydraulic limitations. LMA, leaf and palisade layer thickness, and leaf density were greater at greater light availability but similar heights, highlighting the strong control of light on leaf morphology and cellular structure. Experimental shading decreased both LMA and area-based leaf nitrogen (Narea) and revealed that LMA and Narea were more strongly correlated with height earlier in the growing season and with light later in the growing season. The supply of CO2 to leaves at higher heights appeared to be constrained by stomatal sensitivity to vapor pressure deficit (VPD) or midday leaf water potential, as indicated by increasing δ13C and VPD and decreasing midday Ψl with height. Model simulations showed that daily canopy photosynthesis was biased during the early growing season when seasonality was not accounted for, and was biased throughout the growing season when vertical gradients in LMA and Narea were not accounted for. Overall, our results suggest that leaves acclimate to light soon after leaf expansion, through an accumulation of leaf carbon, thickening of palisade layers and increased LMA, and reduction in stomatal sensitivity to Ψl or VPD. This period of light acclimation in leaves appears to optimize leaf function over time, despite height-related constraints early in the growing season. Our results imply that vertical gradients in leaf functional traits and leaf acclimation to light should be incorporated in canopy function models in order to refine estimates of canopy photosynthesis.

An estimate of heavy quark momentum diffusion coefficient in gluon plasma

We calculate the momentum diffusion coefficient for heavy quarks in SU(3) gluon plasma at temperatures 1-2 times the deconfinement temperature. The momentum diffusion coefficient is extracted from a Monte Carlo calculation of the correlation function of color electric fields, in the leading order of expansion in heavy quark mass. Systematics of the calculation are examined, and compared with perturbtion theory and other estimates.

Structure and superorganization of acetylcholine receptor-rapsyn complexes

The scaffolding protein at the neuromuscular junction, rapsyn, enables clustering of nicotinic acetylcholine receptors in high concentration and is critical for muscle function. Patients with insufficient receptor clustering suffer from muscle weakness. However, the detailed organization of the receptor-rapsyn network is poorly understood: it is unclear whether rapsyn first forms a wide meshwork to which receptors can subsequently dock or whether it only forms short bridges linking receptors together to make a large cluster. Furthermore, the number of rapsyn-binding sites per receptor (a heteropentamer) has been controversial. Here, we show by cryoelectron tomography and subtomogram averaging of Torpedo postsynaptic membrane that receptors are connected by up to three rapsyn bridges, the minimum number required to form a 2D network. Half of the receptors belong to rapsyn-connected groups comprising between two and fourteen receptors. Our results provide a structural basis for explaining the stability and low diffusion of receptors within clusters.

Combined effect of heterogeneity, anisotropy and saturation on steady state flow and transport: Structure recognition and numerical simulation

1 Natural soil profiles may be interpreted as an arrangement of parts which are characterized by properties like hydraulic conductivity and water retention function. These parts form a complicated structure. Characterizing the soil structure is fundamental in subsurface hydrology because it has a crucial influence on flow and transport and defines the patterns of many ecological processes. We applied an image analysis method for recognition and classification of visual soil attributes in order to model flow and transport through a man-made soil profile. Modeled and measured saturation-dependent effective parameters were compared. We found that characterizing and describing conductivity patterns in soils with sharp conductivity contrasts is feasible. Differently, solving flow and transport on the basis of these conductivity maps is difficult and, in general, requires special care for representation of small-scale processes.

A loss-of-function variant of PTPN22 is associated with reduced risk of systemic lupus erythematosus.

A gain-of-function R620W polymorphism in the PTPN22 gene, encoding the lymphoid tyrosine phosphatase LYP, has recently emerged as an important risk factor for human autoimmunity. Here we report that another missense substitution (R263Q) within the catalytic domain of LYP leads to reduced phosphatase activity. High-resolution structural analysis revealed the molecular basis for this loss of function. Furthermore, the Q263 variant conferred protection against human systemic lupus erythematosus, reinforcing the proposal that inhibition of LYP activity could be beneficial in human autoimmunity.

Plasticity of lipid-protein interactions in the function and topogenesis of the membrane protein lactose permease from Escherichia coli.

Phosphatidylcholine (PC) has been widely used in place of naturally occurring phosphatidylethanolamine (PE) in reconstitution of bacterial membrane proteins. However, PC does not support native structure or function for several reconstituted transport proteins. Lactose permease (LacY) of Escherichia coli, when reconstituted in E. coli phospholipids, exhibits energy-dependent uphill and energy-independent downhill transport function and proper conformation of periplasmic domain P7, which is tightly linked to uphill transport function. LacY expressed in cells lacking PE and containing only anionic phospholipids exhibits only downhill transport and lacks native P7 conformation. Reconstitution of LacY in the presence of E. coli-derived PE, but not dioleoyl-PC, results in uphill transport. We now show that LacY exhibits uphill transport and native conformation of P7 when expressed in a mutant of E. coli in which PC completely replaces PE even though the structure is not completely native. E. coli-derived PC and synthetic PC species containing at least one saturated fatty acid also support the native conformation of P7 dependent on the presence of anionic phospholipids. Our results demonstrate that the different effects of PE and PC species on LacY structure and function cannot be explained by differences in the direct interaction of the lipid head groups with specific amino acid residues alone but are due to more complex effects of the physical and chemical properties of the lipid environment on protein structure. This conclusion is supported by the effect of different lipids on the proper folding of domain P7, which indirectly influences uphill transport function.