Forced detachment of the CD2-CD58 complex

Shock-Induced Damage Mechanism of Perineuronal Nets

The perineuronal net (PNN) region of the brain’s extracellular matrix (ECM) surrounds the neural networks within the brain tissue. The PNN is a protective net-like structure regulating neuronal activity such as neurotransmission, charge balance, and action potential generation. Shock-induced damage of this essential component may lead to neuronal cell death and neurodegenerations. The shock generated during a vehicle accident, fall, or improvised device explosion may produce sufficient energy to damage the structure of the PNN. The goal is to investigate the mechanics of the PNN in reaction to shock loading and to understand the mechanical properties of different PNN components such as glycan, GAG, and protein. In this study, we evaluated the mechanical strength of PNN molecules and the interfacial strength between the PNN components. Afterward, we assessed the PNN molecules’ damage efficiency under various conditions such as shock speed, preexisting bubble, and boundary conditions. The secondary structure altercation of the protein molecules of the PNN was analyzed to evaluate damage intensity under varying shock speeds. At a higher shock speed, damage intensity is more elevated, and hyaluronan (glycan molecule) is most likely to break at the rigid junction. The primary structure of the protein molecules is least likely to fail. Instead, the molecules’ secondary bonds will be altered. Our study suggests that the number of hydrogen bonds during the shock wave propagation is reduced, which leads to the change in protein conformations and damage within the PNN structure. As such, we found a direct connection between shock wave intensity and PNN damage.

Structure-based identification of inhibitors disrupting the CD2-CD58 interactions

The immune system has very intricate mechanisms of fighting against the invading infections which are accomplished by a sequential event of molecular interactions in the body. One of the crucial phenomena in this process is the recognition of T-cells by the antigen-presenting cells (APCs), which is initiated by the rapid interaction between both cell surface receptors, i.e., CD2 located on T-cells and CD58 located on APCs. Under various pathological conditions, which involve undesired immune response, inhibiting the CD2-CD58 interactions becomes a therapeutically relevant opportunity. Herein we present an extensive work to identify novel inhibiting agents of the CD2-CD58 interactions. Classical molecular dynamics (MD) simulations of the CD2-CD58 complex highlighted a series of crucial CD58 residues responsible for the interactions with CD2. Based on such results, a pharmacophore map, complementary to the CD2-binding site of CD58, was created and employed for virtual screening of ~ 300,000 available compounds. On the ~ 6000 compounds filtered from pharmacophore mapping, ADME screening leads to ~ 350 molecules. Molecular docking was then performed on these molecules, and fifteen compounds emerged with significant binding energy (< - 50 kcal/mol) for CD58. Finally, short MD simulations were performed in triplicate on each complex (i) to provide a microscopic view of the ligand binding and (ii) to rule out possibly weak binders of CD58 from the identified hits. At last, we suggest eight compounds for in vitro testing that were identified as promising hits to bind CD58 with a high binding affinity.

Steered molecular dynamic simulations of conformational lock of Cu, Zn-superoxide dismutase

The conformational lock was a bio-thermodynamic theory to explain the characteristics of interfaces in oligomeric enzymes and their effects on catalytic activity. The previous studies on superoxide dismutases (Cu, Zn-SODs) showed that the dimeric structure contributed to the high catalytic efficiency and the stability. In this study, steered molecular dynamics simulations were used firstly to study the main interactions between two subunits of Cu, Zn-SODs. The decomposition process study showed that there were not only four pairs of hydrogen bonds but also twenty-five residue pairs participating hydrophobic interactions between A and B chains of SOD, and van der Waals interactions occupied a dominant position among these residue pairs. Moreover, the residue pairs of hydrogen bonds played a major role in maintaining the protein conformation. The analysis of the energy and conformational changes in the SMD simulation showed that there were two groups (two conformational locks) between A and B chains of SOD. The first group consisted of one hydrogen-bond residues pair and seven hydrophobic interactions residues pairs with a total average energy of −30.10 KJ/mol, and the second group of three hydrogen-bond residues pair and eighteen hydrophobic interactions residues pairs formed with a total average energy of −115.23 KJ/mol.

Investigating cyclic peptides inhibiting CD2-CD58 interactions through molecular dynamics and molecular docking methods

