Search

When dispersed and cultured in a multielectrode dish (MED), suprachiasmatic nucleus (SCN) neurons express fast oscillations of firing rate (FOFR; fast relative to the circadian cycle), with burst duration ∼10 min, and interburst interval varying from 20 to 60 min in different cells but remaining nevertheless rather regular in individual cells. In many cases, separate neurons in distant parts of the 1 mm recording area of a MED exhibited correlated FOFR. Neither the mechanism of FOFR nor the mechanism of their synchronization among neurons is known. Based on recent data implicating vasoactive intestinal polypeptide (VIP) as a key intercellular synchronizing agent, we built a model in which VIP acts as both a feedback regulator to generate FOFR in individual neurons, and a diffusible synchronizing agent to produce coherent electrical output of a neuronal network. In our model, VIP binding to its (VPAC2) receptors acts through Gs G-proteins to activate adenylyl cyclase (AC), increase intracellular cAMP, and open cyclic-nucleotide-gated (CNG) cation channels, thus depolarizing the cell and generating neuronal firing to release VIP. In parallel, slowly developing homologous desensitization and internalization of VPAC2 receptors terminates elevation of cAMP and thereby provides an interpulse silent interval. Through mathematical modeling, we show that this VIP/VPAC2/AC/cAMP/CNG-channel mechanism is sufficient for generating reliable FOFR in single neurons. When our model for FOFR is combined with a published model of synchronization of circadian rhythms based on VIP/VPAC2 and Per gene regulation synchronization of circadian rhythms is significantly accelerated. These results suggest that (a) auto/paracrine regulation by VIP/VPAC2 and intracellular AC/cAMP/CNG-channels are sufficient to provide robust FOFR and synchrony among neurons in a heterogeneous network, and (b) this system may also participate in synchronization of circadian rhythms. [ABSTRACT FROM AUTHOR]

The article presents research on the release of inflammatory chemokines from Weibel-Palade bodies (WPBs) during exocytosis, with a focus on differential cargo mobilisation after transient fusion with the plasma membrane. Topics include the selectivity of cargo release by WPBs due to the formation of a narrow fusion pore acting as a molecular size filter; intravascular inflammatory processes with reduced risk of thrombosis; and membrane mobility.

Collective cell migration is a fundamental process, occurring during embryogenesis and cancer metastasis. Neural crest cells exhibit such coordinated migration, where aberrant motion can lead to fatality or dysfunction of the embryo. Migration involves at least two complementary mechanisms: contact inhibition of locomotion (a repulsive interaction corresponding to a directional change of migration upon contact with a reciprocating cell), and co-attraction (a mutual chemoattraction mechanism). Here, we develop and employ a parameterized discrete element model of neural crest cells, to investigate how these mechanisms contribute to long-range directional migration during development. Motion is characterized using a coherence parameter and the time taken to reach, collectively, a target location. The simulated cell group is shown to switch from a diffusive to a persistent state as the response-rate to co-attraction is increased. Furthermore, the model predicts that when co-attraction is inhibited, neural crest cells can migrate into restrictive regions. Indeed, inhibition of co-attraction in vivo and in vitro leads to cell invasion into restrictive areas, confirming the prediction of the model. This suggests that the interplay between the complementary mechanisms may contribute to guidance of the neural crest. We conclude that directional migration is a system property and does not require action of external chemoattractants. [ABSTRACT FROM AUTHOR]

Background: The ever increasing discovery of non-coding RNAs leads to unprecedented demand for the accurate modeling of RNA folding, including the predictions of two-dimensional (base pair) and three-dimensional all-atom structures and folding stabilities. Accurate modeling of RNA structure and stability has far-reaching impact on our understanding of RNA functions in human health and our ability to design RNA-based therapeutic strategies. Results: The Vfold server offers a web interface to predict (a) RNA two-dimensional structure from the nucleotide sequence, (b) three-dimensional structure from the two-dimensional structure and the sequence, and (c) folding thermodynamics (heat capacity melting curve) from the sequence. To predict the two-dimensional structure (base pairs), the server generates an ensemble of structures, including loop structures with the different intra-loop mismatches, and evaluates the free energies using the experimental parameters for the base stacks and the loop entropy parameters given by a coarse-grained RNA folding model (the Vfold model) for the loops. To predict the three-dimensional structure, the server assembles the motif scaffolds using structure templates extracted from the known PDB structures and refines the structure using all-atom energy minimization. Conclusions: The Vfold-based web server provides a user friendly tool for the prediction of RNA structure and stability. The web server and the source codes are freely accessible for public use at “”. [ABSTRACT FROM AUTHOR]

A majority of viruses are composed of long single-stranded genomic RNA molecules encapsulated by protein shells with diameters of just a few tens of nanometers. We examine the extent to which these viral RNAs have evolved to be physically compact molecules to facilitate encapsulation. Measurements of equal-length viral, non-viral, coding and non-coding RNAs show viral RNAs to have among the smallest sizes in solution, i.e., the highest gel-electrophoretic mobilities and the smallest hydrodynamic radii. Using graph-theoretical analyses we demonstrate that their sizes correlate with the compactness of branching patterns in predicted secondary structure ensembles. The density of branching is determined by the number and relative positions of 3-helix junctions, and is highly sensitive to the presence of rare higher-order junctions with 4 or more helices. Compact branching arises from a preponderance of base pairing between nucleotides close to each other in the primary sequence. The density of branching represents a degree of freedom optimized by viral RNA genomes in response to the evolutionary pressure to be packaged reliably. Several families of viruses are analyzed to delineate the effects of capsid geometry, size and charge stabilization on the selective pressure for RNA compactness. Compact branching has important implications for RNA folding and viral assembly. [ABSTRACT FROM AUTHOR]

Late-stage or post-market identification of adverse drug reactions (ADRs) is a significant public health issue and a source of major economic liability for drug development. Thus, reliable in silico screening of drug candidates for possible ADRs would be advantageous. In this work, we introduce a computational approach that predicts ADRs by combining the results of molecular docking and leverages known ADR information from DrugBank and SIDER. We employed a recently parallelized version of AutoDock Vina (VinaLC) to dock 906 small molecule drugs to a virtual panel of 409 DrugBank protein targets. L1-regularized logistic regression models were trained on the resulting docking scores of a 560 compound subset from the initial 906 compounds to predict 85 side effects, grouped into 10 ADR phenotype groups. Only 21% (87 out of 409) of the drug-protein binding features involve known targets of the drug subset, providing a significant probe of off-target effects. As a control, associations of this drug subset with the 555 annotated targets of these compounds, as reported in DrugBank, were used as features to train a separate group of models. The Vina off-target models and the DrugBank on-target models yielded comparable median area-under-the-receiver-operating-characteristic-curves (AUCs) during 10-fold cross-validation (0.60–0.69 and 0.61–0.74, respectively). Evidence was found in the PubMed literature to support several putative ADR-protein associations identified by our analysis. Among them, several associations between neoplasm-related ADRs and known tumor suppressor and tumor invasiveness marker proteins were found. A dual role for interstitial collagenase in both neoplasms and aneurysm formation was also identified. These associations all involve off-target proteins and could not have been found using available drug/on-target interaction data. This study illustrates a path forward to comprehensive ADR virtual screening that can potentially scale with increasing number of CPUs to tens of thousands of protein targets and millions of potential drug candidates. [ABSTRACT FROM AUTHOR]

