Search

Background: The discharge of textile industrial wastewater, a complex effluent containing different chemicals, to the environment causes serious bioenvironmental pollution problems.Therefore, this study aimed to evaluate the use of advanced oxidation methods and find the most efficient method to remove dye from textile effluents.Materials and Methods: This experimental study was conducted on Kashan textile industrial effluent on a laboratory scale and in a closed system. Taking into consideration the fixed amount of O3 dose and H2O2 concentration for the four advanced oxidation methods(UV/O3, UV/H2O2, H2O2/O3, O3/UV/H2O2), in three reaction time (10, 20, 30 minutes) and different pH (4, 6, 8, 10), the effect of the methods on dye removal was evaluated and the results were compared using the complex repetitive measurements.Results: Results showed that, in terms of method used and the contact time, the O3/UV/H2O2 method in 30 minutes had the most efficient result (89.2%) and the UV/H2O2method in 10 minutes time the least one (40.7%). Moreover, in terms of the method used and pH, the O3/UV/H2O2 method in pH=6 had the most efficient result (84.7) and the UV/H2O2method in pH=4 the least one (59.4). The results showed that time had a positive effect on dye removal efficiency (P<0.001).Conclusion: The O3/UV/H2O2 method is the most efficient method in removing dye fromtextile wastewater due to the synchronous presence of multi-pollutants and its synergic effect through producing the highly reactive hydroxyl radicals (OH0). [ABSTRACT FROM AUTHOR]

Limited therapy options due to antibiotic resistance underscore the need for optimization of current diagnostics. In some bacterial species, antimicrobial resistance can be unambiguously predicted based on their genome sequence. In this study, we sequenced the genomes and transcriptomes of 414 drug-resistant clinical Pseudomonas aeruginosa isolates. By training machine learning classifiers on information about the presence or absence of genes, their sequence variation, and expression profiles, we generated predictive models and identified biomarkers of resistance to four commonly administered antimicrobial drugs. Using these data types alone or in combination resulted in high (0.8-0.9) or very high (> 0.9) sensitivity and predictive values. For all drugs except for ciprofloxacin, gene expression information improved diagnostic performance. Our results pave the way for the development of a molecular resistance profiling tool that reliably predicts antimicrobial susceptibility based on genomic and transcriptomic markers. The implementation of a molecular susceptibility test system in routine microbiology diagnostics holds promise to provide earlier and more detailed information on antibiotic resistance profiles of bacterial pathogens and thus could change how physicians treat bacterial infections., (© 2020 The Authors. Published under the terms of the CC BY 4.0 license.)

Protein thermostability engineering is a powerful tool to improve resistance of proteins against high temperatures and thereafter broaden their applications. For efficient protein thermostability engineering, different thermostability-classified data sources including sequences and 3D structures are needed for different protein families. However, no data source is available providing such data easily. It is the first release of ProtDataTherm database for analysis and engineering of protein thermostability which contains more than 14 million protein sequences categorized based on their thermal stability and protein family. This database contains data needed for better understanding protein thermostability and stability engineering. Providing categorized protein sequences and structures as psychrophilic, mesophilic and thermophilic makes this database useful for the development of new tools in protein stability prediction. This database is available at http://profiles.bs.ipm.ir/softwares/protdatatherm. As a proof of concept, the thermostability that improves mutations were suggested for one sample protein belonging to one of protein families with more than 20 mesophilic and thermophilic sequences and with known experimentally measured ΔT of mutations available within ProTherm database.

Export of messenger ribonucleic acids (mRNAs) into the cytoplasm is a fundamental step in gene regulation processes, which is meticulously quality controlled by highly efficient mechanisms in eukaryotic cells. Yet, it remains unclear how the aberrant mRNAs are recognized and retained inside the nucleus. Using a new modelling approach for complex systems, namely the agent-based modelling (ABM) approach, we develop a minimal model of the mRNA quality control (QC) mechanism. Our results demonstrate that regulation of the affinity of RNA-binding proteins (RBPs) to export receptors along with the weak interaction between the nuclear basket protein (Mlp1 or Tpr) and RBPs are the minimum requirements to distinguish and retain aberrant mRNAs. Our results show that the affinity between Tpr and RBPs is optimized to maximize the retention of aberrant mRNAs. In addition, we demonstrate how the length of mRNA affects the QC process. Since longer mRNAs spend more time in the nuclear basket to form a compact conformation and initiate their export, nuclear basket proteins could more easily capture and retain them inside the nucleus.

The genetic information of eukaryotic cells is enclosed within a double-layered nuclear envelope, which comprises an inner and outer nuclear membrane. Several transmembrane proteins locate to the nuclear envelope; however, only two integral protein complexes span the nuclear envelope and connect the inside of the nucleus to the cytoplasm. The nuclear pore complex (NPC) acts as a gateway for molecular exchange between the interior of the nucleus and the cytoplasm, whereas so-called LINC complexes physically link the nucleoskeleton and the cytoskeleton. In this Commentary, we will discuss recent studies that have established direct functional associations between these two complexes. The assembly of NPCs and their even distribution throughout the nuclear envelope is dependent on components of the LINC complex. Additionally, LINC complex formation is dependent on the successful localization of inner nuclear membrane components of LINC complexes and their transport through the NPC. Furthermore, the architecture of the nuclear envelope depends on both protein complexes. Finally, we will present recent evidence showing that LINC complexes can affect nucleo-cytoplasmic transport through the NPC, further highlighting the importance of understanding the associations of these essential complexes at the nuclear envelope., (© 2016. Published by The Company of Biologists Ltd.)

Conformational behavior of intrinsically disordered proteins, such as Phe-Gly repeat domains, alters drastically when they are confined in, and tethered to, nan channels. This has challenged our understanding of how they serve to selectively facilitate translocation of nuclear transport receptor (NTR)-bearing macromolecules. Heterogeneous FG-repeats, tethered to the NPC interior, nonuniformly fill the channel in a diameter-dependent manner and adopt a rapid Brownian motion, thereby forming a porous and highly dynamic polymeric meshwork that percolates in radial and axial directions and features two distinguishable zones: a dense hydrophobic rod-like zone located in the center, and a peripheral low-density shell-like zone. The FG-meshwork is locally disrupted upon interacting with NTR-bearing macromolecules, but immediately reconstructs itself between 0.44 μs and 7.0 μs, depending on cargo size and shape. This confers a perpetually-sealed state to the NPC, and is solely due to rapid Brownian motion of FG-repeats, not FG-repeat hydrophobic bonds. Elongated-shaped macromolecules, both in the presence and absence of NTRs, penetrate more readily into the FG-meshwork compared to their globular counterparts of identical volume and surface chemistry, highlighting the importance of the shape effects in nucleocytoplasmic transport. These results can help our understanding of geometrical effects in, and the design of, intelligent and responsive biopolymer-based materials in nanofiltration and artificial nanopores.

