Search

Solution NMR spectroscopy is a well-established tool with unique advantages for structural studies of RNA molecules. However, for large RNA sequences, the NMR resonances often overlap severely. A reliable way to perform resonance assignment and allow further analysis despite spectral crowding is the use of site-specific isotope labeling in sample preparation. While solid-phase oligonucleotide synthesis has several advantages, RNA length and availability of isotope-labeled building blocks are persistent issues. Purely enzymatic methods represent an alternative and have been presented in the literature. In this study, we report on a method in which we exploit the preference of T7 RNA polymerase for nucleotide monophosphates over triphosphates for the 5' position, which allows 5'-labeling of RNA. Successive ligation to an unlabeled RNA strand generates a site-specifically labeled RNA. We show the successful production of such an RNA sample for NMR studies, report on experimental details and expected yields, and present the surprising finding of a previously hidden set of peaks which reveals conformational exchange in the RNA structure. This study highlights the feasibility of site-specific isotope-labeling of RNA with enzymatic methods. [ABSTRACT FROM AUTHOR]

Despite the progress in prediction of protein complexes over the last decade, recent blind protein complex structure prediction challenges revealed limited success rates (less than 20% models with DockQ score > 0.4) on targets that exhibit significant conformational change upon binding. To overcome limitations in capturing backbone motions, we developed a new, aggressive sampling method that incorporates temperature replica exchange Monte Carlo (T-REMC) and conformational sampling techniques within docking protocols in Rosetta. Our method, ReplicaDock 2.0, mimics induced-fit mechanism of protein binding to sample backbone motions across putative interface residues on-the-fly, thereby recapitulating binding-partner induced conformational changes. Furthermore, ReplicaDock 2.0 clocks in at 150-500 CPU hours per target (protein-size dependent); a runtime that is significantly faster than Molecular Dynamics based approaches. For a benchmark set of 88 proteins with moderate to high flexibility (unbound-to-bound iRMSD over 1.2 Å), ReplicaDock 2.0 successfully docks 61% of moderately flexible complexes and 35% of highly flexible complexes. Additionally, we demonstrate that by biasing backbone sampling particularly towards residues comprising flexible loops or hinge domains, highly flexible targets can be predicted to under 2 Å accuracy. This indicates that additional gains are possible when mobile protein segments are known. Author summary: Proteins bind each other in a highly specific and regulated manner, and these associated dynamics of binding are intimately linked to their function. Conventional techniques of structure determination such as cryo-EM, X-ray crystallography and NMR are time-consuming and arduous. Using a temperature-replica exchange Monte Carlo approach that mimics the kinetic mechanism of "induced fit" binding, we improved prediction of protein complex structures, particularly for targets that exhibit considerable conformational changes upon binding (Interface root mean square deviation (unbound-bound) > 1.2 Å. Capturing these binding-induced conformational changes in proteins can aid us in better understanding biological mechanisms and suggest intervention strategies for disease mechanisms. [ABSTRACT FROM AUTHOR]

Sugars are crucial components in biosystems and industrial applications. In aqueous environments, the natural state of short saccharides or charged glycosaminoglycans is floating and wiggling in solution. Therefore, tools to characterize their structure in a native aqueous environment are crucial but not always available. Here, we show that a combination of Raman/ROA and, on occasions, NMR experiments with Molecular Dynamics (MD) and Quantum Mechanics (QM) is a viable method to gain insights into structural features of sugars in solutions. Combining these methods provides information about accessible ring puckering conformers and their proportions. It also provides information about the conformation of the linkage between the sugar monomers, i.e., glycosidic bonds, allowing for identifying significantly accessible conformers and their relative abundance. For mixtures of sugar moieties, this method enables the deconvolution of the Raman/ROA spectra to find the actual amounts of its molecular constituents, serving as an effective analytical technique. For example, it allows calculating anomeric ratios for reducing sugars and analyzing more complex sugar mixtures to elucidate their real content. Altogether, we show that combining Raman/ROA spectroscopies with simulations is a versatile method applicable to saccharides. It allows for accessing many features with precision comparable to other methods routinely used for this task, making it a viable alternative. Furthermore, we prove that the proposed technique can scale up by studying the complicated Raffinose trisaccharide, and therefore, we expect its wide adoption to characterize sugar structural features in solution. Author summary: Saccharides are an essential part of many biosystems. Not only their identification but also their structural characterization is crucial to understanding their role. For example, animal cells possess a surrounding aqueous layer rich in sugars. Although its sugar content is known, the structures formed by these sugars in such a flexible and mobile environment are lesser understood. A trustworthy characterization of their atomistic structure will foster our knowledge and contribute to drug development. This characterization is also crucial in developing accurate computer models needed to study such environments. Unfortunately, many available structure characterization techniques are inadequate when applied to flexible molecules such as saccharides. Therefore, new methods are welcomed. We use molecular dynamics simulations and quantum chemical calculations to extract structural data of sugars in solution when compared against experimental Raman/ROA/NMR spectra. This combination of techniques allows interpreting these spectra to extract structural features of saccharides. We get, for example, the molecular structure of sugar rings (puckering confirmation) or the preferred orientation of bonds connecting two sugar units (glycosidic bonds). We can determine the presence of different sugar moieties in solutions and their proportions. Overall, this technique should contribute to the understanding of sugars in their native aqueous environment. [ABSTRACT FROM AUTHOR]

Large-scale conformational changes are essential for proteins to function properly. Given that these transition events rarely occur, however, it is challenging to comprehend their underlying mechanisms through experimental and theoretical approaches. In this study, we propose a new computational methodology called internal coordinate normal mode-guided elastic network interpolation (ICONGENI) to predict conformational transition pathways in proteins. Its basic approach is to sample intermediate conformations by interpolating the interatomic distance between two end-point conformations with the degrees of freedom constrained by the low-frequency dynamics afforded by normal mode analysis in internal coordinates. For validation of ICONGENI, it is applied to proteins that undergo open-closed transitions, and the simulation results (i.e., simulated transition pathways) are compared with those of another technique, to demonstrate that ICONGENI can explore highly reliable pathways in terms of thermal and chemical stability. Furthermore, we generate an ensemble of transition pathways through ICONGENI and investigate the possibility of using this method to reveal the transition mechanisms even when there are unknown metastable states on rough energy landscapes. [ABSTRACT FROM AUTHOR]

S-adenosylmethionine (SAM) is one of the most important enzyme substrates. It is vital for the function of various proteins, including large group of methyltransferases (MTs). Intriguingly, some bacterial and eukaryotic MTs, while catalysing the same reaction, possess significantly different topologies, with the former being a knotted one. Here, we conducted a comprehensive analysis of SAM conformational space and factors that affect its vastness. We investigated SAM in two forms: free in water (via NMR studies and explicit solvent simulations) and bound to proteins (based on all data available in the PDB and on all-atom molecular dynamics simulations in water). We identified structural descriptors—angles which show the major differences in SAM conformation between unknotted and knotted methyltransferases. Moreover, we report that this is caused mainly by a characteristic for knotted MTs compact binding site formed by the knot and the presence of adenine-binding loop. Additionally, we elucidate conformational restrictions imposed on SAM molecules by other protein groups in comparison to conformational space in water. Author summary: The topology of a folded polypeptide chain has great impact on the resulting protein function and its interaction with ligands. Interestingly, topological constraints appear to affect binding of one of the most ubiquitous substrates in the cell, S-adenosylmethionine (SAM), to its target proteins. Here, we demonstrate how binding sites of specific proteins restrict SAM conformational freedom in comparison to its unbound state, with a special interest in proteins with non-trivial topology, including an exciting group of knotted methyltransferases. Using a vast array of computational methods combined with NMR experiments, we identify key structural features of knotted methyltransferases that impose unorthodox SAM conformations. We compare them with the characteristics of standard, unknotted SAM binding proteins. These results are significant for understanding differences between analogous, yet topologically different enzymes, as well as for future rational drug design. [ABSTRACT FROM AUTHOR]

A designed disulfide-rich β-hairpin peptide that dimerizes spontaneously served as a hinge-type connection between proteins. Here, we analyze the range of dynamics of this hinge dimer with the aim of proposing new applications for the DNA-encodable peptide and establishing guidelines for the computational analysis of other disulfide hinges. A recent structural analysis based on nuclear magnetic resonance spectroscopy and ion mobility spectrometry revealed an averaged conformation in the hinge region which motivated us to investigate the dynamic behavior using a combination of molecular dynamics simulation, metadynamics and free energy surface analysis to characterize the conformational space available to the hinge. Principal component analysis uncovered two slow modes of the peptide, namely, the opening and closing motion and twisting of the two β-hairpins assembling the hinge. Applying a collective variable (CV) that mimics the first dominating mode, led to a major expansion of the conformational space. The description of the dynamics could be achieved by analysis of the opening angle and the twisting of the β-hairpins and, thus, offers a methodology that can also be transferred to other derivatives. It has been demonstrated that the hinge peptide's lowest energy conformation consists of a large opening angle and strong twist but is separated by small energy barriers and can, thus, adopt a closed and untwisted structure. With the aim of proposing further applications for the hinge peptide, we simulated its behavior in the sterically congested environment of a four-helix bundle. Preliminary investigations show that one helix is pushed out and a three-helix bundle forms. The insights gained into the dynamics of the tetra-disulfide peptide and analytical guidelines developed in this study may contribute to the understanding of the structure and function of more complex hinge-type proteins, such as the IgG antibody family. [ABSTRACT FROM AUTHOR]

Protein structures, usually visualized in various highly idealized forms focusing on the three-dimensional arrangements of secondary structure elements, can also be described as lists of interacting residues or atoms and visualized as two-dimensional distance or contact maps. We show that contact maps provide an ideal tool to describe and analyze differences between structures of proteins in different conformations. Expanding functionality of the PDBFlex server and database developed previously in our group, we describe how analysis of difference contact maps (DCMs) can be used to identify critical interactions stabilizing alternative protein conformations, recognize residues and positions controlling protein functions and build hypotheses as to molecular mechanisms of disease mutations. [ABSTRACT FROM AUTHOR]

Gastrotropin, the intracellular carrier of bile salts in the small intestine, binds two ligand molecules simultaneously in its internal cavity. The molecular rearrangements required for ligand entry are not yet fully clear. To improve our understanding of the binding process we combined molecular dynamics simulations with previously published structural and dynamic NMR parameters. The resulting ensembles reveal two distinct modes of barrel opening with one corresponding to the transition between the apo and holo states, whereas the other affecting different protein regions in both ligation states. Comparison of the calculated structures with NMR-derived parameters reporting on slow conformational exchange processes suggests that the protein undergoes partial unfolding along a path related to the second mode of the identified barrel opening motion. [ABSTRACT FROM AUTHOR]

Combining structural proteomics experimental data with computational methods is a powerful tool for protein structure prediction. Here, we apply a recently-developed approach for de novo protein structure determination based on the incorporation of short-distance crosslinking data as constraints in discrete molecular dynamics simulations (CL-DMD) for the determination of conformational ensemble of the intrinsically disordered protein α-synuclein in the solution. The predicted structures were in agreement with hydrogen-deuterium exchange, circular dichroism, surface modification, and long-distance crosslinking data. We found that α-synuclein is present in solution as an ensemble of rather compact globular conformations with distinct topology and inter-residue contacts, which is well-represented by movements of the large loops and formation of few transient secondary structure elements. Non-amyloid component and C-terminal regions were consistently found to contain β-structure elements and hairpins. [ABSTRACT FROM AUTHOR]

Follicle-stimulating hormone receptor (FSHR) is a G-protein coupled receptor (GPCR) and a prototype of the glycoprotein hormone receptors subfamily of GPCRs. Structural data of the FSHR ectodomain in complex with follicle-stimulating hormone suggests a “pull and lift” activation mechanism that triggers a conformational change on the seven α-helix transmembrane domain (TMD). To analyze the conformational changes of the FSHR TMD resulting from sequence variants associated with reproductive impairment in humans, we set up a computational approach combining helix modeling and molecular simulation methods to generate conformational ensembles of the receptor at room (300 K) and physiological (310 K) temperatures. We examined the receptor dynamics in an explicit membrane environment of polyunsaturated phospholipids and solvent water molecules. The analysis of the conformational dynamics of the functional (N680 and S680) and dysfunctional (mutations at D408) variants of the FSHR allowed us to validate the FSHR-TMD model. Functional variants display a concerted motion of flexible intracellular regions at TMD helices 5 and 6. Disruption of side chain interactions and conformational dynamics were detected upon mutation at D408 when replaced with alanine, arginine, or tyrosine. Dynamical network analysis confirmed that TMD helices 2 and 5 may share communication pathways in the functional FSHR variants, whereas no connectivity was detected in the dysfunctional mutants, indicating that the global dynamics of the FSHR was sensitive to mutations at amino acid residue 408, a key position apparently linked to misfolding and variable cell surface plasma membrane expression of FSHRs with distinct mutations at this position. [ABSTRACT FROM AUTHOR]

Covalent modification of proteins by ubiquitin or ubiquitin chains is one of the most prevalent post-translational modifications in eukaryotes. Different types of ubiquitin chains are assumed to selectively signal respectively modified proteins for different fates. In support of this hypothesis, structural studies have shown that the eight possible ubiquitin dimers adopt different conformations. However, at least in some cases, these structures cannot sufficiently explain the molecular basis of the selective signaling mechanisms. This indicates that the available structures represent only a few distinct conformations within the entire conformational space adopted by a ubiquitin dimer. Here, molecular simulations on different levels of resolution can complement the structural information. We have combined exhaustive coarse grained and atomistic simulations of all eight possible ubiquitin dimers with a suitable dimensionality reduction technique and a new method to characterize protein-protein interfaces and the conformational landscape of protein conjugates. We found that ubiquitin dimers exhibit characteristic linkage type-dependent properties in solution, such as interface stability and the character of contacts between the subunits, which can be directly correlated with experimentally observed linkage-specific properties. [ABSTRACT FROM AUTHOR]

Mutations in superoxide dismutase 1 (SOD1) associated with familial amyotrophic lateral sclerosis (fALS) induce the protein to misfold and aggregate. Missense mutations at more than 80 different amino acid positions have been associated with disease. How these mutations heighten the propensity of SOD1 to misfold and aggregate is unclear. With so many mutations, it is possible that more than one mechanism of aggregation may be involved. Of many possible mechanisms to explain heightened aggregation, one that has been suggested is that mutations that eliminate charged amino acids could diminish repulsive forces that would inhibit aberrant protein:protein interactions. Mutations at twenty-one charged residues in SOD1 have been associated with fALS, but of the 11 Lys residues in the protein, only 1 has been identified as mutated in ALS patients. Here, we examined whether loss of positively charged surface Lys residues in SOD1 would induce misfolding and formation of intracellular inclusions. We mutated four different Lys residues (K30, K36, K75, K91) in SOD1 that are not particularly well conserved, and expressed these variants as fusion proteins with yellow fluorescent protein (YFP) to assess inclusion formation. We also assessed whether these mutations induced binding to a conformation-restricted SOD1 antibody, designated C4F6, which recognizes non-natively folded protein. Although we observed some mutations to cause enhanced C4F6 binding, we did not observe that mutations that reduce charge at these positions caused the protein to form intracellular inclusions. Our findings may have implications for the low frequency of mutations at Lys residues SOD1 in ALS patients. [ABSTRACT FROM AUTHOR]

The structure and function of clinical dosage insulin and its analogues were assessed. This included ‘native insulins’ (human recombinant, bovine, porcine), ‘fast-acting analogues’ (aspart, glulisine, lispro) and ‘slow-acting analogues’ (glargine, detemir, degludec). Analytical ultracentrifugation, both sedimentation velocity and equilibrium experiments, were employed to yield distributions of both molar mass and sedimentation coefficient of all nine insulins. Size exclusion chromatography, coupled to multi-angle light scattering, was also used to explore the function of these analogues. On ultracentrifugation analysis, the insulins under investigation were found to be in numerous conformational states, however the majority of insulins were present in a primarily hexameric conformation. This was true for all native insulins and two fast-acting analogues. However, glargine was present as a dimer, detemir was a multi-hexameric system, degludec was a dodecamer (di-hexamer) and glulisine was present as a dimer-hexamer-dihexamer system. However, size-exclusion chromatography showed that the two hexameric fast-acting analogues (aspart and lispro) dissociated into monomers and dimers due to the lack of zinc in the mobile phase. This comprehensive study is the first time all nine insulins have been characterised in this way, the first time that insulin detemir have been studied using analytical ultracentrifugation and the first time that insulins aspart and glulisine have been studied using sedimentation equilibrium. The structure and function of these clinically administered insulins is of critical importance and this research adds novel data to an otherwise complex functional physiological protein. [ABSTRACT FROM AUTHOR]

The fatty acid-binding proteins L-BABP and Rep1-NCXSQ bind to anionic lipid membranes by electrostatic interactions. According to Molecular Dynamics (MD) simulations, the interaction of the protein macrodipole with the membrane electric field is a driving force for protein binding and orientation in the interface. To further explore this hypothesis, we studied the interactions of these proteins with cationic lipid membranes. As in the case of anionic lipid membranes, we found that both proteins, carrying a negative as well as a positive net charge, were bound to the positively charged membrane. Their major axis, those connecting the bottom of the β-barrel with the α-helix portal domain, were rotated about 180 degrees as compared with their orientations in the anionic lipid membranes. Fourier transform infrared (FTIR) spectroscopy of the proteins showed that the positively charged membranes were also able to induce conformational changes with a reduction of the β-strand proportion and an increase in α-helix secondary structure. Fatty acid-binding proteins (FABPs) are involved in several cell processes, such as maintaining lipid homeostasis in cells. They transport hydrophobic molecules in aqueous medium and deliver them into lipid membranes. Therefore, the interfacial orientation and conformation, both shown herein to be electrostatically determined, have a strong correlation with the specific mechanism by which each particular FABP exerts its biological function. [ABSTRACT FROM AUTHOR]

The N-terminal region of G protein-coupled receptors can be efficiently targeted for the generation of receptor-selective antibodies. These antibodies are useful for the biochemical characterization of the receptors. In this study, we developed a set of criteria to select the optimal epitope and applied them to generate antibodies to the N-terminal region of 34 different G protein-coupled receptors. The antibody characterization revealed that a subset of antibodies exhibited increased recognition of the receptor following agonist treatment and this increase could be blocked by treatment with the receptor antagonist. An analysis of the epitopes showed that those antibodies that exhibit increased recognition are on average twelve residues long, have an overall net negative charge and are enriched in aspartic and glutamic acids. These antibodies are useful since they facilitate studies examining dose dependent increases in recognition of receptors in heterologous cells as well as in native tissue. Another interesting use of these antibodies is that they facilitate measuring changes in receptor recognition in brain following peripheral drug administration; for example, systemic administration of cocaine, a blocker of dopamine transporter that increases local dopamine levels at the synapse, was found to lead to increases in antibody recognition of dopamine receptors in the brain. Taken together these studies, in addition to describing novel tools to study native receptors, provide a framework for the generation of antibodies to G protein-coupled receptors that can detect ligand-induced conformational changes. [ABSTRACT FROM AUTHOR]

Small-angle X-ray scattering is an increasingly popular technique used to detect protein structures and ensembles in solution. However, the refinement of structures and ensembles against SAXS data is often ambiguous due to the low information content of SAXS data, unknown systematic errors, and unknown scattering contributions from the solvent. We offer a solution to such problems by combining Bayesian inference with all-atom molecular dynamics simulations and explicit-solvent SAXS calculations. The Bayesian formulation correctly weights the SAXS data versus prior physical knowledge, it quantifies the precision or ambiguity of fitted structures and ensembles, and it accounts for unknown systematic errors due to poor buffer matching. The method further provides a probabilistic criterion for identifying the number of states required to explain the SAXS data. The method is validated by refining ensembles of a periplasmic binding protein against calculated SAXS curves. Subsequently, we derive the solution ensembles of the eukaryotic chaperone heat shock protein 90 (Hsp90) against experimental SAXS data. We find that the SAXS data of the apo state of Hsp90 is compatible with a single wide-open conformation, whereas the SAXS data of Hsp90 bound to ATP or to an ATP-analogue strongly suggests heterogenous ensembles of a closed and a wide-open state. [ABSTRACT FROM AUTHOR]

SecB, a small tetrameric chaperone in Escherichia coli, plays a crucial role during protein export via the general secretory pathway by binding precursor polypeptides in a nonnative conformation and passing them to SecA, the ATPase of the translocon. The dissociation constants for the interactions are known; however to relate studies in vitro to export in a living cell requires knowledge of the concentrations of the proteins in the cell. Presently in the literature there is no report of a rigorous determination of the intracellular concentration of SecB. The values available vary over 60 fold and the details of the techniques used are not given. Here we use quantitative immunoblotting to determine the level of SecB expressed from the chromosome in E.coli grown in two commonly used media. In rich medium SecB was present at 1.6 ± 0.2 μM and in minimal medium at 2.5 ± 0.6 μM. These values allow studies of SecB carried out in vitro to be applied to the situation in the cell as SecB interacts with its binding partners to move precursor polypeptides through the export pathway. [ABSTRACT FROM AUTHOR]

All organisms have evolved pathways to respond to different forms of cellular stress. The Gcn2 kinase is best known as a regulator of translation initiation in response to starvation for amino acids. Work in budding yeast has showed that the molecular mechanism of GCN2 activation involves the binding of uncharged tRNAs, which results in a conformational change and GCN2 activation. This pathway requires GCN1, which ensures delivery of the uncharged tRNA onto GCN2. However, Gcn2 is activated by a number of other stresses which do not obviously involve accumulation of uncharged tRNAs, raising the question how Gcn2 is activated under these conditions. Here we investigate the requirement for ongoing translation and tRNA binding for Gcn2 activation after different stresses in fission yeast. We find that mutating the tRNA-binding site on Gcn2 or deleting Gcn1 abolishes Gcn2 activation under all the investigated conditions. These results suggest that tRNA binding to Gcn2 is required for Gcn2 activation not only in response to starvation but also after UV irradiation and oxidative stress. [ABSTRACT FROM AUTHOR]

Fluorescence spectroscopy techniques like Förster resonance energy transfer (FRET) and fluorescence correlation spectroscopy (FCS) have become important tools for the in vitro and in vivo investigation of conformational dynamics in biomolecules. These methods rely on the distance-dependent quenching of the fluorescence signal of a donor fluorophore either by a fluorescent acceptor fluorophore (FRET) or a non-fluorescent quencher, as used in FCS with photoinduced electron transfer (PET). The attachment of fluorophores to the molecule of interest can potentially alter the molecular properties and may affect the relevant conformational states and dynamics especially of flexible biomolecules like intrinsically disordered proteins (IDP). Using the intrinsically disordered S-peptide as a model system, we investigate the impact of terminal fluorescence labeling on the molecular properties. We perform extensive molecular dynamics simulations on the labeled and unlabeled peptide and compare the results with in vitro PET-FCS measurements. Experimental and simulated timescales of end-to-end fluctuations were found in excellent agreement. Comparison between simulations with and without labels reveal that the π-stacking interaction between the fluorophore labels traps the conformation of S-peptide in a single dominant state, while the unlabeled peptide undergoes continuous conformational rearrangements. Furthermore, we find that the open to closed transition rate of S-peptide is decreased by at least one order of magnitude by the fluorophore attachment. Our approach combining experimental and in silico methods provides a benchmark for the simulations and reveals the significant effect that fluorescence labeling can have on the conformational dynamics of small biomolecules, at least for inherently flexible short peptides. The presented protocol is not only useful for comparing PET-FCS experiments with simulation results but provides a strategy to minimize the influence on molecular properties when chosing labeling positions for fluorescence experiments. [ABSTRACT FROM AUTHOR]

The critical role of BACE-1 in the formation of neurotoxic ß-amyloid peptides in the brain makes it an attractive target for an efficacious treatment of Alzheimer’s disease. However, the development of clinically useful BACE-1 inhibitors has proven to be extremely challenging. In this study we examine the binding mode of a novel potent inhibitor (compound 1, with IC50 80 nM) designed by synergistic combination of two fragments—huprine and rhein—that individually are endowed with very low activity against BACE-1. Examination of crystal structures reveals no appropriate binding site large enough to accommodate 1. Therefore we have examined the conformational flexibility of BACE-1 through extended molecular dynamics simulations, paying attention to the highly flexible region shaped by loops 8–14, 154–169 and 307–318. The analysis of the protein dynamics, together with studies of pocket druggability, has allowed us to detect the transient formation of a secondary binding site, which contains Arg307 as a key residue for the interaction with small molecules, at the edge of the catalytic cleft. The formation of this druggable “floppy” pocket would enable the binding of multisite inhibitors targeting both catalytic and secondary sites. Molecular dynamics simulations of BACE-1 bound to huprine-rhein hybrid compounds support the feasibility of this hypothesis. The results provide a basis to explain the high inhibitory potency of the two enantiomeric forms of 1, together with the large dependence on the length of the oligomethylenic linker. Furthermore, the multisite hypothesis has allowed us to rationalize the inhibitory potency of a series of tacrine-chromene hybrid compounds, specifically regarding the apparent lack of sensitivity of the inhibition constant to the chemical modifications introduced in the chromene unit. Overall, these findings pave the way for the exploration of novel functionalities in the design of optimized BACE-1 multisite inhibitors. [ABSTRACT FROM AUTHOR]

The eukaryotic group II chaperonin, the chaperonin-containing t-complex polypeptide 1 (CCT), plays an important role in cytosolic proteostasis. It has been estimated that as much as 10% of cytosolic proteins interact with CCT during their folding process. CCT is composed of 8 different paralogous subunits. Due to its complicated structure, molecular and biochemical investigations of CCT have been difficult. In this study, we constructed an expression system for CCT from a thermophilic fungus, Chaetomium thermophilum (CtCCT), by using E. coli as a host. As expected, we obtained recombinant CtCCT with a relatively high yield, and it exhibited fairly high thermal stability. We showed the advantages of the overproduction system by characterizing CtCCT variants containing ATPase-deficient subunits. For diffracted X-ray tracking experiment, we removed all surface exposed cysteine residues, and added cysteine residues at the tip of helical protrusions of selected two subunits. Gold nanocrystals were attached onto CtCCTs via gold-thiol bonds and applied for the analysis by diffracted X-ray tracking. Irrespective of the locations of cysteines, it was shown that ATP binding induces tilting motion followed by rotational motion in the CtCCT molecule, like the archaeal group II chaperonins. When gold nanocrystals were attached onto two subunits in the high ATPase activity hemisphere, the CtCCT complex exhibited a fairly rapid response to the motion. In contrast, the response of CtCCT, which had gold nanocrystals attached to the low-activity hemisphere, was slow. These results clearly support the possibility that ATP-dependent conformational change starts with the high-affinity hemisphere and progresses to the low-affinity hemisphere. [ABSTRACT FROM AUTHOR]

Dynamics of three MET antibody constructs (IgG1, IgG2, and IgG4) and the IgG4-MET antigen complex was investigated by creating their atomic models with an integrative experimental and computational approach. In particular, we used two-dimensional (2D) Electron Microscopy (EM) images, image class averaging, homology modeling, Rapidly exploring Random Tree (RRT) structure sampling, and fitting of models to images, to find the relative orientations of antibody domains that are consistent with the EM images. We revealed that the conformational preferences of the constructs depend on the extent of the hinge flexibility. We also quantified how the MET antigen impacts on the conformational dynamics of IgG4. These observations allow to create testable hypothesis to investigate MET biology. Our protocol may also help describe structural diversity of other antigen systems at approximately 5 Å precision, as quantified by Root-Mean-Square Deviation (RMSD) among good-scoring models. [ABSTRACT FROM AUTHOR]

The non-hemolytic enterotoxin complex (Nhe) is supposed to be the main virulence factor of B. cereus causing a diarrheal outcome of food poisoning. This tripartite toxin consists of the single components NheA, -B and -C all of them being necessary for maximum toxicity. In the past, research activities aiming to elucidate the mode-of-action of Nhe were mostly focused on the B- and C-component. In this study the generation of novel monoclonal antibodies (mAb) and their thorough characterization enabled the determination of key features for NheA. By the means of immunoaffinity chromatography it could be shown that NheA does not interact with -B and -C in solution. Additionally, the establishment of a highly sensitive sandwich-EIA now enables the detection of NheA in B. cereus supernatants down to 20 pg ml-1.Peptide-based epitope mapping in combination with partially deleted recombinant NheA fragments allowed the allocation of the binding regions for the three mAbs under study. Furthermore, by different EIA set-ups the conformational flexibility of NheA could be shown. For two of the antibodies under study different mechanisms of NheA neutralization were proven. Due to prevention of complete pore formation by one of the antibodies, NheA could be detected in an intermediate stage of the tripartite complex on the cell surface. Taken together, the results obtained in the present study allow a refinement of the mode-of-action for the Nhe toxin-complex. [ABSTRACT FROM AUTHOR]

Protein binding often involves conformational changes. Important questions are whether a conformational change occurs prior to a binding event (‘conformational selection’) or after a binding event (‘induced fit’), and how conformational transition rates can be obtained from experiments. In this article, we present general results for the chemical relaxation rates of conformational-selection and induced-fit binding processes that hold for all concentrations of proteins and ligands and, thus, go beyond the standard pseudo-first-order approximation of large ligand concentration. These results allow to distinguish conformational-selection from induced-fit processes—also in cases in which such a distinction is not possible under pseudo-first-order conditions—and to extract conformational transition rates of proteins from chemical relaxation data. [ABSTRACT FROM AUTHOR]

Conformational diversity of the native state plays a central role in modulating protein function. The selection paradigm sustains that different ligands shift the conformational equilibrium through their binding to highest-affinity conformers. Intramolecular vibrational dynamics associated to each conformation should guarantee conformational transitions, which due to its importance, could possibly be associated with evolutionary conserved traits. Normal mode analysis, based on a coarse-grained model of the protein, can provide the required information to explore these features. Herein, we present a novel procedure to identify key positions sustaining the conformational diversity associated to ligand binding. The method is applied to an adequate refined dataset of 188 paired protein structures in their bound and unbound forms. Firstly, normal modes most involved in the conformational change are selected according to their corresponding overlap with structural distortions introduced by ligand binding. The subspace defined by these modes is used to analyze the effect of simulated point mutations on preserving the conformational diversity of the protein. We find a negative correlation between the effects of mutations on these normal mode subspaces associated to ligand-binding and position-specific evolutionary conservations obtained from multiple sequence-structure alignments. Positions whose mutations are found to alter the most these subspaces are defined as key positions, that is, dynamically important residues that mediate the ligand-binding conformational change. These positions are shown to be evolutionary conserved, mostly buried aliphatic residues localized in regular structural regions of the protein like β-sheets and α-helix. [ABSTRACT FROM AUTHOR]

The abundant Fis nucleoid protein selectively binds poorly related DNA sequences with high affinities to regulate diverse DNA reactions. Fis binds DNA primarily through DNA backbone contacts and selects target sites by reading conformational properties of DNA sequences, most prominently intrinsic minor groove widths. High-affinity binding requires Fis-stabilized DNA conformational changes that vary depending on DNA sequence. In order to better understand the molecular basis for high affinity site recognition, we analyzed the effects of DNA sequence within and flanking the core Fis binding site on binding affinity and DNA structure. X-ray crystal structures of Fis-DNA complexes containing variable sequences in the noncontacted center of the binding site or variations within the major groove interfaces show that the DNA can adapt to the Fis dimer surface asymmetrically. We show that the presence and position of pyrimidine-purine base steps within the major groove interfaces affect both local DNA bending and minor groove compression to modulate affinities and lifetimes of Fis-DNA complexes. Sequences flanking the core binding site also modulate complex affinities, lifetimes, and the degree of local and global Fis-induced DNA bending. In particular, a G immediately upstream of the 15 bp core sequence inhibits binding and bending, and A-tracts within the flanking base pairs increase both complex lifetimes and global DNA curvatures. Taken together, our observations support a revised DNA motif specifying high-affinity Fis binding and highlight the range of conformations that Fis-bound DNA can adopt. The affinities and DNA conformations of individual Fis-DNA complexes are likely to be tailored to their context-specific biological functions. [ABSTRACT FROM AUTHOR]

Ribose-5-phosphate isomerase B from Leishmania donovani (LdRpiB) is one of the potential drug targets against visceral leishmaniasis. In the present study, we have targeted several conserved amino acids for mutational analysis (i.e. Cys69, His11, His102, His138, Asp45, Tyr46, Pro47 and Glu149) to gain crucial insights into their role in substrate binding, catalysis and conformational stability of the enzyme. All the eight LdRpiB variants were cloned, sequenced, expressed and purified. C69S, H102N, D45N and E149A mutants exhibited complete loss of enzyme activity indicating that they are indispensable for the enzyme activity. Kinetic parameters were altered in case of H138N, H11N and P47A variants; however Y46F exhibited similar kinetic behaviour as wild type. All the mutants except H138N exhibited altered protein structure as determined by CD and fluorescence spectral analysis. This data was supported by the atomic level details of the conformational changes and substrate binding using molecular dynamic simulations. LdRpiB also exhibited activity with D-form of various aldose substrates in the order of D-ribose > D-talose > D-allose > D-arabinose. Our study provides insights for better understanding of substrate enzyme interactions which can rationalize the process of drug design against parasite RpiB. [ABSTRACT FROM AUTHOR]

The mechanical properties of virus capsids correlate with local conformational dynamics in the capsid structure. They also reflect the required stability needed to withstand high internal pressures generated upon genome loading and contribute to the success of important events in viral infectivity, such as capsid maturation, genome uncoating and receptor binding. The mechanical properties of biological nanoparticles are often determined from monitoring their dynamic deformations in Atomic Force Microscopy nanoindentation experiments; but a comprehensive theory describing the full range of observed deformation behaviors has not previously been described. We present a new theory for modeling dynamic deformations of biological nanoparticles, which considers the non-linear Hertzian deformation, resulting from an indenter-particle physical contact, and the bending of curved elements (beams) modeling the particle structure. The beams’ deformation beyond the critical point triggers a dynamic transition of the particle to the collapsed state. This extreme event is accompanied by a catastrophic force drop as observed in the experimental or simulated force (F)-deformation (X) spectra. The theory interprets fine features of the spectra, including the nonlinear components of the FX-curves, in terms of the Young’s moduli for Hertzian and bending deformations, and the structural damage dependent beams’ survival probability, in terms of the maximum strength and the cooperativity parameter. The theory is exemplified by successfully describing the deformation dynamics of natural nanoparticles through comparing theoretical curves with experimental force-deformation spectra for several virus particles. This approach provides a comprehensive description of the dynamic structural transitions in biological and artificial nanoparticles, which is essential for their optimal use in nanotechnology and nanomedicine applications. [ABSTRACT FROM AUTHOR]

Kinesins are molecular motors which walk along microtubules by moving their heads to different binding sites. The motion of kinesin is realized by a conformational change in the structure of the kinesin molecule and by a diffusion of one of its two heads. In this study, a novel model is developed to account for the 2D diffusion of kinesin heads to several neighboring binding sites (near the surface of microtubules). To determine the direction of the next step of a kinesin molecule, this model considers the extension in the neck linkers of kinesin and the dynamic behavior of the coiled-coil structure of the kinesin neck. Also, the mechanical interference between kinesins and obstacles anchored on the microtubules is characterized. The model predicts that both the kinesin velocity and run length (i.e., the walking distance before detaching from the microtubule) are reduced by static obstacles. The run length is decreased more significantly by static obstacles than the velocity. Moreover, our model is able to predict the motion of kinesin when other (several) motors also move along the same microtubule. Furthermore, it suggests that the effect of mechanical interaction/interference between motors is much weaker than the effect of static obstacles. Our newly developed model can be used to address unanswered questions regarding degraded transport caused by the presence of excessive tau proteins on microtubules. [ABSTRACT FROM AUTHOR]

The Josephin Domain (JD), i.e. the N-terminal domain of Ataxin 3 (At3) protein, is an interesting example of competition between physiological function and aggregation risk. In fact, the fibrillogenesis of Ataxin 3, responsible for the spinocerebbellar ataxia 3, is strictly related to the JD thermodynamic stability. Whereas recent NMR studies have demonstrated that different JD conformations exist, the likelihood of JD achievable conformational states in solution is still an open issue. Marked differences in the available NMR models are located in the hairpin region, supporting the idea that JD has a flexible hairpin in dynamic equilibrium between open and closed states. In this work we have carried out an investigation on the JD conformational arrangement by means of both classical molecular dynamics (MD) and Metadynamics employing essential coordinates as collective variables. We provide a representation of the free energy landscape characterizing the transition pathway from a JD open-like structure to a closed-like conformation. Findings of our in silico study strongly point to the closed-like conformation as the most likely for a Josephin Domain in water. [ABSTRACT FROM AUTHOR]

Suppression of Tumorigenicity 2 (ST2), a member of the interleukin-1 receptor (IL-1R) family, activates type 2 immune responses to pathogens and tissue damage via binding to IL-33. Dysregulated responses contribute to asthma, graft-versus-host and autoinflammatory diseases and disorders. To study ST2 structure for inhibitor development, we performed the principal component (PC) analysis on the crystal structures of IL1-1R1, IL1-1R2, ST2 and the refined ST2 ectodomain (ST2ECD) models, constructed from previously reported small-angle X-ray scattering data. The analysis facilitates mapping of the ST2ECD conformations to PC subspace for characterizing structural changes. Extensive coverage of ST2ECD conformations was then obtained using the accelerated molecular dynamics simulations started with the IL-33 bound ST2ECD structure as instructed by their projected locations on the PC subspace. Cluster analysis of all conformations further determined representative conformations of ST2ECD ensemble in solution. Alignment of the representative conformations with the ST2/IL-33 structure showed that the D3 domain of ST2ECD (containing D1-D3 domains) in most conformations exhibits no clashes with IL-33 in the crystal structure. Our experimental binding data informed that the D1-D2 domain of ST2ECD contributes predominantly to the interaction between ST2ECD and IL-33 underscoring the importance of the D1-D2 domain in binding. Computational binding site assessment revealed one third of the total detected binding sites in the representative conformations may be suitable for binding to potent small molecules. Locations of these sites include the D1-D2 domain ST2ECD and modulation sites conformed to ST2ECD conformations. Our study provides structural models and analyses of ST2ECD that could be useful for inhibitor discovery. [ABSTRACT FROM AUTHOR]

Various neurodegenerative disorders are ascribed to pathogenic molecular processes involving conformational transitions of amyloidogenic proteins into toxic aggregates characterized by their β structures. Accumulating evidence indicates that neuronal cell membranes provide platforms for such conformational transitions of pathogenic proteins as best exemplified by amyloid β (Aβ). Therefore, membrane-bound Aβ species can be promising targets for the development of novel drugs for Alzheimer’s disease. In the present study, solid-state nuclear magnetic resonance spectroscopy has elucidated the membrane-induced conformation of Aβ, in which the disordered N-terminal segment is followed by the stable C-terminal β strand. The data provides an insight into the molecular processes of the conformational transition of Aβ coupled with its assembly into parallel β structures. [ABSTRACT FROM AUTHOR]

The properties of disordered proteins are thought to depend on intrinsic conformational propensities for polyproline II (PPII) structure. While intrinsic PPII propensities have been measured for the common biological amino acids in short peptides, the ability of these experimentally determined propensities to quantitatively reproduce structural behavior in intrinsically disordered proteins (IDPs) has not been established. Presented here are results from molecular simulations of disordered proteins showing that the hydrodynamic radius (Rh) can be predicted from experimental PPII propensities with good agreement, even when charge-based considerations are omitted. The simulations demonstrate that Rh and chain propensity for PPII structure are linked via a simple power-law scaling relationship, which was tested using the experimental Rh of 22 IDPs covering a wide range of peptide lengths, net charge, and sequence composition. Charge effects on Rh were found to be generally weak when compared to PPII effects on Rh. Results from this study indicate that the hydrodynamic dimensions of IDPs are evidence of considerable sequence-dependent backbone propensities for PPII structure that qualitatively, if not quantitatively, match conformational propensities measured in peptides. [ABSTRACT FROM AUTHOR]

One new cyclic hexapeptide named rubipodanin A (1), which is the first identified natural N-desmonomethyl Rubiaceae-type cyclopeptide, together with six known Rubiaceae-type cyclopeptides (2–7) were obtained using the TLC cyclopeptide protosite detection method with ninhydrin from the roots and rhizomes of Rubia podantha. The cyclopeptide structures were elucidated by extensive spectroscopic analysis, including 1D-NMR, 2D-NMR, IR, UV and MS. The solution conformation and biological activities of 1 and RA-V (4) were evaluated, and the results demonstrated that the N9-methyl group plays a vital role in the maintenance of the conformation and bioactivity. [ABSTRACT FROM AUTHOR]

Copper-free click reactions between a dibenzoazocine derivative and azides derived from 5-methyluridine were investigated. The non-catalyzed reaction yielded both regioisomers in an approximately equivalent ratio. The NMR spectra of each regioisomer revealed conformational isomery. The ratio of isomers was dependent on the type of regioisomer and the type of solvent. The synthesis of various analogs, a detailed NMR study and computational modeling provided evidence that the isomery was dependent on the interaction of the azocine and pyrimidine parts. [ABSTRACT FROM AUTHOR]

Chemokine G protein coupled receptors, principally CCR5 or CXCR4, function as co-receptors for HIV-1 entry into CD4+ T cells. Initial binding of the viral envelope glycoprotein (Env) gp120 subunit to the host CD4 receptor induces a cascade of structural conformational changes that lead to the formation of a high-affinity co-receptor-binding site on gp120. Interaction between gp120 and the co-receptor leads to the exposure of epitopes on the viral gp41 that mediates fusion between viral and cell membranes. Soluble CD4 (sCD4) mimetics can act as an activation-based inhibitor of HIV-1 entry in vitro, as it induces similar structural changes in gp120, leading to increased virus infectivity in the short term but to virus Env inactivation in the long term. Despite promising clinical implications, sCD4 displays low efficiency in vivo, and in multiple HIV strains, it does not inhibit viral infection. This has been attributed to the slow kinetics of the sCD4-induced HIV Env inactivation and to the failure to obtain sufficient sCD4 mimetic levels in the serum. Here we present uniquely structured CCR5 co-receptor mimetics. We hypothesized that such mimetics will enhance sCD4-induced HIV Env inactivation and inhibition of HIV entry. Co-receptor mimetics were derived from CCR5 gp120-binding epitopes and functionalized with a palmitoyl group, which mediated their display on the surface of lipid-coated magnetic beads. CCR5-peptidoliposome mimetics bound to soluble gp120 and inhibited HIV-1 infectivity in a sCD4-dependent manner. We concluded that CCR5-peptidoliposomes increase the efficiency of sCD4 to inhibit HIV infection by acting as bait for sCD4-primed virus, catalyzing the premature discharge of its fusion potential. [ABSTRACT FROM AUTHOR]

A sequence of complex conformational changes is required for insulin to bind to the insulin receptor. Recent experimental evidence points to the B chain C-terminal (BC-CT) as the location of these changes in insulin. Here, we present molecular dynamics simulations of insulin that reveal new insights into the structural changes occurring in the BC-CT. We find three key results: 1) The opening of the BC-CT is inherently stochastic and progresses through an open and then a “wide-open” conformation—the wide-open conformation is essential for receptor binding, but occurs only rarely. 2) The BC-CT opens with a zipper-like mechanism, with a hinge at the Phe24 residue, and is maintained in the dominant closed/inactive state by hydrophobic interactions of the neighboring Tyr26, the critical residue where opening of the BC-CT (activation of insulin) is initiated. 3) The mutation Y26N is a potential candidate as a therapeutic insulin analogue. Overall, our results suggest that the binding of insulin to its receptor is a highly dynamic and stochastic process, where initial docking occurs in an open conformation and full binding is facilitated through interactions of insulin receptor residues with insulin in its wide-open conformation. [ABSTRACT FROM AUTHOR]

