Single-Molecule Experiments in Vitro and in Silico

In-Silico Analyses of Molecular Force Sensors for Mechanical Characterization of Biological Systems

Mechanical forces play key roles in biological processes such as cell migration and sensory perception. In recent years molecular force sensors have been developed as tools for in situ force measurements. Here we use all-atom steered molecular dynamics simulations to predict and study the relationship between design parameters and mechanical properties for three types of molecular force sensors commonly used in cellular biological research: two peptide- and one DNA-based. The peptide-based sensors consist of a pair of fluorescent proteins, which can undergo Förster resonance energy transfer (FRET), linked by spider silk (GPGGA)n or synthetic (GGSGGS)n disordered regions. The DNA-based sensor consists of two fluorophore-labeled strands of DNA that can be unzipped or sheared upon force application with a FRET signal as readout of dissociation. We simulated nine sensors, three of each kind. After equilibration, flexible peptide linkers of three different lengths were stretched by applying forces to their N- and C-terminal Cα atoms in opposite directions. Similarly, we equilibrated a DNA-based sensor and pulled on the phosphate atom of the terminal guanine of one strand and a selected phosphate atom on the other strand in the opposite direction. These simulations were performed at constant velocity (0.01 nm/ns - 10 nm/ns) and constant force (10 pN - 500 pN) for all versions of the sensors. Our results show how the force response of these sensors depends on their length, sequence, configuration and loading rate. Mechanistic insights gained from simulations analyses indicate that interpretation of experimental results should consider the influence of transient formation of secondary structure in peptide-based sensors and of overstretching in DNA-based sensors. These predictions can guide optimal fluorophore choice and facilitate the rational design of new sensors for use in protein, DNA, hybrid systems, and molecular devices.

Elasticity and Thermal Stability are Key Determinants of Hearing Rescue by Mini-Protocadherin-15 Proteins

Protocadherin-15 is a core protein component of inner-ear hair-cell tip links pulling on transduction channels essential for hearing and balance. Protocadherin-15 defects can result in non-syndromic deafness or Usher syndrome type 1F (USH1F) with hearing loss, balance deficits, and progressive blindness. Three rationally engineered shortened versions of protocadherin-15 (mini-PCDH15s) amenable for gene therapy have been used to rescue function in USH1F mouse models. Two can successfully or partially rescue hearing, while another one fails. Here we show that despite varying levels of hearing rescue, all three mini-PCDH15 versions can rescue hair-cell mechanotransduction. Negative-stain electron microscopy shows that all three versions form dimers like the wild-type protein, while crystal structures of some engineered fragments show that these can properly fold and bind calcium ions essential for function. In contrast, simulations predict distinct elasticities and nano differential scanning fluorimetry shows differences in melting temperature measurements. Our data suggest that elasticity and thermal stability are key determinants of sustained hearing rescue by mini-PCDH15s.

Mechanical Stability and Unfolding Pathways of Parallel Tetrameric G-Quadruplexes Probed by Pulling Simulations

Guanine quadruplex (GQ) is a noncanonical nucleic acid structure formed by guanine-rich DNA and RNA sequences. Folding of GQs is a complex process, where several aspects remain elusive, despite being important for understanding structure formation and biological functions of GQs. Pulling experiments are a common tool for acquiring insights into the folding landscape of GQs. Herein, we applied a computational pulling strategy─steered molecular dynamics (SMD) simulations─in combination with standard molecular dynamics (MD) simulations to explore the unfolding landscapes of tetrameric parallel GQs. We identified anisotropic properties of elastic conformational changes, unfolding transitions, and GQ mechanical stabilities. Using a special set of structural parameters, we found that the vertical component of pulling force (perpendicular to the average G-quartet plane) plays a significant role in disrupting GQ structures and weakening their mechanical stabilities. We demonstrated that the magnitude of the vertical force component depends on the pulling anchor positions and the number of G-quartets. Typical unfolding transitions for tetrameric parallel GQs involve base unzipping, opening of the G-stem, strand slippage, and rotation to cross-like structures. The unzipping was detected as the first and dominant unfolding event, and it usually started at the 3'-end. Furthermore, results from both SMD and standard MD simulations indicate that partial spiral conformations serve as a transient ensemble during the (un)folding of GQs.

Single-molecule analysis of osmolyte-mediated nanomechanical unfolding behavior of a protein domain

The small organic molecules, known as osmolytes being ubiquitously present in different cell types, affect protein folding, stability and aggregation. However, it is unknown how the osmolytes affect the nanomechanical unfolding behavior of protein domain. Here, we show the osmolyte-dependent mechanical unfolding properties of protein titin immunoglobulin-27 (I27) domain using an atomic force microscopy (AFM)-based single-molecule force spectroscopy. We found that amines and methylamines improved the mechanical stability of I27 domain, whereas polyols had no effect. Interestingly, glycine betaine (GB) or trimethylamine-N-oxide (TMAO) increased the average unfolding force of the protein domain. The kinetic parameters analyzed at single-molecule level reveal that stabilizing effect of osmolytes is due to a decrease in the unfolding rate constant of I27, which was confirmed by molecular dynamics simulations. Our study reveals different effects that diverse osmolytes have on the mechanical properties of the protein, and suggests the potential use of osmolytes in modulating the mechanical stability of proteins required for various nano-biotechnological applications.

Assessing mitochondria-targeted acyl hydroquinones on the mitochondrial platelet function and cytotoxic activity: Role of the linker length

INTRODUCTION

The use of triphenylphosphonium cation (TPP+) linked to phenolic compounds by alkyl chains has a significant relevance as a mitochondrial delivery strategy in biomedicine because it affects mitochondrial bioenergetics in models of noncommunicable diseases such as cancer and cardiovascular-related conditions. Studies indicate that a long alkyl chain (10-12 carbon) increases the mitochondrial accumulation of TPP+-linked drugs. In contrast, other studies show that these compounds are consistently toxic to micromolar concentrations (as observed in platelets). In the present study, we evaluated the in vitro effect of three series of triphenylphosphonium-linked acyl hydroquinones derivates on the metabolism and function of human platelets using 3-9 carbons for the alkyl linker. Those were assessed to determine the role of the length of the alkyl chain linker on platelet toxicity.

METHODS

Human platelets were exposed in vitro to different concentrations (2-40 μM) of every compound; cellular viability, phosphatidylserine exposition, mitochondrial membrane potential (ΔΨm), intracellular calcium release, and intracellular ROS generation were assessed by flow cytometry. An in silico energetic profile was generated with Umbrella sampling molecular dynamics (MD).