α-Actinin is an essential actin cross-linker involved in cytoskeletal organization and dynamics. The molecular conformation of α-actinin's actin-binding domain (ABD) regulates its association with actin and thus mutations in this domain can lead to severe pathogenic conditions. A point mutation at lysine 255 in human α-actinin-4 to glutamate increases the binding affinity resulting in stiffer cytoskeletal structures. The role of different ABD conformations and the effect of K255E mutation on ABD conformations remain elusive. To evaluate the impact of K255E mutation on ABD binding to actin we use all-atom molecular dynamics and free energy calculation methods and study the molecular mechanism of actin association in both wild-type α-actinin and in the K225E mutant. Our models illustrate that the strength of actin association is indeed sensitive to the ABD conformation, predict the effect of K255E mutation--based on simulations with the K237E mutant chicken α-actinin--and evaluate the mechanism of α-actinin binding to actin. Furthermore, our simulations showed that the calmodulin domain binding to the linker region was important for regulating the distance between actin and ABD. Our results provide valuable insights into the molecular details of this critical cellular phenomenon and further contribute to an understanding of cytoskeletal dynamics in health and disease., (Copyright © 2016 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Nucleocytoplasmic transport has been the subject of a large body of research in the past few decades. Recently, the focus of investigations in this field has shifted from studies of the overall function of the nuclear pore complex (NPC) to the examination of the role of different domains of phenylalanine-glycine nucleoporin (FG Nup) sequences on the NPC function. In our recent bioinformatics study, we showed that FG Nups have some evolutionarily conserved sequence-based features that might govern their physical behavior inside the NPC. We proposed the 'like charge regions' (LCRs), sequences of charged residues with only one type of charge, as one of the features that play a significant role in the formation of FG network inside the central channel. In this study, we further explore the role of LCRs in the distribution of FG Nups, using a recently developed coarse-grained molecular dynamics model. Our results demonstrate how LCRs affect the formation of two transport pathways. While some FG Nups locate their FG network at the center of the NPC forming a homogeneous meshwork of FG repeats, other FG Nups cover the space adjacent to the NPC wall. LCRs in the former group, i.e. FG Nups that form an FG domain at the center, tend to regulate the size of the highly dense, doughnut-shaped FG meshwork and leave a small low FG density area at the center of the pore for passive diffusion. On the other hand, LCRs in the latter group of FG Nups enable them to maximize their interactions and cover a larger space inside the NPC to increase its capability to transport numerous cargos at the same time. Finally, a new viewpoint is proposed that reconciles different models for the nuclear pore selective barrier function.

We introduce a new representation and feature extraction method for biological sequences. Named bio-vectors (BioVec) to refer to biological sequences in general with protein-vectors (ProtVec) for proteins (amino-acid sequences) and gene-vectors (GeneVec) for gene sequences, this representation can be widely used in applications of deep learning in proteomics and genomics. In the present paper, we focus on protein-vectors that can be utilized in a wide array of bioinformatics investigations such as family classification, protein visualization, structure prediction, disordered protein identification, and protein-protein interaction prediction. In this method, we adopt artificial neural network approaches and represent a protein sequence with a single dense n-dimensional vector. To evaluate this method, we apply it in classification of 324,018 protein sequences obtained from Swiss-Prot belonging to 7,027 protein families, where an average family classification accuracy of 93%±0.06% is obtained, outperforming existing family classification methods. In addition, we use ProtVec representation to predict disordered proteins from structured proteins. Two databases of disordered sequences are used: the DisProt database as well as a database featuring the disordered regions of nucleoporins rich with phenylalanine-glycine repeats (FG-Nups). Using support vector machine classifiers, FG-Nup sequences are distinguished from structured protein sequences found in Protein Data Bank (PDB) with a 99.8% accuracy, and unstructured DisProt sequences are differentiated from structured DisProt sequences with 100.0% accuracy. These results indicate that by only providing sequence data for various proteins into this model, accurate information about protein structure can be determined. Importantly, this model needs to be trained only once and can then be applied to extract a comprehensive set of information regarding proteins of interest. Moreover, this representation can be considered as pre-training for various applications of deep learning in bioinformatics. The related data is available at Life Language Processing Website: http://llp.berkeley.edu and Harvard Dataverse: http://dx.doi.org/10.7910/DVN/JMFHTN.

The nuclear pore complex (NPC) is the portal for bidirectional transportation of cargos between the nucleus and the cytoplasm. While most of the structural elements of the NPC, i.e. nucleoporins (Nups), are well characterized, the exact transport mechanism is still under much debate. Many of the functional Nups are rich in phenylalanine-glycine (FG) repeats and are believed to play the key role in nucleocytoplasmic transport. We present a bioinformatics study conducted on more than a thousand FG Nups across 252 species. Our results reveal the regulatory role of polar residues and specific sequences of charged residues, named 'like charge regions' (LCRs), in the formation of the FG network at the center of the NPC. Positively charged LCRs prepare the environment for negatively charged cargo complexes and regulate the size of the FG network. The low number density of charged residues in these regions prevents FG domains from forming a relaxed coil structure. Our results highlight the significant role of polar interactions in FG network formation at the center of the NPC and demonstrate that the specific localization of LCRs, FG motifs, charged, and polar residues regulate the formation of the FG network at the center of the NPC.

von Willebrand factor (VWF) is a naturally collapsed protein that participates in primary haemostasis and coagulation events. The clotting process is triggered by the adsorption and conformational changes of the plasma VWFs localized to the collagen fibres found near the site of injury. We develop coarse-grained models to simulate the adsorption dynamics of VWF flowing near the adhesive collagen fibres at different shear rates and investigate the effect of factors such as interaction and cooperativity of VWFs on the success of adsorption events. The adsorption probability of a flowing VWF confined to the receptor field is enhanced when it encounters an adhered VWF in proximity to the collagen receptors. This enhancement is observed within a wide range of shear rates and is mostly controlled by the attractive van der Waals interactions rather than the hydrodynamic interactions among VWF monomers. The cooperativity between the VWFs acts as an effective mechanism for enhancing VWF adsorption to the collagen fibres. Additionally, this implies that the adsorption of such molecules is nonlinearly dependent on the density of flowing VWFs. These findings are important for studies of primary haemostasis as well as general adsorption dynamics processes in polymer physics., (© 2015 The Author(s).)

