Estimation of Hydrogen-Exchange Protection Factors from MD Simulation Based on Amide Hydrogen Bonding Analysis.

Hydrogen exchange (HX) studies have provided critical insight into our understanding of protein folding, structure, and dynamics. More recently, hydrogen exchange mass spectrometry (HX-MS) has become a widely applicable tool for HX studies. The interpretation of the wealth of data generated by HX-MS experiments as well as other HX methods would greatly benefit from the availability of exchange predictions derived from structures or models for comparison with experiment. Most reported computational HX modeling studies have employed solvent-accessible-surface-area based metrics in attempts to interpret HX data on the basis of structures or models. In this study, a computational HX-MS prediction method based on classification of the amide hydrogen bonding modes mimicking the local unfolding model is demonstrated. Analysis of the NH bonding configurations from molecular dynamics (MD) simulation snapshots is used to determine partitioning over bonded and nonbonded NH states and is directly mapped into a protection factor (PF) using a logistics growth function. Predicted PFs are then used for calculating deuteration values of peptides and compared with experimental data. Hydrogen exchange MS data for fatty acid synthase thioesterase (FAS-TE) collected for a range of pHs and temperatures was used for detailed evaluation of the approach. High correlation between prediction and experiment for observable fragment peptides is observed in the FAS-TE and additional benchmarking systems that included various apo/holo proteins for which literature data were available. In addition, it is shown that HX modeling can improve experimental resolution through decomposition of in-exchange curves into rate classes, which correlate with prediction from MD. Successful rate class decompositions provide further evidence that the presented approach captures the underlying physical processes correctly at the single residue level. This assessment is further strengthened in a comparison of residue resolved protection factor predictions for staphylococcal nuclease with NMR data, which was also used to compare prediction performance with other algorithms described in the literature. The demonstrated transferable and scalable MD based HX prediction approach adds significantly to the available tools for HX-MS data interpretation based on available structures and models.

Structural and functional characterization of BaiA, an enzyme involved in secondary bile acid synthesis in human gut microbe.

Despite significant influence of secondary bile acids on human health and disease, limited structural and biochemical information is available for the key gut microbial enzymes catalyzing its synthesis. Herein, we report apo- and cofactor bound crystal structures of BaiA2, a short chain dehydrogenase/reductase from Clostridium scindens VPI 12708 that represent the first protein structure of this pathway. The structures elucidated the basis of cofactor specificity and mechanism of proton relay. A conformational restriction involving Glu42 located in the cofactor binding site seems crucial in determining cofactor specificity. Limited flexibility of Glu42 results in imminent steric and electrostatic hindrance with 2'-phosphate group of NADP(H). Consistent with crystal structures, steady state kinetic characterization performed with both BaiA2 and BaiA1, a close homolog with 92% sequence identity, revealed specificity constant (kcat /KM ) of NADP(+) at least an order of magnitude lower than NAD(+) . Substitution of Glu42 with Ala improved specificity toward NADP(+) by 10-fold compared to wild type. The cofactor bound structure uncovered a novel nicotinamide-hydroxyl ion (NAD(+) -OH(-) ) adduct contraposing previously reported adducts. The OH(-) of the adduct in BaiA2 is distal to C4 atom of nicotinamide and proximal to 2'-hydroxyl group of the ribose moiety. Moreover, it is located at intermediary distances between terminal functional groups of active site residues Tyr157 (2.7 Å) and Lys161 (4.5 Å). Based on these observations, we propose an involvement of NAD(+) -OH(-) adduct in proton relay instead of hydride transfer as noted for previous adducts.

Mapping conformational dynamics of proteins using torsional dynamics simulations.

All-atom molecular dynamics simulations are widely used to study the flexibility of protein conformations. However, enhanced sampling techniques are required for simulating protein dynamics that occur on the millisecond timescale. In this work, we show that torsional molecular dynamics simulations enhance protein conformational sampling by performing conformational search in the low-frequency torsional degrees of freedom. In this article, we use our recently developed torsional-dynamics method called Generalized Newton-Euler Inverse Mass Operator (GNEIMO) to study the conformational dynamics of four proteins. We investigate the use of the GNEIMO method in simulations of the conformationally flexible proteins fasciculin and calmodulin, as well as the less flexible crambin and bovine pancreatic trypsin inhibitor. For the latter two proteins, the GNEIMO simulations with an implicit-solvent model reproduced the average protein structural fluctuations and sample conformations similar to those from Cartesian simulations with explicit solvent. The application of GNEIMO with replica exchange to the study of fasciculin conformational dynamics produced sampling of two of this protein's experimentally established conformational substates. Conformational transition of calmodulin from the Ca(2+)-bound to the Ca(2+)-free conformation occurred readily with GNEIMO simulations. Moreover, the GNEIMO method generated an ensemble of conformations that satisfy about half of both short- and long-range interresidue distances obtained from NMR structures of holo to apo transitions in calmodulin. Although unconstrained all-atom Cartesian simulations have failed to sample transitions between the substates of fasciculin and calmodulin, GNEIMO simulations show the transitions in both systems. The relatively short simulation times required to capture these long-timescale conformational dynamics indicate that GNEIMO is a promising molecular-dynamics technique for studying domain motion in proteins.