Allosteric effects of mutations, ligand binding, or post-translational modifications on protein function occur through changes to the protein’s shape, or conformation. In a cell, there are many copies of the same protein, all experiencing these perturbations in a dynamic fashion and fluctuating through different conformations and activity states. According to the “conformational selection and population shift” theory, ligand binding selects a particular conformation. This perturbs the ensemble and induces a population shift. In a new PLOS Biology paper, Melacini and colleagues describe a novel model of protein regulation, the “Double-Conformational Selection Model”, which demonstrates how two tandem ligand-binding domains interact to regulate protein function. Here we explain how tandem domains with tuned interactions—but not single domains—can provide a blueprint for sensitive activation sensors within a narrow window of ligand concentration, thereby promoting signaling control. [ABSTRACT FROM AUTHOR]

Due to its inhibition of the Abl kinase domain in the BCR-ABL fusion protein, imatinib is strikingly effective in the initial stage of chronic myeloid leukemia with more than 90% of the patients showing complete remission. However, as in the case of most targeted anti-cancer therapies, the emergence of drug resistance is a serious concern. Several drug-resistant mutations affecting the catalytic domain of Abl and other tyrosine kinases are now known. But, despite their importance and the adverse effect that they have on the prognosis of the cancer patients harboring them, the molecular mechanism of these mutations is still debated. Here by using long molecular dynamics simulations and large-scale free energy calculations complemented by in vitro mutagenesis and microcalorimetry experiments, we model the effect of several widespread drug-resistant mutations of Abl. By comparing the conformational free energy landscape of the mutants with those of the wild-type tyrosine kinases we clarify their mode of action. It involves significant and complex changes in the inactive-to-active dynamics and entropy/enthalpy balance of two functional elements: the activation-loop and the conserved DFG motif. What is more the T315I gatekeeper mutant has a significant impact on the binding mechanism itself and on the binding kinetics. [ABSTRACT FROM AUTHOR]

