End-to-end distance distributions and intrachain diffusion constants in unfolded polypeptide chains indicate intramolecular hydrogen bond formation

Artificial intelligence assisted identification of potential tau aggregation inhibitors: ligand- and structure-based virtual screening, in silico ADME, and molecular dynamics study

Alzheimer's disease (AD) is a severe, growing, multifactorial disorder affecting millions of people worldwide characterized by cognitive decline and neurodegeneration. The accumulation of tau protein into paired helical filaments is one of the major pathological hallmarks of AD and has gained the interest of researchers as a potential drug target to treat AD. Lately, Artificial Intelligence (AI) has revolutionized the drug discovery process by speeding it up and reducing the overall cost. As a part of our continuous effort to identify potential tau aggregation inhibitors, and leveraging the power of AI, in this study, we used a fully automated AI-assisted ligand-based virtual screening tool, PyRMD to screen a library of 12 million compounds from the ZINC database to identify potential tau aggregation inhibitors. The preliminary hits from virtual screening were filtered for similar compounds and pan-assay interference compounds (the compounds containing reactive functional groups which can interfere with the assays) using RDKit. Further, the selected compounds were prioritized based on their molecular docking score with the binding pocket of tau where the binding pockets were identified using replica exchange molecular dynamics simulation. Thirty-three compounds showing good docking scores for all the tau clusters were selected and were further subjected to in silico pharmacokinetic prediction. Finally, top 10 compounds were selected for molecular dynamics simulation and MMPBSA binding free energy calculations resulting in the identification of UNK_175, UNK_1027, UNK_1172, UNK_1173, UNK_1237, UNK_1518, and UNK_2181 as potential tau aggregation inhibitors.

A molecular dynamics simulation study of glycine/serine octapeptides labeled with 2,3-diazabicyclo[2.2.2]oct-2-ene fluorophore

The compound 2,3-diazabicyclo[2.2.2]oct-2-ene (DBO) is a versatile fluorophore widely used in Förster resonance energy transfer (FRET) spectroscopy studies due to its remarkable sensitivity, enabling precise donor-acceptor distance measurements, even for short peptides. Integrating time-resolved and FRET spectroscopies with molecular dynamics simulations provides a robust approach to unravel the structure and dynamics of biopolymers in a solution. This study investigates the structural behavior of three octapeptide variants: Trp-(Gly-Ser)3-Dbo, Trp-(GlyGly)3-Dbo, and Trp-(SerSer)3-Dbo, where Dbo represents the DBO-containing modified aspartic acid, using molecular dynamics simulations. Glycine- and serine-rich amino acid fragments, common in flexible protein regions, play essential roles in functional properties. Results show excellent agreement between end-to-end distances, orientational factors from simulations, and the available experimental and theoretical data, validating the reliability of the GROMOS force field model. The end-to-end distribution, modeled using three Gaussian distributions, reveals a complex shape, confirmed by cluster analysis highlighting a limited number of significant conformations dominating the peptide landscape. All peptides predominantly adopt a disordered state in the solvent, yet exhibit a compact shape, aligning with the model of disordered polypeptide chains in poor solvents. Conformations show marginal dependence on chain composition, with Ser-only chains exhibiting slightly more elongation. This study enhances our understanding of peptide behavior, providing valuable insights into their structural dynamics in solution.

From isolated polyelectrolytes to star-like assemblies: the role of sequence heterogeneity on the statistical structure of the intrinsically disordered neurofilament-low tail domain

Intrinsically disordered proteins (IDPs) are a subset of proteins that lack stable secondary structure. Given their polymeric nature, previous mean-field approximations have been used to describe the statistical structure of IDPs. However, the amino-acid sequence heterogeneity and complex intermolecular interaction network have significantly impeded the ability to get proper approximations. One such case is the intrinsically disordered tail domain of neurofilament low (NFLt), which comprises a 50 residue-long uncharged domain followed by a 96 residue-long negatively charged domain. Here, we measure two NFLt variants to identify the impact of the NFLt two main subdomains on its complex interactions and statistical structure. Using synchrotron small-angle x-ray scattering, we find that the uncharged domain of the NFLt induces attractive interactions that cause it to self-assemble into star-like polymer brushes. On the other hand, when the uncharged domain is truncated, the remaining charged N-terminal domains remain isolated in solution with typical polyelectrolyte characteristics. We further discuss how competing long- and short-ranged interactions within the polymer brushes dominate their ensemble structure and, in turn, their implications on previously observed phenomena in NFL native and diseased states.

Structural biases in disordered proteins are prevalent in the cell

Intrinsically disordered proteins and protein regions (IDPs) are prevalent in all proteomes and are essential to cellular function. Unlike folded proteins, IDPs exist in an ensemble of dissimilar conformations. Despite this structural plasticity, intramolecular interactions create sequence-specific structural biases that determine an IDP ensemble's three-dimensional shape. Such structural biases can be key to IDP function and are often measured in vitro, but whether those biases are preserved inside the cell is unclear. Here we show that structural biases in IDP ensembles found in vitro are recapitulated inside human-derived cells. We further reveal that structural biases can change in a sequence-dependent manner due to changes in the intracellular milieu, subcellular localization, and intramolecular interactions with tethered well-folded domains. We propose that the structural sensitivity of IDP ensembles can be leveraged for biological function, can be the underlying cause of IDP-driven pathology or can be used to design disorder-based biosensors and actuators.

Fluorescence-readout as a powerful macromolecular characterisation tool

The last few decades have witnessed significant progress in synthetic macromolecular chemistry, which can provide access to diverse macromolecules with varying structural complexities, topology and functionalities, bringing us closer to the aim of controlling soft matter material properties with molecular precision. To reach this goal, the development of advanced analytical techniques, allowing for micro-, molecular level and real-time investigation, is essential. Due to their appealing features, including high sensitivity, large contrast, fast and real-time response, as well as non-invasive characteristics, fluorescence-based techniques have emerged as a powerful tool for macromolecular characterisation to provide detailed information and give new and deep insights beyond those offered by commonly applied analytical methods. Herein, we critically examine how fluorescence phenomena, principles and techniques can be effectively exploited to characterise macromolecules and soft matter materials and to further unravel their constitution, by highlighting representative examples of recent advances across major areas of polymer and materials science, ranging from polymer molecular weight and conversion, architecture, conformation to polymer self-assembly to surfaces, gels and 3D printing. Finally, we discuss the opportunities for fluorescence-readout to further advance the development of macromolecules, leading to the design of polymers and soft matter materials with pre-determined and adaptable properties.

Probing the inner local density of complex macromolecules by pyrene excimer formation

