pH-induced structural change of a multi-tryptophan protein MPT63 with immunoglobulin-like fold: identification of perturbed tryptophan residue/residues

The structural change of M. tuberculosis MPT63, which is predominantly a β-sheet protein having an immunoglobulin like fold, has been investigated in the pH range 7.5–1.5 using various biophysical techniques along with low-temperature phosphorescence (LTP) spectroscopy. MPT63 contains four Tryptophan (Trp) residues at 26, 48, 82, and 129. Although circular dichroism, steady-state and time-resolved fluorescence, time-resolved anisotropy, 1-aniline-8-naphthalene sulfonic (ANS) acid binding, and analytical ultracentrifuge depict more open largely unfolded structure of MPT63 at pH 1.5 and also more accessible nature of Trp residues to neutral quencher at pH 1.5, it is, however, not possible to assign the specific Trp residue/residues being perturbed. This problem has been resolved using LTP of MPT63, which shows optically resolved four distinct (0, 0) bands corresponding to four Trp residues in the pH range 7.5–3.0. LTP at pH 1.5 clearly reveals that the solvent-exposed Trp 82 and the almost buried Trp 129 are specifically affected compared with Trp 48 and Trp 26. Lys 8 and Lys 27 are predicted to affect Trp 129. Tyrosine residues are found to be silent even at pH 1.5. This type of specific perturbation in a multi-Trp protein has not been addressed before. LTP further indicates that partially exposed Trp 48 is preferentially quenched by acrylamide compared with other Trp residues at both pH 7.5 and 1.5. The solvent-exposed Trp 82 is surprisingly found to be not quenched by acrylamide at pH 7.5. All the results are obtained using micromolar concentration of protein and without using any Trp-substituted mutant.     

Probing the energy landscape of protein folding/unfolding transition states

Previous molecular dynamics (MD) simulations of the thermal denaturation of chymotrypsin inhibitor 2 (CI2) have provided atomic-resolution models of the transition state ensemble that is well supported by experimental studies. Here, we use simulations to further investigate the energy landscape around the transition state region. Nine structures within approximately 35 ps and 3 A C(alpha) RMSD of the transition state ensemble identified in a previous 498 K thermal denaturation simulation were quenched under the quasi-native conditions of 335 K and neutral pH. All of the structures underwent hydrophobically driven collapse in response to the drop in temperature. Structures less denatured than the transition state became structurally more native-like, while structures that were more denatured than the transition state tended to show additional loss of native structure. The structures in the immediate region of the transition state fluctuated between becoming more and less native-like. All of the starting structures had the same native-like topology and were quite similar (within 3.5 A C(alpha) RMSD). That the structures all shared native-like topology, yet diverged into either more or less native-like structures depending on which side of the transition state they occupied on the unfolding trajectory, indicates that topology alone does not dictate protein folding. Instead, our results suggest that a detailed interplay of packing interactions and interactions with water determine whether a partially denatured protein will become more native-like under refolding conditions.     

Accelerated simulation of unfolding and refolding of a large single chain globular protein

We have developed novel strategies for contracting simulation times in protein dynamics that enable us to study a complex protein with molecular weight in excess of 34 kDa. Starting from a crystal structure, we produce unfolded and then refolded states for the protein. We then compare these quantitatively using both established and new metrics for protein structure and quality checking. These include use of the programs Concoord and Darvols. Simulation of protein-folded structure well beyond the molten globule state and then recovery back to the folded state is itself new, and our results throw new light on the protein-folding process. We accomplish this using a novel cooling protocol developed for this work.     

Folding simulations of Trp‐cage mini protein in explicit solvent using biasing potential replica‐exchange molecular dynamics simulations