RESULTS AND CONCLUSIONS

There was an increase in cytotoxic activity directly related to the length of the acyl chain and lipophilicity, as seen by three techniques, and this was consistent with a decrease in ΔΨm. The in silico energetic profiles point out that the permeability of the mitochondrial membrane may be involved in the cytotoxicity of phosphonium salts. This information may be relevant for the design of new TPP+ -based drugs with a safe cardiovascular profile.

Two molecule force spectroscopy on ligand-receptor interactions

Many biological processes involve the rupture of multiple ligand-receptors or multivalent ligand-receptors. It is challenging to study the rupture of such parallelly arranged multiple ligand-receptors due to the difficulties in engineering such systems in a well-controlled fashion. Here we report the use of two-molecule force spectroscopy to investigate the rupture of two parallelly arranged monomeric streptavidin (mSA)-biotin complexes. By using SpyCatcher-SpyTag chemistry, we successfully engineered a molecular twin of biotin, in which two biotins are arranged in parallel. By reacting mSA with twin biotin, we constructed parallelly arranged two mSA-biotin complexes for force spectroscopy experiments. The incorporation of single molecule fingerprint domains into our mSA-biotin dimers allowed us to identify and assign the rupture events of the parallelly arranged mSA-biotin complexes without any ambiguity in the two-molecule force spectroscopy experiments. Our results revealed that the rupture force of the parallel dimer mSA-biotin is 172 pN at a pulling speed of 400 nm s-1, which is about 1.6 times of that of single mSA-biotin (105 pN). Furthermore, our findings indicate that the two mSA-biotin behave as non-interacting, independent ligand-receptors. The strategy we demonstrated here can be extended to other ligand-receptors and may open up an avenue toward rigorously testing the theoretic predictions proposed in various models regarding the rupture of multiple parallel ligand-receptors.

Complexity of Guanine Quadruplex Unfolding Pathways Revealed by Atomistic Pulling Simulations

Guanine quadruplexes (GQs) are non-canonical nucleic acid structures involved in many biological processes. GQs formed in single-stranded regions often need to be unwound by cellular machinery, so their mechanochemical properties are important. Here, we performed steered molecular dynamics simulations of human telomeric GQs to study their unfolding. We examined four pulling regimes, including a very slow setup with pulling velocity and force load accessible to high-speed atomic force microscopy. We identified multiple factors affecting the unfolding mechanism, i.e.,: (i) the more the direction of force was perpendicular to the GQ channel axis (determined by GQ topology), the more the base unzipping mechanism happened, (ii) the more parallel the direction of force was, GQ opening and cross-like GQs were more likely to occur, (iii) strand slippage mechanism was possible for GQs with an all-anti pattern in a strand, and (iv) slower pulling velocity led to richer structural dynamics with sampling of more intermediates and partial refolding events. We also identified that a GQ may eventually unfold after a force drop under forces smaller than those that the GQ withstood before the drop. Finally, we found out that different unfolding intermediates could have very similar chain end-to-end distances, which reveals some limitations of structural interpretations of single-molecule spectroscopic data.

Big Data e Inteligência Artificial para pesquisa translacional na Covid-19: revisão rápida

RESUMO O objetivo deste estudo foi identificar como a Inteligência Artificial (IA) vem sendo utilizada para a pesquisa translacional no contexto da Covid-19. Foi realizada uma revisão rápida para identificar o uso de técnicas de IA na translação de tecnologias para o enfrentamento da Covid-19. Empregou-se estratégia de busca com base em termos MeSH e seus respectivos sinônimos em sete bases de dados. Dos 59 artigos identificados, oito foram incluídos. Foram identificadas 11 experiências que usaram IA para a pesquisa translacional em Covid-19: predição de eficácia medicamentosa; predição de patogenicidade do Sars-CoV-2; diagnóstico de imagem para Covid-19; predição de incidência de Covid-19; estimativas de impacto da Covid-19 na sociedade; automatização de sanitização de ambientes hospitalares e clínicos; rastreio de pessoas infectadas e possivelmente infectadas; monitoramento do uso de máscaras; predição de gravidade de pacientes; estratificação de risco do paciente; e predição de recursos hospitalares. A pesquisa translacional pode ajudar no desenvolvimento produtivo e industrial em saúde, especialmente quando apoiada em métodos de IA, uma ferramenta cada vez mais importante, sobretudo quando se discute a Quarta Revolução Industrial e suas aplicações na saúde.

Methyl Esters of (E)-Cinnamic Acid: Activity against the Plant-Parasitic Nematode Meloidogyne incognita and In Silico Interaction with Histone Deacetylase

(E)-Cinnamaldehyde is very active against Meloidogyne incognita but has low persistence in soil. To circumvent this problem, esters of cinnamic acid were evaluated as a substitute for (E)-cinnamaldehyde. The best results under assays with M. incognita second-stage juveniles (J2) were obtained for the methyl esters of (E)-p-fluoro- (13), (E)-p-chloro- (14), and (E)-p-bromocinnamic acid (15), which showed lethal concentrations to 50% (LC₅₀) J2 of 168, 95, and 216 μg/mL, respectively. Under the same conditions, the LC₅₀ values for the nematicides carbofuran and fluensulfone were 160 and 34 μg/mL, respectively. Substances 13-15 were also active against nematode eggs, which account for most of the M. incognita population in the field. According to an in silico study, substances 13-15 can act against the nematode through inhibition of histone deacetylase. Therefore, esters 13-15 and histone deacetylase are potentially useful for the rational design of new nematicides for the control of M. incognita.

Elastic versus brittle mechanical responses predicted for dimeric cadherin complexes