The direct relationship existing between the average rate constant 〈k〉 for pyrene excimer formation and the local concentration [Py]loc of ground-state pyrenyl labels covalently attached to a macromolecule was established for 55 pyrene-labeled macromolecules (PyLM). These PyLM belonged to three different families of macromolecules with the first representing short monodisperse linear chains end-labeled with pyrene (polystyrene, poly(ethylene oxide), and poly(N-isopropyl acrylamide)), the second representing long polydisperse linear chains randomly labeled with pyrene (poly(methyl acrylate), poly(methyl methacrylate), polystyrene, poly(butyl methacrylate), poly(methoxyethyl methacrylate), and poly(N-isopropyl acrylamide)), and the third being comprised of two series of pyrene end-labeled low generation dendrimers with a bis(hydroxymethyl)propionic acid or a polyamidoamine backbone. The assumption, that the polymeric segments probed by an excited pyrenyl label covalently attached to one of these macromolecules obeyed Gaussian statistics, enabled the calculation of their square root average squared end-to-end distance (LPy), which was applied to calculate [Py]loc. The log-log plots of 〈k〉 as a function of [Py]loc yielded straight lines with a slope of unity for all families of macromolecules studied in four different organic solvents demonstrating the validity and generality of the 〈k〉-vs.-[Py]loc relationship. Since an experimentalist knows how the the pyrenyl labels are covalently attached onto a macromolecule, [Py]loc offers a means to probe the local density of a macromolecule, which can be employed to characterize its conformation in solution. Consequently, the 〈k〉-vs.-[Py]loc relationship provides a novel experimental means to probe the conformation of macromolecules which should establish pyrene excimer formation as an appealing method for conformational studies of macromolecules in solution, which should nicely complement scattering techniques.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. II. Effect of Temperature and Internal Friction

Internal dynamics of proteins are essential for protein folding and function. Dynamics in unfolded proteins are of particular interest since they are the basis for many cellular processes like folding, misfolding, aggregation, and amyloid formation and also determine the properties of intrinsically disordered proteins (IDPs). It is still an open question of what governs motions in unfolded proteins and whether they encounter major energy barriers. Here we use triplet-triplet energy transfer (TTET) in unfolded homopolypeptide chains and IDPs to characterize the barriers for local and long-range loop formation. The results show that the formation of short loops encounters major energy barriers with activation energies (Ea) up to 18 kJ/mol (corrected for effects of temperature on water viscosity) with very little dependence on amino acid sequence. For poly(Gly-Ser) and polySer chains the barrier decreases with increasing loop size and reaches a limiting value of 4.6 ± 0.4 kJ/mol for long and flexible chains. This observation is in accordance with the concept of internal friction encountered by chain motions due to steric effects, which is high for local motions and decreases with increasing loop size. Comparison with the results from the viscosity dependence of loop formation shows a negative correlation between Ea and the sensitivity of the reaction to solvent viscosity (α) in accordance with the Grote-Hynes theory of memory friction. The Arrhenius pre-exponential factor (A) also decreases with increasing loop size, indicating increased entropic costs for loop formation. Long-range loop formation in the investigated sequences derived from IDPs shows increased Ea and A compared with poly(Gly-Ser) and polySer chains. This increase is exclusively due to steric effects that cause additional internal friction, whereas intramolecular hydrogen bonds, dispersion forces, and charge interactions do not affect the activation parameters.

Local and Large-Scale Conformational Dynamics in Unfolded Proteins and IDPs. I. Effect of Solvent Viscosity and Macromolecular Crowding

Protein/solvent interactions largely influence protein dynamics, particularly motions in unfolded and intrinsically disordered proteins (IDPs). Here, we apply triplet-triplet energy transfer (TTET) to investigate the coupling of internal protein motions to solvent motions by determining the effect of solvent viscosity (η) and macromolecular crowding on the rate constants of loop formation (kc) in several unfolded polypeptide chains including IDPs. The results show that the viscosity dependence of loop formation depends on amino acid sequence, loop length, and co-solute size. Below a critical size (rc), co-solutes exert a maximum effect, indicating that under these conditions microviscosity experienced by chain motions matches macroviscosity of the solvent. rc depends on chain stiffness and reflects the length scale of the chain motions, i.e., it is related to the persistence length. Above rc, the effect of solvent viscosity decreases with increasing co-solute size. For co-solutes typically used to mimic cellular environments, a scaling of kc ∝ η-0.1 is observed, suggesting that dynamics in unfolded proteins are only marginally modulated in cells. The effect of solvent viscosity on kc in the small co-solute limit (below rc) increases with increasing chain length and chain flexibility. Formation of long and very flexible loops exhibits a kc ∝ η-1 viscosity dependence, indicating full solvent coupling. Shorter and less flexible loops show weaker solvent coupling with values as low as kc ∝ η-0.75 ± 0.02. Coupling of formation of short loops to solvent motions is very little affected by amino acid sequence, but solvent coupling of long-range loop formation is decreased by side chain sterics.

Diffusion-Enhanced Förster Resonance Energy Transfer in Flexible Peptides: From the Haas-Steinberg Partial Differential Equation to a Closed Analytical Expression

In the huge field of polymer structure and dynamics, including intrinsically disordered peptides, protein folding, and enzyme activity, many questions remain that cannot be answered by methodology based on artificial intelligence, X-ray, or NMR spectroscopy but maybe by fluorescence spectroscopy. The theory of Förster resonance energy transfer (FRET) describes how an optically excited fluorophore transfers its excitation energy through space to an acceptor moiety-with a rate that depends on the distance between donor and acceptor. When the donor and acceptor moiety are conjugated to different sites of a flexible peptide chain or any other linear polymer, the pair could in principle report on chain structure and dynamics, on the site-to-site distance distribution, and on the diffusion coefficient of mutual site-to-site motion of the peptide chain. However, the dependence of FRET on distance distribution and diffusion is not defined by a closed analytical expression but by a partial differential equation (PDE), by the Haas-Steinberg equation (HSE), which can only be solved by time-consuming numerical methods. As a second complication, time-resolved FRET measurements have thus far been deemed necessary. As a third complication, the evaluation requires a computationally demanding but indispensable global analysis of an extended experimental data set. These requirements have made the method accessible to only a few experts. Here, we show how the Haas-Steinberg equation leads to a closed analytical expression (CAE), the Haas-Steinberg-Jacob equation (HSJE), which relates a diffusion-diagnosing parameter, the effective donor-acceptor distance, to the augmented diffusion coefficient, J, composed of the diffusion coefficient, D, and the photophysical parameters that characterize the used FRET method. The effective donor-acceptor distance is easily retrieved either through time-resolved or steady-state fluorescence measurements. Any global fit can now be performed in seconds and minimizes the sum-of-square difference between the experimental values of the effective distance and the values obtained from the HSJE. In summary, the HSJE can give a decisive advantage in applying the speed and sensitivity of FRET spectroscopy to standing questions of polymer structure and dynamics.

Molecular Origin of Internal Friction in Intrinsically Disordered Proteins