Replica exchange molecular dynamics (RexMD) simulations are frequently used for studying structure formation and dynamics of peptides and proteins. A significant drawback of standard temperature RexMD is, however, the rapid increase of the replica number with increasing system size to cover a desired temperature range. A recently developed Hamiltonian RexMD method has been used to study folding of the Trp‐cage protein. It employs a biasing potential that lowers the backbone dihedral barriers and promotes peptide backbone transitions along the replica coordinate. In two independent applications of the biasing potential RexMD method including explicit solvent and starting from a completely unfolded structure the formation of near‐native conformations was observed after 30–40 ns simulation time. The conformation representing the most populated cluster at the final simulation stage had a backbone root mean square deviation of ～1.3 Å from the experimental structure. This was achieved with a very modest number of five replicas making it well suited for peptide and protein folding and refinement studies including explicit solvent. In contrast, during five independent continuous 70 ns molecular dynamics simulations formation of collapsed states but no near native structure formation was observed. The simulations predict a largely collapsed state with a significant helical propensity for the helical domain of the Trp‐cage protein already in the unfolded state. Hydrogen bonded bridging water molecules were identified that could play an active role by stabilizing the arrangement of the helical domain with respect to the rest of the chain already in intermediate states of the protein. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Characterizing the inhibition of α‐synuclein oligomerization by a pharmacological chaperone that prevents prion formation by the protein PrP

Aggregation of the disordered protein α‐synuclein into amyloid fibrils is a central feature of synucleinopathies, neurodegenerative disorders that include Parkinson's disease. Small, pre‐fibrillar oligomers of misfolded α‐synuclein are thought to be the key toxic entities, and α‐synuclein misfolding can propagate in a prion‐like way. We explored whether a compound with anti‐prion activity that can bind to unfolded parts of the protein PrP, the cyclic tetrapyrrole Fe‐TMPyP, was also active against α‐synuclein aggregation. Observing the initial stages of aggregation via fluorescence cross‐correlation spectroscopy, we found that Fe‐TMPyP inhibited small oligomer formation in a dose‐dependent manner. Fe‐TMPyP also inhibited the formation of mature amyloid fibrils in vitro, as detected by thioflavin T fluorescence. Isothermal titration calorimetry indicated Fe‐TMPyP bound to monomeric α‐synuclein with a stoichiometry of 2, and two‐dimensional heteronuclear single quantum coherence NMR spectra revealed significant interactions between Fe‐TMPyP and the C‐terminus of the protein. These results suggest commonalities among aggregation mechanisms for α‐synuclein and the prion protein may exist that can be exploited as therapeutic targets.     

Secondary structure estimation of recombinant <Emphasis Type="Italic">psbH</Emphasis>, encoding a photosynthetic membrane protein of cyanobacterium <Emphasis Type="Italic">Synechocystis</Emphasis> sp. PCC 6803

The PsbH protein of cyanobacterium Synechocystis sp. PCC 6803 was expressed as a fusion protein with glutathione-S transferase (GST) in E. coli grown on a mineral medium enriched in ¹⁵N isotope. After enzymatic cleavage of the fusion protein, the ¹H-¹⁵N-HSQC spectrum of PsbH protein in presence of the detergent β-D-octyl-glucopyranoside (OG) was recorded on a Bruker DRX 500 MHz NMR spectrometer equipped with a 5 mm TXI cryoprobe to enhance the sensitivity and resolution. Non-labelled protein was used for secondary structure estimation by deconvolution from circular dichroism (CD) spectra. Experimental results were compared with our results from a structural model of PsbH using a restraint-based comparative modelling approach combined with molecular dynamics and energetic modelling. We found that PsbH shows 34–38% α-helical structure (Thr36-Ser60), a maximum of around 15% of β-sheet, and 12–19% of β-turn.     

Long-range tertiary interactions in single hammerhead ribozymes bias motional sampling toward catalytically active conformations

Enzymes generally are thought to derive their functional activity from conformational motions. The limited chemical variation in RNA suggests that such structural dynamics may play a particularly important role in RNA function. Minimal hammerhead ribozymes are known to cleave efficiently only in ～ 10-fold higher than physiologic concentrations of Mg(2+) ions. Extended versions containing native loop-loop interactions, however, show greatly enhanced catalytic activity at physiologically relevant Mg(2+) concentrations, for reasons that are still ill-understood. Here, we use Mg(2+) titrations, activity assays, ensemble, and single molecule fluorescence resonance energy transfer (FRET) approaches, combined with molecular dynamics (MD) simulations, to ask what influence the spatially distant tertiary loop-loop interactions of an extended hammerhead ribozyme have on its structural dynamics. By comparing hammerhead variants with wild-type, partially disrupted, and fully disrupted loop-loop interaction sequences we find that the tertiary interactions lead to a dynamic motional sampling that increasingly populates catalytically active conformations. At the global level the wild-type tertiary interactions lead to more frequent, if transient, encounters of the loop-carrying stems, whereas at the local level they lead to an enrichment in favorable in-line attack angles at the cleavage site. These results invoke a linkage between RNA structural dynamics and function and suggest that loop-loop interactions in extended hammerhead ribozymes-and Mg(2+) ions that bind to minimal ribozymes-may generally allow more frequent access to a catalytically relevant conformation(s), rather than simply locking the ribozyme into a single active state.     