Cadherins are a superfamily of adhesion proteins involved in a variety of biological processes that include the formation of intercellular contacts, the maintenance of tissue integrity, and the development of neuronal circuits. These transmembrane proteins are characterized by ectodomains composed of a variable number of extracellular cadherin (EC) repeats that are similar but not identical in sequence and fold. E-cadherin, along with desmoglein and desmocollin proteins, are three classical-type cadherins that have slightly curved ectodomains and engage in homophilic and heterophilic interactions through an exchange of conserved tryptophan residues in their N-terminal EC1 repeat. In contrast, clustered protocadherins are straighter than classical cadherins and interact through an antiparallel homophilic binding interface that involves overlapped EC1 to EC4 repeats. Here we present molecular dynamics simulations that model the adhesive domains of these cadherins using available crystal structures, with systems encompassing up to 2.8 million atoms. Simulations of complete classical cadherin ectodomain dimers predict a two-phased elastic response to force in which these complexes first softly unbend and then stiffen to unbind without unfolding. Simulated α, β, and γ clustered protocadherin homodimers lack a two-phased elastic response, are brittle and stiffer than classical cadherins and exhibit complex unbinding pathways that in some cases involve transient intermediates. We propose that these distinct mechanical responses are important for function, with classical cadherin ectodomains acting as molecular shock absorbers and with stiffer clustered protocadherin ectodomains facilitating overlap that favors binding specificity over mechanical resilience. Overall, our simulations provide insights into the molecular mechanics of single cadherin dimers relevant in the formation of cellular junctions essential for tissue function.

Collective mechanical responses of cadherin-based adhesive junctions as predicted by simulations

Cadherin-based adherens junctions and desmosomes help stabilize cell-cell contacts with additional function in mechano-signaling, while clustered protocadherin junctions are responsible for directing neuronal circuits assembly. Structural models for adherens junctions formed by epithelial cadherin (CDH1) proteins indicate that their long, curved ectodomains arrange to form a periodic, two-dimensional lattice stabilized by tip-to-tip trans interactions (across junction) and lateral cis contacts. Less is known about the exact architecture of desmosomes, but desmoglein (DSG) and desmocollin (DSC) cadherin proteins are also thought to form ordered junctions. In contrast, clustered protocadherin (PCDH)-based cell-cell contacts in neuronal tissues are thought to be responsible for self-recognition and avoidance, and structural models for clustered PCDH junctions show a linear arrangement in which their long and straight ectodomains form antiparallel overlapped trans complexes. Here, we report all-atom molecular dynamics simulations testing the mechanics of minimalistic adhesive junctions formed by CDH1, DSG2 coupled to DSC1, and PCDHγB4, with systems encompassing up to 3.7 million atoms. Simulations generally predict a favored shearing pathway for the adherens junction model and a two-phased elastic response to tensile forces for the adhesive adherens junction and the desmosome models. Complexes within these junctions first unbend at low tensile force and then become stiff to unbind without unfolding. However, cis interactions in both the CDH1 and DSG2-DSC1 systems dictate varied mechanical responses of individual dimers within the junctions. Conversely, the clustered protocadherin PCDHγB4 junction lacks a distinct two-phased elastic response. Instead, applied tensile force strains trans interactions directly, as there is little unbending of monomers within the junction. Transient intermediates, influenced by new cis interactions, are observed after the main rupture event. We suggest that these collective, complex mechanical responses mediated by cis contacts facilitate distinct functions in robust cell-cell adhesion for classical cadherins and in self-avoidance signaling for clustered PCDHs.

Assessing the Stability of Biological Fibrils by Molecular-Scale Simulations

The nanomechanical characterization of several biological fibrils that are the result of protein aggregation via molecular dynamics simulation is nowadays feasible, and together with atomic force microscopy experiments has widened our understanding of the forces in the regime of pN-nN and system sizes of about hundreds of nanometers. Several methodologies have been developed to achieve this target, and they range from the atomistic representation via molecular force fields to coarse-grained strategies that provide comparable results with experiments in a systematic way. In this chapter, we discuss several methodologies for the calculation of mechanical parameters, such as the elastic constants of relevant biological systems. They are presented together with details about parameterization and current limitations. Then, we discuss some of the applications of such methodologies for the description of bacterial filament and β-amyloid systems. Finally, the latest lines of development are discussed.

CHARMM-GUI Drude prepper for molecular dynamics simulation using the classical Drude polarizable force field

Explicit treatment of electronic polarizability in empirical force fields (FFs) represents an extension over a traditional additive or pairwise FF and provides a more realistic model of the variations in electronic structure in condensed phase, macromolecular simulations. To facilitate utilization of the polarizable FF based on the classical Drude oscillator model, Drude Prepper has been developed in CHARMM-GUI. Drude Prepper ingests additive CHARMM protein structures file (PSF) and pre-equilibrated coordinates in CHARMM, PDB, or NAMD format, from which the molecular components of the system are identified. These include all residues and patches connecting those residues along with water, ions, and other solute molecules. This information is then used to construct the Drude FF-based PSF using molecular generation capabilities in CHARMM, followed by minimization and equilibration. In addition, inputs are generated for molecular dynamics (MD) simulations using CHARMM, GROMACS, NAMD, and OpenMM. Validation of the Drude Prepper protocol and inputs is performed through conversion and MD simulations of various heterogeneous systems that include proteins, nucleic acids, lipids, polysaccharides, and atomic ions using the aforementioned simulation packages. Stable simulations are obtained in all studied systems, including 5 μs simulation of ubiquitin, verifying the integrity of the generated Drude PSFs. In addition, the ability of the Drude FF to model variations in electronic structure is shown through dipole moment analysis in selected systems. The capabilities and availability of Drude Prepper in CHARMM-GUI is anticipated to greatly facilitate the application of the Drude FF to a range of condensed phase, macromolecular systems.

Heterophilic and homophilic cadherin interactions in intestinal intermicrovillar links are species dependent

Enterocytes are specialized epithelial cells lining the luminal surface of the small intestine that build densely packed arrays of microvilli known as brush borders. These microvilli drive nutrient absorption and are arranged in a hexagonal pattern maintained by intermicrovillar links formed by 2 nonclassical members of the cadherin superfamily of calcium-dependent cell adhesion proteins: protocadherin-24 (PCDH24, also known as CDHR2) and the mucin-like protocadherin (CDHR5). The extracellular domains of these proteins are involved in heterophilic and homophilic interactions important for intermicrovillar function, yet the structural determinants of these interactions remain unresolved. Here, we present X-ray crystal structures of the PCDH24 and CDHR5 extracellular tips and analyze their species-specific features relevant for adhesive interactions. In parallel, we use binding assays to identify the PCDH24 and CDHR5 domains involved in both heterophilic and homophilic adhesion for human and mouse proteins. Our results suggest that homophilic and heterophilic interactions involving PCDH24 and CDHR5 are species dependent with unique and distinct minimal adhesive units.

Application of in vitro PAMPA technique and in silico computational methods for blood-brain barrier permeability prediction of novel CNS drug candidates