Protein folding and dynamics are controlled by an interplay of thermal and viscosity effects. The effect of viscous drag through the solvent molecules is described by the classic Kramers theory in the high friction limit, which considers the dampening of the reactant molecules in the solution and quantifies the dependence of the reaction rate on the frictional drag. In addition to the external energy dissipation originating from the surrounding solvent molecules, there is an additional mode of internal energy dissipative force operative within the polypeptide chain reflecting the internal resistance of the chain to its conformational alterations. This dry, solvent-independent intrinsic frictional drag is termed internal friction. In the case of natively folded proteins, the physical origin of internal friction is primarily attributed to the intrachain interactions or other nonnative interactions in their compact states. However, the molecular origin of internal friction in intrinsically disordered proteins (IDPs) remains elusive.In this Account, we address this fundamental issue: what are the principal drivers of viscosity-independent (dry) friction in highly solvated, expanded, conformationally flexible, rapidly fluctuating IDPs that do not possess persistent intrachain interactions? IDPs exhibit diffusive conformational dynamics that is predominantly dominated by the segmental motion of the backbone arising due to the dihedral rotations in the Ramachandran Φ-Ψ space. The physical origin of friction in a complex biopolymeric system such as IDPs can be described by classic polymer models, namely, Rouse/Zimm models with internal friction. These one-dimensional models do not invoke torsional fluctuation components. Kuhn's classic description includes the connection between internal friction and microscopic dihedral hopping. Based on our time-resolved fluorescence anisotropy results, we describe that the sequence-dependent, collective, short-range backbone dihedral rotations govern localized internal friction in an archetypal IDP, namely, α-synuclein. The highly sensitive, residue-specific fluorescence depolarization kinetics offers a potent methodology to characterize and quantify the directional decorrelation engendered due to the short-range dihedral relaxation of the polypeptide backbone in the dihedral space. We utilized this characteristic relaxation time scale as our dynamic readout to quantify the site-specific frictional component. Our linear viscosity-dependent model of torsional relaxation time scale furnished a finite nonzero time constant at the zero solvent viscosity representing the solvent-independent internal friction. These results unveil the effect of the degree of dihedral restraining parameter on the internal friction component by showing that a restrained proline residue imparts higher torsional stiffness in the chain segments and, therefore, exhibits higher internal friction. This Account sheds light on the molecular underpinning of the sequence-specific internal friction in IDPs and will be of interest to unmask the role of internal friction in a diverse range of biomolecular processes involving binding-induced folding, allosteric interaction, protein misfolding and aggregation, and biomolecular condensation via phase separation.

Hydration Free Energies of Polypeptides from Popular Implicit Solvent Models versus All-Atom Simulation Results Based on Molecular Quasichemical Theory

Calculating the hydration free energy of a macromolecule in all-atom simulations has long remained a challenge, necessitating the use of models wherein the effect of the solvent is captured without explicit account of solvent degrees of freedom. This situation has changed with developments in the molecular quasi-chemical theory (QCT)─an approach that enables calculation of the hydration free energy of macromolecules within all-atom simulations at the same resolution as is possible for small molecular solutes. The theory also provides a rigorous and physically transparent framework to conceptualize and model interactions in molecular solutions and thus provides a convenient framework to investigate the assumptions in implicit solvent models. In this study, we compare the results using molecular QCT versus predictions from EEF1, ABSINTH, and GB/SA implicit solvent models for polyglycine and polyalanine solutes covering a range of chain lengths and conformations. The hydration free energies or the differences in hydration free energies between conformers obtained from the implicit solvent models do not agree with explicit solvent results, with the deviations being largest for the group additive EEF1 and ABSINTH models. GB/SA does better in capturing the qualitative trends seen in explicit solvent results. Analysis founded on QCT reveals the critical importance of the cooperativity of hydration that is inherent in the hydrophilic and hydrophobic contributions to hydration─physics that is not well captured in additive models but somewhat better accounted for by means of a dielectric in the GB/SA approach.

Heterogeneity in Protein Folding and Unfolding Reactions

Proteins have dynamic structures that undergo chain motions on time scales spanning from picoseconds to seconds. Resolving the resultant conformational heterogeneity is essential for gaining accurate insight into fundamental mechanistic aspects of the protein folding reaction. The use of high-resolution structural probes, sensitive to population distributions, has begun to enable the resolution of site-specific conformational heterogeneity at different stages of the folding reaction. Different states populated during protein folding, including the unfolded state, collapsed intermediate states, and even the native state, are found to possess significant conformational heterogeneity. Heterogeneity in protein folding and unfolding reactions originates from the reduced cooperativity of various kinds of physicochemical interactions between various structural elements of a protein, and between a protein and solvent. Heterogeneity may arise because of functional or evolutionary constraints. Conformational substates within the unfolded state and the collapsed intermediates that exchange at rates slower than the subsequent folding steps give rise to heterogeneity on the protein folding pathways. Multiple folding pathways are likely to represent distinct sequences of structure formation. Insight into the nature of the energy barriers separating different conformational states populated during (un)folding can also be obtained by resolving heterogeneity.

Quantitative Elasticity of Flexible Polymer Chains Using Interferometer-Based AFM

We estimate the elasticity of single polymer chains using atomic force microscope (AFM)-based oscillatory experiments. An accurate estimate of elasticity using AFM is limited by assumptions in describing the dynamics of an oscillating cantilever. Here, we use a home-built fiber-interferometry-based detection system that allows a simple and universal point-mass description of cantilever oscillations. By oscillating the cantilever base and detecting changes in cantilever oscillations with an interferometer, we extracted stiffness versus extension profiles for polymers. For polyethylene glycol (PEG) in a good solvent, stiffness-extension data showed significant deviation from conventional force-extension curves (FECs) measured in constant velocity pulling experiments. Furthermore, modeling stiffness data with an entropic worm-like chain (WLC) model yielded a persistence length of (0.5 ± 0.2 nm) compared to anomaly low value (0.12 nm ± 0.01) in conventional pulling experiments. This value also matched well with equilibrium measurements performed using magnetic tweezers. In contrast, polystyrene (PS) in a poor solvent, like water, showed no deviation between the two experiments. However, the stiffness profile for PS in good solvent (8M Urea) showed significant deviation from conventional force-extension curves. We obtained a persistence length of (0.8 ± 0.2 nm) compared to (0.22 nm ± 0.01) in pulling experiments. Our unambiguous measurements using interferometer yield physically acceptable values of persistence length. It validates the WLC model in good solvents but suggests caution for its use in poor solvents.

Coarse-Grained Simulations of Protein Folding: Bridging Theory and Experiments

Computational coarse-grained models play a fundamental role as a research tool in protein folding, and they are important in bridging theory and experiments. Folding mechanisms are generally discussed using the energy landscape framework, which is well mapped within a class of simplified structure-based models. In this chapter, simplified computer models are discussed with special focus on structure-based ones.

Glycine-Rich Peptides from FUS Have an Intrinsic Ability to Self-Assemble into Fibers and Networked Fibrils

Glycine-rich regions feature prominently in intrinsically disordered regions (IDRs) of proteins that drive phase separation and the regulated formation of membraneless biomolecular condensates. Interestingly, the Gly-rich IDRs seldom feature poly-Gly tracts. The protein fused in sarcoma (FUS) is an exception. This protein includes two 10-residue poly-Gly tracts within the prion-like domain (PLD) and at the interface between the PLD and the RNA binding domain. Poly-Gly tracts are known to be highly insoluble, being potent drivers of self-assembly into solid-like fibrils. Given that the internal concentrations of FUS and FUS-like molecules cross the high micromolar and even millimolar range within condensates, we reasoned that the intrinsic insolubility of poly-Gly tracts might be germane to emergent fluid-to-solid transitions within condensates. To assess this possibility, we characterized the concentration-dependent self-assembly for three non-overlapping 25-residue Gly-rich peptides derived from FUS. Two of the three peptides feature 10-residue poly-Gly tracts. These peptides form either long fibrils based on twisted ribbon-like structures or self-supporting gels based on physical cross-links of fibrils. Conversely, the peptide with similar Gly contents but lacking a poly-Gly tract does not form fibrils or gels. Instead, it remains soluble across a wide range of concentrations. Our findings highlight the ability of poly-Gly tracts within IDRs that drive phase separation to undergo self-assembly. We propose that these tracts are likely to contribute to nucleation of fibrillar solids within dense condensates formed by FUS.