The CD2-CD58 protein-protein interaction is known to favor the recognition of antigen presenting cells by T cells. The structural, energetics, and dynamical properties of three known cyclic CD58 ligands, named P6, P7, and RTD-c, are studied through molecular dynamics (MD) simulations and molecular docking calculations. The ligands are built so as to mimic the C and F β-strands of protein CD2, connected via turn inducers. The MD analyses focus on the location of the ligands with respect to the experimental binding site and on the direct and water-mediated hydrogen bonds (H bonds) they form with CD58. Ligand P6, with a sequence close to the experimental β-strands of CD2, presents characteristics that explain its higher experimental affinity, e.g., the lower mobility and flexibility at the CD58 surface, and the larger number and occurrence frequency of ligand-CD58 H bonds. For the two other ligands, the structural modifications lead to changes in the binding pattern with CD58 and its dynamics. In parallel, a large set of molecular docking calculations, carried out with various search spaces and docking algorithms, are compared to provide a consensus view of the preferred ligand binding modes. The analysis of the ligand side chain locations yields results that are consistent with the CD2-CD58 crystal structure and suggests various binding modes of the experimentally identified hot spot of the ligands, i.e., Tyr86. P6 is shown to form a number of contacts that are also present in the experimental CD2-CD58 structure.

Capstan-like mechanism in hyaluronan-phospholipid systems

Functionality of articular cartilage results from complex interactions between its molecular components. Among many biomolecules, two are of prime importance for lubrication: hyaluronic acid (HA) and phospholipids (PL). The purpose of this study is to discuss a mechanism of interaction between these two components and how their synergies contribute to nanobiolubrication of articular cartilage. Preliminary molecular dynamics simulations have been performed to investigate these interactions by adopting a capstan-like mechanism of action. By applying a constant pulling force to both ends of a HA molecule, wrapped around a PL micelle, we viewed the rotation of the PL micelle. The simulations were performed upon two physicochemical constraints: force- and solvent-dependency. The results show the efficiency of rotation from intermolecular bond creation and annihilation. We found a direct relation between the available surface of the micelle and the magnitude of the force, which varies significantly through the unwinding. The movement of the attached molecules is characterized by a slide-to-roll relation, which is affected by the viscosity of the surrounding medium. As a consequence, two solvents were studied for specific force conditions and the molecular dynamics simulation exhibited double the slide-to-roll coefficient for the viscous solvent as compared to its low-viscosity limit.

Simulation of Reversible Protein-Protein Binding and Calculation of Binding Free Energies Using Perturbed Distance Restraints

Virtually all biological processes depend on the interaction between proteins at some point. The correct prediction of biomolecular binding free-energies has many interesting applications in both basic and applied pharmaceutical research. While recent advances in the field of molecular dynamics (MD) simulations have proven the feasibility of the calculation of protein–protein binding free energies, the large conformational freedom of proteins and complex free energy landscapes of binding processes make such calculations a difficult task. Moreover, convergence and reversibility of resulting free-energy values remain poorly described. In this work, an easy-to-use, yet robust approach for the calculation of standard-state protein–protein binding free energies using perturbed distance restraints is described. In the binding process the conformations of the proteins were restrained, as suggested earlier. Two approaches to avoid end-state problems upon release of the conformational restraints were compared. The method was evaluated by practical application to a small model complex of ubiquitin and the very flexible ubiquitin-binding domain of human DNA polymerase ι (UBM2). All computed free energy differences were closely monitored for convergence, and the calculated binding free energies had a mean unsigned deviation of only 1.4 or 2.5 kJ·mol^–1 from experimental values. Statistical error estimates were in the order of thermal noise. We conclude that the presented method has promising potential for broad applicability to quantitatively describe protein–protein and various other kinds of complex formation.

Investigating the effect of different transducer stiffness values on the contactin complex detachment by steered molecular dynamics

This study investigated the adhesion behavior of Contactin4 (CNTN4), a member of Immunoglobulin Super Family (Ig-SF) of cell adhesion molecules. Contactin4 plays a crucial role in the formation, maintenance, and plasticity of neuronal networks. Contactin in its complex configuration with protein tyrosine phosphatase gamma (PTPRG) was selected for simulation. By utilizing Steered Molecular Dynamics (SMD), the uniaxial force was applied to induce unbinding of the complex, and the force-induced detachment of complex components was probed. Three sets of simulations with three values of transducer stiffness and five pulling speeds were designed. Our results showed the dependence of unbinding force on both accessible parameters of pulling speed and spring stiffness. By increasing the stiffness value and pulling speed the rupture force increased. Accordingly, the dissociation rates due to the Bell's theory based on rupture forces and loading rates were calculated.

Theoretical study of enzymatically catalyzed tautomerization of carbon acids in aqueous solution: quantum calculations and steered molecular dynamics simulations

The thermodynamics and kinetics of enzymatically assisted reactions of carbon acids were studied theoretically in this work. Quantum electronic (QE) structure calculations and steered molecular dynamics (SMD) simulations were carried out. Three 3-butenal tautomerization reactions that proceed from the β,γ-unsaturated reactant (R) to the α,β-unsaturated carbon acid product (P) and occur in two elementary steps through an intermediate (I) were studied, ignoring or including the surrounding aqueous medium in the calculations. The Gibbs free energies of activation of the R ⇆ I enolization and I ⇆ P ketonization steps were found to decrease considerably when residues simulating enzymes were introduced into these processes. Although the processes became slightly more favorable thermodynamically when the solution was included in the simulations, they became less favorable kinetically. The results from SMD simulations of these reactions were qualitatively consistent with the values we obtained using QE as well as those found by other authors in similar studies. Our simulations also allowed us to perform a detailed study of these reactions in solution.