Permeability assessment of small molecules through the blood-brain barrier (BBB) plays a significant role in the development of effective central nervous system (CNS) drug candidates. Since in vivo methods for BBB permeability estimation require a lot of time and resources, in silico and in vitro approaches are becoming increasingly popular nowadays for faster and more economical predictions in early phases of drug discovery. In this work, through application of in vitro parallel artificial membrane permeability assay (PAMPA-BBB) and in silico computational methods we aimed to examine the passive permeability of eighteen compounds, which affect serotonin and dopamine levels in the CNS. The data set was consisted of novel six human dopamine transporter (hDAT) substrates that were previously identified as the most promising lead compounds for further optimisation to achieve neuroprotective effect, twelve approved CNS drugs, and their related compounds. Firstly, PAMPA methods was used to experimentally determine effective BBB permeability (P_e) for all studied compounds and obtained results were further submitted for quantitative structure permeability relationship (QSPR) analysis. QSPR models were built by using three different statistical methods: stepwise multiple linear regression (MLR), partial least square (PLS), and support-vector machine (SVM), while their predictive capability was tested through internal and external validation. Obtained statistical parameters (MLR- R²_pred=-0.10; PLS- R²_pred=0.64, r²_m=0.69, r^/2_m=0.44; SVM- R²_pred=0.57, r²_m=0.72, r^/2_m=0.55) indicated that the SVM model is superior over others. The most important molecular descriptors (H₀p and SolvEMt_3D) were identified and used to propose structural modifications of the examined compounds in order to improve their BBB permeability. Moreover, steered molecular dynamics (SMD) simulation was employed to comprehensively investigate the permeability pathway of compounds through a lipid bilayer. Taken together, the created QSPR model could be used as a reliable and fast pre-screening tool for BBB permeability prediction of structurally related CNS compounds, while performed MD simulations provide a good foundation for future in silico examination.

Calculation of CYP450 protein-ligand binding and dissociation free energy paths

The function of an enzyme depends on its dynamic structure, and the catalytic mechanism has long been an active focus of research. The principle for interpreting protein selectivity and fidelity stems from optimization of the active site upon protein-substrate complexation, i.e., a lock-and-key configuration, on which most protein-substrate molecule binding recognition, and hence drug discovery, relies. Yet another thought has been to incorporate the protein folding interior tunnels for stereo- and regio-selectivity along the protein-substrate or protein-ligand/inhibitor binding process. Free energy calculations provide valuable information for molecular recognition and protein-ligand binding dynamics and kinetics. In this study, we focused on the kinetics of cytochrome P450 proteins (CYP450s) and the protein interior tunnel structure-dynamics relationship in terms of the substrate binding and leaving mechanism. A case in point is given by the prostaglandin H₂ (PGH₂) homologous isomerase of prostacyclin synthase. To calculate the reactant and product traversing the tunnels to and from the heme site, the free energy paths and tunnel potentials of mean force are constructed from steered molecular dynamics simulations and adaptive basing force umbrella sampling simulations. We explore the binding tunnels and critical residue lining characteristics for the ligand traverse and the underlying mechanism of CYP450 activity. Our theoretical analysis provides insights into the decisive role of the substrate tunnel binding process of the CYP450 mechanism and may be useful in drug design and protein engineering contexts.

Understanding the Mechanochemistry of Ladder-Type Cyclobutane Mechanophores by Single Molecule Force Spectroscopy

We have recently reported a series of ladder-type cyclobutane mechanophores, polymers of which can transform from nonconjugated structures to conjugated structures and change many properties at once. These multicyclic mechanophores, namely, exo-ladderane/ene, endo-benzoladderene, and exo-bicyclohexene-peri-naphthalene, have different ring structures fused to the first cyclobutane, significantly different free energy changes for ring-opening, and different stereochemistry. To better understand their mechanochemistry, we used single molecule force spectroscopy (SMFS) to characterize their force-extension behavior and measure the threshold forces. The threshold forces correlate with the activation energy of the first bond, but not with the strain of the fused rings distal to the polymer main chain, suggesting that the activation of these ladder-type mechanophores occurs with similar early transition states, which is supported by force-modified potential energy surface calculations. We further determined the stereochemistry of the mechanically generated dienes and observed significant and variable contour length elongation for these mechanophores both experimentally and computationally. The fundamental understanding of ladder-type mechanophores will facilitate future design of multicyclic mechanophores with amplified force-response and their applications as mechanically responsive materials.

Validity of point-mass model in off-resonance dynamic atomic force microscopy

The quantitative measurement of viscoelasticity of nano-scale entities is an important goal of nanotechnology research and there is considerable progress with advent of dynamic atomic force microscopy. The hydrodynamics of cantilever, the force sensor in AFM measurements, plays a pivotal role in quantitative estimates of nano-scale viscoelasticity. The point-mass (PM) model, wherein the AFM cantilever is approximated as a point-mass with mass-less spring is widely used in dynamic AFM analysis and its validity, particularly in liquid environments, is debated. It is suggested that the cantilever must be treated as a continuous rectangular beam to obtain accurate estimates of nano-scale viscoelasticity of materials it is probing. Here, we derived equations, which relate stiffness and damping coefficient of the material under investigation to measured parameters, by approximating cantilever as a point-mass and also considering the full geometric details. These equations are derived for both tip-excited as well as base-excited cantilevers. We have performed off-resonance dynamic atomic force spectroscopy on a single protein molecule to investigate the validity of widely used PM model. We performed measurements with AFMs equipped with different cantilever excitation methods as well as detection schemes to measure cantilever response. The data was analyzed using both, continuous beam model and the PM model. We found that both models yield same results when the experiments are performed in truly off-resonance regime with small amplitudes and the cantilever stiffness is much higher than the interaction stiffness. Our findings suggest that a simple PM approximation based model is adequate to describe the dynamics, provided care is taken while performing experiments so that the approximations used in these models are valid.

Reconstruction of mechanical unfolding and refolding pathways of proteins with atomic force spectroscopy and computer simulations

