Molecular dynamics simulations of folding processes of a beta-hairpin in an implicit solvent

Computer simulations of beta-hairpin folding are relatively difficult, especially those based on the explicit water model. This greatly limits the complete analysis and understanding of their folding mechanisms. In this paper, we use the generalized Born/solvent accessible implicit solvent model to simulate the folding processes of a nine-residue beta-hairpin. We find that the beta-hairpin can fold into its native structure very easily, even using the traditional molecular dynamics method. This allows us to extract 21 complete folding events and investigate the folding process sufficiently. Our results show that there exist four most stable states on the free energy landscape of the short peptide, one native state and three intermediates. We find that two of the non-native stable states have almost the same potential energy as the native state but with lower entropy. This suggests that the native state can be stabilized entropically. Furthermore, we find that the folding processes of this peptide have common features: to fold into its native state, the peptide undergoes a continuous collapsing-extending-recollapsing process to adjust the positions of the side chains in order to form the native middle inter-strand hydrogen bonds. The formations of these bonds are the key step of the folding process. Once these bonds are formed, the peptide can fold into the native state quickly.     

Entropic elasticity controls nanomechanics of single tropocollagen molecules

We report molecular modeling of stretching single molecules of tropocollagen, the building block of collagen fibrils and fibers that provide mechanical support in connective tissues. For small deformation, we observe a dominance of entropic elasticity. At larger deformation, we find a transition to energetic elasticity, which is characterized by first stretching and breaking of hydrogen bonds, followed by deformation of covalent bonds in the protein backbone, eventually leading to molecular fracture. Our force-displacement curves at small forces show excellent quantitative agreement with optical tweezer experiments. Our model predicts a persistence length xi(p) approximately 16 nm, confirming experimental results suggesting that tropocollagen molecules are very flexible elastic entities. We demonstrate that assembly of single tropocollagen molecules into fibrils significantly decreases their bending flexibility, leading to decreased contributions of entropic effects during deformation. The molecular simulation results are used to develop a simple continuum model capable of describing an entire deformation range of tropocollagen molecules. Our molecular model is capable of describing different regimes of elastic and permanent deformation, without relying on empirical parameters, including a transition from entropic to energetic elasticity.     

Structural basis of compound recognition by adenosine deaminase

Structural snapshots corresponding to various states enable elucidation of the molecular recognition mechanism of enzymes. Adenosine deaminase has two distinct conformations, an open form and a closed form, although it has so far been unclear what factors influence adaptation of the alternative conformations. Herein, we have determined the first nonligated structure as an initial state, which was the open form, and have thereby rationally deduced the molecular recognition mechanism. Inspection of the active site in the nonligated and ligated states indicated that occupancy at one of the water-binding positions in the nonligated state was highly significant in determining alternate conformations. When this position is empty, subsequent movement of Phe65 toward the space induces the closed form. On the other hand, while occupied, the overall conformation remains in the open form. This structural understanding should greatly assist structure-oriented drug design and enable control of the enzymatic activity.     

NO rebinding to myoglobin: a reactive molecular dynamics study

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.     

Fundamental Characteristics of Neuron Adhesion Revealed by Forced Peeling and Time-Dependent Healing

A current bottleneck in the advance of neurophysics is the lack of reliable methods to quantitatively measure the interactions between neural cells and their microenvironment. Here, we present an experimental technique to probe the fundamental characteristics of neuron adhesion through repeated peeling of well-developed neurite branches on a substrate with an atomic force microscopy cantilever. At the same time, a total internal reflection fluorescence microscope is also used to monitor the activities of neural cell adhesion molecules (NCAMs) during detaching. It was found that NCAMs aggregate into clusters at the neurite-substrate interface, resulting in strong local attachment with an adhesion energy of ∼0.1 mJ/m² and sudden force jumps in the recorded force-displacement curve. Furthermore, by introducing a healing period between two forced peelings, we showed that stable neurite-substrate attachment can be re-established in 2–5 min. These findings are rationalized by a stochastic model, accounting for the breakage and rebinding of NCAM-based molecular bonds along the interface, and provide new insights into the mechanics of neuron adhesion as well as many related biological processes including axon outgrowth and nerve regeneration.     