Sirtuins are key regulators of many cellular functions including cell growth, apoptosis, metabolism, and genetic control of age-related diseases. Sirtuins are themselves regulated by their cofactor nicotinamide adenine dinucleotide (NAD+) as well as their reaction product nicotinamide (NAM), the physiological concentrations of which vary during the process of aging. Nicotinamide inhibits sirtuins through the so-called base exchange pathway, wherein rebinding of the reaction product to the enzyme accelerates the reverse reaction. We investigated the mechanism of nicotinamide inhibition of human SIRT3, the major mitochondrial sirtuin deacetylase, in vitro and in silico using experimental kinetic analysis and Molecular Mechanics-Poisson Boltzmann/Generalized Born Surface Area (MM-PB(GB)SA) binding affinity calculations with molecular dynamics sampling. Through experimental kinetic studies, we demonstrate that NAM inhibition of SIRT3 involves apparent competition between the inhibitor and the enzyme cofactor NAD+, contrary to the traditional characterization of base exchange as noncompetitive inhibition. We report a model for base exchange inhibition that relates such kinetic properties to physicochemical properties, including the free energies of enzyme-ligand binding, and estimate the latter through the first reported computational binding affinity calculations for SIRT3:NAD+, SIRT3:NAM, and analogous complexes for Sir2. The computational results support our kinetic model, establishing foundations for quantitative modeling of NAD+/NAM regulation of mammalian sirtuins during aging and the computational design of sirtuin activators that operate through alleviation of base exchange inhibition. [ABSTRACT FROM AUTHOR]

A new coarse-grained model of the E. coli cytoplasm is developed by describing the proteins of the cytoplasm as flexible units consisting of one or more spheres that follow Brownian dynamics (BD), with hydrodynamic interactions (HI) accounted for by a mean-field approach. Extensive BD simulations were performed to calculate the diffusion coefficients of three different proteins in the cellular environment. The results are in close agreement with experimental or previously simulated values, where available. Control simulations without HI showed that use of HI is essential to obtain accurate diffusion coefficients. Anomalous diffusion inside the crowded cellular medium was investigated with Fractional Brownian motion analysis, and found to be present in this model. By running a series of control simulations in which various forces were removed systematically, it was found that repulsive interactions (volume exclusion) are the main cause for anomalous diffusion, with a secondary contribution from HI. [ABSTRACT FROM AUTHOR]

Recent work on the use of mRNA lipoplexes for gene delivery demonstrates the need for a mathematical model that simulates and predicts kinetics and transfection efficiency. The small copy numbers involved make it necessary to use stochastic models and include statistical analysis of the variation observed in the experimental data. The modeling requirements are further complicated by the multi-level nature of the problem, where mRNA molecules are contained in lipoplexes, which are in turn contained in endosomes, where each of these entities displays a behavior of its own. We have created a mathematical model that reproduces both the time courses and the statistical variance observed in recent experiments using single-cell tracking of GFP expression after transfection. By applying a few key simplifications and assumptions, we have limited the number of free parameters to five, which we optimize to match five experimental determinants by means of a simulated annealing algorithm. The models demonstrate the need for modeling of nested species in order to reproduce the shape of the dose-response and expression-level curves. [ABSTRACT FROM AUTHOR]

The success of recombinant virus-like particles (VLPs) for human papillomavirus and hepatitis B demonstrates the potential of VLPs as safe and efficacious vaccines. With new modular designs emerging, the effects of antigen module insertion on the self-assembly and structural integrity of VLPs should be clarified so as to better enabling improved design. Previous work has revealed insights into the molecular energetics of a VLP subunit, capsomere, comparing energetics within various solution conditions known to drive or inhibit self-assembly. In the present study, molecular dynamics (MD) simulations coupled with the molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) method were performed to examine the molecular interactions and energetics in a modular capsomere of a murine polyomavirus (MPV) VLP designed to protect against influenza. Insertion of an influenza antigenic module is found to lower the binding energy within the capsomere, and a more active state is observed in Assembly Buffer as compared with that in Stabilization Buffer, which has been experimentally validated through measurements using differential scanning calorimetry. Further in-depth analysis based on free-energy decomposition indicates that destabilized binding can be attributed to electrostatic interaction induced by the chosen antigen module. These results provide molecular insights into the conformational stability of capsomeres and their abilities to be exploited for antigen presentation, and are expected to be beneficial for the biomolecular engineering of VLP vaccines. [ABSTRACT FROM AUTHOR]

The plant-specific insert is an approximately 100-residue domain found exclusively within the C-terminal lobe of some plant aspartic proteases. Structurally, this domain is a member of the saposin-like protein family, and is involved in plant pathogen defense as well as vacuolar targeting of the parent protease molecule. Similar to other members of the saposin-like protein family, most notably saposins A and C, the recently resolved crystal structure of potato (Solanum tuberosum) plant-specific insert has been shown to exist in a substrate-bound open conformation in which the plant-specific insert oligomerizes to form homodimers. In addition to the open structure, a closed conformation also exists having the classic saposin fold of the saposin-like protein family as observed in the crystal structure of barley (Hordeum vulgare L.) plant-specific insert. In the present study, the mechanisms of tertiary and quaternary conformation changes of potato plant-specific insert were investigated in silico as a function of pH. Umbrella sampling and determination of the free energy change of dissociation of the plant-specific insert homodimer revealed that increasing the pH of the system to near physiological levels reduced the free energy barrier to dissociation. Furthermore, principal component analysis was used to characterize conformational changes at both acidic and neutral pH. The results indicated that the plant-specific insert may adopt a tertiary structure similar to the characteristic saposin fold and suggest a potential new structural motif among saposin-like proteins. To our knowledge, this acidified PSI structure presents the first example of an alternative saposin-fold motif for any member of the large and diverse SAPLIP family. [ABSTRACT FROM AUTHOR]

Nanoscale particles have become promising materials in many fields, such as cancer therapeutics, diagnosis, imaging, drug delivery, catalysis, as well as biosensors. In order to stimulate and facilitate these applications, there is an urgent need for the understanding of the interaction mode between the nano-particles and proteins. In this study, we investigate the orientation and adsorption between several enzymes (cytochrome c, RNase A, lysozyme) and 4 nm/11 nm silica nanoparticles (SNPs) by using molecular dynamics (MD) simulation. Our results show that three enzymes are adsorbed onto the surfaces of both 4 nm and 11 nm SNPs during our MD simulations and the small SNPs induce greater structural stabilization. The active site of cytochrome c is far away from the surface of 4 nm SNPs, while it is adsorbed onto the surface of 11 nm SNPs. We also explore the influences of different groups (-OH, -COOH, -NH2 and CH3) coated onto silica nanoparticles, which show significantly different impacts. Our molecular dynamics results indicate the selective interaction between silicon nanoparticles and enzymes, which is consistent with experimental results. Our study provides useful guides for designing/modifying nanomaterials to interact with proteins for their bio-applications. [ABSTRACT FROM AUTHOR]

Life has effectively colonized most of our planet and extremophilic organisms require specialized enzymes to survive under harsh conditions. Cold-loving organisms (psychrophiles) express heat-labile enzymes that possess a high specific activity and catalytic efficiency at low temperatures. A remarkable universal characteristic of cold-active enzymes is that they show a reduction both in activation enthalpy and entropy, compared to mesophilic orthologs, which makes their reaction rates less sensitive to falling temperature. Despite significant efforts since the early 1970s, the important question of the origin of this effect still largely remains unanswered. Here we use cold- and warm-active trypsins as model systems to investigate the temperature dependence of the reaction rates with extensive molecular dynamics free energy simulations. The calculations quantitatively reproduce the catalytic rates of the two enzymes and further yield high-precision Arrhenius plots, which show the characteristic trends in activation enthalpy and entropy. Detailed structural analysis indicates that the relationship between these parameters and the 3D structure is reflected by significantly different internal protein energy changes during the reaction. The origin of this effect is not localized to the active site, but is found in the outer regions of the protein, where the cold-active enzyme has a higher degree of softness. Several structural mechanisms for softening the protein surface are identified, together with key mutations responsible for this effect. Our simulations further show that single point-mutations can significantly affect the thermodynamic activation parameters, indicating how these can be optimized by evolution. [ABSTRACT FROM AUTHOR]

