Conformation of prion protein repeat peptides probed by FRET measurements and molecular dynamics simulations

Metal Ions Bound to Prion Protein Affect its Interaction with Plasminogen Activation System

Prion diseases are a group of neurodegenerative diseases, which can progress rapidly. Previous data have demonstrated that prion protein (PrP) stimulates activation of plasminogen (Plg) by tissue plasminogen activator (tPA). In this study, using spectroscopic method, we aimed to determine whether PrP’s role in activating Plg is influenced by metal binding. We also investigated the region in PrP involved in binding to tPA and Plg, and whether PrP in fibrillar form behaves the same way as PrP unbound to any metal ion i.e., apo-PrP. We investigated the effect of recombinant mouse PrP (residues 23-231) refolded with nickel, manganese, copper, and a variant devoid of any metal ions, on tPA-catalyzed Plg activation. Using mutant PrP (H95A, H110A), we also investigated whether histidine residues outside the octarepeat region in PrP, which is known to bind tPA and Plg, are also involved in their binding. We demonstrated that apo-PrP is most effective at stimulating Plg. PrP refolded with nickle or manganese behave similar to apo-PrP, and PrP refolded with copper is least effective. The mutant form of PrP did not stimulate Plg activation to the same degree as apo-PrP indicating that the histidine residues outside the octarepeat region are also involved in binding to tPA and Plg. Similarly, the fibrillar form of PrP was ineffective at stimulating Plg activation. Our data suggest that upon loss of copper specifically, a structural rearrangement of PrP occurs that exposes binding sites to Plg and tPA, enhancing the stimulation of Plg activation.

Proteolytic functional amyloid digests pathogenic amyloid

Although amyloids are a well-known pathological structure, functional amyloids are beneficial. Functional amyloids can be engineered to cultivate desired functionality that can destroy malicious amyloids. However, not much is known...

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.

ODiNPred: comprehensive prediction of protein order and disorder

Structural disorder is widespread in eukaryotic proteins and is vital for their function in diverse biological processes. It is therefore highly desirable to be able to predict the degree of order and disorder from amino acid sequence. It is, however, notoriously difficult to predict the degree of local flexibility within structured domains and the presence and nuances of localized rigidity within intrinsically disordered regions. To identify such instances, we used the CheZOD database, which encompasses accurate, balanced, and continuous-valued quantification of protein (dis)order at amino acid resolution based on NMR chemical shifts. To computationally forecast the spectrum of protein disorder in the most comprehensive manner possible, we constructed the sequence-based protein order/disorder predictor ODiNPred, trained on an expanded version of CheZOD. ODiNPred applies a deep neural network comprising 157 unique sequence features to 1325 protein sequences together with the experimental NMR chemical shift data. Cross-validation for 117 protein sequences shows that ODiNPred better predicts the continuous variation in order along the protein sequence, suggesting that contemporary predictors are limited by the quality of training data. The inclusion of evolutionary features reduces the performance gap between ODiNPred and its peers, but analysis shows that it retains greater accuracy for the more challenging prediction of intermediate disorder.

Total degradation of extracellular amyloids by miniature artificial proteases

A miniaturized mimic of the active site of a protease, chymotrypsin, was linked to a target recognition unit to generate "Miniature Artificial Proteases" (mAPs). Time-resolved MALDI-TOF data analyses indicated that mAPs cleaved every amide bond between Lys16-Phe20 of the amyloid β fragment (Aβ12-21) and Aβ1-40, resulting in inhibition of fibrillization and disruption of the preformed amyloid. Such a platform may offer not only new therapeutic options against various amyloidoses but also novel routes for the selective knockdown of specific proteins.

Structural Impact of Phosphorylation and Dielectric Constant Variation on Synaptotagmin's IDR

We used time-resolved Förster resonance energy transfer, circular dichroism, and molecular dynamics simulation to investigate the structural dependence of synaptotagmin 1's intrinsically disordered region (IDR) on phosphorylation and dielectric constant. We found that a peptide corresponding to the full-length IDR sequence, a ∼60-residue strong polyampholyte, can sample structurally collapsed states in aqueous solution, consistent with its κ-predicted behavior, where κ is a sequence-dependent parameter that is used to predict IDR compaction. In implicit solvent simulations of this same sequence, lowering the dielectric constant to more closely mimic the environment near a lipid bilayer surface promoted further sampling of collapsed structures. We then examined the structural tendencies of central region residues of the IDR in isolation. We found that the exocytosis-modulating phosphorylation of Thr¹¹² disrupts a local disorder-to-order transition induced by trifluoroethanol/water mixtures that decrease the solution dielectric constant and stabilize helical structure. Implicit solvent simulations on these same central region residues testing the impact of dielectric constant alone converge on a similar result, showing that helical structure is formed with higher probability at a reduced dielectric. In these helical conformers, lysine-aspartic acid salt bridges contribute to stabilization of transient secondary structure. In contrast, phosphorylation results in formation of salt bridges unsuitable for helix formation. Collectively, these results suggest a model in which phosphorylation and compaction of the IDR sequence regulate structural transitions that in turn modulate neuronal exocytosis.

Quantification of Small Molecule-Protein Interactions using FRET between Tryptophan and the Pacific Blue Fluorophore