Effect of Age and Cytoskeletal Elements on the Indentation-Dependent Mechanical Properties of Chondrocytes

Articular cartilage chondrocytes are responsible for the synthesis, maintenance, and turnover of the extracellular matrix, metabolic processes that contribute to the mechanical properties of these cells. Here, we systematically evaluated the effect of age and cytoskeletal disruptors on the mechanical properties of chondrocytes as a function of deformation. We quantified the indentation-dependent mechanical properties of chondrocytes isolated from neonatal (1-day), adult (5-year) and geriatric (12-year) bovine knees using atomic force microscopy (AFM). We also measured the contribution of the actin and intermediate filaments to the indentation-dependent mechanical properties of chondrocytes. By integrating AFM with confocal fluorescent microscopy, we monitored cytoskeletal and biomechanical deformation in transgenic cells (GFP-vimentin and mCherry-actin) under compression. We found that the elastic modulus of chondrocytes in all age groups decreased with increased indentation (15–2000 nm). The elastic modulus of adult chondrocytes was significantly greater than neonatal cells at indentations greater than 500 nm. Viscoelastic moduli (instantaneous and equilibrium) were comparable in all age groups examined; however, the intrinsic viscosity was lower in geriatric chondrocytes than neonatal. Disrupting the actin or the intermediate filament structures altered the mechanical properties of chondrocytes by decreasing the elastic modulus and viscoelastic properties, resulting in a dramatic loss of indentation-dependent response with treatment. Actin and vimentin cytoskeletal structures were monitored using confocal fluorescent microscopy in transgenic cells treated with disruptors, and both treatments had a profound disruptive effect on the actin filaments. Here we show that disrupting the structure of intermediate filaments indirectly altered the configuration of the actin cytoskeleton. These findings underscore the importance of the cytoskeletal elements in the overall mechanical response of chondrocytes, indicating that intermediate filament integrity is key to the non-linear elastic properties of chondrocytes. This study improves our understanding of the mechanical properties of articular cartilage at the single cell level.     

High Resolution,Large Deformation 3D Traction Force Microscopy

Traction Force Microscopy (TFM) is a powerful approach for quantifying cell-material interactions that over the last two decades has contributed significantly to our understanding of cellular mechanosensing and mechanotransduction. In addition, recent advances in three-dimensional (3D) imaging and traction force analysis (3D TFM) have highlighted the significance of the third dimension in influencing various cellular processes. Yet irrespective of dimensionality, almost all TFM approaches have relied on a linear elastic theory framework to calculate cell surface tractions. Here we present a new high resolution 3D TFM algorithm which utilizes a large deformation formulation to quantify cellular displacement fields with unprecedented resolution. The results feature some of the first experimental evidence that cells are indeed capable of exerting large material deformations, which require the formulation of a new theoretical TFM framework to accurately calculate the traction forces. Based on our previous 3D TFM technique, we reformulate our approach to accurately account for large material deformation and quantitatively contrast and compare both linear and large deformation frameworks as a function of the applied cell deformation. Particular attention is paid in estimating the accuracy penalty associated with utilizing a traditional linear elastic approach in the presence of large deformation gradients.     

Matrix representation of cell structure,functional processes and evolution: Part V of a general theory