Flexibility in biomolecular recognition is essential and critical for many cellular activities. Flexible recognition often leads to moderate affinity but high specificity, in contradiction with the conventional wisdom that high affinity and high specificity are coupled. Furthermore, quantitative understanding of the role of flexibility in biomolecular recognition is still challenging. Here, we meet the challenge by quantifying the intrinsic biomolecular recognition energy landscapes with and without flexibility through the underlying density of states. We quantified the thermodynamic intrinsic specificity by the topography of the intrinsic binding energy landscape and the kinetic specificity by association rate. We found that the thermodynamic and kinetic specificity are strongly correlated. Furthermore, we found that flexibility decreases binding affinity on one hand, but increases binding specificity on the other hand, and the decreasing or increasing proportion of affinity and specificity are strongly correlated with the degree of flexibility. This shows more (less) flexibility leads to weaker (stronger) coupling between affinity and specificity. Our work provides a theoretical foundation and quantitative explanation of the previous qualitative studies on the relationship among flexibility, affinity and specificity. In addition, we found that the folding energy landscapes are more funneled with binding, indicating that binding helps folding during the recognition. Finally, we demonstrated that the whole binding-folding energy landscapes can be integrated by the rigid binding and isolated folding energy landscapes under weak flexibility. Our results provide a novel way to quantify the affinity and specificity in flexible biomolecular recognition. [ABSTRACT FROM AUTHOR]

In vitro cultures of endothelial cells are a widely used model system of the collective behavior of endothelial cells during vasculogenesis and angiogenesis. When seeded in an extracellular matrix, endothelial cells can form blood vessel-like structures, including vascular networks and sprouts. Endothelial morphogenesis depends on a large number of chemical and mechanical factors, including the compliancy of the extracellular matrix, the available growth factors, the adhesion of cells to the extracellular matrix, cell-cell signaling, etc. Although various computational models have been proposed to explain the role of each of these biochemical and biomechanical effects, the understanding of the mechanisms underlying in vitro angiogenesis is still incomplete. Most explanations focus on predicting the whole vascular network or sprout from the underlying cell behavior, and do not check if the same model also correctly captures the intermediate scale: the pairwise cell-cell interactions or single cell responses to ECM mechanics. Here we show, using a hybrid cellular Potts and finite element computational model, that a single set of biologically plausible rules describing (a) the contractile forces that endothelial cells exert on the ECM, (b) the resulting strains in the extracellular matrix, and (c) the cellular response to the strains, suffices for reproducing the behavior of individual endothelial cells and the interactions of endothelial cell pairs in compliant matrices. With the same set of rules, the model also reproduces network formation from scattered cells, and sprouting from endothelial spheroids. Combining the present mechanical model with aspects of previously proposed mechanical and chemical models may lead to a more complete understanding of in vitro angiogenesis. [ABSTRACT FROM AUTHOR]

Protein-ligand recognition plays key roles in many biological processes. One of the most fascinating questions about protein-ligand recognition is to understand its underlying mechanism, which often results from a combination of induced fit and conformational selection. In this study, we have developed a three-pronged approach of Markov State Models, Molecular Dynamics simulations, and flux analysis to determine the contribution of each model. Using this approach, we have quantified the recognition mechanism of the choline binding protein (ChoX) to be ∼90% conformational selection dominant under experimental conditions. This is achieved by recovering all the necessary parameters for the flux analysis in combination with available experimental data. Our results also suggest that ChoX has several metastable conformational states, of which an apo-closed state is dominant, consistent with previous experimental findings. Our methodology holds great potential to be widely applied to understand recognition mechanisms underlining many fundamental biological processes. [ABSTRACT FROM AUTHOR]

Single-molecule tweezers measurements of double-stranded nucleic acids (dsDNA and dsRNA) provide unprecedented opportunities to dissect how these fundamental molecules respond to forces and torques analogous to those applied by topoisomerases, viral capsids, and other biological partners. However, tweezers data are still most commonly interpreted post facto in the framework of simple analytical models. Testing falsifiable predictions of state-of-the-art nucleic acid models would be more illuminating but has not been performed. Here we describe a blind challenge in which numerical predictions of nucleic acid mechanical properties were compared to experimental data obtained recently for dsRNA under applied force and torque. The predictions were enabled by the HelixMC package, first presented in this paper. HelixMC advances crystallography-derived base-pair level models (BPLMs) to simulate kilobase-length dsDNAs and dsRNAs under external forces and torques, including their global linking numbers. These calculations recovered the experimental bending persistence length of dsRNA within the error of the simulations and accurately predicted that dsRNA's “spring-like” conformation would give a two-fold decrease of stretch modulus relative to dsDNA. Further blind predictions of helix torsional properties, however, exposed inaccuracies in current BPLM theory, including three-fold discrepancies in torsional persistence length at the high force limit and the incorrect sign of dsRNA link-extension (twist-stretch) coupling. Beyond these experiments, HelixMC predicted that ‘nucleosome-excluding’ poly(A)/poly(T) is at least two-fold stiffer than random-sequence dsDNA in bending, stretching, and torsional behaviors; Z-DNA to be at least three-fold stiffer than random-sequence dsDNA, with a near-zero link-extension coupling; and non-negligible effects from base pair step correlations. We propose that experimentally testing these predictions should be powerful next steps for understanding the flexibility of dsDNA and dsRNA in sequence contexts and under mechanical stresses relevant to their biology. [ABSTRACT FROM AUTHOR]

Identification of important nodes in complex networks has attracted an increasing attention over the last decade. Various measures have been proposed to characterize the importance of nodes in complex networks, such as the degree, betweenness and PageRank. Different measures consider different aspects of complex networks. Although there are numerous results reported on undirected complex networks, few results have been reported on directed biological networks. Based on network motifs and principal component analysis (PCA), this paper aims at introducing a new measure to characterize node importance in directed biological networks. Investigations on five real-world biological networks indicate that the proposed method can robustly identify actually important nodes in different networks, such as finding command interneurons, global regulators and non-hub but evolutionary conserved actually important nodes in biological networks. Receiver Operating Characteristic (ROC) curves for the five networks indicate remarkable prediction accuracy of the proposed measure. The proposed index provides an alternative complex network metric. Potential implications of the related investigations include identifying network control and regulation targets, biological networks modeling and analysis, as well as networked medicine. [ABSTRACT FROM AUTHOR]

2-Mercaptobenzimidazole (MBI) is widely utilized as a corrosion inhibitor, copper-plating brightener and rubber accelerator. The residue of MBI in the environment is potentially harmful. In the present work, the toxic interaction of MBI with the important antioxidant enzyme copper-zinc superoxide dismutase (Cu/ZnSOD) was investigated using spectroscopic and molecular docking methods. MBI can interact with Cu/ZnSOD to form an MBI-Cu/ZnSOD complex. The binding constant, number of binding sites and thermodynamic parameters were measured, which indicated that MBI could spontaneously bind with Cu/ZnSOD with one binding site through hydrogen bonds and van der Waals forces. MBI bound into the Cu/ZnSOD interface of two subdomains, which caused some microenvironmental and secondary structure changes of Cu/ZnSOD and further resulted in the inhibition of Cu/ZnSOD activity. This work provides direct evidence at a molecular level to show that exposure to MBI could induce changes in the structure and function of the enzyme Cu/ZnSOD. The estimated methods in this work may be applied to probe molecular interactions of biomacromolecules and other pollutants and drugs. [ABSTRACT FROM AUTHOR]