Effect of Organic Solvents on the Structure and Activity of a Minimal Lipase

Lipases are ubiquitously used in chemo-enzymatic synthesis and industrial applications. Nevertheless, the modulation of the activity of lipases by organic solvents still is not fully understood at the molecular level. We systematically investigated the activity and structure of lipase A from Bacillus subtilis in binary water-organic solvent mixtures of dimethyl sulfoxide (DMSO), acetonitrile (ACN), and isopropyl alcohol (IPA) using activity assays, fluorescence spectroscopy, molecular dynamics (MD) simulations, and FRET/MD analysis. The enzymatic activity strongly depended on the type and amount of organic solvent in the reaction media. Whereas IPA and ACN reduced the activity of the enzyme, small concentrations of DMSO led to lipase activation via an uncompetitive mechanism. DMSO molecules did not directly interfere with the binding of the substrate in the active site, contrary to what is known for other solvents and enzymes. We propose that the His156-Asp133 interaction, the binding of organic molecules to the active site, and the water accessibility of the substrate are key factors modulating the catalytic activity. Furthermore, we rationalized the role of solvent descriptors on the regulation of enzymatic activity in mixtures with low concentrations of the organic molecule, with prospective implications for the optimization of biocatalytic processes via solvent tuning.

Systematic transition modeling analysis in the MEF2B-DNA binding interface due to Y69H and K4E variants

Transcriptional coactivator myocyte enhancer factor 2B (MEF2B) mutations are the most common cause of germinal center-derived B-cell non-Hodgkin lymphoma. Despite well-established contributions in lymphomagenesis, the structure-function paradigms of these mutations are largely unknown. Here through in silico approaches, we present structural evaluation of two reported missense variants (K4E and Y69H) in MEF2B to investigate their impact on DNA-binding through molecular dynamics simulation assays. Notably, MEF2B-specific MADs box domain (Lys23, Arg24 and Lys31) and N-terminal loop residues (Gly2, Arg3, Lys4, Lys5, Ile6 and Asn13) contribute in DNA binding, while in MEF2B^K4E, DNA binding is facilitated by Gly2, Arg3 and Arg91 (α3) residues. Conversely, in MEF2B^Y69H, Arg3, Lys5, Ser78, Arg79 and Asn81 residues mediate DNA binding. DNA binding induces pronounced conformational readjustments in MEF2B^WT-specific α1-N-terminal loop region, while MEF2B^Y69H and MEF2B^K4E exhibit fluctuations in both α1 and α3. Hydrogen (H)-bond occupancy analysis reveals a similar DNA binding behavior for MEF2^WT and MEF2B^Y69H, compared to MEF2B^K4E structure. The Anisotropic Network Model analysis depicts α1 and α3 as more fluctuant regions in MEF2B^K4E as compared to other systems. MEF2B^WT and MEF2B^K4E, Tyr69 residue is involved in p300 binding thus possible influence of Y69H variation in the functions other than DNA binding, such as p300 co-activator recruitment may explain the reduced transcriptional activation of MEF2B^Y69H. Thus, present study may provide a structural basis of DNA recognition by pinpointing the underlying conformational changes in the dynamics of MEF2B^K4E, MEF2B^Y69H, and MEF2B^WT structures that may contribute in the identification of novel therapeutic strategies for lymphomagenesis.

Highly Emissive Biological Bilirubin Molecules: Shedding New Light on the Phototherapy Scheme

Bilirubin (BR) is the main end-product of the hemoglobin catabolism. For decades, its photophysics has been mainly discussed in terms of ultrafast deactivation of the excited state in solution, where, indeed, BR shows a very low green emission quantum yield (EQY), 0.03%, resulting from an efficient nonradiative isomerization process. Herein, we present, for the first time, unique and exceptional photophysical properties of solid-state BR, which amend by changing the type of crystal, from a closely packed α crystal to an amorphous loosely packed β crystal. BR α crystals show a very bright red emission with an EQY of ca. 24%, whereas β crystals present, in addition, a low green EQY of ca. 0.5%. By combining density functional theory (DFT) calculations and time-resolved emission spectroscopy, we trace back this dual emission to the presence of two types of BR molecules in the crystal: a "stiff" monomer, M1, distorted by particularly strong internal H-bonds and a "floppy" monomer, M2, having a structure close to that of BR in solution. We assign the red strong emission of BR crystals to M1 present in both the α and β crystals, while the low green emission, only present in the amorphous (β) crystal, is interpreted as M2 emission. Efficient energy-transfer processes from M2 to M1 in the closely packed α crystal are invoked to explain the absence of the green component in its emission spectrum. Interestingly, these unique photophysical properties of BR remain in polar solvents such as water. Based on these unprecedented findings, we propose a new model for the phototherapy scheme of BR inside the human body and highlight the usefulness of BR as a strong biological fluorescent probe.

Screening potential FDA-approved inhibitors of the SARS-CoV-2 major protease 3CLpro through high-throughput virtual screening and molecular dynamics simulation

It has been confirmed that the new coronavirus SARS-CoV-2 caused the global pandemic of coronavirus disease 2019 (COVID-19). Studies have found that 3-chymotrypsin-like protease (3CL^pro) is an essential enzyme for virus replication, and could be used as a potential target to inhibit SARS-CoV-2. In this work, 3CL^pro was used as the target to complete the high-throughput virtual screening of the FDA-approved drugs, and Indinavir and other 10 drugs with high docking scores for 3CL^pro were obtained. Studies on the binding pattern of 3CL^pro and Indinavir found that Indinavir could form the stable hydrogen bond (H-bond) interactions with the catalytic dyad residues His41-Cys145. Binding free energy study found that Indinavir had high binding affinity with 3CL^pro. Subsequently, molecular dynamics simulations were performed on the 3CL^pro and 3CL^pro-Indinavir systems, respectively. The post-dynamic analyses showed that the conformational state of the 3CLpro-Indinavir system transformed significantly and the system tended to be more stable. Moreover, analyses of the residue interaction network (RIN) and H-bond occupancy revealed that the residue-residue interaction at the catalytic site of 3CL^pro was significantly enhanced after binding with Indinavir, which in turn inactivated the protein. In short, through this research, we hope to provide more valuable clues against COVID-19.

The nano-scale viscoelasticity using atomic force microscopy in liquid environment

