Abstract: Nucleic acids have been widely recognized as potential targets in drug discovery and aptamer selection. Quantifying the interactions between small molecules and nucleic acids is critical to discover lead compounds and design novel aptamers. Scoring function is normally employed to quantify the interactions in structure-based virtual screening. However, the predictive power of nucleic acid-ligand scoring functions is still a challenge compared to other types of biomolecular recognition. With the rapid growth of experimentally determined nucleic acid-ligand complex structures, in this work, we develop a knowledge-based scoring function of nucleic acid-ligand interactions, namely SPA-LN. SPA-LN is optimized by maximizing both the affinity and specificity of native complex structures. The development strategy is different from those of previous nucleic acid-ligand scoring functions which focus on the affinity only in the optimization. The native conformation is stabilized while non-native conformations are destabilized by our optimization, making the funnel-like binding energy landscape more biased toward the native state. The performance of SPA-LN validates the development strategy and provides a relatively more accurate way to score the nucleic acid-ligand interactions.

Abstract: Alzheimer's disease (AD) is one of the neurodegenerative disease characterized by progressive neuronal loss in the brain. Its two major hallmarks are extracellular senile plaques and intracellular neurofibrillary tangles (NFTs), formed by aggregation of amyloid β-42 (Aβ-42) and Tau protein respectively. Aβ-42 is a transmembrane protein, which is produced after the sequential action of β- and γ-secretases, thus obtained peptide is released extracellularly and gets deposited on the neuron forming senile plaques. NFTs are composed of microtubule-associated protein-Tau (MAPT). Tau protein's major function is to stabilize the microtubule that provides a track on which the cargo proteins are shuttled and the stabilized microtubule also maintains shape and integrity of the neuronal cell. Tau protein is subjected to various modifications such as phosphorylation, ubiquitination, glycation, acetylation, truncation, glycosylation, deamination, and oxidation; these modifications ultimately lead to its aggregation. Phosphorylation is the major modification and is extensively studied with respect to Tau protein. Tau protein, however, undergoes certain level of phosphorylation and dephosphorylation, which regulates its affinity for microtubule and ultimately leading to microtubule assembly and disassembly. Our main aim was to study the native state of longest isoform of Tau (hTau40WT-4R2N) and its shortest isoform, (hTau23WT-3R0N), at various temperatures such as 10, 25, and 37 °C. Raman spectroscopic results suggested that the proportion of random coils or unordered structure depends on the temperature of the protein environment. Upon increase in the temperature from 10 to 37 °C, the proportion of random coils or unordered structures increased in the case of hTau40WT. However, we did not find a significant effect of temperature on the structure of hTau23WT. This current approach enables one to analyze the global conformation of soluble Tau in solution.

Abstract: Little structural information is available so far on amyloid fibrils consisting of immunoglobulin light chains. It is not understood which features of the primary sequence of the protein result in fibril formation. We report here MAS solid-state NMR studies to identify the structured core of κ-type variable domain light chain fibrils. The core contains residues of the CDR2 and the β-strands D, E, F and G of the native immunoglobulin fold. The assigned core region of the fibril is distinct in comparison to the core identified in a previous solid-state NMR study on AL-09 by Piehl at. al, suggesting that VL fibrils can adopt different topologies. In addition, we investigated a soluble oligomeric intermediate state, previously termed the alternatively folded state (AFS), using NMR and FTIR spectroscopy. The NMR oligomer spectra display a high degree of similarity when compared to the fibril spectra, indicating a high structural similarity of the two aggregation states. Based on comparison to the native state NMR chemical shifts, we suggest that fibril formation via domain-swapping seems unlikely. Moreover, we used our results to test the quality of different amyloid prediction algorithms.

