Real Time Observation of Ultrafast Peptide Conformational Dynamics: Molecular Dynamics Simulation vs Infrared Experiment

Photoinduced reconfiguration to control the protein-binding affinity of azobenzene-cyclized peptides

The impact of next-generation biorecognition elements (ligands) will be determined by the ability to remotely control their binding activity for a target biomolecule in complex environments. Compared to conventional mechanisms for regulating binding affinity (pH, ionic strength, or chaotropic agents), light provides higher accuracy and rapidity, and is particularly suited for labile targets. In this study, we demonstrate a general method to develop azobenzene-cyclized peptide ligands with light-controlled affinity for target proteins. Light triggers a cis/trans isomerization of the azobenzene, which results in a major structural rearrangement of the cyclic peptide from a non-binding to a binding configuration. Critical to this goal are the ability to achieve efficient photo-isomerization under low light dosage and the temporal stability of both cis and trans isomers. We demonstrated our method by designing photo-switchable peptides targeting vascular cell adhesion marker 1 (VCAM1), a cell marker implicated in stem cell function. Starting from a known VCAM1-binding linear peptide, an ensemble of azobenzene-cyclized variants with selective light-controlled binding were identified by combining in silico design with experimental characterization via spectroscopy and surface plasmon resonance. Variant cycloAZOB[G-VHAKQHRN-K] featured rapid, light-controlled binding of VCAM1 (KD,trans/KD,cis ∼ 130). Biotin-cycloAZOB[G-VHAKQHRN-K] was utilized to label brain microvascular endothelial cells (BMECs), showing co-localization with anti-VCAM1 antibodies in cis configuration and negligible binding in trans configuration.

Peptides and pseudopeptide ligands: a powerful toolbox for the affinity purification of current and next-generation biotherapeutics

Following the consolidation of therapeutic proteins in the fight against cancer, autoimmune, and neurodegenerative diseases, recent advancements in biochemistry and biotechnology have introduced a host of next-generation biotherapeutics, such as CRISPR-Cas nucleases, stem and car-T cells, and viral vectors for gene therapy. With these drugs entering the clinical pipeline, a new challenge lies ahead: how to manufacture large quantities of high-purity biotherapeutics that meet the growing demand by clinics and biotech companies worldwide. The protein ligands employed by the industry are inadequate to confront this challenge: while featuring high binding affinity and selectivity, these ligands require laborious engineering and expensive manufacturing, are prone to biochemical degradation, and pose safety concerns related to their bacterial origin. Peptides and pseudopeptides make excellent candidates to form a new cohort of ligands for the purification of next-generation biotherapeutics. Peptide-based ligands feature excellent target biorecognition, low or no toxicity and immunogenicity, and can be manufactured affordably at large scale. This work presents a comprehensive and systematic review of the literature on peptide-based ligands and their use in the affinity purification of established and upcoming biological drugs. A comparative analysis is first presented on peptide engineering principles, the development of ligands targeting different biomolecular targets, and the promises and challenges connected to the industrial implementation of peptide ligands. The reviewed literature is organized in (i) conventional (α-)peptides targeting antibodies and other therapeutic proteins, gene therapy products, and therapeutic cells; (ii) cyclic peptides and pseudo-peptides for protein purification and capture of viral and bacterial pathogens; and (iii) the forefront of peptide mimetics, such as β-/γ-peptides, peptoids, foldamers, and stimuli-responsive peptides for advanced processing of biologics.

Conformational Changes of Trialanine in Water Induced by Vibrational Relaxation of the Amide I Mode

Most of the protein-based diseases are caused by anomalies in the functionality and stability of these molecules. Experimental and theoretical studies of the conformational dynamics of proteins are becoming in this respect essential to understand the origin of these anomalies. However, a description of the conformational dynamics of proteins based on mechano-energetic principles still remains elusive because of the intrinsic high flexibility of the peptide chains, the participation of weak noncovalent interactions, and the role of the ubiquitous water solvent. In this work, the conformational dynamics of trialanine dissolved in water (D2O) is investigated through Molecular Dynamics (MD) simulations combined with instantaneous normal modes (INMs) analysis both at equilibrium and after the vibrational excitation of the C-terminal amide I mode. The conformational equilibrium between α and pPII conformers is found to be altered by the intramolecular relaxation of the amide I mode as a consequence of the different relaxation pathways of each conformer which modify the amount of vibrational energy stored in the torsional motions of the tripeptide, so the α → pPII and pPII → α conversion rates are increased differently. The selectivity of the process comes from the shifts of the vibrational frequencies with the conformational changes that modify the resonance conditions driving the intramolecular energy flows.

Scaling Rules for Vibrational Energy Transport in Globular Proteins