We measured viscoelasticity of two nanoscale systems, single protein molecules and molecular layers of water confined between solid walls. In order to quantify the viscoelastic response of these nanoscale systems in liquid environment, the measurements are performed using two types of atomic force microscopes (AFMs), which employ different detection schemes to measure the cantilever response. We used a deflection detection scheme, available in commercial AFMs, that measures cantilever bending and a fibre-interferometer based detection which measures cantilever displacement. The hydrodynamics of the cantilever is modelled using Euler-Bernoulli equation with appropriate boundary conditions which accommodate both detection schemes. In a direct contradiction with many reports in the literature, the dissipation coefficient of a single octomer of titin I27₈ is found to be immeasurably low. The upper bound on the dissipation coefficient is 5 × 10^-7 kg s^-1, which is much lower than the reported values. The entropic stiffness of single unfolded domains of protein measured using both methods is in the range of 10 mN m^-1. We show that in a conventional deflection detection measurement, the phase of the bending signal can be a primary source of artefacts in the dissipation estimates. It is recognized that the measurement of cantilever displacement, which has negligibly small phase lag due to hydrodynamics of the cantilever at low excitation frequencies, is better suited for ensuring artefact-free measurement of viscoelasticity compared to the measurement of the cantilever bending. Further, it was possible to measure dissipation in molecular layers of water confined between the tip and the substrate using fibre interferometer based AFM with similar experimental parameters. It confirms that the dissipation coefficient of a single I27₈ is below the detection limit of AFM. The results shed light on the discrepancy observed in the measured diffusional dynamics of protein collapse measured using Force spectroscopic techniques and single-molecule optical techniques.

Automated and optimally FRET-assisted structural modeling

FRET experiments can provide state-specific structural information of complex dynamic biomolecular assemblies. However, to overcome the sparsity of FRET experiments, they need to be combined with computer simulations. We introduce a program suite with (i) an automated design tool for FRET experiments, which determines how many and which FRET pairs should be used to minimize the uncertainty and maximize the accuracy of an integrative structure, (ii) an efficient approach for FRET-assisted coarse-grained structural modeling, and all-atom molecular dynamics simulations-based refinement, and (iii) a quantitative quality estimate for judging the accuracy of FRET-derived structures as opposed to precision. We benchmark our tools against simulated and experimental data of proteins with multiple conformational states and demonstrate an accuracy of ~3 Å RMSD_Cα against X-ray structures for sets of 15 to 23 FRET pairs. Free and open-source software for the introduced workflow is available at https://github.com/Fluorescence-Tools. A web server for FRET-assisted structural modeling of proteins is available at http://nmsim.de.

To overcome the limitation of FRET data being too sparse to cover all structural details, FRET experiments need to be carefully designed and complemented with simulations. Here the authors present a toolkit for automated design of FRET experiments, which determines how many and which FRET pairs should be used to maximize the accuracy, and for FRET-assisted structural modeling and refinement at the atomistic level.

Kinetics of Loop Closure in Disordered Proteins: Theory vs Simulations vs Experiments

We study intrachain dynamics of intrinsically disordered proteins, as manifested by the time scales of loop formation, using atomistic simulations, experiment-parametrized coarse-grained models, and one-dimensional theories assuming Markov or non-Markov dynamics along the reaction coordinate. Despite the generally non-Markov character of monomer dynamics in polymers, we find that the simplest model of one-dimensional diffusion along the reaction coordinate (equated to the distance between the loop-forming monomers) well captures the mean first passage times to loop closure measured in coarse-grained and atomistic simulations, which, in turn, agree with the experimental values. This justifies use of the one-dimensional diffusion model in interpretation of experimental data. At the same time, the transition path times for loop closure in longer polypeptide chains show significant non-Markov effects; at intermediate times, these effects are better captured by the generalized Langevin equation model. At long times, however, atomistic simulations predict long tails in the distributions of transition path times, which are at odds with both the one-dimensional diffusion model and the generalized Langevin equation model.

Intrinsically disordered proteins of viruses: Involvement in the mechanism of cell regulation and pathogenesis

Intrinsically disordered proteins (IDPs) possess the property of inherent flexibility and can be distinguished from other proteins in terms of lack of any fixed structure. Such dynamic behavior of IDPs earned the name "Dancing Proteins." The exploration of these dancing proteins in viruses has just started and crucial details such as correlation of rapid evolution, high rate of mutation and accumulation of disordered contents in viral proteome at least understood partially. In order to gain a complete understanding of this correlation, there is a need to decipher the complexity of viral mediated cell hijacking and pathogenesis in the host organism. Further there is necessity to identify the specific patterns within viral and host IDPs such as aggregation; Molecular recognition features (MoRFs) and their association to virulence, host range and rate of evolution of viruses in order to tackle the viral-mediated diseases. The current book chapter summarizes the aforementioned details and suggests the novel opportunities for further research of IDPs senses in viruses.

Protein structural changes characterized by high-pressure, pulsed field gradient diffusion NMR spectroscopy

Pulsed-field gradient NMR spectroscopy is widely used to measure the translational diffusion and hydrodynamic radius (R_h) of biomolecules in solution. For unfolded proteins, the R_h provides a sensitive reporter on the ensemble-averaged conformation and the extent of polypeptide chain expansion as a function of added denaturant. Hydrostatic pressure is a convenient and reversible alternative to chemical denaturants for the study of protein folding, and enables NMR measurements to be performed on a single sample. While the impact of pressure on the viscosity of water is well known, and our water diffusivity measurements agree closely with theoretical expectations, we find that elevated pressures increase the R_h of dioxane and other small molecules by amounts that correlate with their hydrophobicity, with parallel increases in rotational friction indicated by ¹³C longitudinal relaxation times. These data point to a tighter coupling with water for hydrophobic surfaces at elevated pressures. Translational diffusion measurement of the unfolded state of a pressure-sensitized ubiquitin mutant (VA2-ubiquitin) as a function of hydrostatic pressure or urea concentration shows that R_h values of both the folded and the unfolded states remain nearly invariant. At ca 23 Å, the R_h of the fully pressure-denatured state is essentially indistinguishable from the urea-denatured state, and close to the value expected for an idealized random coil of 76 residues. The intrinsically disordered protein (IDP) α-synuclein shows slight compaction at pressures above 2 kbar. Diffusion of unfolded ubiquitin and α-synuclein is significantly impacted by sample concentration, indicating that quantitative measurements need to be carried out under dilute conditions.

Emerging consensus on the collapse of unfolded and intrinsically disordered proteins in water

Establishing the degree of collapse of unfolded or disordered proteins is a fundamental problem in biophysics, because of its relation to protein folding and to the function of intrinsically disordered proteins. However, until recently, different experiments gave qualitatively different results on collapse and there were large discrepancies between experiments and all-atom simulations. New methodology introduced in the past three years has helped to resolve the differences between experiments, and improvements in simulations have closed the gap between experiment and simulation. These advances have led to an emerging consensus on the collapse of disordered proteins in water.

Memory effects in single-molecule force spectroscopy measurements of biomolecular folding

Folding is generally assumed to be a Markov process, without memory. When the molecular motion is coupled to that of a probe as in single-molecule force spectroscopy (SMFS) experiments, however, theory predicts that the coupling to a second Markov process should induce memory when monitoring a projection of the full multi-dimensional motion onto a reduced coordinate. We developed a method to evaluate the time constant of the induced memory from its effects on the autocorrelation function, which can be readily determined from experimental data. Applying this method to both simulated SMFS measurements and experimental trajectories of DNA hairpin folding measured by optical tweezers as a model system, we validated the prediction that the linker induces memory. For these measurements, the timescale of the induced memory was found to be similar to the time required for the force probe to respond to changes in the molecule, and in the regime where the experimentally observed dynamics were not significantly perturbed by probe-molecule coupling artifacts. Memory effects are thus a general feature of SMFS measurements induced by the mechanical connection between the molecule and force probe that should be considered when interpreting experimental data.