Abstract: Genetic, biochemical, and pharmacologic evidence supports the hypothesis that conformationally altered or misfolded protein states enable aggregation and cytotoxicity in the systemic amyloid diseases. Reversible structural fluctuations of natively folded proteins are involved in the aggregation of many degenerative disease associated proteins. Herein, we use antibody light chains (LCs) that form amyloid fibrils in AL amyloidosis to consider an alternative hypothesis of amyloidogenesis: that transient unfolding and incomplete extracellular refolding of secreted proteins can lead to metastable, alternatively folded states that are more susceptible to aggregation or to endoproteolysis that can release aggregation-prone fragments. Refolding of full-length λ6a LC dimers comprising an interchain disulfide bond from heat- or chaotrope-denatured ensembles in buffers yields the native dimeric state as well as alternatively folded dimers and aggregates. LC variants lacking an interchain disulfide bond appear to refold fully reversibly to the native state. The conformation of a backbone peptidyl-proline amide in the LC constant domain, which is cis in the native state, may determine whether the LC refolds back to the native state. A proline to alanine (P147A) LC variant, which cannot form the native cis-amide conformation, forms amyloid fibrils from the alternatively folded ensemble, whereas all the full-length λ6a LCs we have studied to date do not form amyloid under analogous conditions. P147A LC variants are susceptible to endoproteolysis by thrombin, enabling amyloidogenesis of the fragments released. Thus, non-native LC structural ensembles containing a tyrosine 146-proline 147 trans-amide bond can initiate and propagate amyloid formation, either directly or after aberrant endoproteolysis.

Abstract: Amyloid formation of natively folded proteins involves global and/or local unfolding of the native state to form aggregation-prone intermediates. Here we report solid-state NMR structural studies of amyloid derived from wild-type (WT) and more aggressive mutant forms of transthyretin (TTR) to investigate the structural changes associated with effective TTR aggregation. We employed selective 13C-labeling schemes to investigate structural features of β-structured core regions in amyloid states of WT and two mutant forms (V30M and L55P) of TTR. Analyses of the 13C-13C correlation solid-state NMR spectra revealed that WT TTR aggregates contain an amyloid core consisting of native-like CBEF and DAGH β-sheet structures and the mutant TTR amyloids adopt a similar amyloid core structure with native-like CBEF and AGH β-structures. However, the V30M mutant amyloid was shown to have a different DA β-structure. In addition, strand D is more disordered even in the native state of L55P TTR, indicating that the pathogen...

Abstract: For proteins with a single well-defined native state, protein 3Dstructure is a major determinant of sequence evolution. On the other hand, many proteins adopt multiple, distinct native structures under different conditions ("conformational switches"), yet the impact of such native state switching on protein evolution is not fully understood. Here, we performed a proteome-wide analysis of how protein structure impacts sequence evolution for protein conformational switches in Saccharomyces cerevisiae using pooled analysis of sites with similar packing or burial. We observed a strong linear relationship between residue evolutionary rate and residue burial for conformational switches. In addition, we found that conformational switches evolve significantly and consistently more slowly than proteins with a single native state, even after controlling for degree of residue burial or packing. Next, we focused on proteins that switch conformations upon molecular binding. We found that interfacial residues in these conformational switches evolve more slowly than interfacial residues in proteins with a single native state, and that the bound conformation is a better predictor for residue evolutionary rate than the unbound conformation. Our findings suggest that for conformational switches, the necessity to encode multiple distinct native structures under different conditions imposes strong evolutionary constraints on the entire protein, rather than just a few key residues. Our results provide new insights into the structure-evolution relationship of protein conformational switches.

Abstract: Although RNA interactions with K+ and Mg2+ have been studied extensively, much less is known about the third most abundant cation in bacterial cells, putrescine2+, and how RNA folding might be influenced by the three ions in combination. In a new approach, we have observed the competition between Mg2+ and putrescine2+ (in a background of K+) with native, partially unfolded and highly extended conformations of an adenine riboswitch aptamer. With the native state, putrescine2+ is a weak competitor when the ratio of the excess Mg2+ (which neutralizes phosphate charge) to RNA is very low, but becomes much more effective at replacing Mg2+ as the excess Mg2+ in the RNA ion atmosphere increases. Putrescine2+ is even more effective in competing Mg2+ from the extended conformation, independent of the Mg2+ excess. To account for these and other results, we propose that both ions closely approach the surface of RNA secondary structure, but the completely folded RNA tertiary structure develops small pockets of very negative electrostatic potential that are more accessible to the compact charge of Mg2+. The sensitivity of RNA folding to the combination of Mg2+ and putrescine2+ found in vivo depends on the architectures of both the unfolded and native conformations.