We report a new method to quantify the affinity of small molecules for proteins. This method is based on Förster resonance energy transfer (FRET) between endogenous tryptophan (Trp) residues and the coumarin-derived fluorophore Pacific Blue (PB). Tryptophan residues are frequently found in proteins near ligand-binding sites, making this approach potentially applicable to a wide range of systems. To improve access to PB, we developed a scalable multigram synthesis of this fluorophore, starting with inexpensive 2,3,4,5-tetrafluorobenzoic acid. This route was used to synthesize fluorescent derivatives of biotin, as well as lower affinity thiobiotin, iminobiotin, and imidazolidinethione analogues that bind the protein streptavidin. Compared with previously published FRET acceptors for tryptophan, PB proved to be superior in both sensitivity and efficiency. These unique properties of PB enabled direct quantification of dissociation constants (K_d) as well as competitive inhibition constants (K_i) in the micromolar to nanomolar range. In comparison to analogous binding studies using fluorescence polarization, fluorescence quenching, or fluorescence enhancement, affinities determined using Trp-FRET were more precise and accurate as validated using independent isothermal titration calorimetry studies. FRET between tryptophan and PB represents a new tool for the characterization of protein-ligand complexes.

The Emergence of NMDA Receptor Metabotropic Function: Insights from Imaging

The NMDA receptor (R) participates in many important physiological and pathological processes. For example, its activation is required for both long-term potentiation (LTP) and long-term depression (LTD) of synaptic transmission, cellular models of learning and memory. Furthermore, it may play a role in the actions of amyloid-beta on synapses as well as in the signaling leading to cell death following stroke. Until recently, these processes were thought to be mediated by ion-flux through the receptor. Using a combination of imaging and electrophysiological approaches, ion-flux independent functions of the NMDAR were recently examined. In this review, we will discuss the role of metabotropic NMDAR function in LTD and synaptic dysfunction.

Synthesis of double-fluorescent labeled prion protein for FRET analysis

An abnormal form of prion protein (PrP) is considered to be the pathogen in prion diseases. However, the structural details of this abnormal form are not known. To characterize the non-native structure of PrP, we synthesized position-specific double-fluorescent labeled PrP for a fluorescence resonance energy transfer (FRET) experiment. Using FRET, we observed a conformational change in the labeled PrP associated with amyloid fibril formation. The FRET analysis indicated that the distance between fluorescent labeled N- and C-terminal sites of PrP increased upon the formation of amyloid fibrils compared with that of the native state. This approach using FRET analysis is useful for elucidating the structure of abnormal PrP.

Detailed microscopic unfolding pathways of an α-helix and a β-hairpin: direct observation and molecular dynamics

We present a combined experimental and computational study of unfolding pathways of a model 21-residue α-helical heteropeptide (W1H5-21) and a 16-residue β-hairpin (GB41-56). Experimentally, we measured fluorescence energy transfer efficiency as a function of temperature, employing natural tryptophans as donors and dansylated lysines as acceptors. Secondary structural analysis was performed with circular dichroism and Fourier transform infrared spectroscopy. Our studies present markedly different unfolding pathways of the two elementary secondary structural elements. During thermal denaturation, the helical peptide exhibits an initial decrease in length, followed by an increase, while the hairpin undergoes a systematic increase in length. In the complementary computational part of the project, we performed microsecond length replica-exchange molecular dynamics simulations of the peptides in explicit solvent, yielding a detailed microscopic picture of the unfolding processes. For the α-helical peptide, we found a large heterogeneous population of intermediates that are primarily frayed single helices or helix-turn-helix motifs. Unfolding starts at the termini and proceeds through a stable helical region in the interior of the peptide but shifted off-center toward the C-terminus. The simulations explain the experimentally observed non-monotonic variation of helix length with temperature as due primarily to the presence of frayed-end single-helix intermediate structures. For the β-hairpin peptide, our simulations indicate that folding is initiated at the turn, followed by formation of the hairpin in zipper-like fashion, with Cα···Cα contacts propagating from the turn to termini and hairpin hydrogen bonds forming in parallel with these contacts. In the early stages of hairpin formation, the hydrophobic side-chain contacts are only partly populated. Intermediate structures with low numbers of β-hairpin hydrogen bonds have very low populations. This is in accord with the "broken zipper" model of Scheraga. The monotonic increase in length with temperature may be explained by the zipper-like breaking of the hairpin hydrogen bonds and backbone contacts.

A facile route to tailoring peptide-stabilized gold nanoparticles using glutathione as a synthon

The preparation of gold nanoparticles (AuNPs) of high purity and stability remains a major challenge for biological applications. This paper reports a simple synthetic strategy to prepare water-soluble peptide-stabilized AuNPs. Reduced glutathione, a natural tripeptide, was used as a synthon for the growth of two peptide chains directly on the AuNP surface. Both nonpolar (tryptophan and methionine) and polar basic (histidine and dansylated arginine) amino acids were conjugated to the GSH-capped AuNPs. Ultracentrifugation concentrators with polyethersulfone (PES) membranes were used to purify precursor materials in each stage of the multi-step synthesis to minimize side reactions. Thin layer chromatography, transmission electron microscopy, UV-Visible, 1H-NMR, and fluorescence spectroscopies demonstrated that ultracentrifugation produces high purity AuNPs, with narrow polydispersity, and minimal aggregation. More importantly, it allows for more control over the composition of the final ligand structure. Studies under conditions of varying pH and ionic strength revealed that peptide length, charge, and hydrophobicity influence the stability as well as solubility of the peptide-capped AuNPs. The synthetic and purification strategies used provide a facile route for developing a library of tailored biocompatible peptide-stabilized AuNPs for biomedical applications.