The purpose of this paper is to formulate a compact analytical representation of cell structure, functional processes and evolution. In this formulation, the individual molecular structures are represented by their force-field surfaces. Complementary active sites on these surfaces permit molecular interactions. In cells, these interactions are further regulated by barrier systems in time, space, specificity, and energy. In terms of these parameters, evolution can be represented (modeled) as a random walk in a multi-dimensional space, subject to constraints. In this paper, the various parameters are integrated into a single compact matrix (stack) representation (a three index array). Cell life cycle and functional processes can be represented as a sequence of quantized, time-dependent changes in the representation matrix, subject to specified constraints. Cell evolution can be modeled by generating allowed matrix combinations. This theoretical approach has applications in: (1) ordering and interpreting experimental findings into the matrix representation. Missing matrix elements can be predicted, to be confirmed experimentally; (2) theoretical analysis and prediction of cell regulatory processes and the possible pathological failures; (3) theoretical derivation of the possible biological structures and functional processes, modeling possible pathways of cell and molecular evolution in terms of the matrix representation.     

An alternative model for the structural mechanics of hagfish slime threads

Hagfish slime threads are an interesting intermediate filament (IF) system in that, in contrast to hair, they appear to be a series system of IFs, consisting of rods and terminal domains (TDs), without matrix. The protein composition also differs from hair. Published data show that the wet stress-strain curve consists of four regions: I, increasing low modulus; II, plateau; III increasing higher modulus; IV decreasing modulus to break. Beyond 34% extension, there is plastic deformation, which contrasts with the elastic deformation of most biological fibres. The paper considers two explanations, one published by Fudge et al. [D.S. Fudge, K.H. Gardner, V.T. Forsyth, C. Riekel, J.M. Gosline, Biophys. J. 85 (2003) 2015] and one new, of deformation in the four regions. The former regards the TDs as elastic, entropic coils, which extend in region I. The latter regards the TDs as having a plastic, energy-dependent extension in regions III and IV, which is comparable to the drawing of polyester fibres. A rough theoretical model of jumps over energy barriers gives a similar prediction. Twist is suggested as a mechanism for overall thread cohesion. Experimental and theoretical ways forward to a greater understanding of the structural mechanics of hagfish threads are suggested. The behaviour of the total thread/mucus/water system is discussed and some speculations on the defensive mechanisms that have evolved are presented.     

Ion conductance vs. pore gating and selectivity in KcsA channel: Modeling achievements and perspectives

KcsA potassium channel belongs to a wide family of allosteric proteins that switch between closed and open states conformations in response to a stimulus, and act as a regulator of cation activity in living cells. The gating mechanism and cation selectivity of such channels have been extensively studied in the literature, with a revival emphasis these latter years, due to the publication of the crystallized structure of KcsA. Despite the increasing number of research and review papers on these topics, quantitative interpretation of these processes at the atomic scale is far from achieved. On the basis of available experimental and theoretical data, and by including our recent results, we review the progresses in this field of activity and discuss the weaknesses that should be corrected. In this spirit, we partition the channel into the filter, cavity, extra and intracellular media, in order to analyze separately the specificity of each region. Special emphasis is brought to the study of an open state for the channel and to the different properties generated by the opening. The influence of water as a structural and dynamical component of the channel properties in closed and open states, as well as in the sequential motions of the cations, is analyzed using molecular dynamics simulations and ab initio calculations. The polarization and charge transfer effects on the ions’ dynamics and kinetics are discussed in terms of partial charge models.     

Modeling and simulation of chemomechanics at the cell-matrix interface

Chemomechanical characteristics of the extracellular materials with which cells interact can have a profound impact on cell adhesion and migration. To understand and modulate such complex multiscale processes, a detailed understanding of the feedback between a cell and the adjacent microenvironment is crucial. Here, we use computational modeling and simulation to examine the cell-matrix interaction at both the molecular and continuum lengthscales. Using steered molecular dynamics, we consider how extracellular matrix (ECM) stiffness and extracellular pH influence the interaction between cell surface adhesion receptors and extracellular matrix ligands, and we predict potential consequences for focal adhesion formation and dissolution. Using continuum-level finite element simulations and analytical methods to model cell-induced ECM deformation as a function of ECM stiffness and thickness, we consider the implications toward design of synthetic substrata for cell biology experiments that intend to decouple chemical and mechanical cues.     