Abstract: Proteins work only if folded in their native state, but changes in temperature T and pressure P induce their unfolding. Therefore for each protein there is a stability region (SR) in the T-P thermodynamic plane outside which the biomolecule is denaturated. It is known that the extension and shape of the SR depend on i) the specific protein residue-residue interactions in the native state of the amino acids sequence and ii) the water properties at the hydration interface. Here we analyze by Monte Carlo simulations of different coarse-grained protein models in explicit water how changes in i) and ii) affect the SR. We show that the solvent properties ii) are essential to rationalize the SR shape at low T and high P and that our finding are robust with respect to parameter changes and with respect to different protein models. These results can help in developing new strategies for the design of novel synthetic biopolymers.

Abstract: Kallikrein, a physiologically vital serine protease, was investigated for its functional and conformational transitions during chemical (organic solvents, Gdn-HCl), thermal, and pH induced denaturation using biochemical and biophysical techniques and molecular dynamics (MD) simulations approach. The enzyme was exceptionally stable in isopropanol and ethanol showing 110% and 75% activity, respectively, after 96 h, showed moderate tolerance in acetonitrile (45% activity after 72 h) and much lower stability in methanol (40% activity after 24 h) (all the solvents [90% v/v]). Far UV CD and fluorescence spectra indicated apparent reduction in compactness of KLKp structure in isopropanol system. MD simulation studies of the enzyme in isopropanol revealed (1) minimal deviation of the structure from native state (2) marginal increase in radius of gyration and solvent accessible surface area (SASA) of the protein and the active site, and (3) loss of density barrier at the active site possibly leading to increased accessibility of substrate to catalytic triad as compared to methanol and acetonitrile. Although kallikrein was structurally stable up to 90 °C as indicated by secondary structure monitoring, it was functionally stable only up to 45 °C, implicating thermolabile active site geometry. In GdnHCl [1.0 M], 75% of the activity of KLKp was retained after incubation for 4 h, indicating its denaturant tolerance. A molten globule-like structure of KLKp formed at pH 1.0 was more thermostable and exhibited interesting structural transitions in organic solvents. The above results provide deeper understanding of functional and structural stability of the serine proteases at molecular level.

Abstract: MPT63, a major secreted protein from Mycobacterium tuberculosis, has been shown to have immunogenic properties and has been implicated in virulence. MPT63 is a β-sandwich protein containing 11 β strands and a very short stretch of 310 helix. The detailed experimental and computational study reported here investigates the equilibrium unfolding transition of MPT63. It is shown that in spite of being a complete β-sheet protein, MPT63 has a strong propensity towards helix structures in its early intermediates. Far UV-CD and FTIR spectra clearly suggest that the low pH intermediate of MTP63 has enhanced helical content, while fluorescence correlation spectroscopy (FCS) suggests a significant contraction. Molecular dynamics simulation complements the experimental results indicating that the unfolded state of MPT63 traverses through intermediate forms with increased helical characteristics. It is found that this early intermediate contains exposed hydrophobic surface, and is aggregation prone. Although MPT63 is a complete β-sheet protein in its native form, the present findings suggest that the secondary structure preferences of the local interactions in early folding pathway may not always follow the native conformation. Furthermore, the Gly25Ala mutant supports the proposed hypothesis by increasing the non-native helical propensity of the protein structure. This article is protected by copyright. All rights reserved.

Abstract: In order to directly observe the refolding kinetics from a partially misfolded state to a native state in the bottom of the protein-folding funnel, we used a "caging" strategy to trap the β-sheet structure of ubiquitin in a misfolded conformation. We used molecular dynamics simulation to generate the cage-induced, misfolded structure and compared the structure of the misfolded ubiquitin with native ubiquitin. Using laser flash irradiation, the cage can be cleaved from the misfolded structure within one nanosecond, and we monitored the refolding kinetics of ubiquitin from this misfolded state to the native state by photoacoustic calorimetry and photothermal beam deflection techniques on nanosecond to millisecond timescales. Our results showed two refolding events in this refolding process. The fast event is shorter than 20 ns and corresponds to the instant collapse of ubiquitin upon cage release initiated by laser irradiation. The slow event is ~60 μs, derived from a structural rearrangement in β-sheet refolding. The event lasts 10 times longer than the timescale of β-hairpin formation for short peptides as monitored by temperature jump, suggesting that rearrangement of a β-sheet structure from a misfolded state to its native state requires more time than ab initio folding of a β-sheet.