Most proteins in proteomes are large, typically consist of more than one domain and are structurally complex. This often makes studying their mechanical unfolding pathways challenging. Proteins composed of tandem repeat domains are a subgroup of multi-domain proteins that, when stretched, display a saw-tooth pattern in their mechanical unfolding force extension profiles due to their repetitive structure. However, the assignment of force peaks to specific repeats undergoing mechanical unraveling is complicated because all repeats are similar and they interact with their neighbors and form a contiguous tertiary structure. Here, we describe in detail a combination of experimental and computational single-molecule force spectroscopy methods that proved useful for examining the mechanical unfolding and refolding pathways of ankyrin repeat proteins. Specifically, we explain and delineate the use of atomic force microscope-based single molecule force spectroscopy (SMFS) to record the mechanical unfolding behavior of ankyrin repeat proteins and capture their unusually strong refolding propensity that is responsible for generating impressive refolding force peaks. We also describe Coarse Grain Steered Molecular Dynamic (CG-SMD) simulations which complement the experimental observations and provide insights in understanding the unfolding and refolding of these proteins. In addition, we advocate the use of novel coiled-coils-based mechanical polypeptide probes which we developed to demonstrate the vectorial character of folding and refolding of these repeat proteins. The combination of AFM-based SMFS on native and CC-equipped proteins with CG-SMD simulations is powerful not only for ankyrin repeat polypeptides, but also for other repeat proteins and more generally to various multidomain, non-repetitive proteins with complex topologies.

Force probe simulations using an adaptive resolution scheme

Molecular simulations of the forced unfolding and refolding of biomolecules or molecular complexes allow to gain important kinetic, structural and thermodynamic information about the folding process and the underlying energy landscape. In force probe molecular dynamics (FPMD) simulations, one pulls one end of the molecule with a constant velocity in order to induce the relevant conformational transitions. Since the extended configuration of the system has to fit into the simulation box together with the solvent such simulations are very time consuming. Here, we apply a hybrid scheme in which the solute is treated with atomistic resolution and the solvent molecules far away from the solute are described in a coarse-grained manner. We use the adaptive resolution scheme (AdResS) that has very successfully been applied to various examples of equilibrium simulations. We perform FPMD simulations using AdResS on a well studied system, a dimer formed from mechanically interlocked calixarene capsules. The results of the multiscale simulations are compared to all-atom simulations of the identical system and we observe that the size of the region in which atomistic resolution is required depends on the pulling velocity, i.e. the particular non-equilibrium situation. For large pulling velocities a larger all atom region is required. Our results show that multiscale simulations can be applied also in the strong non-equilibrium situations that the system experiences in FPMD simulations.

Reconstruction of ARNT PAS-B Unfolding Pathways by Steered Molecular Dynamics and Artificial Neural Networks

Several experimental studies indicated that large conformational changes, including partial domain unfolding, have a role in the functional mechanisms of the basic helix loop helix Per/ARNT/SIM (bHLH-PAS) transcription factors. Recently, single-molecule atomic force microscopy (AFM) revealed two distinct pathways for the mechanical unfolding of the ARNT PAS-B. In this work we used steered molecular dynamics simulations to gain new insights into this process at an atomistic level. To reconstruct and classify pathways sampled in multiple simulations, we designed an original approach based on the use of self-organizing maps (SOMs). This led us to identify two types of unfolding pathways for the ARNT PAS-B, which are in good agreement with the AFM findings. Analysis of average forces mapped on the SOM revealed a stable conformation of the PAS-B along one pathway, which represents a possible structural model for the intermediate state detected by AFM. The approach here proposed will facilitate the study of other signal transmission mechanisms involving the folding/unfolding of PAS domains.

Atomic Force Microscopy-Based Force Spectroscopy and Multiparametric Imaging of Biomolecular and Cellular Systems

During the last three decades, a series of key technological improvements turned atomic force microscopy (AFM) into a nanoscopic laboratory to directly observe and chemically characterize molecular and cell biological systems under physiological conditions. Here, we review key technological improvements that have established AFM as an analytical tool to observe and quantify native biological systems from the micro- to the nanoscale. Native biological systems include living tissues, cells, and cellular components such as single or complexed proteins, nucleic acids, lipids, or sugars. We showcase the procedures to customize nanoscopic chemical laboratories by functionalizing AFM tips and outline the advantages and limitations in applying different AFM modes to chemically image, sense, and manipulate biosystems at (sub)nanometer spatial and millisecond temporal resolution. We further discuss theoretical approaches to extract the kinetic and thermodynamic parameters of specific biomolecular interactions detected by AFM for single bonds and extend the discussion to multiple bonds. Finally, we highlight the potential of combining AFM with optical microscopy and spectroscopy to address the full complexity of biological systems and to tackle fundamental challenges in life sciences.

Structural insight into toxin secretion by contact-dependent growth inhibition transporters

Bacterial contact-dependent growth inhibition (CDI) systems use a type Vb secretion mechanism to export large CdiA toxins across the outer membrane by dedicated outer membrane transporters called CdiB. Here, we report the first crystal structures of two CdiB transporters from Acinetobacter baumannii and Escherichia coli. CdiB transporters adopt a TpsB fold, containing a 16-stranded transmembrane β-barrel connected to two periplasmic domains. The lumen of the CdiB pore is occluded by an N-terminal α-helix and the conserved extracellular loop 6; these two elements adopt different conformations in the structures. We identified a conserved DxxG motif located on strand β1 that connects loop 6 through different networks of interactions. Structural modifications of DxxG induce rearrangement of extracellular loops and alter interactions with the N-terminal α-helix, preparing the system for α-helix ejection. Using structural biology, functional assays, and molecular dynamics simulations, we show how the barrel pore is primed for CdiA toxin secretion.

[1-(4-chloro-3-nitrobenzenesulfonyl)-1H-indol-3-yl]-methanol, an indole-3-carbinol derivative, inhibits glutamate release in rat cerebrocortical nerve terminals by suppressing the P/Q-type Ca2+ channels and Ca2+/calmodulin/protein kinase A pathway

Indole-3-carbinol (I3C), found in cruciferous vegetables, has been proposed to exhibit neuroprotective effects. This study aimed to investigate the effect of the I3C derivative [1(4-chloro-3-nitrobenzenesulfonyl)-1H-indol-3-yl]-methanol (CIM), which has superior pharmacokinetic properties to I3C, on glutamate release in rat cerebrocortical nerve terminals (synaptosomes). We observed that CIM dose-dependently inhibited glutamate release evoked by the potassium channel blocker 4-aminopyridine (4-AP). CIM-mediated inhibition of glutamate release was attributed to reduced exocytosis, as it correlated with the removal of extracellular calcium and blocking of the vesicular glutamate transporter but not the glutamate transporter. In addition, CIM decreased 4-AP-evoked intrasynaptosomal Ca²⁺ elevation; however, it did not alter the synaptosomal membrane potential. The inhibition of P/Q-typeCa²⁺ channels abolished the effect of CIM on 4-AP-evoked glutamate release, and the effect was not prevented by intracellular Ca²⁺ release inhibitors. Moreover, the molecular docking study showed that CIM exhibited the highest binding affinity with the P/Q-type Ca²⁺channels. Finally, the CIM-mediated inhibition of glutamate release was sensitive to calmodulin, adenylate cyclase (AC), and protein kinase A (PKA) inhibitors. Based on these results, we propose that CIM, through the direct suppression of P/Q-type Ca²⁺ channels, decreases Ca²⁺ influx and the activation of Ca²⁺/calmodulin/AC/PKA signaling, thereby inhibiting glutamate release. This finding is crucial for understanding the role of CIM in the central nervous system and for exploiting its potential in therapeutic interventions.