Exploring peptide/MHC detachment processes using hierarchical natural move Monte Carlo

MOTIVATION

The binding between a peptide and a major histocompatibility complex (MHC) is one of the most important processes for the induction of an adaptive immune response. Many algorithms have been developed to predict peptide/MHC (pMHC) binding. However, no approach has yet been able to give structural insight into how peptides detach from the MHC.

RESULTS

In this study, we used a combination of coarse graining, hierarchical natural move Monte Carlo and stochastic conformational optimization to explore the detachment processes of 32 different peptides from HLA-A*02:01. We performed 100 independent repeats of each stochastic simulation and found that the presence of experimentally known anchor amino acids affects the detachment trajectories of our peptides. Comparison with experimental binding affinity data indicates the reliability of our approach (area under the receiver operating characteristic curve 0.85). We also compared to a 1000 ns molecular dynamics simulation of a non-binding peptide (AAAKTPVIV) and HLA-A*02:01. Even in this simulation, the longest published for pMHC, the peptide does not fully detach. Our approach is orders of magnitude faster and as such allows us to explore pMHC detachment processes in a way not possible with all-atom molecular dynamics simulations.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available for download at http://www.cs.ox.ac.uk/mosaics/.

CONTACT

bernhard.knapp@stats.ox.ac.uk

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Signaling factor interactions with polysaccharide aggregates of bacterial biofilms

Biofilms are surface-attached colonies of bacteria embedded in an extracellular polymeric substance (EPS). Inside the eukaryotic hosts, bacterial biofilms interact with the host cells through signaling factors (SFs). These signaling processes play important roles in the interaction between bacteria and host cells and the outcome of infections and symbiosis. However, how host immune factors diffuse through biofilms is not well understood. Here, we describe synergistic molecular dynamics and experimental approaches for studying the translocation of signaling factors through polysaccharide chain aggregates present in the extracellular matrix of bacterial biofilms. The effect of polysaccharide chain degradation on the energetics of SF-EPS interactions was examined by simulating an EPS consisting of various polysaccharide chain lengths. It is shown that the SF stabilization energy, defined as the average potential of mean force difference between the environments outside and within the matrix, increases linearly with decreasing chain length. This effect has been explained based on the changes in the polysaccharide configurations around the SF. Specifically, shorter chains are packed tightly around the SF, promoting favorable SF-EPS interactions, while longer chains are packed loosely resulting in screening of interactions with neighboring chains. We further investigated the translocation of SFs through the host cell membrane using molecular dynamics simulations. Further, simulations predict the existence of energy barriers greater than 1000 kJ mol(-1) associated with the translocation of the signaling factors necrosis factor-alpha (TNF-α) and granulocyte macrophage colony stimulating factor (GM-CSF) across the lipid bilayer. The agreement of computational and experimental findings motivates future computational studies using a more detailed description of the EPS aimed at understanding the role of the extracellular matrix on biofilm drug resistance.

A computational model of the glycine tautomerization reaction in aqueous solution

A theoretical study of the internal proton transfer reaction of glycine (Gly) in aqueous solution was performed by means of steered molecular dynamics (SMD) simulation with solute-solvent interaction potentials derived from AMBER van der Waals parameters and QM/MM electrostatic charges in solution. Thermodynamic calculations and the analysis of the solvation structure, dynamic properties, and vibrational spectra associated with the species involved in the tautomerization process were performed. The results obtained for the Gibbs free energy activation and reaction (ΔG(≠) =5.28 kcal mol⁻¹ and ΔG(R)=-6.65 kcal mol⁻¹), the solute-solvent interaction energy of the different glycine structures, and the hydrogen-bond lifetimes are in agreement with previous studies. These hydrations drive an increase in the diffusion coefficient and a decrease in the time of reorientation when the process takes place in the direction Z-Gly → TS-Gly → N-Gly. The vibrational spectrum associated with the normal modes of the bridge hydrogen atoms shows the N-H stretching at ν(s)=3,470 cm⁻¹ and ν(as)=3,470 cm⁻¹, the O-H stretching at 3,205 cm⁻¹, and the NHO bending at about 1,400 cm⁻¹, in agreement with previously reported data.

Amide-imide tautomerism of acetohydroxamic acid in aqueous solution: quantum calculation and SMD simulations

The kinetics and the thermodynamics of the amide-imide tautomerizations of acetohydroxamic acid in aqueous solution are studied from a theoretical point of view.

Antigen-antibody interactions and structural flexibility of a femtomolar-affinity antibody