Abstract: The design of artificial metalloenzymes has achieved tremendous progress, although few designs can achieve catalytic performances comparable to that of native enzymes. Moreover, the structure and function of artificial metalloenzymes in non-native states has rarely been explored. Herein, we found that a c-type cytochrome b5 (Cyt b5), N57C/S71C Cyt b5, with heme covalently attached to the protein matrix through two Cys–heme linkages, adopts a non-native state with an open heme site after guanidine hydrochloride (Gdn⋅HCl)-induced unfolding, which facilitates H2O2 activation and substrate binding. Stopped-flow kinetic studies further revealed that c-type Cyt b5 in the non-native state exhibited impressive peroxidase activity comparable to that of native peroxidases, such as the most efficient horseradish peroxidase. This study presents an alternative approach to the design of functional artificial metalloenzymes by exploring enzymatic functions in non-native states.

Abstract: Many proteins fold in apparent two-state behavior, as partially folded intermediates only transiently accumulate and easily escape detection. Besides a native form and a mainly unfolded form, we captured a partially unfolded form of an acyl carrier protein from Micromonospora echinospora (meACP) in the folding/unfolding equilibrium using chemical exchange saturation transfer NMR experiments. The C-terminal region of the partially unfolded form is mainly folded and the N-terminal is unfolded. Furthermore, to understand how the folding process of meACP is influenced by solvent environments, we compared the folding dynamics of meACP in D2O, H2O and low concentration of urea. As the environment becomes more denaturing from D2O to H2O and then to urea, the unfolded state becomes increasingly populated, and the folding rate decreases. Adding a small amount of urea, which does not change solvent viscosity, has little effects on the unfolding rates, while changing H2O to D2O reduces the unfolding rates possibly due to the increase of solvent viscosity. The quantified solvent effects on the protein folding Gibbs energy and activation energy suggest that the transition state of folding may have a similar structure to the native state of the protein.

Abstract: Amyloid light-chain (AL) amyloidosis is a protein-misfolding disease characterized by accumulation of immunoglobulin light chains (LCs) into amyloid fibrils. Dimerization of a full length or variable domain (VL) of LC serves to stabilize the native state and prevent the formation of amyloid fibrils. We here analyzed the thermodynamic properties of dimerization and unfolding reactions by nonamyloidogenic VL from REI LC or its monomeric Y96K mutant using sedimentation velocity and circular dichroism. The data indicate that the equilibrium shifts to native dimerization for wild-type REI VL by elevating temperature due to the negative enthalpy change for dimer dissociation (−81.2 kJ·mol−1). The Y96K mutation did not affect the stability of the monomeric native state but increased amyloidogenicity. These results suggest that the heat-induced native homodimerization is the major factor preventing amyloid formation by wild-type REI VL. Heat-induced native oligomerization may be an efficient strategy to avoid the formation of misfolded aggregates particularly for thermostable proteins that are used at elevated temperatures under conditions where other proteins tend to misfold. Database Structural data are available in the Protein Data Bank under the accession numbers 5XP1 and 5XQY.

Abstract: Single-molecule microscopic experiments can measure the mechanical response of proteins to pulling forces applied externally along different directions (inducing different residue pairs in the proteins by uniaxial tension). This response to external forces away from equilibrium should in principle, correlate with the flexibility or stiffness of proteins in their folded states. Here, a simple topology-based atomistic anisotropic network model (ANM) is shown which captures the protein flexibility as a fundamental property that determines the collective dynamics and hence, the protein conformations in native state. An all-atom ANM is used to define two measures of protein flexibility in the native state. One measure quantifies overall stiffness of the protein and the other one quantifies protein stiffness along a particular direction which is effectively the mechanical resistance of the protein towards external pulling force exerted along that direction. These measures are sensitive to the protein sequence and yields reliable values through computations of normal modes of the protein. ANM at an atomistic level (heavy atoms) explains the experimental (atomic force microscopy) observations viz., different mechanical stability of structurally similar but sequentially distinct proteins which, otherwise were implied to possess similar mechanical properties from analytical/theoretical coarse-grained (backbone only) models. The results are exclusively demonstrated for human fibronectin (FN) protein domains. The topology of interatomic contacts in the folded states of proteins essentially determines the native flexibility. The mechanical differences of topologically similar proteins are captured from a high-resolution (atomic level) ANM at a low computational cost. The relative trend in flexibility of such proteins is reflected in their stability differences that they exhibit while unfolding in atomic force microscopic (AFM) experiments.