Computational studies of the structure, dynamics and native content of amyloid-like fibrils of ribonuclease A

It is a common belief that some residues of a protein are more important than others. In some cases, point mutations of some residues make butterfly effect on the protein structure and function, but in other cases they do not. In addition, the residues important for the protein function tend to be not only conserved but also coevolved with other interacting residues in a protein. Motivated by these observations, the authors propose that there is a network composed of the residues, the residue-residue coevolution network (RRCN), where nodes are residues and links are set when the coevolutionary interaction strengths between residues are sufficiently large. The authors build the RRCN for the 44 diverse protein families. The interaction strengths are calculated by using McBASC algorithm. After constructing the RRCN, the authors identify residues that have high degree of connectivity (hub nodes), and residues that play a central role in network flow of information (C(I) nodes). The authors show that these residues are likely to be functionally important residues. Moreover, the C(I) nodes appear to be more relevant to the function than the hub nodes. Unlike other similar methods, the method described in this study is solely based on sequences. Therefore, the method can be applied to the function annotation of a wider range of proteins.     

Diffusing and colliding: the atomic level folding/unfolding pathway of a small helical protein

Proteins with ultra-fast folding/unfolding kinetics are excellent candidates for study by molecular dynamics. Here, we describe such simulations of a three helix bundle protein, the engrailed homeodomain (En-HD), which folds via the diffusion-collision model. The unfolding pathway of En-HD was characterized by seven simulations of the protein and 12 simulations of its helical fragments yielding over 1.1 micros of simulation time in water. Various conformational states along the unfolding pathway were identified. There is the compact native-like transition state, a U-shaped helical intermediate and an unfolded state with dynamic helical segments. Each of these states is in good agreement with experimental data. Examining these states as well as the transitions between them, we find the role of long-range tertiary contacts, specifically salt-bridges, important in the folding/unfolding pathway. In the folding direction, charged residues form long-range tertiary contacts before the hydrophobic core is formed. The formation of HII is assisted by a specific salt-bridge and by non-specific (fluctuating) tertiary contacts, which we call contact-assisted helix formation. Salt-bridges persist as the protein approaches the transition state, stabilizing HII until the hydrophobic core is formed. To complement this information, simulations of fragments of En-HD illustrate the helical propensities of the individual segments. By thermal denaturation, HII proved to be the least stable helix, unfolding in less than 450 ps at high temperature. We observed the low helical propensity of C-terminal residues from HIII in fragment simulations which, when compared to En-HD unfolding simulations, link the unraveling of HIII to the initial event that drives the unfolding of En-HD.     

Alpha helical crossovers favor right‐handed supersecondary structures by kinetic trapping: The phone cord effect in protein folding

The remarkable predominance of right‐handedness in beta‐alpha‐beta helical crossovers has been previously explained in terms of thermodynamic stability and kinetic accessibility, but a different kinetic trapping mechanism may also play a role. If the beta‐sheet contacts are made before the crossover helix is fully formed, and if the backbone angles of the folding helix follows the energetic pathway of least resistance, then the helix would impart a torque on the ends of the two strands. Such a torque would tear apart a left‐handed conformation but hold together a right‐handed one. Right‐handed helical crossovers predominate even in all‐alpha proteins, where previous explanations based on the preferred twist of the beta sheet do not apply. Using simple molecular simulations, we can reproduce the right‐handed preference in beta‐alpha‐beta units, without imposing specific beta‐strand geometry. The new kinetic trapping mechanism is dubbed the “phone cord effect” because it is reminiscent of the way a helical phone cord forms superhelices to relieve torsional stress. Kinetic trapping explains the presence of a right‐handed superhelical preference in alpha helical crossovers and provides a possible folding mechanism for knotted proteins.     