Antibody-based therapeutics provides novel and efficacious treatments for a number of diseases. Traditional experimental approaches for designing therapeutic antibodies rely on raising antibodies against a target antigen in an immunized animal or directed evolution of antibodies with low affinity for the desired antigen. However, these methods remain time consuming, cannot target a specific epitope and do not lead to broad design principles informing other studies. Computational design methods can overcome some of these limitations by using biophysics models to rationally select antibody parts that maximize affinity for a target antigen epitope. This has been addressed to some extend by OptCDR for the design of complementary determining regions. Here, we extend this earlier contribution by addressing the de novo design of a model of the entire antibody variable region against a given antigen epitope while safeguarding for immunogenicity (Optimal Method for Antibody Variable region Engineering, OptMAVEn). OptMAVEn simulates in silico the in vivo steps of antibody generation and evolution, and is capable of capturing the critical structural features responsible for affinity maturation of antibodies. In addition, a humanization procedure was developed and incorporated into OptMAVEn to minimize the potential immunogenicity of the designed antibody models. As case studies, OptMAVEn was applied to design models of neutralizing antibodies targeting influenza hemagglutinin and HIV gp120. For both HA and gp120, novel computational antibody models with numerous interactions with their target epitopes were generated. The observed rates of mutations and types of amino acid changes during in silico affinity maturation are consistent with what has been observed during in vivo affinity maturation. The results demonstrate that OptMAVEn can efficiently generate diverse computational antibody models with both optimized binding affinity to antigens and reduced immunogenicity. [ABSTRACT FROM AUTHOR]

A focused theme in systems biology is to uncover design principles of biological networks, that is, how specific network structures yield specific systems properties. For this purpose, we have previously developed a reverse engineering procedure to identify network topologies with high likelihood in generating desired systems properties. Our method searches the continuous parameter space of an assembly of network topologies, without enumerating individual network topologies separately as traditionally done in other reverse engineering procedures. Here we tested this CPSS (continuous parameter space search) method on a previously studied problem: the resettable bistability of an Rb-E2F gene network in regulating the quiescence-to-proliferation transition of mammalian cells. From a simplified Rb-E2F gene network, we identified network topologies responsible for generating resettable bistability. The CPSS-identified topologies are consistent with those reported in the previous study based on individual topology search (ITS), demonstrating the effectiveness of the CPSS approach. Since the CPSS and ITS searches are based on different mathematical formulations and different algorithms, the consistency of the results also helps cross-validate both approaches. A unique advantage of the CPSS approach lies in its applicability to biological networks with large numbers of nodes. To aid the application of the CPSS approach to the study of other biological systems, we have developed a computer package that is available in Information S1. [ABSTRACT FROM AUTHOR]

Ultraweak biophoton emission originates from the generation of reactive oxygen species (ROS) that are produced in mitochondria as by-products of cellular respiration. In healthy cells, the concentration of ROS is minimized by a system of biological antioxidants. However, heat shock changes the equilibrium between oxidative stress and antioxidant activity, that is, a rapid rise in temperature induces biophoton emission from ROS. Although the rate and intensity of biophoton emission was observed to increase in response to elevated temperatures, pretreatment at lower high temperatures inhibited photon emission at higher temperatures. Biophoton measurements are useful for observing and evaluating heat shock. [ABSTRACT FROM AUTHOR]

The DAP12-NKG2C activating immunoreceptor complex is one of the multisubunit transmembrane protein complexes in which ligand-binding receptor chains assemble with dimeric signal-transducing modules through non-covalent associations in their transmembrane (TM) domains. In this work, both coarse grained and atomistic molecular dynamic simulation methods were applied to investigate the self-assembly dynamics of the transmembrane domains of the DAP12-NKG2C activating immunoreceptor complex. Through simulating the dynamics of DAP12-NKG2C TM heterotrimer and point mutations, we demonstrated that a five-polar-residue motif including: 2 Asps and 2 Thrs in DAP12 dimer, as well as 1 Lys in NKG2C TM plays an important role in the assembly structure of the DAP12-NKG2C TM heterotrimer. Furthermore, we provided clear evidences to exclude the possibility that another NKG2C could stably associate with the DAP12-NKG2C heterotrimer. Based on the simulation results, we proposed a revised model for the self-assembly of DAP12-NKG2C activating immunoreceptor complex, along with a plausible explanation for the association of only one NKG2C with a DAP12 dimer. [ABSTRACT FROM AUTHOR]

Olfactory adaptation is a fundamental process for the functioning of the olfactory system, but the underlying mechanisms regulating its occurrence in intact olfactory sensory neurons (OSNs) are not fully understood. In this work, we have combined stochastic computational modeling and a systematic pharmacological study of different signaling pathways to investigate their impact during short-term adaptation (STA). We used odorant stimulation and electroolfactogram (EOG) recordings of the olfactory epithelium treated with pharmacological blockers to study the molecular mechanisms regulating the occurrence of adaptation in OSNs. EOG responses to paired-pulses of odorants showed that inhibition of phosphodiesterases (PDEs) and phosphatases enhanced the levels of STA in the olfactory epithelium, and this effect was mimicked by blocking vesicle exocytosis and reduced by blocking cyclic adenosine monophosphate (cAMP)-dependent protein kinase (PKA) and vesicle endocytosis. These results suggest that G-coupled receptors (GPCRs) cycling is involved with the occurrence of STA. To gain insights on the dynamical aspects of this process, we developed a stochastic computational model. The model consists of the olfactory transduction currents mediated by the cyclic nucleotide gated (CNG) channels and calcium ion (Ca2+)-activated chloride (CAC) channels, and the dynamics of their respective ligands, cAMP and Ca2+, and it simulates the EOG results obtained under different experimental conditions through changes in the amplitude and duration of cAMP and Ca2+ response, two second messengers implicated with STA occurrence. The model reproduced the experimental data for each pharmacological treatment and provided a mechanistic explanation for the action of GPCR cycling in the levels of second messengers modulating the levels of STA. All together, these experimental and theoretical results indicate the existence of a mechanism of regulation of STA by signaling pathways that control GPCR cycling and tune the levels of second messengers in OSNs, and not only by CNG channel desensitization as previously thought. [ABSTRACT FROM AUTHOR]

Homology models of mammalian voltage-gated sodium (NaV) channels based on the crystal structures of the bacterial counterparts are needed to interpret the functional data on sodium channels and understand how they operate. Such models would also be invaluable in structure-based design of therapeutics for diseases involving sodium channels such as chronic pain and heart diseases. Here we construct a homology model for the pore domain of the NaV1.4 channel and use the functional data for the binding of µ-conotoxin GIIIA to NaV1.4 to validate the model. The initial poses for the NaV1.4–GIIIA complex are obtained using the HADDOCK protein docking program, which are then refined in molecular dynamics simulations. The binding mode for the final complex is shown to be in broad agreement with the available mutagenesis data. The standard binding free energy, determined from the potential of mean force calculations, is also in good agreement with the experimental value. Because the pore domains of NaV1 channels are highly homologous, the model constructed for NaV1.4 will provide an excellent template for other NaV1 channels. [ABSTRACT FROM AUTHOR]