The Influence of Molecular Reach and Diffusivity on the Efficacy of Membrane-Confined Reactions

Signaling by surface receptors often relies on tethered reactions whereby an enzyme bound to the cytoplasmic tail of a receptor catalyzes reactions on substrates within reach. The overall length and stiffness of the receptor tail, the enzyme, and the substrate determine a biophysical parameter termed the molecular reach of the reaction. This parameter determines the probability that the receptor-tethered enzyme will contact the substrate in the volume proximal to the membrane when separated by different distances within the membrane plane. In this work, we develop particle-based stochastic reaction-diffusion models to study the interplay between molecular reach and diffusion. We find that increasing the molecular reach can increase reaction efficacy for slowly diffusing receptors, whereas for rapidly diffusing receptors, increasing molecular reach reduces reaction efficacy. In contrast, if reactions are forced to take place within the two-dimensional plasma membrane instead of the three-dimensional volume proximal to it or if molecules diffuse in three dimensions, increasing molecular reach increases reaction efficacy for all diffusivities. We show results in the context of immune checkpoint receptors (PD-1 dephosphorylating CD28), a standard opposing kinase-phosphatase reaction, and a minimal two-particle model. The work highlights the importance of the three-dimensional nature of many two-dimensional membrane-confined interactions, illustrating a role for molecular reach in controlling biochemical reactions.

Exploring the effect of inhibitor AKB-9778 on VE-PTP by molecular docking and molecular dynamics simulation

Diabetic macular edema, also known as diabetic eye disease, is mainly caused by the overexpression of vascular endothelial protein tyrosine phosphatase (VE-PTP) at hypoxia/ischemic. AKB-9778 is a known VE-PTP inhibitor that can effectively interact with the active site of VE-PTP to inhibit the activity of VE-PTP. However, the binding pattern of VE-PTP with AKB-9778 and the dynamic implications of AKB-9778 on VE-PTP system at the molecular level are poorly understood. Through molecular docking, it was found that the AKB-9778 was docked well in the binding pocket of VE-PTP by the interactions of hydrogen bond and Van der Waals. Furthermore, after molecular dynamic simulations on VE-PTP system and VE-PTP ^AKB-9778 system, a series of postdynamic analyses found that the flexibility and conformation of the active site undergone an obvious transition after VE-PTP binding with AKB-9778. Moreover, by constructing the RIN, it was found that the different interactions in the active site were the detailed reasons for the conformational differences between these two systems. Thus, the finding here might provide a deeper understanding of AKB-9778 as VE-PTP Inhibitor.

Microsecond sub-domain motions and the folding and misfolding of the mouse prion protein

Protein aggregation appears to originate from partially unfolded conformations that are sampled through stochastic fluctuations of the native protein. It has been a challenge to characterize these fluctuations, under native like conditions. Here, the conformational dynamics of the full-length (23-231) mouse prion protein were studied under native conditions, using photoinduced electron transfer coupled to fluorescence correlation spectroscopy (PET-FCS). The slowest fluctuations could be associated with the folding of the unfolded state to an intermediate state, by the use of microsecond mixing experiments. The two faster fluctuations observed by PET-FCS, could be attributed to fluctuations within the native state ensemble. The addition of salt, which is known to initiate the aggregation of the protein, resulted in an enhancement in the time scale of fluctuations in the core of the protein. The results indicate the importance of native state dynamics in initiating the aggregation of proteins.

Probing Position-Dependent Diffusion in Folding Reactions Using Single-Molecule Force Spectroscopy

Folding of proteins and nucleic acids involves a diffusive search over a multidimensional conformational energy landscape for the minimal-energy structure. When examining the projection of conformational motions onto a one-dimensional reaction coordinate, as done in most experiments, the diffusion coefficient D is generally position dependent. However, it has proven challenging to measure such position-dependence experimentally. We investigated the position-dependence of D in the folding of DNA hairpins as a simple model system in two ways: first, by analyzing the round-trip time to return to a given extension in constant-force extension trajectories measured by force spectroscopy, and second, by analyzing the fall time required to reach a given extension in force jump measurements. These methods yielded conflicting results: the fall time implied a fairly constant D, but the round-trip time implied variations of over an order of magnitude. Comparison of experiments with computational simulations revealed that both methods were strongly affected by experimental artifacts inherent to force spectroscopy measurements, which obscured the intrinsic position-dependence of D. Lastly, we applied Kramers's theory to the kinetics of hairpins with energy barriers located at different positions along the hairpin stem, as a crude probe of D at different stem positions, and we found that D did not vary much as the barrier was moved along the reaction coordinate. This work underlines the difficulties faced when trying to deduce position-dependent diffusion coefficients from experimental folding trajectories.

The role of solvent quality and chain stiffness on the end-to-end contact kinetics of semiflexible polymers

Motivated by loop closure during protein folding and DNA packing, we systemically studied the effects of the solvent quality and chain stiffness on the thermodynamics and kinetics of the end-to-end contact formation for semiflexible polymer chains with reactive ends by Langevin dynamics simulations. In thermodynamics, a rich variety of products of the end-to-end contact have been discovered, such as loop, hairpin, toroid, and rodlike bundle, the populations of which are dependent on the solvent quality and chain stiffness. In kinetics, the overall pathways to form the end-to-end contact have been identified. The change of solvent quality and chain stiffness can tune the roughness of energy landscape and modulate the kinetic partitioning of the end-to-end contact formation pathways, leading to differing kinetic behaviors. In good or poor solvents, the first end-to-end contact rate k _c decreases with increasing the strength of bending stiffness k _θ monotonically. In very poor solvents, however, the dependence of the logarithm of the first end-to-end contact rate ln k _c on k _θ exhibits erratic behavior, which stems from more rugged energy landscape due to the polymer chain getting trapped into the intermediate state composed of the rodlike bundles with two ends in separation. For semiflexible chains, with increasing chain length N, the rate k _c increases initially and then decreases: in good solvents, the rate k _c exhibits a power-law relationship to chain length N with an exponent of ∼-1.50 in the region of long chains, which is in good agreement with the value derived from the experiment in the asymptotic limit of large N; and in poor solvents, the rate k _c exhibits a significantly stronger chain length dependence than those observed in good solvents in the region of long chains due to frustration to form the end-to-end contact along a specific path, especially the scaling exponent between the rate k _c and chain length N is ∼-3.62 for the case of polymer chains with k _θ = 4 at the solvent quality ε _ij = 1, in accord with the value obtained from the experiments.

The Conformational Ensemble of Polyglutamine-14 Chains: Specific Influences of Solubility Tail and Chromophores

We address polyglutamine-14 in aqueous solution with specific chromophores and a solubility chain by means of a multiscale simulation approach, combining atomistic molecular dynamics simulations and coarse-grained Monte-Carlo conformational sampling. Despite the intrinsically disordered nature of the amyloidogenic polyglutamine, we observe transient characteristic structural motifs which exhibit a specific hydrogen bonding pattern. We illustrate the relationship between structure pattern and the distance distribution of a pair of chromophores attached to the peptide termini, in light of specific influence of a short solubility tail and the chromophores themselves on the conformational ensemble.