The tp53 gene is found to be mutated in 50% of all the cancers. The p53 protein, a product of tp53 gene, is a multi-domain protein. It consists of a core DNA binding domain (DBD) which is responsible for its binding and transcription of downstream target genes. The mutations in p53 protein are responsible for creating cancerous conditions and are found to be occurring at a high frequency in the DBD region of p53. Some of these mutations are also known to be temperature sensitive (ts) in nature. They are known to exhibit partial or strong binding with DNA in the temperature range (298–306 K). Whereas, at 310 K and above they show complete loss in binding. We have analyzed the changes in binding and conformational behavior at 300 K and 310 K for three of the ts-mutants viz., V143A, R249S and R175H. QM-MM simulations have been performed on the wild type and the above mentioned ts-mutants for 30 ns each. The optimal estimate of free energy of binding for a particular number of interface hydrogen bonds was calculated using the maximum likelihood method as described by Chodera et. al (2007). This parameter has been observed to be able to mimic the binding affinity of the p53 ts-mutants at 300 K and 310 K. Thus the correlation between MM-GBSA free energy of binding and hydrogen bonds formed by the interface residues between p53 and DNA has revealed the temperature dependent nature of these mutants. The role of main chain dihedrals was obtained by performing dihedral principal component analysis (PCA). This analysis, suggests that the conformational variations in the main chain dihedrals (ϕ and ψ) of the p53 ts-mutants may have caused reduction in the overall stability of the protein. The solvent exposure of the side chains of the interface residues were found to hamper the binding of the p53 to the DNA. Solvent Accessible Surface Area (SASA) also proved to be a crucial property in distinguishing the conformers obtained at 300 K and 310 K for the three ts-mutants from the wild type at 300 K. [ABSTRACT FROM AUTHOR]