Computational studies of vibrational energy flow in biomolecules have to date mapped out transport pathways on a case-by-case basis. To provide a more general approach, we derive scaling rules for vibrational energy transport in a globular protein, which are identified from extensive nonequilibrium molecular dynamics simulations of vibrational energy flow in the villin headpiece subdomain HP36. We parametrize a master equation based on inter-residue, residue-solvent, and heater-residue energy-transfer rates, which closely reproduces the results of the all-atom simulations. From that fit, two scaling rules emerge, one for energy transport along the protein backbone which relies on a diffusion model and another for energy transport between tertiary contacts, which is based on a harmonic model. Requiring only the calculation of mean and variance of relatively few atomic distances, the approach holds the potential to predict the pathways and time scales of vibrational energy flow in large proteins.

How Photoisomerization Drives Peptide Folding and Unfolding: Insights from QM/MM and MM Dynamics Simulations

Photoswitchable azobenzene cross-linkers can control the folding and unfolding of peptides by photoisomerization and can thus regulate peptide affinities and enzyme activities. Using quantum mechanics/molecular mechanics (QM/MM) methods and classical MM force fields, we report the first molecular dynamics simulations of the photoinduced folding and unfolding processes in the azobenzene cross-linked FK-11 peptide. We find that the interactions between the peptide and the azobenzene cross-linker are crucial for controlling the evolution of the secondary structure of the peptide and responsible for accelerating the folding and unfolding events. They also modify the photoisomerization mechanism of the azobenzene cross-linker compared with the situation in vacuo or in solution.

A biomimetic molecular switch at work: coupling photoisomerization dynamics to peptide structural rearrangement

A combined experimental and computational study of a peptide-linked retinal-like molecular switch shows the effects on photoreactivity and the α-helix structure.

Long-range conformational transition of a photoswitchable allosteric protein: molecular dynamics simulation study

A local perturbation of a protein may lead to functional changes at some distal site. An example is the PDZ2 domain of human tyrosine phosphatase 1E, which shows an allosteric transition upon binding to a peptide ligand. Recently Buchli et al. presented a time-resolved study of this transition by covalently linking an azobenzene photoswitch across the binding groove and using a femtosecond laser pulse that triggers the cis-trans photoisomerization of azobenzene. To aid the interpretation of these experiments, in this work seven microsecond runs of all-atom molecular dynamics simulations each for the wild-type PDZ2 in the ligand-bound and -free state, as well as the photoswitchable protein (PDZ2S) in the cis and trans states of the photoswitch, in explicit water were conducted. First the theoretical model is validated by recalculating the available NMR data from the simulations. By comparing the results for PDZ2 and PDZ2S, it is analyzed to what extent the photoswitch indeed mimics the free-bound transition. A detailed description of the conformational rearrangement following the cis-trans photoisomerization of PDZ2S reveals a series of photoinduced structural changes that propagate from the anchor residues of the photoswitch via intermediate secondary structure segments to the C-terminus of PDZ2S. The changes of the conformational distribution of the C-terminal region is considered as the distal response of the isolated allosteric protein.

Azobenzene photoisomerization-induced destabilization of B-DNA

Molecular photoswitches provide a promising way for selective regulation of nanoscaled biological systems. It has been shown that conformational changes of azobenzene, one of the widely used photoswitches, can be used to reversibly control DNA duplex formation. Here, we investigate the conformational response of DNA upon azobenzene binding and isomerization, using a threoninol linker that has been experimentally investigated recently. To this end, nonequilibrium molecular dynamics simulations are carried out using a switching potential describing the photoinduced isomerization. Attachment of azobenzene leads to a distortion of the DNA helical conformation that is similar for the trans and cis forms. However, the trans form is stabilized by favorable stacking interactions whereas the cis form is found to remain flipped out of the basepair-stacked position. Multiple azobenzene attachment augments the distortion in DNA helical conformation. The distorted DNA retains nativelike pairing of bases at ambient temperatures, but shows weaker basepairing compared to native DNA at an elevated temperature.

Simulation of transient infrared spectra of a photoswitchable peptide

In transient infrared (IR) experiments, a molecular system may be photoexcited in a nonstationary conformational state, whose time evolution is monitored via IR spectroscopy with high temporal and structural resolution. As a theoretical formulation of these experiments, this work derives explicit expressions for transient one- and two-dimensional IR spectra and discusses various levels of approximation and sampling strategies. Adopting a photoswitchable octapeptide in water as a representative example, nonequilibrium molecular dynamics simulations are performed and the photoinduced conformational dynamics and associated IR spectra are discussed in detail. Interestingly, it is found that the time scales of dynamics and spectra may differ from residue to residue by up to an order of magnitude. Considering merely the cumulative spectrum of all residues, the contributions of the individual residues largely compensate each other, which may explain the surprisingly small frequency shifts and short photoproduct rise times found in experiment. Even when a localized amide I mode is probed (e.g., via isotope labeling), the vibrational frequency shift is shown to depend in a complicated way on the conformation of the entire peptide as well as on the interaction with the solvent. In this context, various issues concerning the interpretation of transient IR spectra and conformational dynamics in terms of a few exponential time scales are discussed.