Two Orders of Magnitude Variation of Diffusion-Enhanced Förster Resonance Energy Transfer in Polypeptide Chains

A flexible peptide chain displays structural and dynamic properties that correspond to its folding and biological activity. These properties are mirrored in intrachain site-to-site distances and diffusion coefficients of mutual site-to-site motion. Both distance distribution and diffusion determine the extent of Förster resonance energy transfer (FRET) between two sites labeled with a FRET donor and acceptor. The relatively large Förster radii of traditional FRET methods (R₀ > 20 Å) lead to a fairly low contribution of diffusion. We introduced short-distance FRET (sdFRET) where Dbo, an asparagine residue conjugated to 2,3-diazabicyclo[2.2.2]octane, acts as acceptor paired with donors, such as naphtylalanine (NAla), tryptophan, 5-l-fluorotryptophan, or tyrosine. The Förster radii are always close to 10 Å, which makes sdFRET highly sensitive to diffusional motion. We recently found indications that the FRET enhancement caused by diffusion depends symmetrically on the product of the radiative fluorescence lifetime of the donor and the diffusion coefficient. In this study, we varied this product by two orders of magnitude, using both donors of different lifetime, NAla and FTrp, as well as a varying viscogen concentration, to corroborate this statement. We demonstrate the consequences of this relationship in evaluating the impact of viscogenic coadditives on peptide dimensions.

Modulating charge patterning and ionic strength as a strategy to induce conformational changes in intrinsically disordered proteins

We present an analytical theory to describe conformational changes as a function of salt for polymers with a given sequence of charges. We apply this model to describe Intrinsically Disordered Proteins (IDPs) by explicitly accounting for charged residues and their exact placement in the primary sequence while approximating the effect of non-electrostatic interactions at a mean-field level by effective short-range (two body and three-body) interaction parameters. The effect of ions is introduced by treating electrostatic interactions within Debye-Huckle approximation. Using typical values of the short-range mean-field parameters derived from all-atom Monte Carlo simulations (at zero salt), we predict the conformational changes as a function of salt concentration. We notice that conformational transitions in response to changes in ionic strength strongly depend on sequence specific charge patterning. For example, globule to coil transition can be observed upon increasing salt concentration, in stark contrast to uniformly charged polyelectrolyte theories based on net charge only. In addition, it is possible to observe non-monotonic behavior with salt as well. Drastic differences in salt-induced conformational transitions is also evident between two doubly phosphorylated sequences-derived from the same wild type sequence-that only differ in the site of phosphorylation. Similar effects are also predicted between two sequences derived from the same parent sequence differing by a single site mutation where a negative charge is replaced by a positive charge. These effects are purely a result of charge decoration and can only be understood in terms of metrics based on specific placement of charges, and cannot be explained by models based on charge composition alone. Identifying sequences and hot spots within sequences-for post translational modification or charge mutation-using our high-throughput theory will yield fundamental insights into design and biological regulation mediated by phosphorylation and/or local changes in salt concentration.

Dynamics of Disordered Proteins under Confinement: Memory Effects and Internal Friction

Many proteins are disordered under physiological conditions. How efficiently they can search for their cellular targets and how fast they can fold upon target binding is determined by their intrinsic dynamics, which have thus attracted much recent attention. Experiments and molecular simulations show that the inherent reconfiguration timescale for unfolded proteins has a solvent friction component and an internal friction component, and the microscopic origin of the latter, along with its proper mathematical description, has been a topic of considerable debate. Internal friction varies across different proteins of comparable length and increases with decreasing denaturant concentration, showing that it depends on how compact the protein is. Here we report on a systematic atomistic simulation study of how confinement, which induces a more compact unfolded state, affects dynamics and friction in disordered peptides. We find that the average reconfiguration timescales increase exponentially as the peptide's spatial dimensions are reduced; at the same time, confinement broadens the spectrum of relaxation timescales exhibited by the peptide. There are two important implications of this broadening: First, it limits applicability of the common Rouse and Zimm models with internal friction, as those models attempt to capture internal friction effects by introducing a single internal friction timescale. Second, the long-tailed distribution of relaxation times leads to anomalous diffusion effects in the dynamics of intramolecular distances. Analysis and interpretation of experimental signals from various measurements that probe intramolecular protein dynamics (such as single-molecule fluorescence correlation spectroscopy and single-molecule force spectroscopy) rely on the assumption of diffusive dynamics along the distances being probed; hence, our results suggest the need for more general models allowing for anomalous diffusion effects.

Disordered peptide chains in an α-C-based coarse-grained model

We construct a one-bead-per-residue coarse-grained dynamical model to describe intrinsically disordered proteins at significantly longer timescales than in the all-atom models. In this model, inter-residue contacts form and disappear during the course of the time evolution. The contacts may arise between the sidechains, the backbones or the sidechains and backbones of the interacting residues. The model yields results that are consistent with many all-atom and experimental data on these systems. We demonstrate that the geometrical properties of various homopeptides differ substantially in this model. In particular, the average radius of gyration scales with the sequence length in a residue-dependent manner.

Conformational response to charge clustering in synthetic intrinsically disordered proteins

BACKGROUND

Recent theoretical and computational studies have shown that the charge content and, most importantly, the linear distribution of opposite charges are major determinants of conformational properties of intrinsically disordered proteins (IDPs). Charge segregation in a sequence can be measured through κ, which represents a normalized measure of charge asymmetry. A strong inverse correlation between κ and radius of gyration has been previously demonstrated for two independent sets of permutated IDP sequences.

METHODS

We used two well-characterized IDPs, namely measles virus N_TAIL and Hendra virus PNT4, sharing a very similar fraction of charged residues and net charge per residue, but differing in proline (Pro) content. For each protein, we have rationally designed a low- and a high-κ variant endowed with the highest and the lowest κ values compatible with their natural amino acid composition. Then, the conformational properties of wild-type and κ-variants have been assessed by biochemical and biophysical techniques.

RESULTS

We confirmed a direct correlation between κ and protein compaction. The analysis of our original data along with those available from the literature suggests that Pro content may affects the responsiveness to charge clustering.

CONCLUSIONS

Charge clustering promotes IDP compaction, but the extent of its effects depends on the sequence context. Proline residues seem to play a role contrasting compaction.

GENERAL SIGNIFICANCE

These results contribute to the identification of sequence determinants of IDP conformational properties. They may also serve as an asset for rational design of non-natural IDPs with tunable degree of compactness.

Identification and nanomechanical characterization of the fundamental single-strand protofilaments of amyloid α-synuclein fibrils