Marine bacterium Vibrio alginolyticus uses a single polar flagellum to navigate in an aqueous environment. Similar to Escherichia coli cells, the polar flagellar motor has two states; when the motor is counter-clockwise, the cell swims forward and when the motor is clockwise, the cell swims backward. V. alginolyticus also incorporates a direction randomization step at the start of the forward swimming interval by flicking its flagellum. To gain an understanding on how the polar flagellar motor switch is regulated, distributions of the forward Δf and backward Δb intervals are investigated herein. We found that the steady-state probability density functions, P(Δf) and P(Δb), of freely swimming bacteria are strongly peaked at a finite time, suggesting that the motor switch is not Poissonian. The short-time inhibition is sufficiently strong and long lasting, i.e., several hundred milliseconds for both intervals, which is readily observed and characterized. Treating motor reversal dynamics as a first-passage problem, which results from conformation fluctuations of the motor switch, we calculated P(Δf) and P(Δb) and found good agreement with the measurements. [ABSTRACT FROM AUTHOR]

Proteins are flexible, and this flexibility has an essential functional role. Flexibility can be observed in loop regions, rearrangements between secondary structure elements, and conformational changes between entire domains. However, most protein structure alignment methods treat protein structures as rigid bodies. Thus, these methods fail to identify the equivalences of residue pairs in regions with flexibility. In this study, we considered that the evolutionary relationship between proteins corresponds directly to the residue–residue physical contacts rather than the three-dimensional (3D) coordinates of proteins. Thus, we developed a new protein structure alignment method, contact area-based alignment (CAB-align), which uses the residue–residue contact area to identify regions of similarity. The main purpose of CAB-align is to identify homologous relationships at the residue level between related protein structures. The CAB-align procedure comprises two main steps: First, a rigid-body alignment method based on local and global 3D structure superposition is employed to generate a sufficient number of initial alignments. Then, iterative dynamic programming is executed to find the optimal alignment. We evaluated the performance and advantages of CAB-align based on four main points: (1) agreement with the gold standard alignment, (2) alignment quality based on an evolutionary relationship without 3D coordinate superposition, (3) consistency of the multiple alignments, and (4) classification agreement with the gold standard classification. Comparisons of CAB-align with other state-of-the-art protein structure alignment methods (TM-align, FATCAT, and DaliLite) using our benchmark dataset showed that CAB-align performed robustly in obtaining high-quality alignments and generating consistent multiple alignments with high coverage and accuracy rates, and it performed extremely well when discriminating between homologous and nonhomologous pairs of proteins in both single and multi-domain comparisons. The CAB-align software is freely available to academic users as stand-alone software at . [ABSTRACT FROM AUTHOR]