The femtomolar-affinity mutant antibody (4M5.3) generated by directed evolution is interesting because of the potential of antibody engineering. In this study, the mutant and its wild type (4-4-20) were compared in terms of antigen-antibody interactions and structural flexibility to elucidate the effects of directed evolution. For this purpose, multiple steered molecular dynamics (SMD) simulations were performed. The pulling forces of SMD simulations elucidated the regions that form strong attractive interactions in the binding pocket. Structural analysis in these regions showed two important mutations for improving attractive interactions. First, mutation of Tyr102(H) to Ser (sequence numbering of Protein Data Bank entry 1FLR ) played a role in resolving the steric hindrance on the pathway of the antigen in the binding pocket. Second, mutation of Asp31(H) to His played a role in resolving electrostatic repulsion. Potentials of mean force (PMFs) of both the wild type and the mutant showed landscapes that do not include obvious intermediate states and go directly to the bound state. These landscapes were regarded as funnel-like binding free energy landscapes. Furthermore, the structural flexibility based on the fluctuations of the positions of atoms was analyzed. It was shown that the fluctuations in the positions of the antigen and residues in contact with antigen tend to be smaller in the mutant than in the wild type. This result suggested that structural flexibility decreases as affinity is improved by directed evolution. This suggestion is similar to the relationship between affinity and flexibility for in vivo affinity maturation, which was suggested by Romesberg and co-workers [Jimenez, R., et al. (2003) Proc. Natl. Acad. Sci. U.S.A.100, 92-97]. Consequently, the relationship was found to be applicable up to femotomolar affinity levels.

Transient-state fluctuationlike relation for the driving force on a biomolecule