Computational approaches to study protein unfolding: Hen egg white lysozyme as a case study

Four methods are compared to drive the unfolding of a protein: (1) high temperature (T-run), (2) high pressure (P-run), (3) by imposing a gradual increase in the mean radius of the protein using a penalty function added to the physical interaction function (F-run, radial force driven unfolding), and (4) by weak coupling of the difference between the temperature of the radially outward moving atoms and the radially inward moving atoms to an external temperature bath (K-run, kinetic energy driven unfolding). The characteristic features of the four unfolding pathways are analyzed in order to detect distortions due to the size or the type of the applied perturbation, as well as the features that are common to all of them. Hen egg white lysozyme is used as a test system. The simulations are analyzed and compared to experimental data like ¹H-NMR amide proton exchange-folding competition, heat capacity, and compressibility measurements. © 1995 Wiley-Liss, Inc.     

Molecular dynamics simulations and free energy analyses on the dimer formation of an amyloidogenic heptapeptide from human beta2-microglobulin: implication for the protofibril structure

Amyloid formation is associated with many neurodegenerative diseases. Recent findings suggest that early oligomeric aggregates could be major sources of toxicity. We present a computational investigation of the first step of amyloid initiation-dimer formation of a seven residue peptide (NHVTLSQ) from human beta2-microglobulin at pH 2.0, which renders +2.0 units charges to each peptide. A total of over 1.2 micros of simulations with explicit solvent and 1.0 micros of simulations with implicit solvent were conducted. Main-chain conformational restraint was applied to facilitate the formation of ordered dimers. An antiparallel beta-sheet with six main-chain hydrogen bonds was dominant in the implicit solvent simulations. In contrast, no stable dimers were observed in the two negative controls, the mouse heptapeptide (KHDSMAE, +3.0 units charges) and the scrambled human heptapeptide (QVLHTSN). Explicit solvent simulations presented a more complex scenario. The wild-type human heptapeptide formed predominantly antiparallel beta-sheets ( approximately 38%) although parallel ones ( approximately 12%) were also observed. Hydrophobic contacts preceded hydrogen bond saturation in the majority of the association events in the explicit solvent simulations, highlighting the important role of hydrophobic interaction in amyloid initiation. The fact that the mouse dimer dissociated immediately after the removal of conformational restraint suggests that the higher conformational entropy barrier, along with the stronger charge repulsion and weaker hydrophobic interaction, contributed to its inability to form amyloid fibril. The closeness of positive charge pairs in the dimers of the scrambled human heptapeptide may prohibit further beta-sheet extension and fibril growth. Combining the results from simulations and free energy analyses, we propose that the building block for this amyloid fibril is an antiparallel dimer with a two-residue register shift and six main-chain hydrogen bonds. A double-layer protofibril structure is also proposed in which two antiparallel beta-sheets face each other and are held together by hydrophobic staples and hydrogen bonds of the polar side-chains.     

A cytotoxicity,optical spectroscopy and computational binding analysis of 4‐[3‐acetyl‐5‐(acetylamino)‐2‐methyl‐2,3‐dihydro‐1,3,4‐thiadiazole‐2‐yl]phenyl benzoate in calf thymus DNA

In this study the interaction mechanism between newly synthesized 4‐(3‐acetyl‐5‐(acetylamino)‐2‐methyl‐2, 3‐dihydro‐1,3,4‐thiadiazole‐2‐yl) phenyl benzoate (thiadiazole derivative) anticancer active drug with calf thymus DNA was investigated by using various optical spectroscopy techniques along with computational technique. The absorption spectrum shows a clear shift in the lower wavelength region, which may be due to strong hypochromic effect in the ctDNA and the drug. The results of steady state fluorescence spectroscopy show that there is static quenching occurring while increasing the thiadiazole drug concentration in the ethidium bromide‐ctDNA system. Also the binding constant (K), thermo dynamical parameters of enthalpy change (ΔH°), entropy change (ΔS°) Gibbs free energy change (ΔG°) were calculated at different temperature (293 K, 298 K) and the results are in good agreement with theoretically calculated MMGBSA binding analysis. Time resolved emission spectroscopy analysis clearly explains the thiadiazole derivative competitive intercalation in the ethidium bromide‐ctDNA system. Further, molecular docking studies was carried out to understand the hydrogen bonding and hydrophobic interaction between ctDNA and thiadiazole derivative molecule. In addition the docking and molecular dynamics charge distribution analysis was done to understand the internal stability of thiadiazole derivative drug binding sites of ctDNA. The global reactivity of thiadiazole derivative such as electronegativity, electrophilicity and chemical hardness has been calculated.     