Rational design of a new mutant of tobacco etch virus protease in order to increase the in vitro solubility

Background and purpose:

Tobacco etch virus (TEV) protease is a protease with high sequence specificity which is useful for the cleavage of fusion proteins. A major limitation of this enzyme is its relatively poor solubility. This study aimed to investigate the effects of some suggested mutations by online tools and molecular dynamics simulation to improve the solubility of TEV protease in vitro.

Experimental approach:

We designed a rational multi-stage process to determine the solubilizing mutations of TEV protease. At the first stage, all the possible mutations were predicted using online tools such as PoPMuSiC and Eris servers, in which five mutations include N23F, N23L, Q74L, Q74V, and Q74I were suggested for further studies. In the next step, the three dimensional structure of the wild type (WT) and the best mutations were subjected to molecular dynamic simulations to evaluate the dynamic behaviour of the obtained structures. The selected mutation was introduced into the structure using site-directed mutagenesis and expressed in Escherichia coli BL21DE3. After purification, solubility and activity of the purified mutant and WT-TEV proteases were assayed.

Findings /Results:

By considering the analysis of various factors such as structural and solubility properties, one mutant, N23F, was selected for in vitro studies which led to a 1.5 times increase in the solubility compared to the WT while its activity was decreased somewhat.

Conclusion and implications:

We propose N23F mutation, according to computational and experimental analyses for TEV proteases which resulted in a 150% increase in solubility compared to the WT.

In situ formation of H-bonding imidazole chains in break-junction experiments

As a small molecule possessing both strong H-bond donor and acceptor functions, 1H-imidazole can participate in extensive homo- or heteromolecular H-bonding networks. These properties are important in Nature, as imidazole moieties are incorporated in many biologically-relevant compounds. Imidazole also finds applications ranging from corrosion inhibition to fire retardants and photography. We have found a peculiar behaviour of imidazole during scanning tunnelling microscopy-break junction (STM-BJ) experiments, in which oligomeric chains connect the two electrodes and allow efficient charge transport. We attributed this behaviour to the formation of hydrogen-bonding networks, as no evidence of such behaviour was found in 1-methylimidazole (incapable of participating in intramolecular hydrogen bonding). The results are supported by DFT calculations, which confirmed our hypothesis. These findings pave the road to the use of hydrogen-bonding networks for the fabrication of dynamic junctions based on supramolecular interactions.

Probing mechanical properties and failure mechanisms of fibrils of self-assembling peptides

Self-assembling peptides (SAPs) are a promising class of biomaterials amenable to easy molecular design and functionalization. Despite their increasing usage in regenerative medicine, a detailed analysis of their biomechanics at the nanoscale level is still missing. In this work, we propose and validate, in all-atom dynamics, a coarse-grained model to elucidate strain distribution, failure mechanisms and biomechanical effects of functionalization of two SAPs when subjected to both axial stretching and bending forces. We highlight different failure mechanisms for fibril seeds and fibrils, as well as the negligible contribution of the chosen functional motif to the overall system rupture. This approach could lay the basis for the development of "more" coarse-grained models in the long pathway connecting SAP sequences and hydrogel mechanical properties.

Mechanical Unfolding of Spectrin Repeats Induces Water-Molecule Ordering

Mechanical processes are involved at many stages of the development of living cells, and often external forces applied to a biomolecule result in its unfolding. Although our knowledge of the unfolding mechanisms and the magnitude of the forces involved has evolved, the role that water molecules play in the mechanical unfolding of biomolecules has not yet been fully elucidated. To this end, we investigated with steered molecular dynamics simulations the mechanical unfolding of dystrophin's spectrin repeat 1 and related the changes in the protein's structure to the ordering of the surrounding water molecules. Our results indicate that upon mechanically induced unfolding of the protein, the solvent molecules become more ordered and increase their average number of hydrogen bonds. In addition, the unfolded structures originating from mechanical pulling expose an increasing amount of the hydrophobic residues to the solvent molecules, and the uncoiled regions adapt a convex surface with a small radius of curvature. As a result, the solvent molecules reorganize around the protein's small protrusions in structurally ordered waters that are characteristic of the so-called "small-molecule regime," which allows water to maintain a high hydrogen bond count at the expense of an increased structural order. We also determined that the response of water to structural changes in the protein is localized to the specific regions of the protein that undergo unfolding. These results indicate that water plays an important role in the mechanically induced unfolding of biomolecules. Our findings may prove relevant to the ever-growing interest in understanding macromolecular crowding in living cells and their effects on protein folding, and suggest that the hydration layer may be exploited as a means for short-range allosteric communication.

Analogies in 3D molecular visualisations: development of a cell biology animation 'How cells move - a new interpretation of old data'

Cell biology and imaging technology have vastly improved over the past decades, enabling scientists to dissect the inner workings of a cell. In addition to technical limits on spatial and temporal resolution, which obscure the picture at the molecular level, the sheer density and complexity of information impede clear understanding. 3D molecular visualisation has therefore blossomed as a way to translate molecular data in a more tangible form. Whilst the molecular machinery involved in cell locomotion has been extensively studied, existing narratives describing how cells generate the forces that drive movement remain unclear. Polymerisation of a protein called actin is clearly essential. The general belief in the cell migration field is that actin polymerisation's main role is to push the leading edge of the cell forwards, while the rest of the cell follows passively. The cell migration & chemotaxis group at the CRUK Beatson Institute propose an alternative hypothesis, in which actin filaments constitute cables. Motor proteins pull on these cables, causing them to behave like the treads of a tank and drive cell movement. This article describes the development of a 3D animation that uses analogical reasoning to contrast the 'tank' hypothesis for cell locomotion with the current dogma.

Cytotoxicity effects and ionic diffusion of single-wall carbon nanotubes in cell membrane