The formation and spreading of amyloid aggregates from the presynaptic protein α-synuclein in the brain play central roles in the pathogenesis of Parkinson's disease. Here, we use high-resolution atomic force microscopy to investigate the early oligomerization events of α-synuclein with single monomer angstrom resolution. We identify, visualize, and characterize directly the smallest elementary unit in the hierarchical assembly of amyloid fibrils, termed here single-strand protofilaments. We show that protofilaments form from the direct molecular assembly of unfolded monomeric α-synuclein polypeptide chains. To unravel protofilaments' internal structure and elastic properties, we manipulated nanomechanically these species by atomic force spectroscopy. The single-molecule scale identification and characterization of the fundamental unit of amyloid assemblies provide insights into early events underlying their formation and shed light on opportunities for therapeutic intervention at the early stages of aberrant protein self-assembly.

Method-Unifying View of Loop-Formation Kinetics in Peptide and Protein Folding

Protein folding can be described as a probabilistic succession of events in which the peptide chain forms loops closed by specific amino acid residue contacts, herein referred to as loop nodes. To measure loop rates, several photophysical methods have been introduced where a pair of optically active probes is incorporated at selected chain positions and the excited probe undergoes contact quenching (CQ) upon collision with the second probe. The quenching mechanisms involved triplet-triplet energy transfer, photoinduced electron transfer, and collision-induced fluorescence quenching, where the fluorescence of Dbo, an asparagine residue conjugated to 2,3-diazabicyclo[2.2.2]octane, is quenched by tryptophan. The discrepancy between the loop rates afforded from these three CQ techniques has, however, remained unresolved. In analyzing this discrepancy, we now report two short-distance FRET methods where Dbo acts as an energy acceptor in combination with tryptophan and naphtylalanine, two donors with largely different fluorescence lifetimes of 1.3 and 33 ns, respectively. Despite the different quenching mechanisms, the rates from FRET and CQ methods were, surprisingly, of comparable magnitude. This combination of FRET and CQ data led to a unifying physical model and to the conclusion that the rate of loop formation in folding reactions varies not only with the kind and number of residues that constitute the chain but also in particular with the size and properties of the residues that constitute the loop node.

Hypothesis: structural heterogeneity of the unfolded proteins originating from the coupling of the local clusters and the long-range distance distribution

We propose a hypothesis that explains two apparently contradicting observations for the heterogeneity of the unfolded proteins. First, the line confocal method of the single-molecule Förster resonance energy transfer (sm-FRET) spectroscopy revealed that the unfolded proteins possess broad peaks in the FRET efficiency plot, implying the significant heterogeneity that lasts longer than milliseconds. Second, the fluorescence correlation method demonstrated that the unfolded proteins fluctuate in the time scale shorter than 100 ns. To formulate the hypothesis, we first summarize the recent consensus for the structure and dynamics of the unfolded proteins. We next discuss the conventional method of the sm-FRET spectroscopy and its limitations for the analysis of the unfolded proteins, followed by the advantages of the line confocal method that revealed the heterogeneity. Finally, we propose that the structural heterogeneity formed by the local clustering of hydrophobic residues modulates the distribution of the long-range distance between the labeled chromophores, resulting in the broadening of the peak in the FRET efficiency plot. A clustering of hydrophobic residues around the chromophore might further contribute to the broadening. The proposed clusters are important for the understanding of protein folding mechanism.

Tuning the Flexibility of Glycine-Serine Linkers To Allow Rational Design of Multidomain Proteins

Flexible polypeptide linkers composed of glycine and serine are important components of engineered multidomain proteins. We have previously shown that the conformational properties of Gly-Gly-Ser repeat linkers can be quantitatively understood by comparing experimentally determined Förster resonance energy transfer (FRET) efficiencies of ECFP-linker-EYFP proteins to theoretical FRET efficiencies calculated using wormlike chain and Gaussian chain models. Here we extend this analysis to include linkers with different glycine contents. We determined the FRET efficiencies of ECFP-linker-EYFP proteins with linkers ranging in length from 25 to 73 amino acids and with glycine contents of 33.3% (GSSGSS), 16.7% (GSSSSSS), and 0% (SSSSSSS). The FRET efficiency decreased with an increasing linker length and was overall lower for linkers with less glycine. Modeling the linkers using the WLC model revealed that the experimentally observed FRET efficiencies were consistent with persistence lengths of 4.5, 4.8, and 6.2 Å for the GSSGSS, GSSSSS, and SSSSSS linkers, respectively. The observed increase in linker stiffness with reduced glycine content is much less pronounced than that predicted by a classical model developed by Flory and co-workers. We discuss possible reasons for this discrepancy as well as implications for using the stiffer linkers to control the effective concentrations of connected domains in engineered multidomain proteins.

Conformational Heterogeneity and FRET Data Interpretation for Dimensions of Unfolded Proteins

A mathematico-physically valid formulation is required to infer properties of disordered protein conformations from single-molecule Förster resonance energy transfer (smFRET). Conformational dimensions inferred by conventional approaches that presume a homogeneous conformational ensemble can be unphysical. When all possible-heterogeneous as well as homogeneous-conformational distributions are taken into account without prejudgment, a single value of average transfer efficiency 〈E〉 between dyes at two chain ends is generally consistent with highly diverse, multiple values of the average radius of gyration 〈R_g〉. Here we utilize unbiased conformational statistics from a coarse-grained explicit-chain model to establish a general logical framework to quantify this fundamental ambiguity in smFRET inference. As an application, we address the long-standing controversy regarding the denaturant dependence of 〈R_g〉 of unfolded proteins, focusing on Protein L as an example. Conventional smFRET inference concluded that 〈R_g〉 of unfolded Protein L is highly sensitive to [GuHCl], but data from SAXS suggested a near-constant 〈R_g〉 irrespective of [GuHCl]. Strikingly, our analysis indicates that although the reported 〈E〉 values for Protein L at [GuHCl] = 1 and 7 M are very different at 0.75 and 0.45, respectively, the Bayesian R_g² distributions consistent with these two 〈E〉 values overlap by as much as 75%. Our findings suggest, in general, that the smFRET-SAXS discrepancy regarding unfolded protein dimensions likely arise from highly heterogeneous conformational ensembles at low or zero denaturant, and that additional experimental probes are needed to ascertain the nature of this heterogeneity.

Cucurbit[7]uril Induced Formation of FRET-Enabled Unilamellar Lipid Vesicles

A unique fluorescence resonance energy transfer (FRET) process is found to be operational in a unilamellar lipid self-assembly in the aqueous phase. A newly synthesized naphthyl based long chain lipid derivative [N-(naphthalene-1-ylmethyl)tetradecane-1-ammonium chloride, 14NA⁺] forms various self-assembled architectures in the aqueous phase. Controlled changes in lipid concentration lead to a transition of the self-assemblies from micelles to vesicles to rods. In the presence of cucurbit[7]uril (CB7), 14NA⁺ forms a host-guest [2]pseudorotaxane complex (CB7∋14NA⁺) and secondary interactions lead to the formation of a lipid bilayer with hydrophobic pockets situated in between the layers. The change in the structure of 14NA⁺ assemblies, interaction with CB7 and formation of supramolecular assemblies of CB7∋14NA⁺ were examined using light scattering, spectroscopic, and microscopic techniques. Entrapment of a luminescent dye, anthracene within the hydrophobic bilayer of the supramolecular assembly CB7∋14NA⁺ favors a modified luminescent response due to an efficient FRET process. Further, the FRET process could be controlled by thermal and chemical stimuli that induce transformation of unilamellar vesicles.