Design, synthesis and biological studies of survivin dimerization modulators that prolong mitotic cycle.

Survivin, a member of the inhibitor of apoptosis protein (IAP) family proteins, has essential roles in cell division and inhibition of apoptosis. Several clinical studies in cancer patients have shown that the elevated levels of survivin correlate with aggressiveness of the disease and resistance to radiation and chemotherapeutic treatments. Survivin is an integral component of chromosomal passenger complex (CPC) where it binds to borealin and INCENP through its dimerization interface. Thus, disruption of functional survivin along its dimer interface with a small molecule is hypothesized to inhibit the proliferation of cancer cells and sensitize them to therapeutic agents and radiation. Recently, a small molecule (Abbott8) was reported to bind at the dimerization interface of survivin. Further development of this compound was accomplished by computational modeling of the molecular interactions along the dimerization interface, which has led to the design of promising survivin dimerization modulators. Two of the most potent survivin modulators, LLP3 and LLP9 at concentrations between 50 and 100nM, caused delay in mitotic progression and major mitotic defects in proliferating human umbilical vein endothelial cells (HUVEC) and prostate cancer cells (PC3).

Simultaneous binding of two peptidyl ligands by a SRC homology 2 domain.

Src homology 2 (SH2) domains mediate protein-protein interactions by recognizing phosphotyrosine (pY)-containing sequences of target proteins. In all of the SH2 domain-pY peptide interactions described to date, the SH2 domain binds to a single pY peptide. Here, determination of the cocrystal structure of the N-terminal SH2 domain of phosphatase SHP-2 bound to a class IV peptide (VIpYFVP) revealed a noncanonical 1:2 (protein-peptide) complex. The first peptide binds in a canonical manner with its pY side chain inserted in the usual binding pocket, while the second pairs up with the first to form two antiparallel ß-strands that extend the central ß-sheet of the SH2 domain. This unprecedented binding mode was confirmed in the solution phase by NMR experiments and shown to be adopted by pY peptides derived from cellular proteins. Site-directed mutagenesis and surface plasmon resonance studies revealed that the binding of the first peptide is pY-dependent, but phosphorylation is not required for the second peptide. Our findings suggest a potential new function for the SH2 domain as a molecular clamp to promote dimerization of signaling proteins.

Characterization of molecular recognition of STAT3 SH2 domain inhibitors through molecular simulation.

Signal transducer and activator of transcription 3 (STAT3) is an anti-cancer target protein due to its over-activation in tumor cells. The Tyr705-phosphorylated (pTyr) STAT3 binds to the pTyr-recognition site of its Src Homology 2 (SH2) domain of another STAT3 monomer to form a homo-dimer, which then causes cellular anti-apoptosis, proliferation, and tumor invasion. Recently, many STAT3 SH2 dimerization inhibitors have been discovered via both computational and experimental methods. To systematically assess their binding affinities and specificities, for eight representative inhibitors, we utilized molecular docking, molecular dynamics simulation, and ensuing energetic analysis to compare their binding characteristics. The inhibitors' binding free energies were calculated via MMPB(GB)SA, and the STAT3 SH2 binding "hot spots" were evaluated through binding energy decomposition and hydrogen bond (H-bond) distribution analysis. Several conclusions can be drawn: (1) the overall enthalpy-entropy compensation paradigm is preserved for the STAT3 SH2/ligand binding thermodynamics; (2) at one end of the binding spectrum, two compounds bind to SH2 due to their minimum entropic penalties that result from their relative rigidities and increased dynamics of SH2 upon their binding; at the other end of the binding spectrum, one compound shows a typical weak binder behavior due to its loose binding in the SH2's strongest enthalpy-contributing binding subsite; (3) hydrogen bonding seems a strong indicator to evaluate the SH2/ligand binding potency, which echoes a finding that CH/π non-classical H-bond is responsible for some pTyr peptides binding to their corresponding SH2 domains; (4) STAT3 SH2 domain possesses three binding "hot spots": pTyr705-binding pocket with polar residues and contributing the largest binding enthalpy (two-thirds); Leu706 subsite which is the most dynamic and hardest to target; a hydrophobic side pocket which is unique to STAT3 and very targetable, which may offer unique opportunity to design STAT3-specific inhibitors, particularly with fragment-based approach.