Abstract: Hepatitis B virus X protein (HBx) is a multifunctional protein that interacts directly with many host proteins. For example, HBx interacts with anti-apoptotic proteins, Bcl-2 and Bcl-xL, through its BH3-like motif, which leads to elevated cytosolic calcium levels, efficient viral DNA replication and the induction of apoptosis. To facilitate sample preparation and perform detailed structural characterization of the complex between HBx and Bcl-xL, we designed and purified a recombinant HBx BH3-like motif-linker-Bcl-xL fusion protein produced in E. coli. The fusion protein was characterized by size exclusion chromatography, circular dichroism and nuclear magnetic resonance experiments. Our results show that the fusion protein is a monomer in aqueous solution, forms a stable intramolecular complex, and likely retains the native conformation of the complex between Bcl-xL and the HBx BH3-like motif. Furthermore, the HBx BH3-like motif of the intramolecular complex forms an α-helix. These observations indicate that the fusion protein should facilitate structural studies aimed at understanding the interaction between HBx and Bcl-xL at the atomic level.

Abstract: Protein homeostasis, namely the ensemble of cellular mechanisms collectively controlling the activity, stability and conformational states of proteins, depends on energy-consuming processes. De novo protein synthesis requires ATP hydrolysis for peptide bond formation. Controlled degradation by the chaperone-gated proteases requires ATP hydrolysis to unfold target proteins and render their peptide bonds accessible to hydrolysis. During and following translation, different classes of molecular chaperones require ATP hydrolysis to control the conformational state of proteins, favor their folding into their active conformation and avoid, under stress, their conversion into potentially harmful aggregates. Furthermore, specific ATP-fueled unfolding chaperones can dynamically revert aggregation itself. We used here various biochemical assays and physical modeling to show that both bacterial chaperones GroEL (HSP60) and DnaK (HSP70) can use the energy liberated by ATP hydrolysis to maintain proteins in their active state even under conditions that do not favor, thermodynamically, the native state. The energy from ATP hydrolysis is thus injected by the chaperones in the system and converted into an enhanced, non-equilibrium steady-state stabilization of the native state of their substrates. Upon ATP consumption, the chaperone substrates spontaneously revert to their equilibrium non-native state.

Abstract: Defining the role of non-native interactions in directing the course of protein folding or unfolding reactions has been a difficult challenge. In particular, the extent to which such interactions play a productive role by stabilizing the structures of transition states (TSs) found on the folding and unfolding pathways of proteins is not known. On the contrary, it is thought that the TSs are expanded forms of the N state stabilized by native interactions, and it is not known whether non-native interactions can modulate TS structure. In this study of the unfolding of the SH3 domain of PI3 kinase using a microsecond mixing methodology, partial non-native structure formation is shown to occur initially during unfolding. The TS of this partial “folding during unfolding” reaction is more compact than the N state: the apparent rate constant of Trp53 burial during this reaction decreases with an increase in denaturant concentration. Kinetic studies of the unfolding of mutant variants suggest that the unusually co...

Abstract: Rational mutation of proteins based on their structural and dynamic characteristics is a useful strategy for amplifying specific fluctuations in proteins. Here, we show the effects of mutation on the conformational fluctuations and thermodynamic stability of ubiquitin. In particular, we focus on the salt bridge between K11 and E34 and the hydrogen bond between I36 and Q41, which are predicted to control the fluctuation between the basic folded state, N1, and the alternatively folded state, N2, of the protein, using high-pressure NMR spectroscopy. The E34A mutation, which disrupts the salt bridge, did not alter picosecond–to–nanosecond, microsecond–to–millisecond dynamic motions, and stability of the protein, while the Q41N mutation, which destabilizes the hydrogen bond, specifically amplified the N1–N2 conformational fluctuation and decreased stability. Based on the observed thermodynamic stabilities of the various conformational states, we showed that in the Q41N mutant, the N1 state is more significantly destabilized than the N2 state, resulting in an increase in the relative population of N2. Identifying the interactions controlling specific motions of a protein will facilitate molecular design to achieve functional dynamics beyond native state dynamics.