Proteins form large macromolecular assemblies with RNA that govern essential molecular processes. RNA-binding proteins have often been associated with conformational flexibility, yet the extent and functional implications of their intrinsic disorder have never been fully assessed. Here, through large-scale analysis of comprehensive protein sequence and structure datasets we demonstrate the prevalence of intrinsic structural disorder in RNA-binding proteins and domains. We addressed their functionality through a quantitative description of the evolutionary conservation of disordered segments involved in binding, and investigated the structural implications of flexibility in terms of conformational stability and interface formation. We conclude that the functional role of intrinsically disordered protein segments in RNA-binding is two-fold: first, these regions establish extended, conserved electrostatic interfaces with RNAs via induced fit. Second, conformational flexibility enables them to target different RNA partners, providing multi-functionality, while also ensuring specificity. These findings emphasize the functional importance of intrinsically disordered regions in RNA-binding proteins. [ABSTRACT FROM AUTHOR]

The transcription factor p53 is the main tumour suppressor in cells and many cancer types have p53 mutations resulting in a loss of its function. In tumours that retain wild-type p53 function, p53 activity is down-regulated by MDM2 (human murine double minute 2) via a direct protein—protein interaction. We have designed and synthesised two series of 2,5-diketopiperazines as inhibitors of the MDM2-p53 interaction. The first set was designed to directly mimic the α-helical region of the p53 peptide, containing key residues in the i, i+4 and i+7 positions of a natural α-helix. Conformational analysis indicated that 1,3,6-trisubstituted 2,5-diketopiperazines were able to place substituents in the same spatial orientation as an α-helix template. The key step of the synthesis involved the cyclisation of substituted dipeptides. The other set of tetrasubstituted 2,5-diketopiperazines were designed based on structure-based docking studies and the Ugi multicomponent reaction was used for the synthesis. This latter set comprised the most potent inhibitors which displayed micromolar IC50-values in a biochemical fluorescence polarisation assay. [ABSTRACT FROM AUTHOR]