While carbon nanotubes have been put into massive practical industrial, environmental and biomedicine applications, the cytotoxicity effects or the effect to the ionic channels they bring into the living cells need to be thoroughly investigated. In this work, molecular dynamic simulations have been carried out to investigate the ionic diffusion through the single wall armchair carbon nanotube embedded right inside the cell membrane. By modeling a two-membrane system, we build a virtual cytoplasm environment including a cell chamber and an extracellular space, in which a certain amount of solute is dissolved. The system is first brought to its equilibrium by deployment of minimization and then simulated. The results suggested that carbon nanotubes (CNTs) with size less than (12, 12) shall be less cytotoxic since it does not bring any ionic diffusion through the CNT channel, so as to maintain active cytoplasm environment. Another phenomenon we observed is a notable shifting angle of the carbon nanotube which was normal to the surface of cell membrane initially. In general, the inclination angle of the carbon nanotube increases with its radius.

Understanding enhanced mechanical stability of DNA in the presence of intercalated anticancer drug: Implications for DNA associated processes

Most of the anticancer drugs bind to double-stranded DNA (dsDNA) by intercalative-binding mode. Although experimental studies have become available recently, a molecular-level understanding of the interactions between the drug and dsDNA that lead to the stability of the intercalated drug is lacking. Of particular interest are the modifications of the mechanical properties of dsDNA observed in experiments. The latter could affect many biological functions, such as DNA transcription and replication. Here, we probe, via all-atom molecular dynamics (MD) simulations, the change in the mechanical properties of intercalated drug-DNA complexes for two intercalators, daunomycin and ethidium. We find that, upon drug intercalation, the stretch modulus of DNA increases significantly, whereas its persistence length and bending modulus decrease. Steered MD simulations reveal that it requires higher forces to stretch the intercalated dsDNA complexes than the normal dsDNA. Adopting various pulling protocols to study force-induced DNA melting, we find that the dissociation of dsDNA becomes difficult in the presence of intercalators. The results obtained here provide a plausible mechanism of function of the anticancer drugs, i.e., via altering the mechanical properties of DNA. We also discuss long-time consequences of using these drugs, which require further in vivo investigations.

High-speed force spectroscopy: microsecond force measurements using ultrashort cantilevers

Complete understanding of the role of mechanical forces in biological processes requires knowledge of the mechanical properties of individual proteins and living cells. Moreover, the dynamic response of biological systems at the nano- and microscales span over several orders of magnitude in time, from sub-microseconds to several minutes. Thus, access to force measurements over a wide range of length and time scales is required. High-speed atomic force microscopy (HS-AFM) using ultrashort cantilevers has emerged as a tool to study the dynamics of biomolecules and cells at video rates. The adaptation of HS-AFM to perform high-speed force spectroscopy (HS-FS) allows probing protein unfolding and receptor/ligand unbinding up to the velocity of molecular dynamics (MD) simulations with sub-microsecond time resolution. Moreover, application of HS-FS on living cells allows probing the viscoelastic response at short time scales providing deep understanding of cytoskeleton dynamics. In this mini-review, we assess the principles and recent developments and applications of HS-FS using ultrashort cantilevers to probe molecular and cellular mechanics.

The Work of Titin Protein Folding as a Major Driver in Muscle Contraction

Single-molecule atomic force microscopy and magnetic tweezers experiments have demonstrated that titin immunoglobulin (Ig) domains are capable of folding against a pulling force, generating mechanical work that exceeds that produced by a myosin motor. We hypothesize that upon muscle activation, formation of actomyosin cross bridges reduces the force on titin, causing entropic recoil of the titin polymer and triggering the folding of the titin Ig domains. In the physiological force range of 4-15 pN under which titin operates in muscle, the folding contraction of a single Ig domain can generate 200% of the work of entropic recoil and occurs at forces that exceed the maximum stalling force of single myosin motors. Thus, titin operates like a mechanical battery, storing elastic energy efficiently by unfolding Ig domains and delivering the charge back by folding when the motors are activated during a contraction. We advance the hypothesis that titin folding and myosin activation act as inextricable partners during muscle contraction.

Steered molecular dynamics simulations reveal the role of Ca2+ in regulating mechanostability of cellulose-binding proteins

The conversion of cellulosic biomass into biofuels requires degradation of the biomass into fermentable sugars. The most efficient natural cellulase system for carrying out this conversion is an extracellular multi-enzymatic complex named the cellulosome. In addition to temperature and pH stability, mechanical stability is important for functioning of cellulosome domains, and experimental techniques such as Single Molecule Force Spectroscopy (SMFS) have been used to measure the mechanical strength of several cellulosomal proteins. Molecular dynamics computer simulations provide complementary atomic-resolution quantitative maps of domain mechanical stability for identification of experimental leads for protein stabilization. In this study, we used multi-scale steered molecular dynamics computer simulations, benchmarked against new SMFS measurements, to measure the intermolecular contacts that confer high mechanical stability to a family 3 Carbohydrate Binding Module protein (CBM3) derived from the archetypal Clostridium thermocellum cellulosome. Our data predicts that electrostatic interactions in the calcium binding pocket modulate the mechanostability of the cellulose-binding module, which provides an additional design rule for the rational re-engineering of designer cellulosomes for biotechnology. Our data offers new molecular insights into the origins of mechanostability in cellulose binding domains and gives leads for synthesis of more robust cellulose-binding protein modules. On the other hand, simulations predict that insertion of a flexible strand can promote alternative unfolding pathways and dramatically reduce the mechanostability of the carbohydrate binding module, which gives routes to rational design of tailormade fingerprint complexes for force spectroscopy experiments.

Combined molecular/continuum modeling reveals the role of friction during fast unfolding of coiled-coil proteins

Coiled-coils are filamentous proteins that form the basic building block of important force-bearing cellular elements, such as intermediate filaments and myosin motors. In addition to their biological importance, coiled-coil proteins are increasingly used in new biomaterials including fibers, nanotubes, or hydrogels. Coiled-coils undergo a structural transition from an α-helical coil to an unfolded state upon extension, which allows them to sustain large strains and is critical for their biological function. By performing equilibrium and out-of-equilibrium all-atom molecular dynamics (MD) simulations of coiled-coils in explicit solvent, we show that two-state models based on Kramers' or Bell's theories fail to predict the rate of unfolding at high pulling rates. We further show that an atomistically informed continuum rod model accounting for phase transformations and for the hydrodynamic interactions with the solvent can reconcile two-state models with our MD results. Our results show that frictional forces, usually neglected in theories of fibrous protein unfolding, reduce the thermodynamic force acting on the interface, and thus control the dynamics of unfolding at different pulling rates. Our results may help interpret MD simulations at high pulling rates, and could be pertinent to cytoskeletal networks or protein-based artificial materials subjected to shocks or blasts.

