Elucidating solvent contributions to solution reactions with ab initio QM/MM methods.

Computer simulations of reaction processes in solution in general rely on the definition of a reaction coordinate and the determination of the thermodynamic changes of the system along the reaction coordinate. The reaction coordinate often is constituted of characteristic geometrical properties of the reactive solute species, while the contributions of solvent molecules are implicitly included in the thermodynamics of the solute degrees of freedoms. However, solvent dynamics can provide the driving force for the reaction process, and in such cases explicit description of the solvent contribution in the free energy of the reaction process becomes necessary. We report here a method that can be used to analyze the solvent contributions to the reaction activation free energies from the combined QM/MM minimum free-energy path simulations. The method was applied to the self-exchange S(N)2 reaction of CH(3)Cl + Cl(-), showing that the importance of solvent-solute interactions to the reaction process. The results were further discussed in the context of coupling between solvent and solute molecules in reaction processes.

Field methods for capture and measurement of three monkey species in Costa Rica.

A total of 54 free-ranging monkeys were captured and marked in Santa Rosa National Park, Costa Rica, during May 1985, and an additional 17 were captured during March 1986. The animals were darted using a blowpipe or a CO2 gun. The drugs used were Ketaset, Sernylan and Telazol. Ketaset was effective for Cebus capucinus but unsuccessful for Alouatta palliata and Ateles geoffroyi. Sernylan was successful for A. geoffroyi and A. palliata but is no longer commercially available. Telazol proved to be an excellent alternative capture drug for both A. palliata and A. geoffroyi.

Diet-induced obesity alters the differentiation potential of stem cells isolated from bone marrow, adipose tissue and infrapatellar fat pad: the effects of free fatty acids.

INTRODUCTION: Obesity is a major risk factor for several musculoskeletal conditions that are characterized by an imbalance of tissue remodeling. Adult stem cells are closely associated with the remodeling and potential repair of several mesodermally derived tissues such as fat, bone and cartilage. We hypothesized that obesity would alter the frequency, proliferation, multipotency and immunophenotype of adult stem cells from a variety of tissues. MATERIALS AND METHODS: Bone marrow-derived mesenchymal stem cells (MSCs), subcutaneous adipose-derived stem cells (sqASCs) and infrapatellar fat pad-derived stem cells (IFP cells) were isolated from lean and high-fat diet-induced obese mice, and their cellular properties were examined. To test the hypothesis that changes in stem cell properties were due to the increased systemic levels of free fatty acids (FFAs), we further investigated the effects of FFAs on lean stem cells in vitro. RESULTS: Obese mice showed a trend toward increased prevalence of MSCs and sqASCs in the stromal tissues. While no significant differences in cell proliferation were observed in vitro, the differentiation potential of all types of stem cells was altered by obesity. MSCs from obese mice demonstrated decreased adipogenic, osteogenic and chondrogenic potential. Obese sqASCs and IFP cells showed increased adipogenic and osteogenic differentiation, but decreased chondrogenic ability. Obese MSCs also showed decreased CD105 and increased platelet-derived growth factor receptor α expression, consistent with decreased chondrogenic potential. FFA treatment of lean stem cells significantly altered their multipotency but did not completely recapitulate the properties of obese stem cells. CONCLUSIONS: These findings support the hypothesis that obesity alters the properties of adult stem cells in a manner that depends on the cell source. These effects may be regulated in part by increased levels of FFAs, but may involve other obesity-associated cytokines. These findings contribute to our understanding of mesenchymal tissue remodeling with obesity, as well as the development of autologous stem cell therapies for obese patients.

Comparison of extramedullary versus intramedullary referencing for tibial component alignment in total ankle arthroplasty.