The extent of glycation and conformational changes of horse myoglobin (Mb) upon glycation with N-acetyl-glucosamine (GlcNAc), glucose (Glc) and glucosamine (GlcN) were investigated. Among tested sugars, the rate of glycation with GlcN was the most rapid as shown by MALDI and ESI mass spectrometries. Protein oxidation, as evaluated by the amount of carbonyl groups present on Mb, was found to increase exponentially in Mb-Glc conjugates over time, whereas in Mb-GlcN mixtures the carbonyl groups decreased significantly after maximum at 3 days of the reaction. The reaction between GlcN and Mb resulted in a significantly higher amount of α-dicarbonyl compounds, mostly glucosone and 3-deoxyglucosone, ranging from and 27 to 332 mg/L and from 14 to 304 mg/L, respectively. Already at 0.5 days, tertiary structural changes of Mb-GlcN conjugate were observed by altered tryptophan fluorescence. A reduction of metmyoglobin to deoxy-and oxymyoglobin forms was observed on the first day of reaction, coinciding with the greatest amount of glucosone produced. In contrast to native α-helical myoglobin, 41% of the glycated protein sequence was transformed into a β-sheet conformation, as determined by circular dichroism spectropolarimetry. Transmission electron microscopy demonstrated that Mb glycation with GlcN causes the formation of amorphous or fibrous aggregates, started already at 3 reaction days. These aggregates bind to an amyloid-specific dye thioflavin T. With the aid of α-dicarbonyl compounds and advanced products of reaction, this study suggests that the Mb glycation with GlcN induces the unfolding of an initially globular protein structure into amyloid fibrils comprised of a β-sheet structure. [ABSTRACT FROM AUTHOR]

