Binding pathway of retinal to bacterio-opsin: a prediction by molecular dynamics simulations

Formation of bacteriorhodopsin (bR) from apoprotein and retinal has been studied experimentally, but the actual pathway, including the point of entry, is little understood. Molecular dynamics simulations provide a surprisingly clear prediction. A window between bR helices E and F in the transmembrane part of the protein can be identified as an entry point for retinal. Steered molecular dynamics, performed by applying a series of external forces in the range of 200–1000 pN over a period of 0.2 ns to retinal, allows one to extract this chromophore from bR once the Schiff base bond to Lys216 is cleaved. Extraction proceeds until the retinal tail forms a hydrogen bond network with Ala144, Met145, and Ser183 side groups lining the exit/entry window. The manipulation induces a distortion with a fitted root mean square deviation of coordinates (ignoring retinal, water, and hydrogen atoms) of less than 1.9 A by the time the retinal carbonyl reaches the protein surface. The forces needed to extract retinal are due to friction and do not indicate significant potential barriers. The simulations therefore suggest a pathway for the binding of retinal. Water molecules are found to play a crucial role in the binding process

MoMA-LigPath: A web server to simulate protein-ligand unbinding

Protein-ligand interactions taking place far away from the active site, during ligand binding or release, may determine molecular specificity and activity. However, obtaining information about these interactions with experimental or computational methods remains difficult. The computational tool presented in this paper, MoMA-LigPath, is based on a mechanistic representation of the molecular system, considering partial flexibility, and on the application of a robotics-inspired algorithm to explore the conformational space. Such a purely geometric approach, together with the efficiency of the exploration algorithm, enables the simulation of ligand unbinding within very short computing time. Ligand unbinding pathways generated by MoMA-LigPath are a first approximation that can provide very useful information about protein-ligand interactions. When needed, this approximation can be subsequently refined and analyzed using state-of-the-art energy models and molecular modeling methods. MoMA-LigPath is available at http://moma.laas.fr. The web server is free and open to all users, with no login requirement

DNA unzipped under a constant force exhibits multiple metastable intermediates

Single molecule studies, at constant force, of the separation of double-stranded DNA into two separated single strands may provide information relevant to the dynamics of DNA replication. At constant applied force, theory predicts that the unzipped length as a function of time is characterized by jumps during which the strands separate rapidly, followed by long pauses where the number of separated base pairs remains constant. Here, we report previously uncharacterized observations of this striking behavior carried out on a number of identical single molecules simultaneously. When several single lphage molecules are subject to the same applied force, the pause positions are reproducible in each. This reproducibility shows that the positions and durations of the pauses in unzipping provide a sequence-dependent molecular fingerprint. For small forces, the DNA remains in a partially unzipped state for at least several hours. For larger forces, the separation is still characterized by jumps and pauses, but the double-stranded DNA will completely unzip in less than 30 min

Single-molecule pulling: phenomenology and interpretation

Single-molecule pulling techniques have emerged as versatile tools for probing the noncovalent forces holding together the secondary and tertiary structure of macromolecules. They also constitute a way to study at the single-molecule level processes that are familiar from our macroscopic thermodynamic experience. In this Chapter, we summarize the essential phenomenology that is typically observed during single-molecule pulling, provide a general statistical mechanical framework for the interpretation of the equilibrium force spectroscopy and illustrate how to simulate single-molecule pulling experiments using molecular dynamics.Comment: arXiv admin note: text overlap with arXiv:0908.220

Mapping the energy landscape of biomolecules using single molecule force correlation spectroscopy (FCS): Theory and applications

In the current AFM experiments the distribution of unfolding times, P(t), is measured by applying a constant stretching force f_s from which the apparent unfolding rate is obtained. To describe the complexity of the underlying energy landscape requires additional probes that can incorporate the dynamics of tension propagation and relaxation of the polypeptide chain upon force quench. We introduce a theory of force correlation spectroscopy (FCS) to map the parameters of the energy landscape of proteins. In the FCS the joint distribution, P(T,t) of folding and unfolding times is constructed by repeated application of cycles of stretching at constant fs, separated by release periods T during which the force is quenched to f_q<f_s. During the release period, the protein can collapse to a manifold of compact states or refold. We show that P(T,t) can be used to resolve the kinetics of unfolding as well as formation of native contacts and to extract the parameters of the energy landscape using chain extension as the reaction coordinate and P(T,t). We illustrate the utility of the proposed formalism by analyzing simulations of unfolding-refolding trajectories of a coarse-grained protein S1 with beta-sheet architecture for several values of f_s, T and f_q=0. The simulations of stretch-relax trajectories are used to map many of the parameters that characterize the energy landscape of S1.Comment: 23 pages, 9 figures; accepted to Biophysical Journa

Atomic detail visualization of photosynthetic membranes with GPU-accelerated ray tracing

