In situ stress tensor measured in an Alaskan glacier

An experimental program at Worthington Glacier, Alaska, U.S.A., has yielded the first in situ measurements of the full stress tensor in glacier ice. Measurements were made with an array of stiff (low-compliance) normal-force sensors frozen into a borehole at 120 m depth. Freezing in temperate ice was accomplished by a down-hole heat exchanger which extracted heat at the rate of 15 W. Under slowly varying stress conditions, relaxation of stress anomalies by viscous creep following drilling of the hole and installation of the sensors allows for equilibration of measured stresses with far-field stresses. Equilibration of local and far-field stresses was confirmed and pressure sensor calibrated in laboratory experiments prior to the field program. Results of the stress measurements show principal axes of the stress tensor oriented in directions consistent with the geometry of the glacier and broadly consistent with measured englacial strain rate. The magnitudes of stress-tensor components are more error-prone and more sensitive to uncertainty in sensor magnitude than uncertainty in sensor orientation. Mean stress determined by pressure measurements agrees with estimated lithostatic overburden with within approximately 15%. Unexpected results include a stress perturbation lasted about 5 days that caused a rotation of the orientations of the principal stress axes of approximately 5 degrees about an axis pointing in the down-flow direction

Influence of conformational fluctuations on enzymatic activity: modelling the functional motion of beta-secretase