In experiments and simulations the force acting on a single biomolecular system has been observed as a fluctuating quantity if the system is driven under constant velocity. We ask the question that is analogous to transient state entropy production and work fluctuation relations whether the force fluctuations observed in the single biomolecular system satisfy a transient state fluctuationlike relation, and the answer is in the affirmative. Using a constant velocity pulling steered molecular dynamics simulation study for protein unfolding, we confirm that the force fluctuations of this single biomolecular system satisfy a transient-state fluctuationlike relation 1/γ(T,v) ln[P(v)(+f)/P(v)(-f)] = f. P(v)(±f) is the probability of positive and negative values of forces f = f · for a given unfolding velocity of magnitude v and the pulling direction n, nis the unit vector of n, and γ(T,v) is a factor that depends on initial equilibrium temperature T and the unfolding velocity. For different unfolding velocities we find that the system in the nonequilibrium pulling region displays substantial negative fluctuation in its unfolding force when velocity decreases. A negative value of force may indicate the emergence of refolding behavior during protein unfolding. We also find that γ(T,v) ~ T(-δ)v(α) and the system relaxation time τ(T,v) ~ T(δ)v(-(1+α), where α and δ are scaling exponents.

DISSECTING INDIVIDUAL LIGAND–RECEPTOR BONDS WITH A LAMINAR FLOW CHAMBER

The most important function of proteins may well be to bind to other biomolecules. It has long been felt that kinetic rates of bond formation and dissociation between soluble receptors and ligands might account for most features of the binding process. Only theoretical considerations allowed to predict the behaviour of surface-attached receptors from the properties of soluble forms. During the last decade, experimental progress essentially based on flow chambers, atomic force microscopes or biomembrane force probes allowed direct analysis of biomolecule interaction at the single bond level and gave new insight into previously ignored features such as bond mechanical properties or energy landscapes. The aim of this review is (i) to describe the main advances brought by laminar flow chambers, including information on bond response to forces, multiplicity of binding states, kinetics of bond formation between attached structures, effect of molecular environment on receptor efficiency and behaviour of multivalent attachment, (ii) to compare results obtain by this and other techniques on a few well defined molecular systems, and (iii) to discuss the limitations of the flow chamber method. It is concluded that a new framework may be needed to account for the effective behaviour of biomolecule association.

Insights into Ligand-Protein Binding from Local Mechanical Response

Computational studies of ligand–protein binding are crucial for properly designing novel compounds of potential pharmacological interest. In this respect, researchers are increasingly interested in steered molecular dynamics for ligand–protein binding and unbinding studies. In particular, it has been suggested that analyzing the work profiles along the ligand–protein undocking paths could be fruitful. Here, we propose that small portions of work profiles, termed “local mechanical responses” of the system to a steering force, could serve as a universal measure for capturing relevant information about the system under investigation. Specifically, we first collected a high number of steering trajectories using two biological systems of increasing complexity (i.e., alanine dipeptide and (R)-roscovitine/CDK5 complex). Then, we devised a novel postprocessing tool to be applied to the local mechanical responses, to extract structural information related to the biological processes under investigation. Despite the out-of-equilibrium character of the trajectories, the analysis carried out on the work profiles provided pivotal information about the investigated biological processes. This could eventually be applied to drug design.

Effect of tensile force on the mechanical behavior of actin filaments

Actin filaments are the most abundant components of the cellular cytoskeleton, and play critical roles in various cellular functions such as migration, division and shape control. In these activities, mechanical tension causes structural changes in the double-helical structure of the actin filament, which is a key modulator of cytoskeletal reorganization. This study performed large-scale molecular dynamics (MD) and steered MD simulations to quantitatively analyze the effects of tensile force on the mechanical behavior of actin filaments. The results revealed that when a tensile force of 200pN was applied to a filament consisting of 14 actin subunits, the twist angle of the filament decreased by approximately 20°, corresponding to a rotation of approximately -2° per subunit, representing a critical structural change in actin filaments. Based on these structural changes, the variance in filament length and twist angle was found to decrease, leading to increases in extensional and torsional stiffness. Torsional stiffness increased significantly under the tensile condition, and the ratio of filament stiffness under tensile force to that under no external force increased significantly on longer temporal scales. The results obtained from this study contribute to the understanding of mechano-chemical interactions concerning actin dynamics, showing that increased tensile force in the filament prevents actin regulatory proteins from binding to the filament.

Design rules for biomolecular adhesion: lessons from force measurements

Cell adhesion to matrix, other cells, or pathogens plays a pivotal role in many processes in biomolecular engineering. Early macroscopic methods of quantifying adhesion led to the development of quantitative models of cell adhesion and migration. The more recent use of sensitive probes to quantify the forces that alter or manipulate adhesion proteins has revealed much greater functional diversity than was apparent from population average measurements of cell adhesion. This review highlights theoretical and experimental methods that identified force-dependent molecular properties that are central to the biological activity of adhesion proteins. Experimental and theoretical methods emphasized in this review include the surface force apparatus, atomic force microscopy, and vesicle-based probes. Specific examples given illustrate how these tools have revealed unique properties of adhesion proteins and their structural origins.

Elasticity and rupture of a multi-domain neural cell adhesion molecule complex

The neural cell adhesion molecule (NCAM) plays an important role in nervous system development. NCAM forms a complex between its terminal domains Ig1 and Ig2. When NCAM of cell A and of cell B connect to each other through complexes Ig12(A)/Ig12(B), the relative mobility of cells A and B and membrane tension exerts a force on the Ig12(A)/Ig12(B) complex. In this study, we investigated the response of the complex to force, using steered molecular dynamics. Starting from the structure of the complex from the Ig1-Ig2-Ig3 fragment, we first demonstrated that the complex, which differs in dimensions from a previous structure from the Ig1-Ig2 fragment in the crystal environment, assumes the same extension when equilibrated in solvent. We then showed that, when the Ig12(A)/Ig12(B) complex is pulled apart with forces 30-70 pN, it exhibits elastic behavior (with a spring constant of approximately 0.03 N/m) because of the relative reorientation of domains Ig1 and Ig2. At higher forces, the complex ruptures; i.e., Ig12(A) and Ig12(B) separate. The interfacial interactions between Ig12(A) and Ig12(B), monitored throughout elastic extension and rupture, identify E16, F19, K98, and L175 as key side chains stabilizing the complex.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Analyzing different parameters of steered molecular dynamics for small membrane interacting molecules

The aim of this report is two-fold: First, to show the applicability of the Steered Molecular Dynamics (SMD) methodology for analyzing non-specific interactions governing the membrane affinity process of small biological molecules. Second, to point out a correlation between the system response and certain combinations of the SMD parameters (spring-elastic-constant and pulling-group). For these purposes, a simplified membrane model was used, modeled as a non-polar region limited by two polar aqueous media in a continuous dielectric representation. Polarization-induced effects at both interfaces were taken into account by the "electrostatic images" method. To perform SMD simulations a harmonic external force, representing a spring acting on a selected atom, forces the molecule to "break" its interaction with the surrounding environment by extracting it out of the membrane. With this approach, small molecules and peptides, with known affinity for the membrane environment, were studied: the zwitterionic tryptophan residue and a pentapeptide AcWLKLL. The SMD parameters, spring-elastic-constant and pulling-group, were varied and combined in order to analyze the systems responses in each case. It was observed that, the spring stiffness was crucial to reveal specific events that occur during the molecule behavior; hence, it was directly responsible for the sensitivity of this methodology. The pulling-group selected highly influenced on the reaction pathway, a fact that it was not observed with other parameters; consequently, force profiles are like the "fingerprints" of these induced pathways. The potential profile for the tryptophan was recovered from the SMD simulations being in good agreement with that estimated by an approximation method. With this rather simple model approach, SMD methodology has proven to be suitable for revealing the main interactions that govern the membrane affinity processes of small molecules and peptides.

Nonlinear mechanical modeling of cell adhesion

Cell adhesion, which is mediated by the receptor-ligand bonds, plays an essential role in various biological processes. Previous studies often described the force-extension relationship of receptor-ligand bond with linear assumption. However, the force-extension relationship of the bond is intrinsically nonlinear, which should have significant influence on the mechanical behavior of cell adhesion. In this work, a nonlinear mechanical model for cell adhesion is developed, and the adhesive strength was studied at various bond distributions. We find that the nonlinear mechanical behavior of the receptor-ligand bonds is crucial to the adhesive strength and stability. This nonlinear behavior allows more bonds to achieve large bond force simultaneously, and therefore the adhesive strength becomes less sensitive to the change of bond density at the outmost periphery of the adhesive area. In this way, the strength and stability of cell adhesion are soundly enhanced. The nonlinear model describes the cell detachment behavior better than the linear model.

Cytochrome c(2) Exit Strategy: Dissociation Studies and Evolutionary Implications

Small, water-soluble, type c cytochromes form a transient network connecting major bioenergetic membrane protein complexes in both photosynthesis and respiration. In the photosynthesis cycle of Rhodobacter sphaeroides, cytochrome c2 (cyt c2) docks to the reaction center (RC), undergoes electron transfer, and exits for the cytochrome bc1 complex. Translations of cyt c2 about the RC-cyt c2 docking interface and surrounding membrane reveal possible exit pathways. A pathway at a minimal elevation allowed by the architecture of the RC is analyzed using both an all-atom steered molecular dynamics simulation of the RC-cyt c2 complex and a bioinformatic analysis of the structures and sequences of cyt c. The structure-based phylogenetic analysis allows for the identification of structural elements that have evolved to satisfy the requirements of having multiple functional partners. The patterns of evolutionary variation obtained from the phylogenetic analysis of both docking partners of cyt c2 reveal conservation of key residues involved in the interaction interfaces that would be candidates for further experimental studies. Additionally, using the molecular mechanics Poisson-Boltzmann surface area method we calculate that the binding free energy of reduced cyt c2 to the RC is nearly 6 kcal/mol more favorable than with oxidized cyt c2. The redox-dependent variations lead to changes in structural flexibility, behavior of the interfacial water molecules, and eventually changes in the binding free energy of the complex.

Impact of salt bridges on the equilibrium binding and adhesion of human CD2 and CD58

This study describes quantitative investigations of the impact of single charge mutations on equilibrium binding, kinetics, and the adhesion strength of the CD2-CD58 interaction. Previously steered molecular dynamics simulations guided the selection of the charge mutants investigated, which include the CD2 mutants D31A, K41A, K51A, and K91A. This set includes mutations in which the previous cell aggregation and binding data either agreed or disagreed with the steered molecular dynamics predictions. Surface plasmon resonance measurements quantified the solution binding properties. Adhesion was quantified with the surface force apparatus, which was used previously to study the closely related CD2-CD48 interaction. The results reveal roles that these salt bridges play in equilibrium binding and adhesion. We discuss both the molecular basis of this behavior and its implications for cell adhesion.

Free-energy reconstruction from experiments performed under different biasing programs

Recently developed nonequilibrium statistical physics relationships, including Jarzynski's equality and the Crooks fluctuation theorem, have been used to calculate equilibrium thermodynamic properties using data from both laboratory and computational experiments. Although Jarzynski's derivation does not include an explicit time dependency, prior work utilizing the relationship to reconstruct free-energy surfaces has combined data from experiments performed under identical conditions. Here, a formalism is developed for combining data from a variety of biasing protocols, as in dynamic force spectroscopy experiments. The method is then demonstrated on data from simulations conducted under a wide range of pulling velocities and with a random biasing protocol.

A structure-based sliding-rebinding mechanism for catch bonds

Catch bonds, whose lifetimes are prolonged by force, have been observed in selectin-ligand interactions and other systems. Several biophysical models have been proposed to explain this counterintuitive phenomenon, but none was based on the structure of the interacting molecules and the noncovalent interactions at the binding interface. Here we used molecular dynamics simulations to study changes in structure and atomic-level interactions during forced unbinding of P-selectin from P-selectin glycoprotein ligand-1. A mechanistic model for catch bonds was developed based on these observations. In the model, "catch" results from forced opening of an interdomain hinge that tilts the binding interface to allow two sides of the contact to slide against each other. Sliding promotes formation of new interactions and even rebinding to the original state, thereby slowing dissociation and prolonging bond lifetimes. Properties of this sliding-rebinding mechanism were explored using a pseudoatom representation and Monte Carlo simulations. The model has been supported by its ability to fit experimental data and can be related to previously proposed two-pathway models.

Partitioning of anesthetics into a lipid bilayer and their interaction with membrane-bound peptide bundles

Molecular dynamics simulations have been performed to investigate the partitioning of the volatile anesthetic halothane from an aqueous phase into a coexisting hydrated bilayer, composed of 1,2-dioleoyl-sn-glycero-3-phosphocholine (DOPC) lipids, with embedded alpha-helical peptide bundles based on the membrane-bound portions of the alpha- and delta-subunits, respectively, of nicotinic acetylcholine receptor. In the molecular dynamics simulations halothane molecules spontaneously partitioned into the DOPC bilayer and then preferentially occupied regions close to lipid headgroups. A single halothane molecule was observed to bind to tyrosine (Tyr-277) residue in the alpha-subunit, an experimentally identified specific binding site. The binding of halothane attenuated the local loop dynamics of alpha-subunit and significantly influenced global concerted motions suggesting anesthetic action in modulating protein function. Steered molecular dynamics calculations on a single halothane molecule partitioned into a DOPC lipid bilayer were performed to probe the free energy profile of halothane across the lipid-water interface and rationalize the observed spontaneous partitioning. Partitioned halothane molecules affect the hydrocarbon chains of the DOPC lipid, by lowering of the hydrocarbon tilt angles. The anesthetic molecules also caused a decrease in the number of peptide-lipid contacts. The observed local and global effects of anesthetic binding on protein motions demonstrated in this study may underlie the mechanism of action of anesthetics at a molecular level.

Sequential unfolding of individual helices of bacterioopsin observed in molecular dynamics simulations of extraction from the purple membrane

Multiple molecular dynamics simulations of bacterioopsin pulling from its C-terminus show that its alpha-helices unfold individually. In the first metastable state observed in the simulations, helix G is unfolded at its C-terminal segment while the rest of helix G (residues 200-216) is folded and opposes resistance because of a salt-bridge network consisting of Asp-212 and Lys-216 on helix G and Arg-82 and Asp-85 on helix C. Helix G unfolds inside the bundle because the external force is applied to its C-terminal end in a direction perpendicular to the surface of the membrane. Inversely, helix F has to flip by 180 degrees to exit from the membrane because the applied force and the helical N-C axis point in opposite directions. At the highest peak of the force, which cannot be interpreted in single-molecule force spectroscopy experiments, helix F has a pronounced kink at Pro-186. Mutation of Pro-186 and/or the charged side chains mentioned above, which are involved in very favorable electrostatic interactions in the low-dielectric region of the membrane, are expected to reduce the highest peak of the force. Helices E and D unfold in a similar way to helices G and F, respectively. Hence, the force-distance profile and sequence of events during forced unfolding of bacterioopsin are influenced by the up-and-down topology of the seven-helix bundle. The sequential extraction of individual helices from the membrane suggests that the spontaneous (un)folding of bacterioopsin proceeds through metastable bundles of fewer than seven helices. The metastable states observed in the simulations provide atomic level evidence that corroborates the interpretation of very recent force spectroscopy experiments of bacteriorhodopsin refolding.

Calculating potentials of mean force from steered molecular dynamics simulations

Steered molecular dynamics (SMD) permits efficient investigations of molecular processes by focusing on selected degrees of freedom. We explain how one can, in the framework of SMD, employ Jarzynski's equality (also known as the nonequilibrium work relation) to calculate potentials of mean force (PMF). We outline the theory that serves this purpose and connects nonequilibrium processes (such as SMD simulations) with equilibrium properties (such as the PMF). We review the derivation of Jarzynski's equality, generalize it to isobaric--isothermal processes, and discuss its implications in relation to the second law of thermodynamics and computer simulations. In the relevant regime of steering by means of stiff springs, we demonstrate that the work on the system is Gaussian-distributed regardless of the speed of the process simulated. In this case, the cumulant expansion of Jarzynski's equality can be safely terminated at second order. We illustrate the PMF calculation method for an exemplary simulation and demonstrate the Gaussian nature of the resulting work distribution.

Molecular mechanisms of cellular mechanics

Mechanical forces play an essential role in cellular processes as input, output, and signals. Various protein complexes in the cell are designed to handle, transform and use such forces. For instance, proteins of muscle and the extracellular matrix can withstand considerable stretching forces, hearing-related and mechanosensory proteins can transform weak mechanical stimuli into electrical signals, and regulatory proteins are suited to forcing DNA into loops to control gene expression. Here we review the structure-function relationship of four protein complexes with well defined and representative mechanical functions. The first example is titin, a protein that confers passive elasticity on muscle. The second system is the elastic extracellular matrix protein, fibronectin, and its cellular receptor integrin. The third protein system is the transduction apparatus in hearing and other mechanical senses, likely containing cadherin and ankyrin repeats. The last system is the lac repressor protein, which regulates gene expression by looping DNA. This review focuses on atomic level descriptions of the physical mechanisms underlying the various mechanical functions of the stated proteins.

Mechanical Unfolding of Segment-Swapped Protein G Dimer: Results from Replica Exchange Molecular Dynamics Simulations

The protein G dimer (pdb code 1Q10) is a mutated dimeric form of the immunoglobulin-binding domain B1 of streptococcal protein G, in which the two monomeric units have swapped elements of their secondary structure. We have used replica exchange molecular dynamics simulations to study how this dimer responds to a mechanical force that pulls the N-terminus of one unit and the C-terminus of the other apart. We have further compared the mechanical response of the dimer to that of the protein G monomer. When the pulling force is low enough, the mechanical unfolding can be viewed as a thermally activated barrier crossing process. For each protein, we have computed the corresponding free energy barrier and its dependence on the pulling force. While the dimer is found to be less resistant to mechanical unfolding than its monomeric counterpart, the two proteins exhibit essentially the same mechanical unfolding mechanism involving separation of the terminal parallel strands. On the basis of our results, we speculate that the mechanical properties of natural adhesives, composites, fibers, and other materials may be optimized not only at a single molecule level but also at the mesoscopic level through the interactions among individual chains.

Molecular simulation of the binding of nerve growth factor peptide mimics to the receptor tyrosine kinase A

Nerve growth factor (NGF) mimics play an important role for therapies that target the receptor tyrosine kinase A (trkA). The N-terminal fragment of the NGF (N-term@NGF) was previously demonstrated to be an important determinant for affinity and specificity in the binding to trkA. Here we use a variety of computational tools (contact surface analysis and free energy predictions) to identify residues playing a key role for the binding to the receptor. Molecular dynamics simulations are then used to investigate the stability of complexes between trkA and peptides mimicking N-term@NGF. Steered molecular dynamics calculations are finally performed to investigate the process of detaching the peptide from the receptor. Three disruptive events are observed, the first involving the breaking of all intermolecular interactions except two salt bridges, which break subsequently.

Force meets chemistry: Analysis of mechanochemical conversion in focal adhesions using fluorescence recovery after photobleaching

Mechanotransduction--the process by which mechanical forces are converted into changes of intracellular biochemistry--is critical for normal cell and tissue function. Integrins facilitate mechanochemical conversion by transferring physical forces from the extracellular matrix, across the cell surface, and to cytoskeletal-associated proteins within focal adhesions. It is likely that force alters biochemistry at these sites by altering molecular binding affinities of a subset of focal adhesion proteins, but this has been difficult to quantify within living cells. Here, we describe how the fluorescence recovery after photobleaching (FRAP) technique can be adapted and used in conjunction with mathematical models to directly measure force-dependent alterations in molecular binding and unbinding rate constants of individual focal adhesion proteins in situ. We review these recent findings, and discuss the strengths and limitations of this approach for analysis of mechanochemical signaling in focal adhesions and other cellular structures. The ability to quantify molecular binding rate constants in the physical context of the living cytoplasm should provide new insight into the molecular basis of cellular mechanotransduction. It also may facilitate future efforts to bridge biological experimentation and mathematical modeling in our quest for a systems biology level description of cell regulation.

Change of the unbinding mechanism upon a mutation: a molecular dynamics study of an antibody-hapten complex

We study forced unbinding of fluorescein from the wild type (WT) and a mutant [H(H58)A] of the single-chain variable-fragment (scFv) anti-fluorescein antibody FITC-E2 by molecular dynamics simulations using various pulling techniques. A large number of long simulations were needed to obtain statistically meaningful results as both the wild type and the H(H58)A mutant unbinding occurs through multiple pathways, often with metastable intermediates. For the wild type, the rate-limiting step in the unbinding process corresponds to the breaking of the non-native interactions characteristic of a specific intermediate. The H(H58)A mutation disfavors the occurrence of this intermediate. Two events where the hapten partially unbinds in the absence of pulling force are observed in extensive equilibrium simulations of the wild type, and their analysis indicates that forced unbinding and spontaneous unbinding proceed along similar pathways. The different unbinding mechanisms observed in the simulations suggest a possible reason for the difference in the experimental off-rate between the two antibodies. We predict mutations that are expected to modulate the occurrence of the unbinding intermediate. For two such new mutants [H(H58)A and S(H52)A], our predictions are validated in silico by additional simulations. The accompanying paper in this issue by Honegger et al. reports the X-ray structure of FITC-E2 with a derivative of fluorescein, which was used as the starting conformation for the work presented here.

Scalable molecular dynamics with NAMD

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

Nanomechanics of adhesion proteins

Recent advances in molecular force measurements have resulted in the quantification of the nanomechanical properties of single molecular bonds, and elucidated novel relationships between molecular architecture and biomolecular adhesion. The measured forces to rupture single intermolecular bonds revealed novel and unexpected ways that proteins respond to mechanical force. Measurement of the magnitude of interprotein forces and the distances over which they act further determined how protein architecture may contribute to both the stability and structural organization of adhesive junctions.