A dimeric equilibrium intermediate nucleates Drp1 reassembly on mitochondrial membranes for fission

Drp1 catalyzes mitochondrial division, but the mechanisms remain elusive. The mitochondrial lipid cardiolipin stimulates Drp1 activity and supports membrane constriction. In addition, Drp1 populates two polymeric states that equilibrate via a dimeric intermediate. Dimers nucleate Drp1 reassembly on mitochondria for fission.

The GTPase dynamin-related protein 1 (Drp1) catalyzes mitochondrial division, but the mechanisms remain poorly understood. Much of what is attributed to Drp1’s mechanism of action in mitochondrial membrane fission parallels that of prototypical dynamin in endocytic vesicle scission. Unlike the case for dynamin, however, no lipid target for Drp1 activation at the mitochondria has been identified. In addition, the oligomerization properties of Drp1 have not been well established. We show that the mitochondria-specific lipid cardiolipin is a potent stimulator of Drp1 GTPase activity, as well as of membrane tubulation. We establish further that under physiological conditions, Drp1 coexists as two morphologically distinct polymeric species, one nucleotide bound in solution and the other membrane associated, which equilibrate via a dimeric assembly intermediate. With two mutations, C300A and C505A, that shift Drp1 polymerization equilibria in opposite directions, we demonstrate that dimers, and not multimers, potentiate the reassembly and reorganization of Drp1 for mitochondrial membrane remodeling both in vitro and in vivo.

Experimental and computational study of BODIPY dye-labeled cavitand dynamics

Understanding the distance distribution and dynamics between moieties attached to the walls of a resorcin[4]arene cavitand, which is switchable between an expanded kite and a contracted vase form, might enable the use of this molecular system for the study of fundamental distance-dependent interactions. Toward this goal, a combined experimental and molecular dynamics (MD) simulation study on donor/acceptor borondipyrromethene (BODIPY) dye-labeled cavitands present in the vase and kite forms was performed. Direct comparison between anisotropy decays calculated from MD simulations with experimental fluorescence anisotropy data showed excellent agreement, indicating that the simulations provide an accurate representation of the dynamics of the system. Distance distributions between the BODIPY dyes were established by comparing time-resolved Förster resonance energy transfer experiments and MD simulations. Fluorescence intensity decay curves emulated on the basis of the MD trajectories showed good agreement with the experimental data, suggesting that the simulations present an accurate picture of the distance distributions and dynamics in this molecular system and provide an important tool for understanding the behavior of extended molecular systems and designing future applications.

Alternate pleckstrin homology domain orientations regulate dynamin-catalyzed membrane fission

The isolated dynamin PH domain is an assembly-independent sensor of membrane curvature but not a curvature generator. In full-length dynamin, the PH alternates between two different orientations on the membrane surface during the GTP hydrolysis cycle, causing dramatic fluctuations in the diameter of dynamin polymers.

The self-assembling GTPase dynamin catalyzes endocytic vesicle scission via membrane insertion of its pleckstrin homology (PH) domain. However, the molecular mechanisms underlying PH domain–dependent membrane fission remain obscure. Membrane-curvature–sensing and membrane-curvature–generating properties have been attributed, but it remains to be seen whether the PH domain is involved in either process independent of dynamin self-assembly. Here, using multiple fluorescence spectroscopic and microscopic techniques, we demonstrate that the isolated PH domain does not act to bend membranes but instead senses high membrane curvature through hydrophobic insertion into the membrane bilayer. Furthermore, we use a complementary set of short- and long-distance Förster resonance energy transfer approaches to distinguish PH-domain orientation from proximity at the membrane surface in full-length dynamin. We reveal, in addition to the GTP-sensitive “hydrophobic mode,” the presence of an alternate, GTP-insensitive “electrostatic mode” of PH domain–membrane interactions that retains dynamin on the membrane surface during the GTP hydrolysis cycle. Stabilization of this alternate orientation produces dramatic variations in the morphology of membrane-bound dynamin spirals, indicating that the PH domain regulates membrane fission through the control of dynamin polymer dynamics.

Minimalist Approaches to Protein Labelling: Getting the Most Fluorescent Bang for Your Steric Buck

Fluorescence methods allow one to monitor protein conformational changes, protein–protein associations, and proteolysis in real time, at the single molecule level and in living cells. The information gained in such experiments is a function of the spectroscopic techniques used and the strategic placement of fluorophore labels within the protein structure. There is often a trade-off between size and utility for fluorophores, whereby large size can be disruptive to the protein’s fold or function, but valuable characteristics, such as visible wavelength absorption and emission or brightness, require sizable chromophores. Three major types of fluorophore readouts are commonly used: (1) Förster resonance energy transfer (FRET); (2) photoinduced electron transfer (PET); and (3) environmental sensitivity. This review focuses on those probes small enough to be incorporated into proteins during ribosomal translation, which allows the probes to be placed on the interiors of proteins as they are folded during synthesis. The most broadly useful method for doing so is site-specific unnatural amino acid (UAA) mutagenesis. We discuss the use of UAA probes in applications relying on FRET, PET, and environmental sensitivity. We also briefly review other methods of protein labelling and compare their relative merits to UAA mutagenesis. Finally, we discuss small probes that have thus far been used only in synthetic peptides, but which have unusual value and may be candidates for incorporation using UAA methods.