Alzheimer's disease (AD) is a complex neurodegenerative disorder involving multiple cellular and molecular processes. The discovery of drug molecules capable of targeting multiple factors involved in AD pathogenesis would greatly facilitate in improving therapeutic strategies. The repositioning of existing non-toxic drugs could dramatically reduce the time and costs involved in developmental and clinical trial stages. In this study, preliminary screening of 140 FDA approved nervous system drugs by docking suggested the viability of the tricyclic group of antidepressants against three major AD targets, viz. Acetylcholinesterase (AChE), β-secretase (BACE-1), and amyloid β (Aβ) aggregation, with one member, protriptyline, showing highest inhibitory activity. Detailed biophysical assays, together with isothermal calorimetry, fluorescence quenching experiments, kinetic studies and atomic force microscopy established the strong inhibitory activity of protriptyline against all three major targets. The molecular basis of inhibition was supported with comprehensive molecular dynamics simulations. Further, the drug inhibited glycation induced amyloid aggregation, another important causal factor in AD progression. This study has led to the discovery of protriptyline as a potent multi target directed ligand and established its viability as a promising candidate for AD treatment. [ABSTRACT FROM AUTHOR]

We present a coarse-grained two dimensional mechanical model for the microtubule-tau bundles in neuronal axons in which we remove taus, as can happen in various neurodegenerative conditions such as Alzheimers disease, tauopathies, and chronic traumatic encephalopathy. Our simplified model includes (i) taus modeled as entropic springs between microtubules, (ii) removal of taus from the bundles due to phosphorylation, and (iii) a possible depletion force between microtubules due to these dissociated phosphorylated taus. We equilibrate upon tau removal using steepest descent relaxation. In the absence of the depletion force, the transverse rigidity to radial compression of the bundles falls to zero at about 60% tau occupancy, in agreement with standard percolation theory results. However, with the attractive depletion force, spring removal leads to a first order collapse of the bundles over a wide range of tau occupancies for physiologically realizable conditions. While our simplest calculations assume a constant concentration of microtubule intercalants to mediate the depletion force, including a dependence that is linear in the detached taus yields the same collapse. Applying percolation theory to removal of taus at microtubule tips, which are likely to be the protective sites against dynamic instability, we argue that the microtubule instability can only obtain at low tau occupancy, from 0.06–0.30 depending upon the tau coordination at the microtubule tips. Hence, the collapse we discover is likely to be more robust over a wide range of tau occupancies than the dynamic instability. We suggest in vitro tests of our predicted collapse. [ABSTRACT FROM AUTHOR]

The hierarchical structure of spider dragline silk is composed of two major constituents, the amorphous phase and crystalline units, and its mechanical response has been attributed to these prime constituents. Silk mechanics, however, might also be influenced by the resistance against sliding of these two phases relative to each other under load. We here used atomistic molecular dynamics (MD) simulations to obtain friction forces for the relative sliding of the amorphous phase and crystalline units of Araneus diadematus spider silk. We computed the coefficient of viscosity of this interface to be in the order of 102 Ns/m2 by extrapolating our simulation data to the viscous limit. Interestingly, this value is two orders of magnitude smaller than the coefficient of viscosity within the amorphous phase. This suggests that sliding along a planar and homogeneous surface of straight polyalanine chains is much less hindered than within entangled disordered chains. Finally, in a simple finite element model, which is based on parameters determined from MD simulations including the newly deduced coefficient of viscosity, we assessed the frictional behavior between these two components for the experimental range of relative pulling velocities. We found that a perfectly relative horizontal motion has no significant resistance against sliding, however, slightly inclined loading causes measurable resistance. Our analysis paves the way towards a finite element model of silk fibers in which crystalline units can slide, move and rearrange themselves in the fiber during loading. [ABSTRACT FROM AUTHOR]

Many signaling events require the binding of cytoplasmic proteins to cell membranes by recognition of specific charged lipids, such as phosphoinositol-phosphates. As a model for a protein-membrane binding site, we consider one charged phosphoinositol phosphate (PtdIns(3)P) embedded in a phosphatidylcholine bilayer. As the protein-membrane binding is driven by electrostatic interactions, continuum solvent models require an accurate representation of the electrostatic potential of the phosphoinositol phosphate-containing membrane. We computed and analyzed the electrostatic potentials of snapshots taken at regular intervals from molecular dynamics simulations of the bilayer. We observe considerable variation in the electrostatic potential of the bilayer both along a single simulation and between simulations performed with the GAFF or CHARMM c36 force fields. However, we find that the choice of GAFF or CHARMM c36 parameters has little effect on the electrostatic potential of a given configuration of the bilayer with a PtdIns(3)P embedded in it. From our results, we propose a remedian averaging method for calculating the electrostatic potential of a membrane system that is suitable for simulations of protein-membrane binding with a continuum solvent model. [ABSTRACT FROM AUTHOR]

Although the crystal structures of the HIV-1 gp120 core bound and pre-bound by CD4 are known, the details of dynamics involved in conformational equilibrium and transition in relation to gp120 function have remained elusive. The homology models of gp120 comprising the N- and C-termini and loops V3 and V4 in the CD4-bound and CD4-unbound states were built and subjected to molecular dynamics (MD) simulations to investigate the differences in dynamic properties and molecular motions between them. The results indicate that the CD4-bound gp120 adopted a more compact and stable conformation than the unbound form during simulations. For both the unbound and bound gp120, the large concerted motions derived from essential dynamics (ED) analyses can influence the size/shape of the ligand-binding channel/cavity of gp120 and, therefore, were related to its functional properties. The differences in motion direction between certain structural components of these two forms of gp120 were related to the conformational interconversion between them. The free energy calculations based on the metadynamics simulations reveal a more rugged and complex free energy landscape (FEL) for the unbound than for the bound gp120, implying that gp120 has a richer conformational diversity in the unbound form. The estimated free energy difference of ∼−6.0 kJ/mol between the global minimum free energy states of the unbound and bound gp120 indicates that gp120 can transform spontaneously from the unbound to bound states, revealing that the bound state represents a high-probability “ground state” for gp120 and explaining why the unbound state resists crystallization. Our results provide insight into the dynamics-and-function relationship of gp120, and facilitate understandings of the thermodynamics, kinetics and conformational control mechanism of HIV-1 gp120. [ABSTRACT FROM AUTHOR]

The motility of spermatozoa of both Lytechinus pictus and Strongylocentrotus purpuratus sea urchin species is modulated by the egg-derived decapeptide speract via an oscillatory [Ca2+]-dependent signaling pathway. Comprehension of this pathway is hence directly related to the understanding of regulated sperm swimming. Niflumic acid (NFA), a nonsteroidal anti-inflammatory drug alters several ion channels. Though unspecific, NFA profoundly affects how sea urchin sperm respond to speract, increasing the [Ca2+]i oscillation period, amplitude, peak and average level values of the responses in immobilized and swimming cells. A previous logical network model we developed for the [Ca2+] dynamics of speract signaling cascade in sea urchin sperm allows integrated dissection of individual and multiple actions of NFA. Among the channels affected by NFA are: hyperpolarization-activated and cyclic nucleotide gated Na+ channels (HCN), [Ca2+]-dependent Cl− channels (CaCC) and [Ca2+]-dependent K+ channels (CaKC), all present in the sea urchin genome. Here, using our model we investigated the effect of blocking in silico HCN and CaCC channels suggested by experiments. Regarding CaKC channels, arguments can be provided for either their blockage or activation by NFA. Our study yielded two scenarios compliant with experimental observations: i) under CaKC inhibition, this [Ca2+]-dependent K+ channel should be different from the Slo1 channel and ii) under activation of the CaKC channel, another [Ca2+] channel not considered previously in the network is required, such as the pH-dependent CatSper channel. Additionally, our findings predict cause-effect relations resulting from a selective inhibition of those channels. Knowledge of these relations may be of consequence for a variety of electrophysiological studies and have an impact on drug related investigations. Our study contributes to a better grasp of the network dynamics and suggests further experimental work. [ABSTRACT FROM AUTHOR]