Abstract: The sensitivity of the stability of folded states and free energy landscapes to the differences in the hydrophobic content of the core residues has been studied for the set of 16-residue trpzips, namely, Trpzip4, Trpzip5 and Trpzip6. The combination of principal component analysis and different secondary structure order metrics as reaction coordinates has been used to characterize and identify all the underlying attractive basins corresponding to the folded and the unfolded states for each trpzip at 300 K. Our results reveal that even a single mutation in the hydrophobic core perturbs the stability of the folded peptide and the conformational preferences for the partially folded and unfolded states significantly, leading to concomitant alterations in the free energy landscape of trpzips. Trpzip4 is observed to have the most rugged and variegated free energy landscape with occurrence of four metastable unfolded states in addition to the folded native state. In contrast, Trpzip5 and Trpzip6 are characterized by two such metastable states. The order metrics pertaining to the rigidity of the turn residues and the distances between the side chains of the hydrophobic core residues have been found to be most revealing to understand the degree of discrimination among the folded states of different peptides in addition to the unfolded states. Our results suggest that both turn propensity and hydrophobic interactions influence the thermodynamics of the folding pathways of trpzips. The implications of the sequence dependent response of amino acids, effect of aromatic stacking interactions and packing of protein's interior for shaping the free energy landscape of the peptides have been highlighted.

Abstract: According to the thermodynamic hypothesis, the native state of proteins is uniquely defined by their amino acid sequence. On the other hand, according to Levinthal, the native state is just a local minimum of the free energy and a given amino acid sequence, in the same thermodynamic conditions, can assume many, very different structures that are as thermodynamically stable as the native state. This is the Levinthal limit explored in this work. Using computer simulations, we compare the interactions that stabilize the native state of four different proteins with those that stabilize three non-native states of each protein and find that the nature of the interactions is very similar for all such 16 conformers. Furthermore, an enhancement of the degree of fluctuation of the non-native conformers can be explained by an insufficient relaxation to their local free energy minimum. These results favor Levinthal’s hypothesis that protein folding is a kinetic non-equilibrium process.

Abstract: Structural account of collagen in its native environment is limited, and a water-mediated hydrogen bonding network as its stability analogue needs to be contemplated site specifically. We present in this Article that, through natural abundance 13C chemical shift anisotropy (CSA) measurement of collagen in its native state (inside bone matrix), structural and mechanistic insight of the water-mediated hydrogen bonding network can be achieved. The 13C CSA of backbone and side chain residues of collagen, perturbed under dehydration and H/D exchange conditions, can be measured. The changes in 13C CSA values due to perturbation of water content have resulted in determining the site-specific accessibility of water molecules. A further hydration topology map (HTM) representing water accessibility of native collagen has been generated on the basis of changes in 13C CSA values. Our results signify the water accessibility of the proline/hydroxyproline carbonyl region is larger compared to the hydroxyproline Cγ site ...

Abstract: The temperature trends of the standard thermodynamic functions of the native and denatured protein in solution are considered within the concept of excess mixing functions. It is assumed that some protein molecules adopt an intermediate state between native and denatured forms within the temperature range between cold and thermal denaturation and form metastable microphases as a result of a specific interaction with water. A phase diagram in the temperature–standard entropy coordinate plane representing an isobar family is proposed. Two limiting isobars are characterized by an entropy jump, which reflects the first-order phase transition between the native and denatured states. The isobars in the intermediate temperature range are represented as van der Waals curves, which reflect the equilibrium between the main phase of the molecules in native state and microphases. The difference between the phases disappears at critical points. It is assumed that the supercritical range is a macroscopically homogeneous single phase zone of reduced stability, which is represented by a dynamic system of monomers and oligomers of the native protein, monomers and clusters of the protein with partially unfolded structure. The phase diagram is collated with the elliptic phase diagram in the temperature–osmotic pressure plane.