Ensemble and single-molecule detected time-resolved FRET methods in studies of protein conformations and dynamics

Most proteins are nanomachines that are selected to execute specific functions and therefore should have some degree of flexibility. The driving force that excites specific motions of domains and smaller chain elements is the thermal fluctuations of the solvent bath which are channeled to selected modes of motions by the structural constraints. Consequently characterization of the ensembles of conformers of proteins and their dynamics should be expressed in statistical terms, i.e., determination of probability distributions of the various conformers. This can be achieved by measurements of time-resolved dynamic non-radiative excitation energy transfer (trFRET) within ensembles of site specifically labeled protein molecules. Distributions of intramolecular segmental end-to-end distances and their fast fluctuations can be determined, and fast and slow conformational transitions within selected sections of the molecule can be monitored and analyzed. Both ensemble and single-molecule detection methods can be applied for data collection. In combination with synchronization methods, time-resolved FRET was also used for studies of fast conformational transitions, in particular the folding/unfolding transitions.

Förster resonance energy transfer microscopy and spectroscopy for localizing protein-protein interactions in living cells

The fundamental theory of Förster resonance energy transfer (FRET) was established in the 1940s. Its great power was only realized in the past 20 years after different techniques were developed and applied to biological experiments. This success was made possible by the availability of suitable fluorescent probes, advanced optics, detectors, microscopy instrumentation, and analytical tools. Combined with state-of-the-art microscopy and spectroscopy, FRET imaging allows scientists to study a variety of phenomena that produce changes in molecular proximity, thereby leading to many significant findings in the life sciences. In this review, we outline various FRET imaging techniques and their strengths and limitations; we also provide a biological model to demonstrate how to investigate protein-protein interactions in living cells using both intensity- and fluorescence lifetime-based FRET microscopy methods.

On the foldability of tryptophan-containing tetra- and pentapeptides: an exhaustive molecular dynamics study

Short peptides serve as minimal model systems to decipher the determinants of foldability due to their simplicity arising from their smaller size, their ability to echo protein-like structural characteristics, and their direct implication in force field validation. Here, we describe an effort to identify small peptides that can still form stable structures in aqueous solutions. We followed the in silico folding of a selected set of 8640 tryptophan-containing tetra- and pentapeptides through 15 210 molecular dynamics simulations amounting to a total of 272.46 μs using explicit representation of the solute and full treatment of the electrostatics. The evaluation and sorting of peptides is achieved through scoring functions, which include terms based on interatomic vector distances, atomic fluctuations, and rmsd matrices between successive frames of a trajectory. Highly scored peptides are studied further via successive simulation rounds of increasing simulation length and using different empirical force fields. Our method suggested only a handful of peptides with strong foldability prognosis. The discrepancies between the predictions of the various force fields for such short sequences are also extensively discussed. We conclude that the vast majority of such short peptides do not adopt significantly stable structures in water solutions, at least based on our computational predictions. The present work can be utilized in the rational design and engineering of bioactive peptides with desired molecular properties.

Ensemble FRET methods in studies of intrinsically disordered proteins

The main structural characteristic of intrinsically disordered proteins (IDPs) or intrinsically disordered regions of globular proteins is that they exist as ensembles of multiple conformers which can continuously interconvert, and at times, form ensembles of a more restricted number of conformers. Characterization of the disordered state and transitions to partially or fully ordered states of such ensembles must be expressed in statistical terms, i.e., determination of probability distributions of the various conformers. This can be achieved by measurements of time-resolved dynamic non-radiative excitation energy transfer within ensembles of site-specifically labeled IDP molecules. Distributions of intramolecular segmental end-to-end distances and their fast fluctuations can be determined and fast and slow conformational transitions within selected sections of the molecule can be monitored and analyzed.

Aggregation kinetics of interrupted polyglutamine peptides

Abnormally expanded polyglutamine domains are associated with at least nine neurodegenerative diseases, including Huntington's disease. Expansion of the glutamine region facilitates aggregation of the impacted protein, and aggregation has been linked to neurotoxicity. Studies of synthetic peptides have contributed substantially to our understanding of the mechanism of aggregation because the underlying biophysics of polyglutamine-mediated association can be probed independent of their context within a larger protein. In this report, interrupting residues were inserted into polyglutamine peptides (Q20), and the impact on conformational and aggregation properties was examined. A peptide with two alanine residues formed laterally aligned fibrillar aggregates that were similar to the uninterrupted Q20 peptide. Insertion of two proline residues resulted in soluble, nonfibrillar aggregates, which did not mature into insoluble aggregates. In contrast, insertion of a β-turn template (D)PG rapidly accelerated aggregation and resulted in a fibrillar aggregate morphology with little lateral alignment between fibrils. These results are interpreted to indicate that (a) long-range nonspecific interactions lead to the formation of soluble oligomers, while maturation of oligomers into fibrils requires conformational conversion and (b) that soluble oligomers dynamically interact with each other, while insoluble aggregates are relatively inert. Kinetic analysis revealed that the increase in aggregation caused by the (D)PG insert is inconsistent with the nucleation-elongation mechanism of aggregation featuring a monomeric β-sheet nucleus. Rather, the data support a mechanism of polyglutamine aggregation by which monomers associate into soluble oligomers, which then undergo slow structural rearrangement to form sedimentable aggregates.