BACKGROUND: The majority of total ankle arthroplasty (TAA) systems use extramedullary alignment guides for tibial component placement. However, at least 1 system offers intramedullary referencing. In total knee arthroplasty, studies suggest that tibial component placement is more accurate with intramedullary referencing. The purpose of this study was to compare the accuracy of extramedullary referencing with intramedullary referencing for tibial component placement in total ankle arthroplasty. METHODS: The coronal and sagittal tibial component alignment was evaluated on the postoperative weight-bearing anteroposterior (AP) and lateral radiographs of 236 consecutive fixed-bearing TAAs. Radiographs were measured blindly by 2 investigators. The postoperative alignment of the prosthesis was compared with the surgeon's intended alignment in both planes. The accuracy of tibial component alignment was compared between the extramedullary and intramedullary referencing techniques using unpaired t tests. Interrater and intrarater reliabilities were assessed with intraclass correlation coefficients (ICCs). RESULTS: Eighty-three tibial components placed with an extramedullary referencing technique were compared with 153 implants placed with an intramedullary referencing technique. The accuracy of the extramedullary referencing was within a mean of 1.5 ± 1.4 degrees and 4.1 ± 2.9 degrees in the coronal and sagittal planes, respectively. The accuracy of intramedullary referencing was within a mean of 1.4 ± 1.1 degrees and 2.5 ± 1.8 degrees in the coronal and sagittal planes, respectively. There was a significant difference (P < .001) between the 2 techniques with respect to the sagittal plane alignment. Interrater ICCs for coronal and sagittal alignment were high (0.81 and 0.94, respectively). Intrarater ICCs for coronal and sagittal alignment were high for both investigators. CONCLUSIONS: Initial sagittal plane tibial component alignment was notably more accurate when intramedullary referencing was used. Further studies are needed to determine the effect of this difference on clinical outcomes and long-term survivability of the implants. LEVEL OF EVIDENCE: Level III, retrospective comparative study.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

High-speed label-free functional photoacoustic microscopy of mouse brain in action.

We present fast functional photoacoustic microscopy (PAM) for three-dimensional high-resolution, high-speed imaging of the mouse brain, complementary to other imaging modalities. We implemented a single-wavelength pulse-width-based method with a one-dimensional imaging rate of 100 kHz to image blood oxygenation with capillary-level resolution. We applied PAM to image the vascular morphology, blood oxygenation, blood flow and oxygen metabolism in both resting and stimulated states in the mouse brain.

A New Method for Modeling Free Surface Flows and Fluid-structure Interaction with Ocean Applications

The computational modeling of ocean waves and ocean-faring devices poses numerous challenges. Among these are the need to stably and accurately represent both the fluid-fluid interface between water and air as well as the fluid-structure interfaces arising between solid devices and one or more fluids. As techniques are developed to stably and accurately balance the interactions between fluid and structural solvers at these boundaries, a similarly pressing challenge is the development of algorithms that are massively scalable and capable of performing large-scale three-dimensional simulations on reasonable time scales. This dissertation introduces two separate methods for approaching this problem, with the first focusing on the development of sophisticated fluid-fluid interface representations and the second focusing primarily on scalability and extensibility to higher-order methods.

We begin by introducing the narrow-band gradient-augmented level set method (GALSM) for incompressible multiphase Navier-Stokes flow. This is the first use of the high-order GALSM for a fluid flow application, and its reliability and accuracy in modeling ocean environments is tested extensively. The method demonstrates numerous advantages over the traditional level set method, among these a heightened conservation of fluid volume and the representation of subgrid structures.

Next, we present a finite-volume algorithm for solving the incompressible Euler equations in two and three dimensions in the presence of a flow-driven free surface and a dynamic rigid body. In this development, the chief concerns are efficiency, scalability, and extensibility (to higher-order and truly conservative methods). These priorities informed a number of important choices: The air phase is substituted by a pressure boundary condition in order to greatly reduce the size of the computational domain, a cut-cell finite-volume approach is chosen in order to minimize fluid volume loss and open the door to higher-order methods, and adaptive mesh refinement (AMR) is employed to focus computational effort and make large-scale 3D simulations possible. This algorithm is shown to produce robust and accurate results that are well-suited for the study of ocean waves and the development of wave energy conversion (WEC) devices.

Path Optimization in Free Energy Calculations

Free energy calculations are a computational method for determining thermodynamic quantities, such as free energies of binding, via simulation.

Currently, due to computational and algorithmic limitations, free energy calculations are limited in scope.

In this work, we propose two methods for improving the efficiency of free energy calculations.

First, we expand the state space of alchemical intermediates, and show that this expansion enables us to calculate free energies along lower variance paths.

We use Q-learning, a reinforcement learning technique, to discover and optimize paths at low computational cost.

Second, we reduce the cost of sampling along a given path by using sequential Monte Carlo samplers.

We develop a new free energy estimator, pCrooks (pairwise Crooks), a variant on the Crooks fluctuation theorem (CFT), which enables decomposition of the variance of the free energy estimate for discrete paths, while retaining beneficial characteristics of CFT.

Combining these two advancements, we show that for some test models, optimal expanded-space paths have a nearly 80% reduction in variance relative to the standard path.

Additionally, our free energy estimator converges at a more consistent rate and on average 1.8 times faster when we enable path searching, even when the cost of path discovery and refinement is considered.