First Report of Two Colletotrichum Species Associated with Bitter Rot on Apple Fruit in Korea - C. fructicola and C. siamense.

Bitter rot caused by the fungal genus Colletotrichum is a well-known, common disease of apple and causes significant yield loss. In 2013, six fungal strains were isolated from Fuji apple fruits exhibiting symptoms of bitter rot from Andong, Korea. These strains were identified as Colletotrichum fructicola and C. siamense based on morphological characteristics and multilocus sequence analysis of the internal transcribed spacer rDNA, actin, calmodulin, chitin synthase, and glyceraldehyde-3-phosphate dehydrogenase Pathogenicity tests confirmed the involvement of C. fructicola and C. siamense in the development of disease symptoms on apple fruits. This is the first report of C. fructicola and C. siamense causing bitter rot on apple fruit in Korea.

Phenomenological Study on Burnout Experience of Clinical Nurses Who have Turnover Intention.

PURPOSE: This study was aimed to understand the meaning and essentials of the experience of burnout for hospital nurses with turnover intention. METHODS: The design was a qualitative research of phenomenological study. Participants: Seven hospital nurses who had worked over three years and had experiences of turnover intention in a hospital with over 400 beds were included. RESULTS: Nine meaningful themes related to burnout experiences and four theme clusters of 1) battery warning sounds almost out; 2) the player who hit the drum and double-headed drum; 3) the target flying arrow without a break; and 4) the pendulum swaying to turn over. Registered nurses (RNs) felt burnout with an overload of work and by the thought that it was illegal action for registered nurses to receive insufficient rewards for their work. RNs also experienced there were no problem solving strategies to verbal violence by patient and medical team. CONCLUSION: The findings show that burnout experiences for those who had turnover intention was developed from the insight that insufficient training to do work independently with over-load for nurses was not ethical. It suggests that it is necessary to rethink training systems for nursing and hospitals to relieve turnover intention.

A collaborative visual analytics suite for protein folding research.

Molecular dynamics (MD) simulation is a crucial tool for understanding principles behind important biochemical processes such as protein folding and molecular interaction. With the rapidly increasing power of modern computers, large-scale MD simulation experiments can be performed regularly, generating huge amounts of MD data. An important question is how to analyze and interpret such massive and complex data. One of the (many) challenges involved in analyzing MD simulation data computationally is the high-dimensionality of such data. Given a massive collection of molecular conformations, researchers typically need to rely on their expertise and prior domain knowledge in order to retrieve certain conformations of interest. It is not easy to make and test hypotheses as the data set as a whole is somewhat "invisible" due to its high dimensionality. In other words, it is hard to directly access and examine individual conformations from a sea of molecular structures, and to further explore the entire data set. There is also no easy and convenient way to obtain a global view of the data or its various modalities of biochemical information. To this end, we present an interactive, collaborative visual analytics tool for exploring massive, high-dimensional molecular dynamics simulation data sets. The most important utility of our tool is to provide a platform where researchers can easily and effectively navigate through the otherwise "invisible" simulation data sets, exploring and examining molecular conformations both as a whole and at individual levels. The visualization is based on the concept of a topological landscape, which is a 2D terrain metaphor preserving certain topological and geometric properties of the high dimensional protein energy landscape. In addition to facilitating easy exploration of conformations, this 2D terrain metaphor also provides a platform where researchers can visualize and analyze various properties (such as contact density) overlayed on the top of the 2D terrain. Finally, the software provides a collaborative environment where multiple researchers can assemble observations and biochemical events into storyboards and share them in real time over the Internet via a client-server architecture. The software is written in Scala and runs on the cross-platform Java Virtual Machine. Binaries and source code are available at http://www.aylasoftware.org and have been released under the GNU General Public License.

Impairment of glioma stem cell survival and growth by a novel inhibitor for Survivin-Ran protein complex.