Effect of acidic pH on the stability of α‐synuclein dimers

Environmental factors, such as acidic pH, facilitate the assembly of α‐synuclein (α‐Syn) in aggregates, but the impact of pH on the very first step of α‐Syn aggregation remains elusive. Recently, we developed a single‐molecule approach that enabled us to measure directly the stability of α‐Syn dimers. Unlabeled α‐Syn monomers were immobilized on a substrate, and fluorophore‐labeled monomers were added to the solution to allow them to form dimers with immobilized α‐Syn monomers. The dimer lifetimes were measured directly from the fluorescence bursts on the time trajectories. Herein, we applied the single‐molecule tethered approach for probing of intermolecular interaction to characterize the effect of acidic pH on the lifetimes of α‐Syn dimers. The experiments were performed at pH 5 and 7 for wild‐type α?Syn and for two mutants containing familial type mutations E46K and A53T. We demonstrate that a decrease of pH resulted in more than threefold increase in the α‐Syn dimers lifetimes with some variability between the α‐Syn species. We hypothesize that the stabilization effect is explained by neutralization of residues 96–140 of α‐Syn and this electrostatic effect facilitates the association of the two monomers. Given that dimerization is the first step of α‐Syn aggregation, we posit that the electrostatic effect thereby contributes to accelerating α‐Syn aggregation at acidic pH. © 2016 Wiley Periodicals, Inc. Biopolymers 105: 715–724, 2016.     

Enhanced sampling of peptide and protein conformations using replica exchange simulations with a peptide backbone biasing-potential

During replica exchange molecular dynamics (RexMD) simulations, several replicas of a system are simulated at different temperatures in parallel allowing for exchange between replicas at frequent intervals. This technique allows significantly improved sampling of conformational space and is increasingly being used for structure prediction of peptides and proteins. A drawback of the standard temperature RexMD is the rapid increase of the replica number with increasing system size to cover a desired temperature range. In an effort to limit the number of replicas, a new Hamiltonian-RexMD method has been developed that is specifically designed to enhance the sampling of peptide and protein conformations by applying various levels of a backbone biasing potential for each replica run. The biasing potential lowers the barrier for backbone dihedral transitions and promotes enhanced peptide backbone transitions along the replica coordinate. The application on several peptide cases including in all cases explicit solvent indicates significantly improved conformational sampling when compared with standard MD simulations. This was achieved with a very modest number of 5-7 replicas for each simulation system making it ideally suited for peptide and protein folding simulations as well as refinement of protein model structures in the presence of explicit solvent.     

Unbinding and unfolding of adhesion protein complexes through stretching: Interplay between shear and tensile mechanical clamps

Using coarse‐grained molecular dynamics simulations, we analyze mechanically induced dissociation and unfolding of the protein complex CD48‐2B4. This heterodimer is an indispensable component of the immunological system: 2B4 is a receptor on natural killer cells whereas CD48 is expressed on surfaces of various immune cells. So far, its mechanostability has not been assessed either experimentally or theoretically. We find that the dissociation processes strongly depend on the direction of pulling and may take place in several pathways. Interestingly, the CD48‐2B4 interface can be divided into three distinct patches that act as units when resisting the pulling forces. At experimentally accessible pulling speeds, the characteristic mechanostability forces are in the range between 100 and 200 pN, depending on the pulling direction. These characteristic forces need not be associated with tensile forces involved in the act of separation of the complex because prior shear‐involving unraveling within individual proteins may give rise to a higher force peak. Proteins 2014; 82:3144–3153. © 2014 Wiley Periodicals, Inc.     

分子动力学模拟尿素对水溶液中蛋白质构象转变的影响