Multivalent Binding of a Ligand-Coated Particle: Role of Shape, Size, and Ligand Heterogeneity

We utilize a multiscale modeling framework to study the effect of shape, size, and ligand composition on the efficacy of binding of a ligand-coated particle to a substrate functionalized with the target receptors. First, we show how molecular dynamics along with steered molecular dynamics calculations can be used to accurately parameterize the molecular-binding free energy and the effective spring constant for a receptor-ligand pair. We demonstrate this for two ligands that bind to the α₅β₁-domain of integrin. Next, we show how these effective potentials can be used to build computational models at the meso- and continuum-scales. These models incorporate the molecular nature of the receptor-ligand interactions and yet provide an inexpensive route to study the multivalent interaction of receptors and ligands through the construction of Bell potentials customized to the molecular identities. We quantify the binding efficacy of the ligand-coated-particle in terms of its multivalency, binding free-energy landscape, and the losses in the configurational entropies. We show that 1) the binding avidity for particle sizes less than 350 nm is set by the competition between the enthalpic and entropic contributions, whereas that for sizes above 350 nm is dominated by the enthalpy of binding; 2) anisotropic particles display higher levels of multivalent binding compared to those of spherical particles; and 3) variations in ligand composition can alter binding avidity without altering the average multivalency. The methods and results presented here have wide applications in the rational design of functionalized carriers and also in understanding cell adhesion.

Enhancement of protein mechanical stability: Correlated deformations are handcuffed by ligand binding

As revealed by previous experiments, protein mechanical stability can be effectively regulated by ligand binding with the binding site distant from the force-bearing region. However, the mechanism for such long-range allosteric control of protein mechanics is still largely unknown. In this work, we use protein topology-based elastic network model (ENM) and all-atomic steered molecular dynamics (SMD) simulations to study the impact of ligand binding on protein mechanical stability in two systems, i.e., GB1 and CheY-binding P2-domain of CheA (CBDCheA). Both ENM and SMD results show that the ligand binding has considerable and negligible effects on the mechanical stability of these two proteins, respectively. These results are consistent with the experimental observations. A physical mechanism for the enhancement of protein mechanical stability was then proposed: the correlated deformations of the force-bearing region and the binding site are handcuffed by the binding of ligand. The handcuff effect suppresses the propagation of internal force in the force-bearing region, thus improving the resistance to the loading force. Our study indicates that ENM method can effectively identify the structure motifs allosterically related to the deformation in the force bearing region, as well as the force propagation pathway within the structure of the studied proteins. Hence, it should be helpful to understand the molecular origin of the different mechanical properties in response to ligand binding for GB1 and CBDCheA.

Generalization of the elastic network model for the study of large conformational changes in biomolecules

The elastic network (EN) is a prime model that describes the long-time dynamics of biomolecules. However, the use of harmonic potentials renders this model insufficient for studying large conformational changes of proteins (e.g. stretching of proteins, folding and thermal unfolding). Here, we extend the capabilities of the EN model by using a harmonic approximation described by Lennard-Jones (LJ) interactions for far contacts and native contacts obtained from the standard overlap criterion as in the case of Gō-like models. While our model is validated against the EN model by reproducing the equilibrium properties for a number of proteins, we also show that the model is suitable for the study of large conformation changes by providing various examples. In particular, this is illustrated on the basis of pulling simulations that predict with high accuracy the experimental data on the rupture force of the studied proteins. Furthermore, in the case of DDFLN4 protein, our pulling simulations highlight the advantages of our model with respect to Gō-like approaches, where the latter fail to reproduce previous results obtained by all-atom simulations that predict an additional characteristic peak for this protein. In addition, folding simulations of small peptides yield different folding times for α-helix and β-hairpin, in agreement with experiment, in this way providing further opportunities for the application of our model in studying large conformational changes of proteins. In contrast to the EN model, our model is suitable for both normal mode analysis and molecular dynamics simulation. We anticipate that the proposed model will find applications in a broad range of problems in biology, including, among others, protein folding and thermal unfolding.

TPP riboswitch aptamer: Role of Mg2+ ions, ligand unbinding, and allostery

Riboswitches are non-coding RNAs that regulate gene expression in response to the binding of metabolites. Their abundance in bacteria makes them ideal drug targets. The prokaryotic thiamine pyrophosphate (TPP) riboswitch regulates gene expression in a wide range of bacteria by undergoing conformational changes in response to the binding of TPP. Although an experimental structure for the aptamer domain of the riboswitch is now available, details of the conformational changes that occur during the binding of the ligand, and the factors that govern these conformational changes, are still not clear. This study employs microsecond-scale molecular dynamics simulations to provide insights into the functioning of the riboswitch aptamer in atomistic detail. A mechanism for the transmission of conformational changes from the ligand-binding site to the P1 switch helix is proposed. Mg²⁺ ions in the binding site play a critical role in anchoring the ligand to the riboswitch. Finally, modeling the egress of TPP from the binding site reveals a two-step mechanism for TPP unbinding. Findings from this study can motivate the design of future studies aimed at modulating the activity of this drug target.

A Bayesian Nonparametric Approach to Single Molecule Förster Resonance Energy Transfer

We develop a Bayesian nonparametric framework to analyze single molecule FRET (smFRET) data. This framework, a variation on infinite hidden Markov models, goes beyond traditional hidden Markov analysis, which already treats photon shot noise, in three critical ways: (1) it learns the number of molecular states present in a smFRET time trace (a hallmark of nonparametric approaches), (2) it accounts, simultaneously and self-consistently, for photophysical features of donor and acceptor fluorophores (blinking kinetics, spectral cross-talk, detector quantum efficiency), and (3) it treats background photons. Point 2 is essential in reducing the tendency of nonparametric approaches to overinterpret noisy single molecule time traces and so to estimate states and transition kinetics robust to photophysical artifacts. As a result, with the proposed framework, we obtain accurate estimates of single molecule properties even when the supplied traces are excessively noisy, subject to photoartifacts, and of short duration. We validate our method using synthetic data sets and demonstrate its applicability to real data sets from single molecule experiments on Holliday junctions labeled with conventional fluorescent dyes.