PURPOSE: Glioblastoma multiforme (GBM) is a devastating disease. Recent studies suggest that the stem cell properties of GBM contribute to the development of therapy resistance. EXPERIMENTAL DESIGN: The expression of Survivin and Ran was evaluated by immunohistochemistry with GBM tissues, and quantitative reverse transcriptase (qRT)-PCR and immunocytochemistry with patient-derived GBM sphere cultures. With a computational structure-based drug design, 11 small-molecule compounds were designed, synthesized, and evaluated as inhibitor candidates for the molecular interaction of Survivin protein. The molecular mechanism of the lead compound, LLP-3, was determined by Western blot, ELISA, in situ proximity ligation assay, and immunocytochemistry. The effects of LLP-3 treatment on GSCs were evaluated both in vitro and in vivo. Quantitative immunohistochemistry was carried out to compare Survivin expression in tissues from 44 newly diagnosed and 31 recurrent post-chemoradiation GBM patients. Lastly, the sensitivities of temozolomide-resistant GBM spheres to LLP-3 were evaluated in vitro. RESULTS: Survivin and Ran were strongly expressed in GBM tissues, particularly in the perivasculature, and also in patient-derived GSC cultures. LLP-3 treatment disrupted the Survivin-Ran protein complex in cancer cells and abolished the growth of patient-derived GBM spheres in vitro and in vivo. This inhibition was dependent on caspase activity and associated with p53 status of cells. Immunohistochemistry showed that Survivin expression is significantly increased in recurrent GBM compared with newly diagnosed tumors, and temozolomide-resistant GBM spheres exhibited high sensitivities to LLP-3 treatment. CONCLUSIONS: Disruption of the Survivin-Ran complex by LLP-3 abolishes survival and growth of GSCs both in vitro and in vivo, indicating an attractive novel therapeutic approach for GBM.

Equipartition Principle for Internal Coordinate Molecular Dynamics.

The principle of equipartition of (kinetic) energy for all-atom Cartesian molecular dynamics states that each momentum phase space coordinate on the average has ½kT of kinetic energy in a canonical ensemble. This principle is used in molecular dynamics simulations to initialize velocities, and to calculate statistical properties such as entropy. Internal coordinate molecular dynamics (ICMD) models differ from Cartesian models in that the overall kinetic energy depends on the generalized coordinates and includes cross-terms. Due to this coupled structure, no such equipartition principle holds for ICMD models. In this paper we introduce non-canonical modal coordinates to recover some of the structural simplicity of Cartesian models and develop a new equipartition principle for ICMD models. We derive low-order recursive computational algorithms for transforming between the modal and physical coordinates. The equipartition principle in modal coordinates provides a rigorous method for initializing velocities in ICMD simulations thus replacing the ad hoc methods used until now. It also sets the basis for calculating conformational entropy using internal coordinates.

Structure refinement of protein low resolution models using the GNEIMO constrained dynamics method.

The challenge in protein structure prediction using homology modeling is the lack of reliable methods to refine the low resolution homology models. Unconstrained all-atom molecular dynamics (MD) does not serve well for structure refinement due to its limited conformational search. We have developed and tested the constrained MD method, based on the generalized Newton-Euler inverse mass operator (GNEIMO) algorithm for protein structure refinement. In this method, the high-frequency degrees of freedom are replaced with hard holonomic constraints and a protein is modeled as a collection of rigid body clusters connected by flexible torsional hinges. This allows larger integration time steps and enhances the conformational search space. In this work, we have demonstrated the use of torsional GNEIMO method without using any experimental data as constraints, for protein structure refinement starting from low-resolution decoy sets derived from homology methods. In the eight proteins with three decoys for each, we observed an improvement of ~2 Å in the rmsd in coordinates to the known experimental structures of these proteins. The GNEIMO trajectories also showed enrichment in the population density of native-like conformations. In addition, we demonstrated structural refinement using a "freeze and thaw" clustering scheme with the GNEIMO framework as a viable tool for enhancing localized conformational search. We have derived a robust protocol based on the GNEIMO replica exchange method for protein structure refinement that can be readily extended to other proteins and possibly applicable for high throughput protein structure refinement.

Outcome predictors of initial treatment with topical lubricant and parafunctional habit control in burning mouth syndrome (BMS).

Burning mouth syndrome (BMS) is a common chronic pain condition which mainly affects elderly women. The concomitant prescription of topical lubricants with oral parafunctional habit control has been reported as an effective initial approach for patients with BMS. In this study, we have investigated outcome predictors of this initial treatment in patients with BMS. One hundred forty patients with BMS (12 men and 128 women, mean age 59.9±10.7 years) were instructed to avoid oral parafunctions and to use topical lubricant for 2 weeks. The patients were sub-grouped according to psychological status, salivary flow rate, presence of psychiatric medications, symptom area and duration, symptom severity, presence of oral parafunctions, and accompanying oral complaints. The changes in symptoms were analyzed and compared between sub-groups. Subjects with T-scores≤50 for each psychological symptom dimension, a flow rate of stimulated whole saliva (SWS)>0.5 ml/min, no psychiatric medications, and a greater degree of initial symptoms (VAS≥5) displayed greater decreases in symptoms compared with their counterparts. In conclusion, psychological status, psychiatric medications, flow rate SWS, and initial symptom severity can be outcome predictors of the initial treatment approach for patients with BMS.