In this study we estimate the subcutaneous tissue counter pressure during drug infusion from a series of injections of insulin in type 2 diabetic patients using a non-invasive method. We construct a model for the pressure evolution in subcutaneous tissue based on mass continuity and the flow laws of a porous medium. For equivalent injection forces we measure the change in the infusion rate between injections in air at atmospheric pressure and in tissue. From a best fit with our model, we then determine the flow permeability as well as the bulk modulus of the tissue, estimated to be of the order 10−11–10−10 m2 and 105 Pa, respectively. The permeability is in good agreement with reported values for adipose porcine tissue. We suggest our model as a general way to estimate the pressure build-up in tissue during subcutaneous injection. [ABSTRACT FROM AUTHOR]

Voltage Dependent Anion-selective Channels (VDACs) are pore-forming proteins located in the outer mitochondrial membrane. They are responsible for the access of ions and energetic metabolites into the inner membrane transport systems. Three VDAC isoforms exist in mammalian, but their specific role is unknown. In this work we have performed extensive (overall ∼5 µs) Molecular Dynamics (MD) simulations of the human VDAC isoforms to detect structural and conformational variations among them, possibly related to specific functional roles of these proteins. Secondary structure analysis of the N-terminal domain shows a high similarity among the three human isoforms of VDAC but with a different plasticity. In particular, the N-terminal domain of the hVDAC1 is characterized by a higher plasticity, with a ∼20% occurrence for the ‘unstructured’ conformation throughout the folded segment, while hVDAC2, containing a peculiar extension of 11 amino acids at the N-terminal end, presents an additional 310-helical folded portion comprising residues 10′ to 3, adhering to the barrel wall. The N-terminal sequences of hVDAC isoforms are predicted to have a low flexibility, with possible consequences in the dynamics of the human VDACs. Clear differences were found between hVDAC1 and hVDAC3 against hVDAC2: a significantly modified dynamics with possible important consequence on the voltage-gating mechanism. Charge distribution inside and at the mouth of the pore is responsible for a different preferential localization of ions with opposite charge and provide a valuable rationale for hVDAC1 and hVDAC3 having a Cl−/K+ selectivity ratio of 1.8, whereas hVDAC2 of 1.4. Our conclusion is that hVDAC isoforms, despite sharing a similar scaffold, have modified working features and a biological work is now requested to give evidence to the described dissimilarities. [ABSTRACT FROM AUTHOR]

Bacterial cholesterol oxidase is commonly used as an experimental tool to reduce cellular cholesterol content. That the treatment also generates the poorly degradable metabolite 4-cholesten-3-one (cholestenone) has received less attention. Here, we investigated the membrane partitioning of cholestenone using simulations and cell biological experiments and assessed the functional effects of cholestenone in human cells. Atomistic simulations predicted that cholestenone reduces membrane order, undergoes faster flip-flop and desorbs more readily from membranes than cholesterol. In primary human fibroblasts, cholestenone was released from membranes to physiological extracellular acceptors more avidly than cholesterol, but without acceptors it remained in cells over a day. To address the functional effects of cholestenone, we studied fibroblast migration during wound healing. When cells were either cholesterol oxidase treated or part of cellular cholesterol was exchanged for cholestenone with cyclodextrin, cell migration during 22 h was markedly inhibited. Instead, when a similar fraction of cholesterol was removed using cyclodextrin, cells replenished their cholesterol content in 3 h and migrated similarly to control cells. Thus, cholesterol oxidation produces long-term functional effects in cells and these are in part due to the generated membrane active cholestenone. [ABSTRACT FROM AUTHOR]

The movement of juvenile loggerhead turtles (n = 42) out-fitted with satellite tags and released in oceanic waters off New Caledonia was examined and compared with ocean circulation data. Merging of the daily turtle movement data with drifter buoy movements, OSCAR (Ocean Surface Current Analyses - Real time) circulation data, and three different vertical strata (0–5 m, 0–40 m, 0–100 m) of HYCOM (HYbrid Coordinate Ocean Model) circulation data indicated the turtles were swimming against the prevailing current in a statistically significant pattern. This was not an artifact of prevailing directions of current and swimming, nor was it an artifact of frictional slippage. Generalized additive modeling was used to decompose the pattern of swimming into spatial and temporal components. The findings are indicative of a positive rheotaxis whereby an organism is able to detect the current flow and orient itself to swim into the current flow direction or otherwise slow down its movement. Potential mechanisms for the means and adaptive significance of rheotaxis in oceanic juvenile loggerhead turtles are discussed. [ABSTRACT FROM AUTHOR]

Cutaneous pain sensations are mediated largely by C-nociceptors consisting of both mechano-sensitive (CM) and mechano-insensitive (CMi) fibres that can be distinguished from one another according to their characteristic axonal properties. In healthy skin and relative to CMi fibres, CM fibres show a higher initial conduction velocity, less activity-dependent conduction velocity slowing, and less prominent post-spike supernormality. However, after sensitization with nerve growth factor, the electrical signature of CMi fibres changes towards a profile similar to that of CM fibres. Here we take a combined experimental and modelling approach to examine the molecular basis of such alterations to the excitation thresholds. Changes in electrical activation thresholds and activity-dependent slowing were examined in vivo using single-fibre recordings of CM and CMi fibres in domestic pigs following NGF application. Using computational modelling, we investigated which axonal mechanisms contribute most to the electrophysiological differences between the fibre classes. Simulations of axonal conduction suggest that the differences between CMi and CM fibres are strongly influenced by the densities of the delayed rectifier potassium channel (Kdr), the voltage-gated sodium channels NaV1.7 and NaV1.8, and the Na+/K+-ATPase. Specifically, the CM fibre profile required less Kdr and NaV1.8 in combination with more NaV1.7 and Na+/K+-ATPase. The difference between CM and CMi fibres is thus likely to reflect a relative rather than an absolute difference in protein expression. In support of this, it was possible to replicate the experimental reduction of the ADS pattern of CMi nociceptors towards a CM-like pattern following intradermal injection of nerve growth factor by decreasing the contribution of Kdr (by 50%), increasing the Na+/K+-ATPase (by 10%), and reducing the branch length from 2 cm to 1 cm. The findings highlight key molecules that potentially contribute to the NGF-induced switch in nociceptors phenotype, in particular NaV1.7 which has already been identified clinically as a principal contributor to chronic pain states such as inherited erythromelalgia. [ABSTRACT FROM AUTHOR]

Protein disulfide isomerase (PDI) composed of four thioredoxin-like domains a, b, b', and a', is a key enzyme catalyzing oxidative protein folding in the endoplasmic reticulum. Large scale molecular dynamics simulations starting from the crystal structures of human PDI (hPDI) in the oxidized and reduced states were performed. The results indicate that hPDI adopts more compact conformations in solution than in the crystal structures, which are stabilized primarily by inter-domain interactions, including the salt bridges between domains a and b' observed for the first time. A prominent feature of the compact conformations is that the two catalytic domains a and a' can locate close enough for intra-molecular electron transfer, which was confirmed by the characterization of an intermediate with a disulfide between the two domains. Mutations, which disrupt the inter-domain interactions, lead to decreased reductase activity of hPDI. Our molecular dynamics simulations and biochemical experiments reveal the intrinsic conformational dynamics of hPDI and its biological impact. [ABSTRACT FROM AUTHOR]