Fast Force Loading Disrupts Molecular Binding Stability in Human and Mouse Cell Adhesions

Force plays critical roles in cell adhesion and mechano-signaling, partially by regulating the dissociation rate, i.e., off-rate, of receptor-ligand bonds. However, the mechanism of such regulation still remains elusive. As a controversial topic of the field, when measuring the “off-rate vs. force” relation of the same molecular system, different dynamic force spectroscopy (DFS) assays (namely, force-clamp and force-ramp assays) often yield contradictive results. Such discrepancies hurdled our further understanding of molecular binding, and casted doubt on the existing theoretical models. In this work, we used a live-cell DFS technique, biomembrane force probe, to measure the single-bond dissociation in three receptor-ligand systems which respectively have important functions in vascular and immune systems: human platelet GPIbα-VWF, mouse T cell receptor-OVA peptide:MHC, and mouse platelet integrin αIIbβ3-fibrinogen. Using force-clamp and force-ramp assays in parallel, we identified that the force loading disrupted the stability of molecular bonds in a rate-dependent manner. This disruptive effect was achieved by the transitioning of bonds between two dissociation states: faster force loading induces more bonds to adopt the fast-dissociating state (and less to adopt the slow-dissociating state). Based on this mechanism, a new biophysical model of bond dissociation was established which took into account the effects of both force magnitude and loading rate. Remarkably, this model reconciled the results from the two assays in all three molecular systems under study. Our discoveries provided a new paradigm for understanding how force regulates receptor-ligand interactions and a guideline for the proper use of DFS technologies. Furthermore, our work highlighted the opportunity of using different DFS assays to answer specific biological questions in the field of cell adhesion and mechano-signaling     

Strength of hydrogen bond network takes crucial roles in the dissociation process of inhibitors from the HIV-1 protease binding pocket

To understand the underlying mechanisms of significant differences in dissociation rate constant among different inhibitors for HIV-1 protease, we performed steered molecular dynamics (SMD) simulations to analyze the entire dissociation processes of inhibitors from the binding pocket of protease at atomistic details. We found that the strength of hydrogen bond network between inhibitor and the protease takes crucial roles in the dissociation process. We showed that the hydrogen bond network in the cyclic urea inhibitors AHA001/XK263 is less stable than that of the approved inhibitor ABT538 because of their large differences in the structures of the networks. In the cyclic urea inhibitor bound complex, the hydrogen bonds often distribute at the flap tips and the active site. In contrast, there are additional accessorial hydrogen bonds formed at the lateral sides of the flaps and the active site in the ABT538 bound complex, which take crucial roles in stabilizing the hydrogen bond network. In addition, the water molecule W301 also plays important roles in stabilizing the hydrogen bond network through its flexible movement by acting as a collision buffer and helping the rebinding of hydrogen bonds at the flap tips. Because of its high stability, the hydrogen bond network of ABT538 complex can work together with the hydrophobic clusters to resist the dissociation, resulting in much lower dissociation rate constant than those of cyclic urea inhibitor complexes. This study may provide useful guidelines for design of novel potent inhibitors with optimized interactions.     

Geminate recombination of nitric oxide to endothelial nitric-oxide synthase and mechanistic implications.

The nitric-oxide synthase (NOS) catalyzes the oxidation of L-arginine to L-citrulline and NO through consumption of oxygen bound to the heme. Because NO is produced close to the heme and may bind to it, its subsequent role in a regulatory mechanism should be scrutinized. We therefore examined the kinetics of NO rebinding after photodissociation in the heme pocket of human endothelial NOS by means of time-resolved absorption spectroscopy. We show that geminate recombination of NO indeed occurs and that this process is strongly modulated by L-Arg. This NO rebinding occurs in a multiphasic fashion and spans over 3 orders of magnitude. In both ferric and ferrous states of the heme, a fast nonexponential picosecond geminate rebinding first takes place followed by a slower nanosecond phase. The rates of both phases decreased, whereas their relative amplitudes are changed by the presence of L-Arg; the overall effect is a slow down of NO rebinding. For the isolated oxygenase domain, the picosecond rate is unchanged, but the relative amplitude of the nanosecond binding decreased. We assigned the nanosecond kinetic component to the rebinding of NO that is still located in the protein core but not in the heme pocket. The implications for a mechanism of regulation involving NO binding are discussed.     