Folding of small proteins using constrained molecular dynamics.

The focus of this paper is to examine whether conformational search using constrained molecular dynamics (MD) method is more enhanced and enriched toward "native-like" structures compared to all-atom MD for the protein folding as a model problem. Constrained MD methods provide an alternate MD tool for protein structure prediction and structure refinement. It is computationally expensive to perform all-atom simulations of protein folding because the processes occur on a time scale of microseconds. Compared to the all-atom MD simulation, constrained MD methods have the advantage that stable dynamics can be achieved for larger time steps and the number of degrees of freedom is an order of magnitude smaller, leading to a decrease in computational cost. We have developed a generalized constrained MD method that allows the user to "freeze and thaw" torsional degrees of freedom as fit for the problem studied. We have used this method to perform all-torsion constrained MD in implicit solvent coupled with the replica exchange method to study folding of small proteins with various secondary structural motifs such as, α-helix (polyalanine, WALP16), ß-turn (1E0Q), and a mixed motif protein (Trp-cage). We demonstrate that constrained MD replica exchange method exhibits a wider conformational search than all-atom MD with increased enrichment of near-native structures. "Hierarchical" constrained MD simulations, where the partially formed helical regions in the initial stretch of the all-torsion folding simulation trajectory of Trp-cage were frozen, showed a better sampling of near-native structures than all-torsion constrained MD simulations. This is in agreement with the zipping-and-assembly folding model put forth by Dill and co-workers for folding proteins. The use of hierarchical "freeze and thaw" clustering schemes in constrained MD simulation can be used to sample conformations that contribute significantly to folding of proteins.

Dynamic ligand-induced-fit simulation via enhanced conformational samplings and ensemble dockings: a survivin example.

Survivin is an anticancer drug target due to its overexpression in tumor cells in a homodimer form. Abbott Laboratories has identified a small molecule binding site near the dimerization interface in a high-throughput-screening (HTS)-NMR experiment. A benchmarking of the binding mode of the compound Abbott8 aided in the search for the ligand-induced-fit receptor structure by exploring the conformational space of the survivin dimer. We performed ensemble dockings with Abbott8 against a large set of conformations sampled via replica exchange molecular dynamics (REMD). This enhanced sampling allowed the reproduction of the holo-NMR experimental binding mode. Surprisingly, the major structural change in the best-REMD snapshot corresponding to the small molecule induced-fit happens in the so-called "survivin mitosis/apoptosis switch loop", consistent with the X-ray crystal structure of survivin-monomer/borealin/INCENP chromosomal passenger complex (CPC), as the distance between Phe93 and Phe101 increased. To verify this hypothetical pathway for the induced-fit conformational change, we utilized morphed intermediate structures that combined the X-ray data and the best-REMD snapshot, and the potential of mean force (PMF) of the survivin dimer was constructed with umbrella sampling (US) followed by a multiple Bennett acceptance ratio estimator (MBAR). It revealed a 3-4 kcal/mol free energy barrier along the reaction coordinate, and the complex is stabilized by the gain of the binding energies of Abbott8. This free energy barrier might prohibit the reproduction of the experimental binding mode from the regular NTP-MD ensemble docking that we had tried. The combination of REMD generalized ensemble sampling with ensemble docking and free energy pathway analysis may provide a novel research protocol for the simulation of protein-ligand induced-fit recognition.

Are gas-phase models of interstellar chemistry tenable? The case of methanol.

We consider the case of methanol production in cold dark clouds, also known as quiescent cores, for which recent work shows that a purely gas-phase synthesis is unlikely to produce a sufficient amount to explain the observational fractional abundance of approximately 10(-9). Moreover, recent experiments appear to confirm a previous hypothesis that methanol can be formed on cold grain surfaces by the hydrogenation of CO via successive reactions with hydrogen atoms. In this paper we consider two ways of including the surface formation of methanol into chemical models of cold dark clouds. First, we use a gas-phase model and artificially include the surface formation of methanol in the same manner that the formation of molecular hydrogen is included. Secondly, we utilize a gas-grain code with a new mechanism for desorption following exothermic chemical reactions on grain surfaces. The latter method can reproduce the observed fractional abundance of gas-phase methanol and many other gas-phase species in the well-studied cold dark cloud TMC1-CP but the best fit to the observational data occurs at times significantly later than at ages estimated from gas-phase models.