The development of noninvasive approaches to facilitate the regeneration of post-traumatic nerve injury is important for clinical rehabilitation. In this study, we investigated the effective dose of noninvasive 808-nm low-level laser therapy (LLLT) on sciatic nerve crush rat injury model. Thirty-six male Sprague Dawley rats were divided into 6 experimental groups: a normal group with or without 808-nm LLLT at 8 J/cm2 and a sciatic nerve crush injury group with or without 808-nm LLLT at 3, 8 or 15 J/cm2. Rats were given consecutive transcutaneous LLLT at the crush site and sacrificed 20 days after the crush injury. Functional assessments of nerve regeneration were analyzed using the sciatic functional index (SFI) and hindlimb range of motion (ROM). Nerve regeneration was investigated by measuring the myelin sheath thickness of the sciatic nerve using transmission electron microscopy (TEM) and by analyzing the expression of growth-associated protein 43 (GAP43) in sciatic nerve using western blot and immunofluorescence staining. We found that sciatic-injured rats that were irradiated with LLLT at both 3 and 8 J/cm2 had significantly improved SFI but that a significant improvement of ROM was only found in rats with LLLT at 8 J/cm2. Furthermore, the myelin sheath thickness and GAP43 expression levels were significantly enhanced in sciatic nerve-crushed rats receiving 808-nm LLLT at 3 and 8 J/cm2. Taken together, these results suggest that 808-nm LLLT at a low energy density (3 J/cm2 and 8 J/cm2) is capable of enhancing sciatic nerve regeneration following a crush injury. [ABSTRACT FROM AUTHOR]

Vascular smooth muscle α2C-adrenoceptors (α2C-ARs) mediate vasoconstriction of small blood vessels, especially arterioles. Studies of endogenous receptors in human arteriolar smooth muscle cells (referred to as microVSM) and transiently transfected receptors in heterologous HEK293 cells show that the α2C-ARs are perinuclear receptors that translocate to the cell surface under cellular stress and elicit a biological response. Recent studies in microVSM unraveled a crucial role of Rap1A-Rho-ROCK-F-actin pathways in receptor translocation, and identified protein-protein interaction of α2C-ARs with the actin binding protein filamin-2 as an essential step in the process. To better understand the molecular nature and specificity of this interaction, in this study, we constructed comparative models of human α2C-AR and human filamin-2 proteins. Finally, we performed in silico protein-protein docking to provide a structural platform for the investigation of human α2C-AR and filamin-2 interactions. We found that electrostatic interactions seem to play a key role in this complex formation which manifests in interactions between the C-terminal arginines of α2C-ARs (particularly R454 and R456) and negatively charged residues from filamin-2 region between residues 1979 and 2206. Phylogenetic and sequence analysis showed that these interactions have evolved in warm-blooded animals. [ABSTRACT FROM AUTHOR]

Electroporation-based therapies are powerful biotechnological tools for enhancing the delivery of exogeneous agents or killing tissue with pulsed electric fields (PEFs). Electrochemotherapy (ECT) and gene therapy based on gene electrotransfer (EGT) both use reversible electroporation to deliver chemotherapeutics or plasmid DNA into cells, respectively. In both ECT and EGT, the goal is to permeabilize the cell membrane while maintaining high cell viability in order to facilitate drug or gene transport into the cell cytoplasm and induce a therapeutic response. Irreversible electroporation (IRE) results in cell kill due to exposure to PEFs without drugs and is under clinical evaluation for treating otherwise unresectable tumors. These PEF therapies rely mainly on the electric field distributions and do not require changes in tissue temperature for their effectiveness. However, in immediate vicinity of the electrodes the treatment may results in cell kill due to thermal damage because of the inhomogeneous electric field distribution and high current density during the electroporation-based therapies. Therefore, the main objective of this numerical study is to evaluate the influence of pulse number and electrical conductivity in the predicted cell kill zone due to irreversible electroporation and thermal damage. Specifically, we simulated a typical IRE protocol that employs ninety 100-µs PEFs. Our results confirm that it is possible to achieve predominant cell kill due to electroporation if the PEF parameters are chosen carefully. However, if either the pulse number and/or the tissue conductivity are too high, there is also potential to achieve cell kill due to thermal damage in the immediate vicinity of the electrodes. Therefore, it is critical for physicians to be mindful of placement of electrodes with respect to critical tissue structures and treatment parameters in order to maintain the non-thermal benefits of electroporation and prevent unnecessary damage to surrounding healthy tissue, critical vascular structures, and/or adjacent organs. [ABSTRACT FROM AUTHOR]

Proteins in thermophilic organisms remain stable and function optimally at high temperatures. Owing to their important applicability in many industrial processes, such thermostable proteins have been studied extensively, and several structural factors attributed to their enhanced stability. How these factors render the emergent property of thermostability to proteins, even in situations where no significant changes occur in their three-dimensional structures in comparison to their mesophilic counter-parts, has remained an intriguing question. In this study we treat Lipase A from Bacillus subtilis and its six thermostable mutants in a unified manner and address the problem with a combined complex network-based analysis and molecular dynamic studies to find commonality in their properties. The Protein Contact Networks (PCN) of the wild-type and six mutant Lipase A structures developed at a mesoscopic scale were analyzed at global network and local node (residue) level using network parameters and community structure analysis. The comparative PCN analysis of all proteins pointed towards important role of specific residues in the enhanced thermostability. Network analysis results were corroborated with finer-scale molecular dynamics simulations at both room and high temperatures. Our results show that this combined approach at two scales can uncover small but important changes in the local conformations that add up to stabilize the protein structure in thermostable mutants, even when overall conformation differences among them are negligible. Our analysis not only supports the experimentally determined stabilizing factors, but also unveils the important role of contacts, distributed throughout the protein, that lead to thermostability. We propose that this combined mesoscopic-network and fine-grained molecular dynamics approach is a convenient and useful scheme not only to study allosteric changes leading to protein stability in the face of negligible over-all conformational changes due to mutations, but also in other molecular networks where change in function does not accompany significant change in the network structure. [ABSTRACT FROM AUTHOR]

Cell–matrix and cell–cell mechanosensing are important in many cellular processes, particularly for epithelial cells. A crucial question, which remains unexplored, is how the mechanical microenvironment is altered as a result of changes to multicellular tissue structure during cancer progression. In this study, we investigated the influence of the multicellular tissue architecture on mechanical properties of the epithelial component of the mammary acinus. Using creep compression tests on multicellular breast epithelial structures, we found that pre-malignant acini with no lumen (MCF10AT) were significantly stiffer than normal hollow acini (MCF10A) by 60%. This difference depended on structural changes in the pre-malignant acini, as neither single cells nor normal multicellular acini tested before lumen formation exhibited these differences. To understand these differences, we simulated the deformation of the acini with different multicellular architectures and calculated their mechanical properties; our results suggest that lumen filling alone can explain the experimentally observed stiffness increase. We also simulated a single contracting cell in different multicellular architectures and found that lumen filling led to a 20% increase in the “perceived stiffness” of a single contracting cell independent of any changes to matrix mechanics. Our results suggest that lumen filling in carcinogenesis alters the mechanical microenvironment in multicellular epithelial structures, a phenotype that may cause downstream disruptions to mechanosensing. [ABSTRACT FROM AUTHOR]

Lumen formation is important for morphogenesis; however, an unanswered question is whether it involves the collective migration of epithelial cells. Here, using a collagen gel overlay culture method, we show that Madin-Darby canine kidney cells migrated collectively and formed a luminal structure in a collagen gel. Immediately after the collagen gel overlay, an epithelial sheet folded from the periphery, migrated inwardly, and formed a luminal structure. The inhibition of integrin-β1 or Rac1 activity decreased the migration rate of the peripheral cells after the sheets folded. Moreover, lumen formation was perturbed by disruption of apical-basolateral polarity induced by transforming growth factor-β1. These results indicate that cell migration and cell polarity play an important role in folding. To further explore epithelial sheet folding, we developed a computer-simulated mechanical model based on the rigidity of the extracellular matrix. It indicated a soft substrate is required for the folding movement. [ABSTRACT FROM AUTHOR]