Insights into intramolecular Trp and His side-chain orientation and stereospecific π interactions surrounding metal centers: an investigation using protein metal-site mimicry in solution

Metal-binding scaffolds incorporating a Trp/His-paired epitope are instrumental in giving novel insights into the physicochemical basis of functional and mechanistic versatility conferred by the Trp-His interplay at a metal site. Herein, by coupling biometal site mimicry and (1)H and (13)C NMR spectroscopy experiments, modular constructs EDTA-(L-Trp, L-His) (EWH; EDTA=ethylenediamino tetraacetic acid) and DTPA-(L-Trp, L-His) (DWH; DTPA=diethylenetriamino pentaacetic acid) were employed to dissect the static and transient physicochemical properties of hydrophobic/hydrophilic aromatic interactive modes surrounding biometal centers. The binding feature and identities of the stoichiometric metal-bound complexes in solution were investigated by using (1)H and (13)C NMR spectroscopy, which facilitated a cross-validation of the carboxylate, amide oxygen, and tertiary amino groups as the primary ligands and indole as the secondary ligand, with the imidazole (Im) N3 nitrogen being weakly bound to metals such as Ca(2+) owing to a multivalency effect. Surrounding the metal centers, the stereospecific orientation of aromatic rings in the diastereoisomerism is interpreted with the Ca(2+)-EWH complex. With respect to perturbed Trp side-chain rotamer heterogeneity, drastically restricted Trp side-chain flexibility and thus a dynamically constrained rotamer interconversion due to π interactions is evident from the site-selective (13)C NMR spectroscopic signal broadening of the Trp indolyl C3 atom. Furthermore, effects of Trp side-chain fluctuation on indole/Im orientation were the subject of a 2D NMR spectroscopy study by using the Ca(2+)-bound state; a C-H2(indolyl)/C-H5(Im(+)) connectivity observed in the NOESY spectra captured direct evidence that the N-H1 of the Ca(2+)-Im(+) unit interacted with the pyrrole ring of the indole unit in Ca(2+)-bound EWH but not in DWH, which is assignable to a moderately static, anomalous, T-shaped, interplanar π(+)-π stacking alignment. Nevertheless, a comparative (13)C NMR spectroscopy study of the two homologous scaffolds revealed that the overall response of the indole unit arises predominantly from global attractions between the indole ring and the entire positively charged first coordination sphere. The study thus demonstrates the coordination-sphere/geometry dependence of the Trp/His side-chain interplay, and established that π interactions allow (13)C NMR spectroscopy to offer a new window for investigating Trp rotamer heterogeneity near metal-binding centers.

An improved open-channel structure of MscL determined from FRET confocal microscopy and simulation

Mechanosensitive channels act as molecular transducers of mechanical force exerted on the membrane of living cells by opening in response to membrane bilayer deformations occurring in physiological processes such as touch, hearing, blood pressure regulation, and osmoregulation. Here, we determine the likely structure of the open state of the mechanosensitive channel of large conductance using a combination of patch clamp, fluorescence resonance energy transfer (FRET) spectroscopy, data from previous electron paramagnetic resonance experiments, and molecular and Brownian dynamics simulations. We show that structural rearrangements of the protein can be measured in similar conditions as patch clamp recordings while controlling the state of the pore in its natural lipid environment by modifying the lateral pressure distribution via the lipid bilayer. Transition to the open state is less dramatic than previously proposed, while the N terminus remains anchored at the surface of the membrane where it can either guide the tilt of or directly translate membrane tension to the conformation of the pore-lining helix. Combining FRET data obtained in physiological conditions with simulations is likely to be of great value for studying conformational changes in a range of multimeric membrane proteins.

Molecular simulations of protein disorder

Protein disorder is abundant in proteomes throughout all kingdoms of life and serves many biologically important roles. Disordered states of proteins are challenging to study experimentally due to their structural heterogeneity and tendency to aggregate. Computer simulations, which are not impeded by these properties, have recently emerged as a useful tool to characterize the conformational ensembles of intrinsically disordered proteins. In this review, we provide a survey of computational studies of protein disorder with an emphasis on the interdisciplinary nature of these studies. The application of simulation techniques to the study of disordered states is described in the context of experimental and bioinformatics approaches. Experimental data can be incorporated into simulations, and simulations can provide predictions for experiment. In this way, simulations have been integrated into the existing methodologies for the study of disordered state ensembles. We provide recent examples of simulations of disordered states from the literature and our own work. Throughout the review, we emphasize important predictions and biophysical understanding made possible through the use of simulations. This review is intended as both an overview and a guide for structural biologists and theoretical biophysicists seeking accurate, atomic-level descriptions of disordered state ensembles.

Orientational averaging of dye molecules attached to proteins in Förster resonance energy transfer measurements: insights from a simulation study