Voltage-gated potassium (K+) channels are present in all living systems. Despite high structural similarities in the transmembrane domains (TMD), this K+ channel type segregates into at least two main functional categories—hyperpolarization-activated, inward-rectifying (Kin) and depolarization-activated, outward-rectifying (Kout) channels. Voltage-gated K+ channels sense the membrane voltage via a voltage-sensing domain that is connected to the conduction pathway of the channel. It has been shown that the voltage-sensing mechanism is the same in Kin and Kout channels, but its performance results in opposite pore conformations. It is not known how the different coupling of voltage-sensor and pore is implemented. Here, we studied sequence and structural data of voltage-gated K+ channels from animals and plants with emphasis on the property of opposite rectification. We identified structural hotspots that alone allow already the distinction between Kin and Kout channels. Among them is a loop between TMD S5 and the pore that is very short in animal Kout, longer in plant and animal Kin and the longest in plant Kout channels. In combination with further structural and phylogenetic analyses this finding suggests that outward-rectification evolved twice and independently in the animal and plant kingdom. [ABSTRACT FROM AUTHOR]

Intramolecular mobility and conformational changes of flexible loops have important roles in the structural and functional integrity of proteins. The Achaetomium sp. Xz8 endo-polygalacturonase (PG8fn) of glycoside hydrolase (GH) family 28 is distinguished for its high catalytic activity (28,000 U/mg). Structure modeling indicated that PG8fn has a flexible T3 loop that folds partly above the substrate in the active site, and forms a hydrogen bond to the substrate by a highly conserved residue Asn94 in the active site cleft. Our research investigates the catalytic roles of Asn94 in T3 loop which is located above the catalytic residues on one side of the substrate. Molecular dynamics simulation performed on the mutant N94A revealed the loss of the hydrogen bond formed by the hydroxyl group at O34 of pentagalacturonic acid and the crucial ND2 of Asn94 and the consequent detachment and rotation of the substrate away from the active site, and that on N94Q caused the substrate to drift away from its place due to the longer side chain. In line with the simulations, site-directed mutagenesis at this site showed that this position is very sensitive to amino acid substitutions. Except for the altered Km values from 0.32 (wild type PG8fn) to 0.75–4.74 mg/ml, all mutants displayed remarkably lowered kcat (~3–20,000 fold) and kcat/Km (~8–187,500 fold) values and significantly increased △(△G) values (5.92–33.47 kJ/mol). Taken together, Asn94 in the GH28 T3 loop has a critical role in positioning the substrate in a correct way close to the catalytic residues. [ABSTRACT FROM AUTHOR]