It is widely recognized that propagation of electrophysiological signals between the soma and dendrites of neurons differs depending on direction, i.e. it is asymmetric. How this asymmetry influences the activation of voltage-gated dendritic channels, and consequent neuronal behavior, remains unclear. Based on the analysis of asymmetry in several types of motoneurons, we extended our previous methodology for reducing a fully reconstructed motoneuron model to a two-compartment representation that preserved asymmetric signal propagation. The reduced models accurately replicated the dendritic excitability and the dynamics of the anatomical model involving a persistent inward current (PIC) dispersed over the dendrites. The relationship between asymmetric signal propagation and dendritic excitability was investigated using the reduced models while varying the asymmetry in signal propagation between the soma and the dendrite with PIC density constant. We found that increases in signal attenuation from soma to dendrites increased the activation threshold of a PIC (hypo-excitability), whereas increases in signal attenuation from dendrites to soma decreased the activation threshold of a PIC (hyper-excitability). These effects were so strong that reversing the asymmetry in the soma-to-dendrite vs. dendrite-to-soma attenuation, reversed the correlation between PIC threshold and distance of this current source from the soma. We propose the tight relation of the asymmetric signal propagation to the input resistance in the dendrites as a mechanism underlying the influence of the asymmetric signal propagation on the dendritic excitability. All these results emphasize the importance of maintaining the physiological asymmetry in dendritic signaling not only for normal function of the cells but also for biophysically realistic simulations of dendritic excitability. [ABSTRACT FROM AUTHOR]

Rapid and selective ion transport is essential for the generation and regulation of electrical signaling pathways in living organisms. Here, we use molecular dynamics (MD) simulations with an applied membrane potential to investigate the ion flux of bacterial sodium channel NaVMs. 5.9 µs simulations with 500 mM NaCl suggest different mechanisms for inward and outward flux. The predicted inward conductance rate of ∼27±3 pS, agrees with experiment. The estimated outward conductance rate is 15±3 pS, which is considerably lower. Comparing inward and outward flux, the mean ion dwell time in the selectivity filter (SF) is prolonged from 13.5±0.6 ns to 20.1±1.1 ns. Analysis of the Na+ distribution revealed distinct patterns for influx and efflux events. In 32.0±5.9% of the simulation time, the E53 side chains adopted a flipped conformation during outward conduction, whereas this conformational change was rarely observed (2.7±0.5%) during influx. Further, simulations with dihedral restraints revealed that influx is less affected by the E53 conformational flexibility. In contrast, during outward conduction, our simulations indicate that the flipped E53 conformation provides direct coordination for Na+. The free energy profile (potential of mean force calculations) indicates that this conformational change lowers the putative barriers between sites SCEN and SHFS during outward conduction. We hypothesize that during an action potential, the increased Na+ outward transition propensities at depolarizing potentials might increase the probability of E53 conformational changes in the SF. Subsequently, this might be a first step towards initiating slow inactivation. [ABSTRACT FROM AUTHOR]

Diffusion barriers are effective means for constraining protein lateral exchange in cellular membranes. In Saccharomyces cerevisiae, they have been shown to sustain parental identity through asymmetric segregation of ageing factors during closed mitosis. Even though barriers have been extensively studied in the plasma membrane, their identity and organization within the nucleus remains poorly understood. Based on different lines of experimental evidence, we present a model of the composition and structural organization of a nuclear diffusion barrier during anaphase. By means of spatial stochastic simulations, we propose how specialised lipid domains, protein rings, and morphological changes of the nucleus may coordinate to restrict protein exchange between mother and daughter nuclear lobes. We explore distinct, plausible configurations of these diffusion barriers and offer testable predictions regarding their protein exclusion properties and the diffusion regimes they generate. Our model predicts that, while a specialised lipid domain and an immobile protein ring at the bud neck can compartmentalize the nucleus during early anaphase; a specialised lipid domain spanning the elongated bridge between lobes would be entirely sufficient during late anaphase. Our work shows how complex nuclear diffusion barriers in closed mitosis may arise from simple nanoscale biophysical interactions. [ABSTRACT FROM AUTHOR]

Advances reported over the last few years and the increasing availability of protein crystal structure data have greatly improved structure-based druggability approaches. However, in practice, nearly all druggability estimation methods are applied to protein crystal structures as rigid proteins, with protein flexibility often not directly addressed. The inclusion of protein flexibility is important in correctly identifying the druggability of pockets that would be missed by methods based solely on the rigid crystal structure. These include cryptic pockets and flexible pockets often found at protein-protein interaction interfaces. Here, we apply an approach that uses protein modeling in concert with druggability estimation to account for light protein backbone movement and protein side-chain flexibility in protein binding sites. We assess the advantages and limitations of this approach on widely-used protein druggability sets. Applying the approach to all mammalian protein crystal structures in the PDB results in identification of 69 proteins with potential druggable cryptic pockets. [ABSTRACT FROM AUTHOR]

A module for fast determination of reduction potentials, E°, of redox-active proteins has been implemented in the CHARMM INterface and Graphics (CHARMMing) web portal (www.charmming.org). The free energy of reduction, which is proportional to E°, is composed of an intrinsic contribution due to the redox site and an environmental contribution due to the protein and solvent. Here, the intrinsic contribution is selected from a library of pre-calculated density functional theory values for each type of redox site and redox couple, while the environmental contribution is calculated from a crystal structure of the protein using Poisson-Boltzmann continuum electrostatics. An accompanying lesson demonstrates a calculation of E°. In this lesson, an ionizable residue in a [4Fe-4S]-protein that causes a pH-dependent E° is identified, and the E° of a mutant that would test the identification is predicted. This demonstration is valuable to both computational chemistry students and researchers interested in predicting sequence determinants of E° for mutagenesis. [ABSTRACT FROM AUTHOR]

A lesson utilizing a coarse-grained (CG) G-like model has been implemented into the CHARMM INterface and Graphics (CHARMMing) web portal (www.charmming.org) to the Chemistry at HARvard Macromolecular Mechanics (CHARMM) molecular simulation package. While widely used to model various biophysical processes, such as protein folding and aggregation, CG models can also serve as an educational tool because they can provide qualitative descriptions of complex biophysical phenomena for a relatively cheap computational cost. As a proof of concept, this lesson demonstrates the construction of a CG model of a small globular protein, its simulation via Langevin dynamics, and the analysis of the resulting data. This lesson makes connections between modern molecular simulation techniques and topics commonly presented in an advanced undergraduate lecture on physical chemistry. It culminates in a straightforward analysis of a short dynamics trajectory of a small fast folding globular protein; we briefly describe the thermodynamic properties that can be calculated from this analysis. The assumptions inherent in the model and the data analysis are laid out in a clear, concise manner, and the techniques used are consistent with those employed by specialists in the field of CG modeling. One of the major tasks in building the G-like model is determining the relative strength of the nonbonded interactions between coarse-grained sites. New functionality has been added to CHARMMing to facilitate this process. The implementation of these features into CHARMMing helps automate many of the tedious aspects of constructing a CG G model. The CG model builder and its accompanying lesson should be a valuable tool to chemistry students, teachers, and modelers in the field. [ABSTRACT FROM AUTHOR]

Correlated inter-domain motions in proteins can mediate fundamental biochemical processes such as signal transduction and allostery. Here we characterize at structural level the inter-domain coupling in a multidomain enzyme, Adenylate Kinase (AK), using computational methods that exploit the shape information encoded in residual dipolar couplings (RDCs) measured under steric alignment by nuclear magnetic resonance (NMR). We find experimental evidence for a multi-state equilibrium distribution along the opening/closing pathway of Adenylate Kinase, previously proposed from computational work, in which inter-domain interactions disfavour states where only the AMP binding domain is closed. In summary, we provide a robust experimental technique for study of allosteric regulation in AK and other enzymes. [ABSTRACT FROM AUTHOR]