Förster resonance energy transfer is an increasingly popular method for studying protein folding at single molecule resolution. By attaching dye molecules to particular residues in a protein molecule and measuring the energy transfer to the acceptor dye on excitation of the donor dye, information about the separation of the dyes can be obtained. Here we use an atomistic coarse-grained molecular model of the protein and dyes to look at the assumption that the dyes rotate freely during the donor decay time. We find that although complete orientational averaging does not always occur, the consequences of this are not extreme. Even in the native state, the errors in efficiency, which result from incorrectly assuming kappa2=2/3, are smaller than the typical experimental error of an efficiency measurement. The orientational freedom of the dyes originates both from the dynamics of the linker and dye molecules and also from the movements of the protein chain itself. In the unfolded state, the movements of the protein chain are sufficient to provide complete, or almost complete, orientational averaging within the donor lifetime. Increasing the rigidity of the dyes therefore has only a very small effect on the measured efficiencies in the unfolded state. In the native state the contribution of the linker and dye dynamics to orientational averaging is larger; nevertheless increasing the rigidity still has only a small effect on the measured efficiency.

Examining polyglutamine peptide length: a connection between collapsed conformations and increased aggregation

Abnormally expanded polyglutamine domains in proteins are associated with several neurodegenerative diseases, of which the best known is Huntington's. Expansion of the polyglutamine domain facilitates aggregation of the affected protein, and several studies directly link aggregation to neurotoxicity. The age of onset of disease is inversely correlated with the length of the polyglutamine domain; this correlation motivates an examination of the role of the length of the domain on aggregation. In this investigation, peptides containing 8 to 24 glutamines were synthesized, and their conformational and aggregation properties were examined. All peptides lacked secondary structure. Fluorescence resonance energy transfer studies revealed that the peptides became increasingly collapsed as the number of glutamine residues increased. The effective persistence length was estimated to decrease from approximately 11 to approximately 7 A as the number of glutamines increased from 8 to 24. A comparison of our data with theoretical results suggests that phosphate-buffered saline is a good solvent for Q8 and Q12, a theta solvent for Q16, and a poor solvent for Q20 and Q24. By dynamic light scattering, we observed that Q16, Q20, and Q24, but not Q8 or Q12, immediately formed soluble aggregates upon dilution into phosphate-buffered saline at 37 degrees C. Thus, Q16 stands at the transition point between good and poor solvent and between stable and aggregation-prone peptide. Examination of aggregates by transmission electron microscopy, along with kinetic assays for sedimentation, provided evidence indicating that soluble aggregates mature into sedimentable aggregates. Together, the data support a mechanism of aggregation in which monomer collapse is accompanied by formation of soluble oligomers; these soluble species lack regular secondary structure but appear morphologically similar to the sedimentable aggregates into which they eventually mature.

Molecular dynamics simulations reveal that AEDANS is an inert fluorescent probe for the study of membrane proteins

Computer simulations were carried out of a number of AEDANS-labeled single cysteine mutants of a small reference membrane protein, M13 major coat protein, covering 60% of its primary sequence. M13 major coat protein is a single membrane-spanning, α-helical membrane protein with a relatively large water-exposed region in the N-terminus. In 10-ns molecular dynamics simulations, we analyze the behavior of the AEDANS label and the native tryptophan, which were used as acceptor and donor in previous FRET experiments. The results indicate that AEDANS is a relatively inert environmental probe that can move unhindered through the lipid membrane when attached to a membrane protein.

Simulation of structure, orientation, and energy transfer between AlexaFluor molecules attached to MscL

Measurements of time-resolved fluorescence anisotropy and fluorescence resonance energy transfer are finding many applications in the study of biological macromolecules as they enable structural properties of the host molecules to be determined in their natural environment. A difficulty in interpreting these experiments is that they both require knowledge of the relative orientation of the fluorophores, a property that is almost impossible to measure. Here we conduct simulations of AlexaFluor488 and AlexaFluor568 attached to two sites on the membrane channel MscL to provide an alternative mechanism for determining the likely configurations and orientational freedom of the fluorophores, as well as the most likely value of the orientation factor kappa(2) for energy transfer between them. The fluorophores are relatively mobile, and are found to be more so when immersed in bulk water than when they interact with the lipid membrane. The fluorophores never insert deeply into the lipid, despite their hydrophobic linkers and aromatic headgroup structures. Properties such as the fluorescence anisotropy decay can be predicted from simulations of the fluorophores in bulk water that closely match experimental data. In contrast, when the fluorophores were attached to the large MscL protein it was difficult to sample all the possible configurations of the fluorophores due to the computational time required. While this approach is likely to provide useful data on solvent-accessible fluorophores attached to small proteins, simulations lasting >50 ns or the use of biasing forces are required to accurately predict orientation factors for use in energy transfer experiments on larger membrane-bound proteins.

Expansion of the octarepeat domain alters the misfolding pathway but not the folding pathway of the prion protein

A misfolded conformation of the prion protein (PrP), PrP (Sc), is the essential component of prions, the infectious agents that cause transmissible neurodegenerative diseases. Insertional mutations that lead to an increase in the number of octarepeats (ORs) in PrP are linked to familial human prion disease. In this study, we investigated how expansion of the OR domain causes PrP to favor a prion-like conformation. Therefore, we compared the conformational and aggregation modulating properties of wild-type versus expanded OR domains, either as a fusion construct with the protein G B1 domain (GB1-OR) or as an integral part of full-length mouse PrP (MoPrP). Using circular dichroism spectroscopy, we first demonstrated that ORs are not unfolded but exist as an ensemble of three distinct conformers: polyproline helix-like, beta-turn, and "Trp-related". Domain expansion had little effect on the conformation of GB1-OR fusion proteins. When part of MoPrP however, OR domain expansion changed PrP's folding landscape, not by hampering the production of native alpha-helical monomers but by greatly reducing the propensity to form amyloid and by altering the assembly of misfolded, beta-rich aggregates. These features may relate to subtle pH-dependent conformational differences between wild-type and mutant monomers. In conclusion, we propose that PrP insertional mutations are pathogenic because they enhance specific misfolding pathways of PrP rather than by undermining native folding. This idea was supported by a trial bioassay in transgenic mice overexpressing wild-type MoPrP, where intracerebral injection of recombinant MoPrP with an expanded OR domain but not wild-type MoPrP caused prion disease.