采用分子动力学方法和全原子模型研究尿素和水分子对模型蛋白S-肽链结构转化的影响。模拟结果显示S-肽链的变性速率常数k值随着尿素浓度的增加而先降低后升高,在尿素浓度为2.9 mol/L时达到最低值。模拟了不同尿素浓度下尿素-肽链、水-肽链以及肽链分子氢键的形成状况。结果表明:尿素浓度较低时,尿素分子与S-肽链的极性氨基酸侧链形成氢键,但不破坏其分子内的骨架氢键,尿素在S-肽链水化层外形成限制性空间,增强了S-肽链的稳定性。随着尿素的升高,尿素分子进入S-肽链内部并与其内部氨基酸残基形成氢键,导致S-肽链的骨架氢键丧失,S-肽链发生去折叠。上述模拟结果与文献报道的实验结果一致,从分子水平上揭示了尿素对蛋白质分子结构变化的影响机制,对于研究和发展蛋白质折叠及稳定化技术具有指导意义。     

Investigating energy‐based pool structure selection in the structure ensemble modeling with experimental distance constraints: The example from a multidomain protein Pub1

The structural variations of multidomain proteins with flexible parts mediate many biological processes, and a structure ensemble can be determined by selecting a weighted combination of representative structures from a simulated structure pool, producing the best fit to experimental constraints such as interatomic distance. In this study, a hybrid structure‐based and physics‐based atomistic force field with an efficient sampling strategy is adopted to simulate a model di‐domain protein against experimental paramagnetic relaxation enhancement (PRE) data that correspond to distance constraints. The molecular dynamics simulations produce a wide range of conformations depicted on a protein energy landscape. Subsequently, a conformational ensemble recovered with low‐energy structures and the minimum‐size restraint is identified in good agreement with experimental PRE rates, and the result is also supported by chemical shift perturbations and small‐angle X‐ray scattering data. It is illustrated that the regularizations of energy and ensemble‐size prevent an arbitrary interpretation of protein conformations. Moreover, energy is found to serve as a critical control to refine the structure pool and prevent data overfitting, because the absence of energy regularization exposes ensemble construction to the noise from high‐energy structures and causes a more ambiguous representation of protein conformations. Finally, we perform structure‐ensemble optimizations with a topology‐based structure pool, to enhance the understanding on the ensemble results from different sources of pool candidates.     

Differential temporal beta‐diversity patterns of native and non‐native arthropod species in a fragmented native forest landscape

An important factor that hinders the management of non‐native species is a general lack of information regarding the biogeography of non‐natives, and, in particular, their rates of turnover. Here, we address this research gap by analysing differences in temporal beta‐diversity (using both pairwise and multiple‐time dissimilarity metrics) between native and non‐native species, using a novel time‐series dataset of arthropods sampled in native forest fragments in the Azores. We use a null model approach to determine whether temporal beta‐diversity was due to deterministic processes or stochastic colonisation and extinction events, and linear modelling selection to assess the factors driving variation in temporal beta‐diversity between plots. In accordance with our predictions, we found that the temporal beta‐diversity was much greater for non‐native species than for native species, and the null model analyses indicated that the turnover of non‐native species was due to stochastic events. No predictor variables were found to explain the turnover of native or non‐native species. We attribute the greater turnover of non‐native species to source‐sink processes and the close proximity of anthropogenic habitats to the fragmented native forest plots sampled in our study. Thus, our findings point to ways in which the study of turnover can be adapted for future applications in habitat island systems. The implications of this for biodiversity conservation and management are significant. The high rate of stochastic turnover of non‐native species indicates that attempts to simply reduce the populations of non‐native species in situ within native habitats may not be successful. A more efficient management strategy would be to interrupt source‐sink dynamics by improving the harsh boundaries between native and adjacent anthropogenic habitats.     

Meeting report: SMART timing—principles of single molecule techniques course at the University of Michigan 2014

Four days after the announcement of the 2014 Nobel Prize in Chemistry for “the development of super‐resolved fluorescence microscopy” based on single molecule detection, the Single Molecule Analysis in Real‐Time (SMART) Center at the University of Michigan hosted a “Principles of Single Molecule Techniques 2014” course. Through a combination of plenary lectures and an Open House at the SMART Center, the course took a snapshot of a technology with an especially broad and rapidly expanding range of applications in the biomedical and materials sciences. Highlighting the continued rapid emergence of technical and scientific advances, the course underscored just how brightly the future of the single molecule field shines. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 296–302, 2015.