Considerable insight into the functional activity of proteins and enzymes can be obtained by studying the low-energy conformational distortions that the biopolymer can sustain. We carry out the characterization of these large scale structural changes for a protein of considerable pharmaceutical interest, the human $\beta$-secretase. Starting from the crystallographic structure of the protein, we use the recently introduced beta-Gaussian model to identify, with negligible computational expenditure, the most significant distortion occurring in thermal equilibrium and the associated time scales. The application of this strategy allows to gain considerable insight into the putative functional movements and, furthermore, helps to identify a handful of key regions in the protein which have an important mechanical influence on the enzymatic activity despite being spatially distant from the active site. The results obtained within the Gaussian model are validated through an extensive comparison against an all-atom Molecular Dynamics simulation.Comment: To be published in a special issue of J. Phys.: Cond. Mat. (Bedlewo Workshop

Using an Augmented Lagrangian Method and block fracturing in the DDA method

This paper presents two extensions to the Discontinuous Deformation Analysis (DDA) method orginally proposed by Shi for modeling the response of blocky rock masses to mechanical loading. The first extension consists of improving the block contact algorithm. An Augmented Lagrangian Method is used to replace the Penalty Method orginally proposed. It allows Lagrange multipliers to be introduced without increasing the number of equations that need to be solved and thus, block contract forces can be calculated more accurately. A block fracturing capability based on a three-parameter Mohr-Coulomb criterion represents the second extension. It allows for shear or tensile fracturing of intact blocks and the formation of smaller blocks

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

A modified empirical criterion for strength of transversely anisotropic rocks with metamorphic origin

A modified empirical criterion is proposed to determine the strength of transversely anisotropic rocks. In this regard, mechanical properties of intact anisotropic slate obtained from three different districts of Iran were taken into consideration. Afterward, triaxial rock strength criterion introduced by Rafiai was modified for transversely anisotropic rocks. The criterion was modified by adding a new parameter α for taking the influence of strength anisotropy into consideration. The results obtained have shown that the parameter α can be considered as the strength reduction parameter due to rock anisotropy. The modified criterion was compared to the modified Hoek–Brown (Saroglou and Tsiambaos) and Ramamurthy criteria for different anisotropic rocks. It was concluded that the criterion proposed in this paper is a more accurate and precise criterion in predicting the strength of anisotropic rocks

The Use of Experimental Structures to Model Protein Dynamics

The number of solved protein structures submitted in the Protein Data Bank (PDB) has increased dramatically in recent years. For some specific proteins, this number is very high—for example, there are over 550 solved structures for HIV-1 protease, one protein that is essential for the life cycle of human immunodeficiency virus (HIV) which causes acquired immunodeficiency syndrome (AIDS) in humans. The large number of structures for the same protein and its variants include a sample of different conformational states of the protein. A rich set of structures solved experimentally for the same protein has information buried within the dataset that can explain the functional dynamics and structural mechanism of the protein. To extract the dynamics information and functional mechanism from the experimental structures, this chapter focuses on two methods—Principal Component Analysis (PCA) and Elastic Network Models (ENM). PCA is a widely used statistical dimensionality reduction technique to classify and visualize high-dimensional data. On the other hand, ENMs are well-established simple biophysical method for modeling the functionally important global motions of proteins. This chapter covers the basics of these two. Moreover, an improved ENM version that utilizes the variations found within a given set of structures for a protein is described. As a practical example, we have extracted the functional dynamics and mechanism of HIV-1 protease dimeric structure by using a set of 329 PDB structures of this protein. We have described, step by step, how to select a set of protein structures, how to extract the needed information from the PDB files for PCA, how to extract the dynamics information using PCA, how to calculate ENM modes, how to measure the congruency between the dynamics computed from the principal components (PCs) and the ENM modes, and how to compute entropies using the PCs. We provide the computer programs or references to software tools to accomplish each step and show how to use these programs and tools. We also include computer programs to generate movies based on PCs and ENM modes and describe how to visualize them

Specialized dynamical properties of promiscuous residues revealed by simulated conformational ensembles

The ability to interact with different partners is one of the most important features in proteins. Proteins that bind a large number of partners (hubs) have been often associated with intrinsic disorder. However, many examples exist of hubs with an ordered structure, and evidence of a general mechanism promoting promiscuity in ordered proteins is still elusive. An intriguing hypothesis is that promiscuous binding sites have specific dynamical properties, distinct from the rest of the interface and pre-existing in the protein isolated state. Here, we present the first comprehensive study of the intrinsic dynamics of promiscuous residues in a large protein data set. Different computational methods, from coarse-grained elastic models to geometry-based sampling methods and to full-atom Molecular Dynamics simulations, were used to generate conformational ensembles for the isolated proteins. The flexibility and dynamic correlations of interface residues with a different degree of binding promiscuity were calculated and compared considering side chain and backbone motions, the latter both on a local and on a global scale. The study revealed that (a) promiscuous residues tend to be more flexible than nonpromiscuous ones, (b) this additional flexibility has a higher degree of organization, and (c) evolutionary conservation and binding promiscuity have opposite effects on intrinsic dynamics. Findings on simulated ensembles were also validated on ensembles of experimental structures extracted from the Protein Data Bank (PDB). Additionally, the low occurrence of single nucleotide polymorphisms observed for promiscuous residues indicated a tendency to preserve binding diversity at these positions. A case study on two ubiquitin-like proteins exemplifies how binding promiscuity in evolutionary related proteins can be modulated by the fine-tuning of the interface dynamics. The interplay between promiscuity and flexibility highlighted here can inspire new directions in protein-protein interaction prediction and design methods. © 2013 American Chemical Society

Engaging Undergraduates to Solve Global Health Challenges: A New Approach Based on Bioengineering Design

Recent reports have highlighted the need for educational programs to prepare students for careers developing and disseminating new interventions that improve global public health. Because of its multi-disciplinary, design-centered nature, the field of Biomedical Engineering can play an important role in meeting this challenge. This article describes a new program at Rice University to give undergraduate students from all disciplines a broad background in bioengineering and global health and provides an initial assessment of program impact. Working in partnership with health care providers in developing countries, students in the Beyond Traditional Borders (BTB) initiative learn about health challenges of the poor and put this knowledge to work immediately, using the engineering design process as a framework to formulate solutions to complex global health challenges. Beginning with a freshman design project and continuing through a capstone senior design course, the BTB curriculum uses challenges provided by partners in the developing world to teach students to integrate perspectives from multiple disciplines, and to develop leadership, communication, and teamwork skills. Exceptional students implement their designs under the guidance of clinicians through summer international internships. Since 2006, 333 students have designed more than 40 technologies and educational programs; 28 have been implemented in sub-Saharan Africa, Latin America, the Caribbean, southeast Asia, and the United States. More than 18,000 people have benefited from these designs. 95% of alumni who completed an international internship reported that participation in the program changed or strengthened their career plans to include a focus on global health medicine, research, and/or policy. Empowering students to use bioengineering design to address real problems is an effective way to teach the new generation of leaders needed to solve global health challenges

Recurrence quantification analysis as a tool for the characterization of molecular dynamics simulations

A molecular dynamics simulation of a Lennard-Jones fluid, and a trajectory of the B1 immunoglobulin G-binding domain of streptococcal protein G (B1-IgG) simulated in water are analyzed by recurrence quantification, which is noteworthy for its independence from stationarity constraints, as well as its ability to detect transients, and both linear and nonlinear state changes. The results demonstrate the sensitivity of the technique for the discrimination of phase sensitive dynamics. Physical interpretation of the recurrence measures is also discussed.Comment: 7 pages, 8 figures, revtex; revised for review for Phys. Rev. E (clarifications and expansion of discussion)-- addition of the 8 postscript figures previously omitted, but unchanged from version