The cellular process responsible for providing energy for most life on Earth, namely, photosynthetic light-harvesting, requires the cooperation of hundreds of proteins across an organelle, involving length and time scales spanning several orders of magnitude over quantum and classical regimes. Simulation and visualization of this fundamental energy conversion process pose many unique methodological and computational challenges. We present, in two accompanying movies, light-harvesting in the photosynthetic apparatus found in purple bacteria, the so-called chromatophore. The movies are the culmination of three decades of modeling efforts, featuring the collaboration of theoretical, experimental, and computational scientists. We describe the techniques that were used to build, simulate, analyze, and visualize the structures shown in the movies, and we highlight cases where scientific needs spurred the development of new parallel algorithms that efficiently harness GPU accelerators and petascale computers

Dominant Folding Pathways of a WW Domain

We investigate the folding mechanism of the WW domain Fip35 using a realistic atomistic force field by applying the Dominant Reaction Pathways (DRP) approach. We find evidence for the existence of two folding pathways, which differ by the order of formation of the two hairpins. This result is consistent with the analysis of the experimental data on the folding kinetics of WW domains and with the results obtained from large-scale molecular dynamics (MD) simulations of this system. Free-energy calculations performed in two coarse-grained models support the robustness of our results and suggest that the qualitative structure of the dominant paths are mostly shaped by the native interactions. Computing a folding trajectory in atomistic detail only required about one hour on 48 CPU's. The gain in computational efficiency opens the door to a systematic investigation of the folding pathways of a large number of globular proteins

Actomyosin Interaction: Mechanical and Energetic Properties in Different Nucleotide Binding States

The mechanics of the actomyosin interaction is central in muscle contraction and intracellular trafficking. A better understanding of the events occurring in the actomyosin complex requires the examination of all nucleotide-dependent states and of the energetic features associated with the dynamics of the cross-bridge cycle. The aim of the present study is to estimate the interaction strength between myosin in nucleotide-free, ATP, ADP·Pi and ADP states and actin monomer. The molecular models of the complexes were constructed based on cryo-electron microscopy maps and the interaction properties were estimated by means of a molecular dynamics approach, which simulate the unbinding of the complex applying a virtual spring to the core of myosin protein. Our results suggest that during an ATP hydrolysis cycle the affinity of myosin for actin is modulated by the presence and nature of the nucleotide in the active site of the myosin motor domain. When performing unbinding simulations with a pulling rate of 0.001 nm/ps, the maximum pulling force applied to the myosin during the experiment is about 1nN. Under these conditions the interaction force between myosin and actin monomer decreases from 0.83 nN in the nucleotide-free state to 0.27 nN in the ATP state, and increases to 0.60 nN after ATP hydrolysis and Pi release from the complex (ADP state)

Molecular Basis of Ligand Dissociation in β-Adrenergic Receptors

The important and diverse biological functions of β-adrenergic receptors (βARs) have promoted the search for compounds to stimulate or inhibit their activity. In this regard, unraveling the molecular basis of ligand binding/unbinding events is essential to understand the pharmacological properties of these G protein-coupled receptors. In this study, we use the steered molecular dynamics simulation method to describe, in atomic detail, the unbinding process of two inverse agonists, which have been recently co-crystallized with β1 and β2ARs subtypes, along four different channels. Our results indicate that this type of compounds likely accesses the orthosteric binding site of βARs from the extracellular water environment. Importantly, reconstruction of forces and energies from the simulations of the dissociation process suggests, for the first time, the presence of secondary binding sites located in the extracellular loops 2 and 3 and transmembrane helix 7, where ligands are transiently retained by electrostatic and Van der Waals interactions. Comparison of the residues that form these new transient allosteric binding sites in both βARs subtypes reveals the importance of non-conserved electrostatic interactions as well as conserved aromatic contacts in the early steps of the binding process

Analyzing the forces binding a restriction endonuclease to DNA using a synthetic nanopore

Restriction endonucleases are used prevalently in recombinant DNA technology because they bind so stably to a specific target sequence and, in the presence of cofactors, cleave double-helical DNA specifically at a target sequence at a high rate. Using synthetic nanopores along with molecular dynamics (MD), we have analyzed with atomic resolution how a prototypical restriction endonuclease, EcoRI, binds to the DNA target sequence—GAATTC—in the absence of a Mg2+ ion cofactor. We have previously shown that there is a voltage threshold for permeation of DNA bound to restriction enzymes through a nanopore that is associated with a nanonewton force required to rupture the complex. By introducing mutations in the DNA, we now show that this threshold depends on the recognition sequence and scales linearly with the dissociation energy, independent of the pore geometry. To predict the effect of mutation in a base pair on the free energy of dissociation, MD is used to qualitatively rank the stability of bonds in the EcoRI–DNA complex. We find that the second base in the target sequence exhibits the strongest binding to the protein, followed by the third and first bases, with even the flanking sequence affecting the binding, corroborating our experiments