Molecular dynamics simulations of two tandem octarepeats from the mammalian prion protein: fully Cu2+-bound and metal-free forms

Molecular dynamics simulations have been conducted on a model fragment (Ac-PHGGGWGQPHGGGW-NH(2)) of the prion protein octarepeat domain, both in the Cu(2+)-bound and metal-free forms. The copper-bound models are based on the consensus structure of the core Cu(2+)-binding site of an individual octarepeat, relevant to the fully Cu(2+)-occupied prion protein octarepeat region. The model peptides contain Cu(2+) bound through a His imidazole ring and two deprotonated amide N-atoms in the peptide backbone supplied by the following two Gly residues. Both the copper-bound and metal-free models have been simulated with the OPLS all-atom force field with the GROMACS molecular dynamics package. These simulations, with two tandem copper-binding sites, represent the minimum model necessary to observe potential structuring between the copper-binding sites in the octarepeat region. The GWGQ residues constitute a flexible linker region that predominantly adopts a turn, serving to bring adjacent His residues into close proximity. The consequent formation of stable structures demonstrates that the copper-bound octarepeat region allows the copper-coordinating sites to come into van der Waals contact, packing into particular orientations to further stabilize the bend in the GWGQ linker region.

Fretting about FRET: correlation between kappa and R

Molecular dynamics simulations were used to examine the structural dynamics of two fluorescent probes attached to a typical protein, hen egg-white lysozyme (HEWL). The donor probe (D) was attached via a succinimide group, consistent with the commonly-used maleimide conjugation chemistry, and the acceptor probe (A) was bound into the protein as occurs naturally for HEWL and the dye Eosin Y. The <kappa(2)> is found to deviate significantly from the theoretical value and high correlation between the orientation factor kappa and the distance R is observed. The correlation is quantified using several possible fixed A orientations and correlation as high as 0.80 is found between kappa and R and as high as 0.68 between kappa(2) and R. The presence of this correlation highlights the fact that essentially all fluorescence-detected resonance energy transfer studies have assumed that kappa and R are independent--an assumption that is clearly not justified in the system studied here. The correlation results in the quantities <kappa(2)R(-)(6)> and <kappa(2)> < R(-)(6)> differing by a factor of 1.6. The observed correlation between kappa and R is caused by the succinimide linkage between the D and HEWL, which is found to be relatively inflexible.

Equilibrium unfolding of the poly(glutamic acid)20 helix

The equilibrium structural ensemble of a 20-residue polyglutamic acid peptide (E(20)) was studied with FRET, circular dichroism, and molecular dynamics (MD) simulations. A FRET donor, o-aminobenzamide, and acceptor, 3-nitrotyrosine, were introduced at the N- and C-termini, respectively. Circular dichroism, steady state FRET, and time-resolved FRET measurements were employed to characterize the fraction helix and end-to-end distance under different pH conditions: pH 4 (60% alpha-helix), pH 6 (0% alpha-helix), and pH 9 (0% alpha-helix). At pH 4, the end-to-end distance was measured at 24 A and determined to be considerably less than the 31 A predicted for an alpha-helix of the same length. At pH 6 and 9, the end-to-end distance was measured at > 31 and 39 A respectively, both which are determined to be considerably greater than the 27 A predicted for a freely jointed random coil of the same length. To better understand the physical forces underlying the unusual helix-coil transition in this peptide, three theoretical MD models of E(20) were constructed: (1) a pure alpha-helix, (2) an alpha-helix with equivalent attractive intramolecular contacts, and (3) a weak alpha-helix with termini-weighted intramolecular contacts ("sticky ends"). Using MD simulations, the bent helix structure calculated from Model 3 was found to be the closest in agreement with the experimental data.

FRET study of membrane proteins: determination of the tilt and orientation of the N-terminal domain of M13 major coat protein

A formalism for membrane protein structure determination was developed. This method is based on steady-state FRET data and information about the position of the fluorescence maxima on site-directed fluorescent labeled proteins in combination with global data analysis utilizing simulation-based fitting. The methodology was applied to determine the structural properties of the N-terminal domain of the major coat protein from bacteriophage M13 reconstituted into unilamellar DOPC/DOPG (4:1 mol/mol) vesicles. For our purpose, the cysteine mutants A7C, A9C, N12C, S13C, Q15C, A16C, S17C, and A18C in the N-terminal domain of this protein were produced and specifically labeled with the fluorescence probe AEDANS. The energy transfer data from the natural Trp-26 to AEDANS were analyzed assuming a two-helix protein model. Furthermore, the polarity Stokes shift of the AEDANS fluorescence maxima is taken into account. As a result the orientation and tilt of the N-terminal protein domain with respect to the bilayer interface were obtained, showing for the first time, to our knowledge, an overall alpha-helical protein conformation from amino acid residues 12-46, close to the protein conformation in the intact phage.