Mapping the Dynamics Landscape of Conformational Transitions in Enzyme: The Adenylate Kinase Case

Conformational transition describes the essential dynamics and mechanism of enzymes in pursuing their various functions. The fundamental and practical challenge to researchers is to quantitatively describe the roles of large-scale dynamic transitions for regulating the catalytic processes. In this study, we tackled this challenge by exploring the pathways and free energy landscape of conformational changes in adenylate kinase (AdK), a key ubiquitous enzyme for cellular energy homeostasis. Using explicit long-timescale (up to microseconds) molecular dynamics and bias-exchange metadynamics simulations, we determined at the atomistic level the intermediate conformational states and mapped the transition pathways of AdK in the presence and absence of ligands. There is clearly chronological operation of the functional domains of AdK. Specifically in the ligand-free AdK, there is no significant energy barrier in the free energy landscape separating the open and closed states. Instead there are multiple intermediate conformational states, which facilitate the rapid transitions of AdK. In the ligand-bound AdK, the closed conformation is energetically most favored with a large energy barrier to open it up, and the conformational population prefers to shift to the closed form coupled with transitions. The results suggest a perspective for a hybrid of conformational selection and induced fit operations of ligand binding to AdK. These observations, depicted in the most comprehensive and quantitative way to date, to our knowledge, emphasize the underlying intrinsic dynamics of AdK and reveal the sophisticated conformational transitions of AdK in fulfilling its enzymatic functions. The developed methodology can also apply to other proteins and biomolecular systems.     

Control of nitric oxide dynamics by guanylate cyclase in its activated state

Soluble guanylate cyclase (sGC) is the target of nitric oxide (NO) released by nitric-oxide synthase in endothelial cells, inducing an increase of cGMP synthesis in response. This heterodimeric protein possesses a regulatory subunit carrying a heme where NO binding occurs, while the second subunit harbors the catalytic site. The binding of NO and the subsequent breaking of the bond between the proximal histidine and the heme-Fe(2+) are assumed to induce conformational changes, which are the origin of the catalytic activation. At the molecular level, the activation and deactivation mechanisms are unknown, as is the dynamics of NO once in the heme pocket. Using ultrafast time-resolved absorption spectroscopy, we measured the kinetics of NO rebinding to sGC after photodissociation. The main spectral transient in the Soret band does not match the equilibrium difference spectrum of NO-liganded minus unliganded sGC, and the geminate rebinding was found to be monoexponential and ultrafast (tau = 7.5 ps), with a relative amplitude close to unity (0.97). These characteristics, so far not observed in other hemoproteins, indicate that NO encounters a high energy barrier for escaping from the heme pocket once the His-Fe(2+) bond has been cleaved; this bond does not reform before NO recombination. The deactivation of isolated sGC cannot occur by only simple diffusion of NO from the heme; therefore, several allosteric states may be inferred, including a desensitized one, to induce NO release. Thus, besides the structural change leading to activation, a consequence of the decoupling of the proximal histidine may also be to induce a change of the heme pocket distal geometry, which raises the energy barrier for NO escape, optimizing the efficiency of NO trapping. The non-single exponential character of the NO picosecond rebinding coexists only with the presence of the protein structure surrounding the heme, and the single exponential rate observed in sGC is very likely to be due to a closed conformation of the heme pocket. Our results emphasize the physiological importance of NO geminate recombination in hemoproteins like nitric-oxide synthase and sGC and show that the protein structure controls NO dynamics in a manner adapted to their function. This control of ligand dynamics provides a regulation at molecular level in the function of these enzymes.