An important goal in molecular biology is to understand functional changes upon single-point mutations in proteins. Doing so through a detailed characterization of structure spaces and underlying energy landscapes is desirable but continues to challenge methods based on Molecular Dynamics. In this paper we propose a novel algorithm, SIfTER, which is based instead on stochastic optimization to circumvent the computational challenge of exploring the breadth of a protein’s structure space. SIfTER is a data-driven evolutionary algorithm, leveraging experimentally-available structures of wildtype and variant sequences of a protein to define a reduced search space from where to efficiently draw samples corresponding to novel structures not directly observed in the wet laboratory. The main advantage of SIfTER is its ability to rapidly generate conformational ensembles, thus allowing mapping and juxtaposing landscapes of variant sequences and relating observed differences to functional changes. We apply SIfTER to variant sequences of the H-Ras catalytic domain, due to the prominent role of the Ras protein in signaling pathways that control cell proliferation, its well-studied conformational switching, and abundance of documented mutations in several human tumors. Many Ras mutations are oncogenic, but detailed energy landscapes have not been reported until now. Analysis of SIfTER-computed energy landscapes for the wildtype and two oncogenic variants, G12V and Q61L, suggests that these mutations cause constitutive activation through two different mechanisms. G12V directly affects binding specificity while leaving the energy landscape largely unchanged, whereas Q61L has pronounced, starker effects on the landscape. An implementation of SIfTER is made available at . We believe SIfTER is useful to the community to answer the question of how sequence mutations affect the function of a protein, when there is an abundance of experimental structures that can be exploited to reconstruct an energy landscape that would be computationally impractical to do via Molecular Dynamics. [ABSTRACT FROM AUTHOR]

HIV virus-like particles (VLPs) present the HIV envelope protein in its native conformation, providing an ideal vaccine antigen. To enhance the immunogenicity of the VLP vaccine, we sought to improve upon two components; the route of administration and the additional adjuvant. Using HIV VLPs, we evaluated sub-cheek as a novel route of vaccine administration when combined with other conventional routes of immunization. Of five combinations of distinct prime and boost sequences, which included sub-cheek, intranasal, and intradermal routes of administration, intranasal prime and sub-cheek boost (IN+SC) resulted in the highest HIV-specific IgG titers among the groups tested. Using the IN+SC regimen we tested the adjuvant VesiVax Conjugatable Adjuvant Lipid Vesicles (CALV) + monophosphoryl lipid A (MPLA) at MPLA concentrations of 0, 7.5, 12.5, and 25 μg/dose in combination with our VLPs. Mice that received 12.5 or 25 μg/dose MPLA had the highest concentrations of Env-specific IgG2c (20.7 and 18.4 μg/ml respectively), which represents a Th1 type of immune response in C57BL/6 mice. This was in sharp contrast to mice which received 0 or 7.5 μg MPLA adjuvant (6.05 and 5.68 μg/ml of IgG2c respectively). In contrast to IgG2c, MPLA had minor effects on Env-specific IgG1; therefore, 12.5 and 25 μg/dose of MPLA induced the optimal IgG1/IgG2c ratio of 1.3. Additionally, the percentage of germinal center B cells increased significantly from 15.4% in the control group to 31.9% in the CALV + 25 μg MPLA group. These mice also had significantly more IL-2 and less IL-4 Env-specific CD8+ T cells than controls, correlating with an increased percentage of Env-specific central memory CD4+ and CD8+ T cells. Our study shows the strong potential of IN+SC as an efficacious route of administration and the effectiveness of VLPs combined with MPLA adjuvant to induce Env specific Th1-oriented HIV-specific immune responses. [ABSTRACT FROM AUTHOR]

Human retinal macular pigment (MP) is formed by the carotenoids lutein and zeaxanthin (including the isomer meso-zeaxanthin). MP has several functions in improving visual performance and protecting against the damaging effects of light, and MP levels are used as a proxy for macular health–specifically, to predict the likelihood of developing age-related macular degeneration. While the roles of these carotenoids in retinal health have been the object of intense study in recent years, precise mechanistic details of their protective action remain elusive. We have measured the Raman signals originating from MP carotenoids in ex vivo human retinal tissue, in order to assess their structure and conformation. We show that it is possible to distinguish between lutein and zeaxanthin, by their excitation profile (related to their absorption spectra) and the position of their ν1 Raman mode. In addition, analysis of the ν4 Raman band indicates that these carotenoids are present in a specific, constrained conformation in situ, consistent with their binding to specific proteins as postulated in the literature. We discuss how these conclusions relate to the function of these pigments in macular protection. We also address the possibilities for a more accurate, consistent measurement of MP levels by Raman spectroscopy. [ABSTRACT FROM AUTHOR]