Numerous biological functions of a cell, including polarization, differentiation, division, and migration, rely on its ability to endure mechanical forces generated by the cytoskeleton on the nucleus. Coupling of the cytoskeleton and nucleoskeleton is ultimately mediated by LINC complexes that are formed via a strong interaction between SUN- and KASH-domain-containing proteins in the nuclear envelope. These complexes are mechanosensitive and essential for the transmission of forces between the cytoskeleton and nucleoskeleton, and the progression of cellular mechanotransduction. Herein, using molecular dynamics, we examine the effect of tension on the human SUN2-KASH2 complex and show that it is remarkably stable under physiologically relevant tensile forces and large strains. However, a covalent disulfide bond between two highly conserved cysteine residues of SUN2 and KASH2 is crucial for the stability of this interaction and the transmission of forces through the complex., (Copyright © 2015 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Axonal microtubule (MT) bundles crosslinked by microtubule-associated protein (MAP) tau are responsible for vital biological functions such as maintaining mechanical integrity and shape of the axon as well as facilitating axonal transport. Breaking and twisting of MTs have been previously observed in damaged undulated axons. Such breaking and twisting of MTs is suggested to cause axonal swellings that lead to axonal degeneration, which is known as "diffuse axonal injury". In particular, overstretching and torsion of axons can potentially damage the axonal cytoskeleton. Following our previous studies on mechanical response of axonal MT bundles under uniaxial tension and compression, this work seeks to characterize the mechanical behavior of MT bundles under pure torsion as well as a combination of torsional and tensile loads using a coarse-grained computational model. In the case of pure torsion, a competition between MAP tau tensile and MT bending energies is observed. After three turns, a transition occurs in the mechanical behavior of the bundle that is characterized by its diameter shrinkage. Furthermore, crosslink spacing is shown to considerably influence the mechanical response, with larger MAP tau spacing resulting in a higher rate of turns. Therefore, MAP tau crosslinking of MT filaments protects the bundle from excessive deformation. Simultaneous application of torsion and tension on MT bundles is shown to accelerate bundle failure, compared to pure tension experiments. MAP tau proteins fail in clusters of 10-100 elements located at the discontinuities or the ends of MT filaments. This failure occurs in a stepwise fashion, implying gradual accumulation of elastic tensile energy in crosslinks followed by rupture. Failure of large groups of interconnecting MAP tau proteins leads to detachment of MT filaments from the bundle near discontinuities. This study highlights the importance of torsional loading in axonal damage after traumatic brain injury., (Copyright © 2015 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Molecular traffic between the nucleus and the cytoplasm is regulated by the nuclear pore complex (NPC), which acts as a highly selective channel perforating the nuclear envelope in eukaryotic cells. The human immunodeficiency virus (HIV) exploits the nucleocytoplasmic pathway to export its RNA transcripts across the NPC to the cytoplasm. Despite extensive study on the HIV life cycle and the many drugs developed to target this cycle, no current drugs have been successful in targeting the critical process of viral nuclear export, even though HIV's reliance on a single host protein, CRM1, to export its unspliced and partially spliced RNA transcripts makes it a tempting target. Due to recent findings implicating a DEAD-box helicase, DDX3, in HIV replication and a member of the export complex, it has become an appealing target for anti-HIV drug inhibition. In the present research, we have applied a hybrid computational protocol to analyze protein-protein interactions in the HIV mRNA export cycle. This method is based on molecular docking followed by molecular dynamics simulation and accompanied by approximate free energy calculation (MM/GBSA), computational alanine scanning, clustering, and evolutionary analysis. We highlight here some of the most likely binding modes and interfacial residues between DDX3 and CRM1 both in the absence and presence of RanGTP. This work shows that although DDX3 can bind to free CRM1, addition of RanGTP leads to more concentrated distribution of binding modes and stronger binding between CRM1 and RanGTP.

mRNA export from the nucleus is an essential step in the expression of every protein- coding gene in eukaryotes, but many aspects of this process remain poorly understood. The density of export receptors that must bind an mRNA to ensure export, as well as how receptor distribution affects transport dynamics, is not known. It is also unclear whether the rate-limiting step for transport occurs at the nuclear basket, in the central channel, or on the cytoplasmic face of the nuclear pore complex. Using previously published biophysical and biochemical parameters of mRNA export, we implemented a three-dimensional, coarse-grained, agent-based model of mRNA export in the nanosecond regime to gain insight into these issues. On running the model, we observed that mRNA export is sensitive to the number and distribution of transport receptors coating the mRNA and that there is a rate-limiting step in the nuclear basket that is potentially associated with the mRNA reconfiguring itself to thread into the central channel. Of note, our results also suggest that using a single location-monitoring mRNA label may be insufficient to correctly capture the time regime of mRNA threading through the pore and subsequent transport. This has implications for future experimental design to study mRNA transport dynamics., (© 2014 Azimi, Bulat, et al. This article is distributed by The American Society for Cell Biology under license from the author(s). Two months after publication it is available to the public under an Attribution–Noncommercial–Share Alike 3.0 Unported Creative Commons License (http://creativecommons.org/licenses/by-nc-sa/3.0).)

Metallic nanopillar/nanowires are emerging as promising platforms for biological applications, as they allow for the direct characterization and regulation of cell function. Herein we study the response of cells to a versatile nanopillar platform. Nanopillar arrays of various shape, size, and spacing and different nanopillar-substrate interfacial strengths were fabricated and interfaced with fibroblasts and several unique cell-nanopillar interactions were observed using high resolution scanning electron microscopy. Nanopillar penetration, engulfment, tilting, lift off and membrane thinning, were observed by manipulating nanopillar material, size, shape and spacing. These unique cell responses to various nanostructures can be employed for a wide range of applications including the design of highly sensitive nano-electrodes for single-cell probing., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

The function-blocking, non-RGD-containing, and primate-specific mouse monoclonal antibody 17E6 binds the αV subfamily of integrins. 17E6 is currently in phase II clinical trials for treating cancer. To elucidate the structural basis of recognition and the molecular mechanism of inhibition, we crystallized αVβ3 ectodomain in complex with the Fab fragment of 17E6. Protein crystals grew in presence of the activating cation Mn(2+). The integrin in the complex and in solution assumed the genuflected conformation. 17E6 Fab bound exclusively to the Propeller domain of the αV subunit. At the core of αV-Fab interface were interactions involving Propeller residues Lys-203 and Gln-145, with the latter accounting for primate specificity. The Propeller residue Asp-150, which normally coordinates Arg of the ligand Arg-Gly-Asp motif, formed contacts with Arg-54 of the Fab that were expected to reduce soluble FN10 binding to cellular αVβ3 complexed with 17E6. This was confirmed in direct binding studies, suggesting that 17E6 is an allosteric inhibitor of αV integrins., (© 2014 by The American Society for Biochemistry and Molecular Biology, Inc.)

A broad range of human diseases are associated with bacterial infections, often initiated by specific adhesion of a bacterium to the target environment. Despite the significant role of bacterial adhesion in human infectious diseases, details and mechanisms of bacterial adhesion have remained elusive. Herein, we study the physical interactions between Staphylococcus aureus, a type of micro-organism relevant to infections associated with medical implants, and nanocrystalline (nc) nickel nanostructures with various columnar features, including solid core, hollow, x-shaped and c-shaped pillars. Scanning electron microscopy results show the tendency of these bacterial cells to attach to the nickel nanostructures. Moreover, unique single bacterium attachment characteristics were observed on nickel nanostructures with dimensions comparable to the size of a single bacterium., (Copyright © 2014 Elsevier Ltd. All rights reserved.)

While much has been devoted to the study of transport mechanisms through the nuclear pore complex (NPC), the specifics of interactions and binding between export transport receptors and the NPC periphery have remained elusive. Recent work has demonstrated a binding interaction between the exportin CRM1 and the unstructured carboxylic tail of Tpr, on the nuclear basket. Strong evidence suggests that this interaction is vital to the functions of CRM1. Using molecular dynamics simulations and a newly refined method for determining binding regions, we have identified nine candidate binding sites on CRM1 for C-Tpr. These include two adjacent to RanGTP--from which one is blocked in the absence of RanGTP--and three next to the binding region of the cargo Snurportin. We report two additional interaction sites between C-Tpr and Snurportin, suggesting a possible role for Tpr import into the nucleus. Using bioinformatics tools we have conducted conservation analysis and functional residue prediction investigations to identify which parts of the obtained binding sites are inherently more important and should be highlighted. Also, a novel measure based on the ratio of available solvent accessible surface (RASAS) is proposed for monitoring the ligand/receptor binding process.

Several in vitro studies have shown the presence of an affinity gradient in nuclear pore complex proteins for the import receptor Importinβ, at least partially contributing to nucleocytoplasmic transport, while others have historically argued against the presence of such a gradient. Nonetheless, the existence of an affinity gradient has remained an uncharacterized contributing factor. To shed light on the affinity gradient theory and better characterize how the existence of such an affinity gradient between the nuclear pore and the import receptor may influence the nucleocytoplasmic traffic, we have developed a general-purpose agent based modeling (ABM) framework that features a new method for relating rate constants to molecular binding and unbinding probabilities, and used our ABM approach to quantify the effects of a wide range of forward and reverse nucleoporin-Importinβ affinity gradients. Our results indicate that transport through the nuclear pore complex is maximized with an effective macroscopic affinity gradient of 2000 µM, 200 µM and 10 µM in the cytoplasmic, central channel and nuclear basket respectively. The transport rate at this gradient is approximately 10% higher than the transport rate for a comparable pore lacking any affinity gradient, which has a peak transport rate when all nucleoporins have an affinity of 200 µM for Importinβ. Furthermore, this optimal ratio of affinity gradients is representative of the ratio of affinities reported for the yeast nuclear pore complex--suggesting that the affinity gradient seen in vitro is highly optimized.

Integrin αIIbβ3 is a member of the integrin family of transmembrane proteins present on the plasma membrane of platelets. Integrin αIIbβ3 is widely known to regulate the process of thrombosis via activation at its cytoplasmic side by talin and interaction with the soluble fibrinogen. It is also reported that three groups of interactions restrain integrin family members in the inactive state, including a set of salt bridges on the cytoplasmic side of the transmembrane domain of the integrin α- and β-subunits known as the inner membrane clasp, hydrophobic packing of a few transmembrane residues on the extracellular side between the α- and β-subunits that is known as the outer membrane clasp, and the key interaction group of the βA domain (located on the β-subunit head domain) with the βTD (proximal to the plasma membrane on the β-subunit). However, molecular details of this key interaction group as well as events that lead to detachment of the βTD and βA domains have remained ambiguous. In this study, we use molecular dynamics models to take a comprehensive outside-in and inside-out approach at exploring how integrin αIIbβ3 is activated. First, we show that talin's interaction with the membrane-proximal and membrane-distal regions of integrin cytoplasmic-transmembrane domains significantly loosens the inner membrane clasp. Talin also interacts with an additional salt bridge (R734-E1006), which facilitates integrin activation through the separation of the integrin's α- and β-subunits. The second part of our study classifies three types of interactions between RGD peptides and the extracellular domains of integrin αIIbβ3. Finally, we show that the interaction of the Arg of the RGD sequence may activate integrin via disrupting the key interaction group between K350 on the βA domain and S673/S674 on the βTD., (Copyright © 2013 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Vinculin can interact with F-actin both in recruitment of actin filaments to the growing focal adhesions and also in capping of actin filaments to regulate actin dynamics. Using molecular dynamics, both interactions are simulated using different vinculin conformations. Vinculin is simulated either with only its vinculin tail domain (Vt), with all residues in its closed conformation, with all residues in an open I conformation, and with all residues in an open II conformation. The open I conformation results from movement of domain 1 away from Vt; the open II conformation results from complete dissociation of Vt from the vinculin head domains. Simulation of vinculin binding along the actin filament showed that Vt alone can bind along the actin filaments, that vinculin in its closed conformation cannot bind along the actin filaments, and that vinculin in its open I conformation can bind along the actin filaments. The simulations confirm that movement of domain 1 away from Vt in formation of vinculin 1 is sufficient for allowing Vt to bind along the actin filament. Simulation of Vt capping actin filaments probe six possible bound structures and suggest that vinculin would cap actin filaments by interacting with both S1 and S3 of the barbed-end, using the surface of Vt normally occluded by D4 and nearby vinculin head domain residues. Simulation of D4 separation from Vt after D1 separation formed the open II conformation. Binding of open II vinculin to the barbed-end suggests this conformation allows for vinculin capping. Three binding sites on F-actin are suggested as regions that could link to vinculin. Vinculin is suggested to function as a variable switch at the focal adhesions. The conformation of vinculin and the precise F-actin binding conformation is dependent on the level of mechanical load on the focal adhesion.

Integrin clustering plays a pivotal role in a host of cell functions. Hetero-dimeric integrin adhesion receptors regulate cell migration, survival, and differentiation by communicating signals bidirectionally across the plasma membrane. Thus far, crystallographic structures of integrin components are solved only separately, and for some integrin types. Also, the sequence of interactions that leads to signal transduction remains ambiguous. Particularly, it remains controversial whether the homo-dimerization of integrin transmembrane domains occurs following the integrin activation (i.e. when integrin ectodomain is stretched out) or if it regulates integrin clustering. This study employs molecular dynamics modeling approaches to address these questions in molecular details and sheds light on the crucial effect of the plasma membrane. Conducting a normal mode analysis of the intact αllbβ3 integrin, it is demonstrated that the ectodomain and transmembrane-cytoplasmic domains are connected via a membrane-proximal hinge region, thus merely transmembrane-cytoplasmic domains are modeled. By measuring the free energy change and force required to form integrin homo-oligomers, this study suggests that the β-subunit homo-oligomerization potentially regulates integrin clustering, as opposed to α-subunit, which appears to be a poor regulator for the clustering process. If α-subunits are to regulate the clustering they should overcome a high-energy barrier formed by a stable lipid pack around them. Finally, an outside-in activation-clustering scenario is speculated, explaining how further loading the already-active integrin affects its homo-oligomerization so that focal adhesions grow in size.

The mechanisms by which living cells respond to mechanical stimuli are not yet fully understood. It has been suggested that mechanosensing proteins play an important role in mechanotransduction because their binding affinities are directly affected by the external stress. α-Actinin is an actin cross-linker and may act as a mechanosensor in adhesion sites. Its interaction with vinculin is suggested to be mechanically regulated. In this study, the free energy of activation is explored using the umbrella sampling method. An activation trajectory is generated in which α-actinin's vinculin-binding site swings out of the rod domain, leading to approximately an 8 kcal/mol free energy release. The activation trajectory reveals several local and global conformational changes along the activation pathway accompanied by the breakage of a number of key interactions stabilizing the inhibited structure. These results may shed light on the role of α-actinin in cellular mechanotransduction and focal adhesion formation., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

With a range of desirable mechanical and optical properties, single wall carbon nanotubes (SWCNTs) are a promising material for nanobiotechnologies. SWCNTs also have potential as biomaterials for modulation of cellular structures. Previously, we showed that highly purified, dispersed SWCNTs grossly alter F-actin inside cells. F-actin plays critical roles in the maintenance of cell structure, force transduction, transport and cytokinesis. Thus, quantification of SWCNT-actin interactions ranging from molecular, sub-cellular and cellular levels with both structure and function is critical for developing SWCNT-based biotechnologies. Further, this interaction can be exploited, using SWCNTs as a unique actin-altering material. Here, we utilized molecular dynamics simulations to explore the interactions of SWCNTs with actin filaments. Fluorescence lifetime imaging microscopy confirmed that SWCNTs were located within ~5 nm of F-actin in cells but did not interact with G-actin. SWCNTs did not alter myosin II sub-cellular localization, and SWCNT treatment in cells led to significantly shorter actin filaments. Functionally, cells with internalized SWCNTs had greatly reduced cell traction force. Combined, these results demonstrate direct, specific SWCNT alteration of F-actin structures which can be exploited for SWCNT-based biotechnologies and utilized as a new method to probe fundamental actin-related cellular processes and biophysics.

Vinculin phosphorylation has been implicated as a potential mechanism for focal adhesion growth and maturation. Four vinculin residues-Y100, S1033, S1045, and Y1065-are phosphorylated by kinases during focal adhesion maturation. In this study, phosphorylation at each of these residues is simulated using molecular dynamics models. The simulations demonstrate that once each phosphorylated vinculin structure is at equilibrium, significant local conformational changes result that may impact either vinculin activation or vinculin binding to actin and PIP2. Simulation of vinculin activation after phosphorylation shows that the added phosphoryl groups can prime vinculin for activation. It remains to be seen if vinculin can be phosphorylated at S1033 in vivo, but these simulations highlight that in the event of a S1033 phophorylation vinculin will likely be primed for activation., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Microtubule bundles cross-linked by tau protein serve a variety of neurological functions including maintaining mechanical integrity of the axon, promoting axonal growth, and facilitating cargo transport. It has been observed that axonal damage in traumatic brain injury leads to bundle disorientation, loss of axonal viability, and cognitive impairment. This study investigates the initial mechanical response of axonal microtubule bundles under uniaxial tension using a discrete bead-spring representation. Mechanisms of failure due to traumatic stretch loading and their impact on the mechanical response and stability are also characterized. This study indicates that cross-linked axonal microtubule bundles in tension display stiffening behavior similar to a power-law relationship from nonaffine network deformations. Stretching of cross-links and microtubule bending were the primary deformation modes at low stresses. Microtubule stretch was negligible up to tensile stresses of ∼1 MPa. Bundle failure occurred by failure of cross-links leading to pull-out of microtubules and loss of bundle integrity. This may explain the elongation, undulation, and delayed elasticity of axons following traumatic stretch loading. More extensively cross-linked bundles withstood higher tensile stresses before failing. The bundle mechanical behavior uncovered by these computational techniques should guide future experiments on stretch-injured axons., (Copyright © 2012 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

The dynamic nature of cell adhesion and detachment is critically important to a variety of physiological and pathophysiological phenomena. Much, however, still remains uncertain and controversial about the mechanochemical players and processes involved in cellular adhesion and detachment. This leads to the need for quantitative characterization of the adhesion and detachment of anchorage-dependent cells. Here, cell adhesion and detachment up to subcellular level are examined using gold surfaces modified with a thiol-functionalized arginine-glycine-aspartic acid (RGD) peptide. A thiol self-assembled monolayer (SAM) on top of the gold surfaces is reductively desorbed with activation potential to spatiotemporally manipulate both cell adhesion and detachment. This method maintains cells of interest living and intact during experiments, making it possible to quantify cell adhesion and detachment as close as possible to in vivo conditions. Experimental characterizations for NIH 3T3 fibroblasts are carried out with a focus on the following issues: the effect of the size and geometric shape of gold surfaces on cell adhesion; the effect of cell confluency, cell shape, and activation potential magnitude on cell detachment; changes in the material properties of cells during cell detachment. The findings of this study should lead to better understanding of cellular dynamics in anchorage-dependent cells., (Copyright © 2011 Elsevier Ltd. All rights reserved.)

The nuclear pore complex (NPC) regulates molecular traffic across the nuclear envelope (NE). Selective transport happens on the order of milliseconds and the length scale of tens of nanometers; however, the transport mechanism remains elusive. Central to the transport process is the hydrophobic interactions between karyopherins (kaps) and Phe-Gly (FG) repeat domains. Taking into account the polymeric nature of FG-repeats grafted on the elastic structure of the NPC, and the kap-FG hydrophobic affinity, we have established a coarse-grained model of the NPC structure that mimics nucleocytoplasmic transport. To establish a foundation for future works, the methodology and biophysical rationale behind the model is explained in details. The model predicts that the first-passage time of a 15 nm cargo-complex is about 2.6±0.13 ms with an inverse Gaussian distribution for statistically adequate number of independent Brownian dynamics simulations. Moreover, the cargo-complex is primarily attached to the channel wall where it interacts with the FG-layer as it passes through the central channel. The kap-FG hydrophobic interaction is highly dynamic and fast, which ensures an efficient translocation through the NPC. Further, almost all eight hydrophobic binding spots on kap-β are occupied simultaneously during transport. Finally, as opposed to intact NPCs, cytoplasmic filaments-deficient NPCs show a high degree of permeability to inert cargos, implying the defining role of cytoplasmic filaments in the selectivity barrier.

The nuclear pore complex (NPC) is the gatekeeper of the nucleus, capable of actively discriminating between the active and inert cargo while accommodating a high rate of translocations. The biophysical mechanisms underlying transport, however, remain unclear due to the lack of information about biophysical factors playing role in transport. Based on published experimental data, we have established a coarse-grained model of an intact NPC structure to examine nucleocytoplasmic transport with refined spatial and temporal resolutions. Using our model, we estimate the transport time versus cargo sizes. Our findings suggest that the mean transport time of cargos smaller than 15 nm is independent of size, while beyond this size, there is a sharp increase in the mean transport time. The model confirms that kap-FG hydrophobicity is sufficient for active cargo transport. Moreover, our model predicts that during translocation, small and large cargo-complexes are hydrophobically attached to FG-repeat domains for 86 and 96% of their transport time, respectively. Inside the central channel FG-repeats form a thick layer on the wall leaving an open tube. The cargo-complex is almost always attached to this layer and diffuses back and forth, regardless of the cargo size. Finally, we propose a plausible model for transport in which the NPC can be viewed as a lubricated gate. This model incorporates basic assumptions underlying virtual-gate and reduction-of-dimensionality models with the addition of the FG-layer inside the central channel acting as a lubricant., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Focal adhesions are critical to a number of cellular processes that involve mechanotransduction and mechanical interaction with the cellular environment. The growth and strengthening of these focal adhesions is dependent on the interaction between talin and vinculin. This study investigates said interaction and how vinculin activation influences it. Using molecular dynamics, the interaction between talin's vinculin binding site (VBS) and vinculin's domain 1 (D1) is simulated both before and after vinculin activation. The simulations of VBS binding to vinculin before activation suggest the proximity of the vinculin tail to D1 prevents helical movement in D1 and thus prevents binding of VBS. In contrast, interaction of VBS with activated vinculin shows the possibility of complete VBS insertion into D1. In the simulations of both activated and autoinhibited vinculin where VBS fails to fully bind, VBS demonstrates significant hydrophobic interaction with surface residues in D1. These interactions link VBS to D1 even without its proper insertion into the hydrophobic core. Together these simulations suggest VBS binds to vinculin with the following mechanism: 1), VBS links to D1 via surface hydrophobic interactions; 2), vinculin undergoes activation and D1 is moved away from the vinculin tail; 3), helices in D1 undergo conformational change to allow VBS binding; and 4), VBS inserts itself into the hydrophobic core of D1., (Copyright © 2011 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Quantitative description of cell mechanics has challenged biological scientists for the past two decades. Various structural models have been attempted to analyze the structure of the cytoskeleton. One important aspect that has been largely ignored in all these modeling approaches is related to the flexural and buckling behavior of microtubular filaments. The objective of this paper is to explore the influence of this flexural and buckling behavior in cytoskeletal mechanics.In vitro the microtubules are observed to buckle in the first mode, reminiscent of a free, simply-supported beam. In vivo images of microtubules, however, indicate that the buckling mostly occurs in higher modes. This buckling mode switch takes place mostly because of the lateral support of microtubules via their connections to actin and intermediate filaments. These lateral loads are exerted throughout the microtubule length and yield a considerable bending behavior that, unless properly accounted for, would produce erroneous results in the modeling and analysis of the cytoskeletal mechanics.One of the promising attempts towards mechanical modeling of the cytoskeleton is the tensegrity model, which simplifies the complex network of cytoskeletal filaments into a combination merely of tension-bearing actin filaments and compression-bearing microtubules. Interestingly, this discrete model can qualitatively explain many experimental observations in cell mechanics. However, evidence suggests that the simplicity of this model may undermine the accuracy of its predictions, given the model's underlying assumption that "every single member bears solely either tensile or compressive behavior," i.e. neglecting the flexural behavior of the microtubule filaments. We invoke an anisotropic continuum model for microtubules and compare the bending energy stored in a single microtubule with its axial strain energy at the verge of buckling. Our results suggest that the bending energy can exceed the axial energy of microtubules by 40 folds. A modification to tensegrity model is, therefore, proved necessary in order to take into account the flexural response of microtubules. The concept of "bendo-tensegrity" is proposed as a modification to contemporary cytoskeletal tensegrity models.

Diffusion plays a key role in many biochemical reaction systems seen in nature. Scenarios where diffusion behavior is critical can be seen in the cell and subcellular compartments where molecular crowding limits the interaction between particles. We investigate the application of a computational method for modeling the diffusion of molecules and macromolecules in three-dimensional solutions using agent based modeling. This method allows for realistic modeling of a system of particles with different properties such as size, diffusion coefficients, and affinity as well as the environment properties such as viscosity and geometry. Simulations using these movement probabilities yield behavior that mimics natural diffusion. Using this modeling framework, we simulate the effects of molecular crowding on effective diffusion and have validated the results of our model using Langevin dynamics simulations and note that they are in good agreement with previous experimental data. Furthermore, we investigate an extension of this framework where single discrete cells can contain multiple particles of varying size in an effort to highlight errors that can arise from discretization that lead to the unnatural behavior of particles undergoing diffusion. Subsequently, we explore various algorithms that differ in how they handle the movement of multiple particles per cell and suggest an algorithm that properly accommodates multiple particles of various sizes per cell that can replicate the natural behavior of these particles diffusing. Finally, we use the present modeling framework to investigate the effect of structural geometry on the directionality of diffusion in the cell cytoskeleton with the observation that parallel orientation in the structural geometry of actin filaments of filopodia and the branched structure of lamellipodia can give directionality to diffusion at the filopodia-lamellipodia interface.

Understanding the biomechanical properties and the effect of biomechanical force on epithelial cells is key to understanding how epithelial cells form uniquely shaped structures in two or three-dimensional space. Nevertheless, with the limitations and challenges posed by biological experiments at this scale, it becomes advantageous to use mathematical and 'in silico' (computational) models as an alternate solution. This paper introduces a single-cell-based model representing the cross section of a typical tissue. Each cell in this model is an individual unit containing several sub-cellular elements, such as the elastic plasma membrane, enclosed viscoelastic elements that play the role of cytoskeleton, and the viscoelastic elements of the cell nucleus. The cell membrane is divided into segments where each segment (or point) incorporates the cell's interaction and communication with other cells and its environment. The model is capable of simulating how cells cooperate and contribute to the overall structure and function of a particular tissue; it mimics many aspects of cellular behavior such as cell growth, division, apoptosis and polarization. The model allows for investigation of the biomechanical properties of cells, cell-cell interactions, effect of environment on cellular clusters, and how individual cells work together and contribute to the structure and function of a particular tissue. To evaluate the current approach in modeling different topologies of growing tissues in distinct biochemical conditions of the surrounding media, we model several key cellular phenomena, namely monolayer cell culture, effects of adhesion intensity, growth of epithelial cell through interaction with extra-cellular matrix (ECM), effects of a gap in the ECM, tensegrity and tissue morphogenesis and formation of hollow epithelial acini. The proposed computational model enables one to isolate the effects of biomechanical properties of individual cells and the communication between cells and their microenvironment while simultaneously allowing for the formation of clusters or sheets of cells that act together as one complex tissue.

Vinculin activation plays a critical role in focal adhesion initiation and formation. In its native state, vinculin is in an autoinhibitory conformation in which domain 1 prevents interaction of the vinculin tail domain with actin by steric hindrance. Once activated, vinculin is able to interact with both actin and talin. Several hypotheses have been put forth addressing the mechanisms of vinculin activation. One set of studies suggests that vinculin interaction with talin is sufficient to cause activation, whereas another set of studies suggests that a simultaneous interaction with several binding partners is necessary to achieve vinculin activation. Using molecular-dynamics (MD) simulations, we investigate the mechanisms of vinculin activation and suggest both a trajectory of conformational changes leading to vinculin activation, and key structural features that are likely involved in stabilizing the autoinhibited conformation. Assuming that the simultaneous interaction of vinculin with both actin and talin causes a stretching force on vinculin, and that vinculin activation results from a removal of steric hindrance blocking the actin-binding sites, we simulate with MD the stretching and activation of vinculin. The MD simulations are further confirmed by normal-mode analysis and simulation after residue modification. Taken together, the results of these simulations suggest that bending of the vinculin-binding-site region in vinculin away from the vinculin tail is the likely trajectory of vinculin activation., (2010 Biophysical Society. Published by Elsevier Inc. All rights reserved.)

Filamins are actin binding proteins that contribute to cytoskeletal integrity and biochemical scaffolds during mechanochemical signal transductions. Structurally, human filamins are dimers composed of an actin-binding domain with 24 immunoglobulin (Ig)-like repeats. In this study, we focus on the recently solved high-resolution crystal structure of Ig-like repeats 19-21 of filamin-A (IgFLNa-R19-R21). IgFLNa-R19-21 is of marked importance because it contains the binding site for integrins and facilitates the dynamic ability of filamin-A to communicate with the extracellular environment. However, the structure of filamin-A shows an interesting domain arrangement where the integrin binding site on IgFLNa-R21 is hindered sterically by IgFLNa-R20. Thus, a number of hypotheses on the regulation of filamin-A exist. Using molecular dynamics simulations we evaluated the effects of two primary regulators of filamin-A, force and phosphorylation. We find that a tensile force of 40 pN is sufficient to initiate the partial removal of the autoinhibition on the integrin binding site of IgFLNa-R21. Force coupled to phosphorylation at Ser(2152), however, affords complete dissociation of autoinhibition with a decreased force requirement. Phosphorylation seems to decrease the threshold for removing the IgFLNa-R20 beta-strand inhibitor within 300 ps with 40 pN tensile force. Furthermore, the molecular dynamic trajectories illustrate phosphorylation of Ser(2152) without force is insufficient to remove autoinhibition. We believe the results of this study implicate filamin-A as a tunable mechanosensor, where its sensitivity can be modulated by the degree of phosphorylation.

The process of epithelial morphogenesis is ubiquitous in animal development, but much remains to be learned about the mechanisms that shape epithelial tissues. The follicle cell (FC) epithelium encapsulating the growing germline of Drosophila is an excellent system to study fundamental elements of epithelial development. During stages 8 to 10 of oogenesis, the FC epithelium transitions between simple geometries-cuboidal, columnar and squamous-and redistributes cell populations in processes described as posterior migration, squamous cell flattening and main body cell columnarization. Here we have carried out a quantitative morphometric analysis of these poorly understood events in order to establish the parameters of and delimit the potential processes that regulate the transitions. Our results compel a striking revision of accepted views of these phenomena, by showing that posterior migration does not involve FC movements, that there is no role for columnar cell apical constriction in FC morphogenesis, and that squamous cell flattening may be a compliant response to germline growth. We utilize mechanical modeling involving finite element computational technologies to demonstrate that time-varying viscoelastic properties and growth are sufficient to account for the bulk of the FC morphogenetic changes.

The aortic heart valve undergoes geometric and mechanical changes over time. The cusps of a normal, healthy valve thicken and become less extensible over time. In the disease calcific aortic stenosis (CAS), calcified nodules progressively stiffen the cusps. The local mechanical changes in the cusps, due to either normal aging or pathological processes, affect overall function of the valve. In this paper, we propose a computational model for the aging aortic valve that connects local changes to overall valve function. We extend a previous model for the healthy valve to describe aging. To model normal/uncomplicated aging, leaflet thickness and extensibility are varied versus age according to experimental data. To model calcification, initial sites are defined and a simple growth law is assumed. The nodules then grow over time, so that the area of calcification increases from one model to the next model representing greater age. Overall valve function is recorded for each individual model to yield a single simulation of valve function over time. This simulation is the first theoretical tool to describe the temporal behavior of aortic valve calcification. The ability to better understand and predict disease progression will aid in design and timing of patient treatments for CAS.

alpha-Actinin is an actin crosslinking molecule that can serve as a scaffold and maintain dynamic actin filament networks. As a crosslinker in the stressed cytoskeleton, alpha-actinin can retain conformation, function, and strength. alpha-Actinin has an actin binding domain and a calmodulin homology domain separated by a long rod domain. Using molecular dynamics and normal mode analysis, we suggest that the alpha-actinin rod domain has flexible terminal regions which can twist and extend under mechanical stress, yet has a highly rigid interior region stabilized by aromatic packing within each spectrin repeat, by electrostatic interactions between the spectrin repeats, and by strong salt bridges between its two anti-parallel monomers. By exploring the natural vibrations of the alpha-actinin rod domain and by conducting bending molecular dynamics simulations we also predict that bending of the rod domain is possible with minimal force. We introduce computational methods for analyzing the torsional strain of molecules using rotating constraints. Molecular dynamics extension of the alpha-actinin rod is also performed, demonstrating transduction of the unfolding forces across salt bridges to the associated monomer of the alpha-actinin rod domain.

Cells can sense mechanical force in regulating focal adhesion assembly. One vivid example is the force-induced recruitment of vinculin to reinforce initial contacts between a cell and the extracellular matrix. Crystal structures of the unbound proteins and bound complex between the vinculin head subdomain (Vh1) and the talin vinculin binding site 1 (VBS1) indicate that vinculin undergoes a conformational change upon binding to talin. However, the molecular basis for this event and the precise nature of the binding pathway remain elusive. In this article, molecular dynamics is used to investigate the binding mechanism of Vh1 and VBS1 under minimal constraints to facilitate binding. One simulation demonstrates binding of the two molecules in the complete absence of external force. VBS1 makes early hydrophobic contact with Vh1 by positioning the critical hydrophobic residues (L608, L615, and L622) in the groove formed by helices 1 and 2 of Vh1. The solvent-exposed hydrophobic residues (V619 and L623) then gradually penetrate the hydrophobic core of Vh1, thus further separating helix 1 from helix 2. These critical residues are highly conserved as large hydrophobic side groups in other vinculin binding sites; studies also have demonstrated that these residues are essential in Vh1-VBS1 binding. Similar binding mechanisms are also demonstrated in separate molecular dynamics simulations of Vh1 binding to other vinculin binding sites both in talin and alpha-actinin.

The basic structure of the nuclear pore complex (NPC), conserved across almost all organisms from yeast to humans, persists in featuring an octagonal symmetry involving the nucleoporins that constitute the NPC ring. In this article, we seek to understand and evaluate the potential biomechanical reasons for this eightfold symmetry. Our analytical investigation shows that the eightfold symmetry maximizes the bending stiffness of each of the eight NPC spokes while our computational analyses identify the most likely deformation modes, frequencies, and associated kinetic energies of the NPC. These modes have energies close to other published findings using membrane analysis of the nuclear membrane pore opening, and deformation states in agreement with experimental observations. A better understanding of NPC mechanics is essential for characterizing the nucleocytoplasmic transport, which has a central importance in cell biology.

Rearrangement of the actin cytoskeleton is integral to cell shape and function. Actin-binding proteins, e.g., filamin, can naturally contribute to the mechanics and function of the actin cytoskeleton. The molecular mechanical bases for filamin's function in actin cytoskeletal reorganization are examined here using molecular dynamics simulations. Simulations are performed by applying forces ranging from 25 pN to 125 pN for 2.5 ns to the rod domain of filamin. Applying small loads ( approximately 25 pN) to filamin's rod domain supplies sufficient energy to alter the conformation of the N-terminal regions of the rod. These forces break local hydrogen bond coordination often enough to allow side chains to find new coordination partners, in turn leading to drastic changes in the conformation of filamin, for example, increasing the hydrophobic character of the N-terminal rod region and, alternatively, activating the C-terminal region to become increasingly stiff. These changes in conformation can lead to changes in the affinity of filamin for its binding partners. Therefore, filamin can function to transduce mechanical signals as well as preserve topology of the actin cytoskeleton throughout the rod domain.

Force-induced changes in protein conformation are thought to be responsible for certain cellular responses to mechanical force. Changes in conformation subsequently initiate a biochemical response by alterations in, for example, binding affinity to another protein or enzymatic activity. Here, a model of protein extension under external forcing is created inspired by Kramers' theory for reaction rate kinetics in liquids. The protein is assumed to have two distinct conformational states: a relaxed state, C(1), preferred in the absence of external force, and an extended state, C(2), favored under force application. In the context of mechanotransduction, the extended state is a conformation from which the protein can initiate signaling. Appearance and persistence of C(2) are assumed to lead to transduction of the mechanical signal into a chemical one. The protein energy landscape is represented by two harmonic wells of stiffness kappa(1) and kappa(2), whose minima correspond to conformations C(1) and C(2). First passage time t(f) from C(1) to C(2) is determined from the Fokker-Plank equation employing several different approaches found in the literature. These various approaches exhibit significant differences in behavior as force increases. Although the level of applied force and the energy difference between states largely determine equilibrium, the dominant influence on t(f) is the height of the transition state. Distortions in the energy landscape due to force can also have a significant influence, however, exhibiting a weaker force dependence than exponential as previously reported, approaching a nearly constant value at a level of force that depends on the ratio kappa(1)/kappa(2). Two model systems are used to demonstrate the utility of this approach: a short alpha-helix undergoing a transition between two well-defined states and a simple molecular motor.

Living cells respond to mechanical stimulation in a variety of ways that affect nearly every aspect of their function. Such responses can range from changes in cell morphology to activation of signaling cascades and changes in cell phenotype. Although the biochemical signaling pathways activated by mechanical stimulus have been extensively studied, little is known of the basic mechanisms by which mechanical force is transduced into a biochemical signal, or how the cell changes its behavior or properties in response to external or internal stresses. One hypothesis is that forces transmitted via individual proteins either at the site of cell adhesion to its surroundings or within the stress-bearing members of the cytoskeleton cause conformational changes that alter their binding affinity to other intracellular molecules. This altered equilibrium state can subsequently either initiate a biochemical signaling cascade or produce more immediate and local structural changes. To understand the phenomena related to mechanotransduction, the mechanics and chemistry of single molecules that form the signal transduction pathways must be examined. This paper presents a range of case studies that seek to explore the molecular basis of mechanical signal sensation and transduction, with particular attention to their macroscopic manifestation in the cell properties, e.g. in focal adhesion remodeling due to local application of force or changes in cytoskeletal rheology and remodeling due to cellular deformation.

Atherosclerotic lesion localization to regions of disturbed flow within certain arterial geometries, in humans and experimental animals, suggests an important role for local hemodynamic forces in atherogenesis. To explore how endothelial cells (EC) acquire functional/dysfunctional phenotypes in response to vascular region-specific flow patterns, we have used an in vitro dynamic flow system to accurately reproduce arterial shear stress waveforms on cultured human EC and have examined the effects on EC gene expression by using a high-throughput transcriptional profiling approach. The flow patterns in the carotid artery bifurcations of several normal human subjects were characterized by using 3D flow analysis based on actual vascular geometries and blood flow profiles. Two prototypic arterial waveforms, "athero-prone" and "athero-protective," were defined as representative of the wall shear stresses in two distinct regions of the carotid artery (carotid sinus and distal internal carotid artery) that are typically "susceptible" or "resistant," respectively, to atherosclerotic lesion development. These two waveforms were applied to cultured EC, and cDNA microarrays were used to analyze the differential patterns of EC gene expression. In addition, the differential effects of athero-prone vs. athero-protective waveforms were further characterized on several parameters of EC structure and function, including actin cytoskeletal organization, expression and localization of junctional proteins, activation of the NF-kappaB transcriptional pathway, and expression of proinflammatory cytokines and adhesion molecules. These global gene expression patterns and functional data reveal a distinct phenotypic modulation in response to the wall shear stresses present in atherosclerosis-susceptible vs. atherosclerosis-resistant human arterial geometries.

Vascular endothelial cells rapidly transduce local mechanical forces into biological signals through numerous processes including the activation of focal adhesion sites. To examine the mechanosensing capabilities of these adhesion sites, focal adhesion translocation was monitored over the course of 5 min with GFP-paxillin while applying nN-level magnetic trap shear forces to the cell apex via integrin-linked magnetic beads. A nongraded steady-load threshold for mechanotransduction was established between 0.90 and 1.45 nN. Activation was greatest near the point of forcing (<7.5 microm), indicating that shear forces imposed on the apical cell membrane transmit nonuniformly to the basal cell surface and that focal adhesion sites may function as individual mechanosensors responding to local levels of force. Results from a continuum, viscoelastic finite element model of magnetocytometry that represented experimental focal adhesion attachments provided support for a nonuniform force transmission to basal surface focal adhesion sites. To further understand the role of force transmission on focal adhesion activation and dynamics, sinusoidally varying forces were applied at 0.1, 1.0, 10, and 50 Hz with a 1.45 nN offset and a 2.25 nN maximum. At 10 and 50 Hz, focal adhesion activation did not vary with spatial location, as observed for steady loading, whereas the response was minimized at 1.0 Hz. Furthermore, applying the tyrosine kinase inhibitors genistein and PP2, a specific Src family kinase inhibitor, showed tyrosine kinase signaling has a role in force-induced translocation. These results highlight the mutual importance of force transmission and biochemical signaling in focal adhesion mechanotransduction.

A three-dimensional viscoelastic finite element model is developed for cell micromanipulation by magnetocytometry. The model provides a robust tool for analysis of detailed strain/stress fields induced in the cell monolayer produced by forcing one microbead attached atop a single cell or cell monolayer on a basal substrate. Both the membrane/cortex and the cytoskeleton are modeled as Maxwell viscoelastic materials, but the structural effect of the membrane/cortex was found to be negligible on the timescales corresponding to magnetocytometry. Numerical predictions are validated against experiments performed on NIH 3T3 fibroblasts and previous experimental work. The system proved to be linear with respect to cytoskeleton mechanical properties and bead forcing. Stress and strain patterns were highly localized, suggesting that the effects of magnetocytometry are confined to a region extending <10 microm from the bead. Modulation of cell height has little effect on the results, provided the monolayer is >5 micro m thick. NIH 3T3 fibroblasts exhibited a viscoelastic timescale of approximately 1 s and a shear modulus of approximately 1000 Pa.