FRET study of membrane proteins: simulation-based fitting for analysis of membrane protein embedment and association

A new formalism for the simultaneous determination of the membrane embedment and aggregation of membrane proteins is developed. This method is based on steady-state Förster (or fluorescence) resonance energy transfer (FRET) experiments on site-directed fluorescence labeled proteins in combination with global data analysis utilizing simulation-based fitting. The simulation of FRET was validated by a comparison with a known analytical solution for energy transfer in idealized membrane systems. The applicability of the simulation-based fitting approach was verified on simulated FRET data and then applied to determine the structural properties of the well-known major coat protein from bacteriophage M13 reconstituted into unilamellar DOPC/DOPG (4:1 mol/mol) vesicles. For our purpose, the cysteine mutants Y24C, G38C, and T46C of this protein were produced and specifically labeled with the fluorescence label AEDANS. The energy transfer data from the natural tryptophan at position 26, which is used as a donor, to AEDANS were analyzed assuming a helix model for the transmembrane domain of the protein. As a result of the FRET data analysis, the topology and bilayer embedment of this domain were quantitatively characterized. The resulting tilt of the transmembrane helix of the protein is 18 +/- 2 degrees. The tryptophan is located at a distance of 8.5 +/- 0.5 A from the membrane center. No specific aggregation of the protein was found. The methodology developed here is not limited to M13 major coat protein and can be used in principle to study the bilayer embedment of any small protein with a single transmembrane domain.

Simulating FRET from Tryptophan: Is the Rotamer Model Correct?

We present a computational model study designed to simulate the results of time-resolved fluorescence spectra of tryptophan in proteins. In such measurements, the occurrence of more than one fluorescence lifetime is generally attributed to the existence of several tryptophan rotamers and/or structural conformations of the protein structure. The protein system we chose for this initial study is the tetracycline repressor (TetR), an interesting model system for the investigation of the mechanisms of transcriptional regulation. Fluorescence resonance energy transfer (FRET) from tryptophan to tetracycline is frequently observed in complexes of the TetR with the antibiotic tetracycline. We use a combined classical/quantum mechanical approach to model the structure and the spectroscopic properties of the TetR-tetracycline complex. A classical molecular dynamics simulation provides input geometries for semiempirical quantum mechanical/molecular mechanical (QM/MM) single-point configuration interaction (CI) calculations, which are used to calculate tryptophan vertical absorption and fluorescence energies and intensities as well as relative FRET rate constants. These rate constants together with the Einstein coefficients for spontaneous emission and an assumed rate for nonradiative deactivation allow us to simulate fluorescence decay curves with and without FRET and for the entire ensemble as well as for individual rotamers. Our results indicate that the classical "rotamer model", used to explain the multiexponential fluorescence-decay curves of time-resolved tryptophan emission spectra, can be extended to systems with FRET acceptors present in the protein matrix but that the interpretation of the fitted lifetimes is different to that usually used.

Conformational Distribution of a 14-Residue Peptide in Solution: A Fluorescence Resonance Energy Transfer Study

We demonstrate here that a nitrile-derivatized phenylalanine residue, p-cyanophenylalanine (Phe(CN)), and tryptophan (Trp) constitute a novel donor-acceptor pair for fluorescence resonance energy transfer (FRET). The Förster distance of this FRET pair was determined to be approximately 16 A and hence is well suited for determining relatively short separation distances. To validate the applicability of this FRET pair in conformational studies, we studied the conformational heterogeneity of a 14-residue amphipathic peptide, Mastoparan X (MPx peptide), in water and 7 M urea solution as well as at different temperatures. Specifically, seven nitrile-derivatized mutants of the MPx peptide, each containing a Phe(CN) residue that replaces different positions along the peptide sequence (i.e., from position 5 to 11) and serves as a resonance energy donor to the native Trp residue at position 3, were studied spectroscopically. The FRET efficiencies obtained from these peptides allowed us to gain a global picture regarding the conformational distribution of the MPx peptide in different environments. Our results suggest that the MPx molecules exist in water as an ensemble of rather compact conformations, with a radius of gyration of approximately 4.2 A, whereas in 7 M urea the radius of gyration increases to approximately 6.5 A, indicating that the peptide conformations become more extended under this condition. However, we found that temperature had only a negligible effect on the size of the MPx peptide, underlining the difference between the thermally and chemically denatured states of polypeptides. The application of the Gaussian chain or the wormlike chain model allowed us to further obtain the probability distribution function of the separation distance between any two residues along the peptide sequence. We found that the effective bond length of the MPx peptide, obtained by using the Gaussian chain model, is 2.78 A in water and 4.28 A in 7 M urea.

Imaging protein molecules using FRET and FLIM microscopy

Förster (or fluorescence) resonance energy transfer (FRET) and fluorescence lifetime imaging (FLIM) have moved center stage and are increasingly forming part of multifaceted imaging approaches. They are complementary methodologies that can be applied to advanced quantitative analyses. The widening application of FRET and FLIM has been driven by the availability of suitable fluorophores, increasingly sophisticated microscopy systems, methodologies to correct spectral bleed-through, and the ease with which FRET can be combined with other techniques. FRET and FLIM have recently found use in several applications: in the analysis of protein-protein interactions with high spatial and temporal specificity (e.g. clustering), in the study of conformational changes, in the analysis of binding sequences, and in applications such as high-throughput screening.