An automated approach for defining core atoms and domains in an ensemble of NMR-derived protein structures

Stepwise activation of the pro-apoptotic protein Bid at mitochondrial membranes

Caspase-8-cleaved Bid (cBid) associates with mitochondria and promotes the activation of BAX, leading to mitochondria outer membrane permeabilization (MOMP) and apoptosis. However, current structural models of cBid are largely based on studies using membrane vesicles and detergent micelles. Here we employ spin-label ESR and site-directed PEGylation methods to identify conformations of cBid at real mitochondrial membranes, revealing stepwise mechanisms in the activation process. Upon the binding of cBid to mitochondria, its structure is reorganized to expose the BH3 domain while leaving the structural integrity only slightly altered. The mitochondria-bound cBid is in association with Mtch2 and it remains in the primed state until interacting with BAX. The interaction subsequently triggers the fragmentation of cBid, causes large conformational changes, and promotes BAX-mediated MOMP. Our results reveal structural differences of cBid between mitochondria and other lipid-like environments and, moreover, highlight the role of the membrane binding in modifying cBid structure and assisting the inactive-to-active transition in function.

Patterns in Protein Flexibility: A Comparison of NMR "Ensembles", MD Trajectories, and Crystallographic B-Factors

Proteins are molecular machines requiring flexibility to function. Crystallographic B-factors and Molecular Dynamics (MD) simulations both provide insights into protein flexibility on an atomic scale. Nuclear Magnetic Resonance (NMR) lacks a universally accepted analog of the B-factor. However, a lack of convergence in atomic coordinates in an NMR-based structure calculation also suggests atomic mobility. This paper describes a pattern in the coordinate uncertainties of backbone heavy atoms in NMR-derived structural “ensembles” first noted in the development of FindCore2 (previously called Expanded FindCore: DA Snyder, J Grullon, YJ Huang, R Tejero, GT Montelione, Proteins: Structure, Function, and Bioinformatics 82 (S2), 219–230) and demonstrates that this pattern exists in coordinate variances across MD trajectories but not in crystallographic B-factors. This either suggests that MD trajectories and NMR “ensembles” capture motional behavior of peptide bond units not captured by B-factors or indicates a deficiency common to force fields used in both NMR and MD calculations.

Low Basicity as a Characteristic for Atypical Ligands of Serotonin Receptor 5-HT2

Serotonin receptors are extensively examined by academic and industrial researchers, due to their vital roles, which they play in the organism and constituting therefore important drug targets. Up to very recently, it was assumed that the basic nitrogen in compound structure is a necessary component to make it active within this receptor system. Such nitrogen interacts in its protonated form with the aspartic acid from the third transmembrane helix (D3x32) forming a hydrogen bond tightly fitting the ligand in the protein binding site. However, there are several recent studies that report strong serotonin receptor affinity also for compounds without a basic moiety in their structures. In the study, we carried out a comprehensive in silico analysis of the low-basicity phenomenon of the selected serotonin receptor ligands. We focused on the crystallized representatives of the proteins of 5-HT_1B, 5-HT_2A, 5-HT_2B, and 5-HT_2C receptors, and examined the problem both from the ligand- and structure-based perspectives. The study was performed for the native proteins, and for D3x32A mutants. The investigation resulted in the determination of nonstandard structural requirements for activity towards serotonin receptors, which can be used in the design of new nonbasic ligands.

Modern Technologies of Solution Nuclear Magnetic Resonance Spectroscopy for Three-dimensional Structure Determination of Proteins Open Avenues for Life Scientists

Nuclear magnetic resonance (NMR) spectroscopy is a powerful technique for structural studies of chemical compounds and biomolecules such as DNA and proteins. Since the NMR signal sensitively reflects the chemical environment and the dynamics of a nuclear spin, NMR experiments provide a wealth of structural and dynamic information about the molecule of interest at atomic resolution. In general, structural biology studies using NMR spectroscopy still require a reasonable understanding of the theory behind the technique and experience on how to recorded NMR data. Owing to the remarkable progress in the past decade, we can easily access suitable and popular analytical resources for NMR structure determination of proteins with high accuracy. Here, we describe the practical aspects, workflow and key points of modern NMR techniques used for solution structure determination of proteins. This review should aid NMR specialists aiming to develop new methods that accelerate the structure determination process, and open avenues for non-specialist and life scientists interested in using NMR spectroscopy to solve protein structures.

Strain analysis of protein structures and low dimensionality of mechanical allosteric couplings

In many proteins, especially allosteric proteins that communicate regulatory states from allosteric to active sites, structural deformations are functionally important. To understand these deformations, dynamical experiments are ideal but challenging. Using static structural information, although more limited than dynamical analysis, is much more accessible. Underused for protein analysis, strain is the natural quantity for studying local deformations. We calculate strain tensor fields for proteins deformed by ligands or thermal fluctuations using crystal and NMR structure ensembles. Strains-primarily shears-show deformations around binding sites. These deformations can be induced solely by ligand binding at distant allosteric sites. Shears reveal quasi-2D paths of mechanical coupling between allosteric and active sites that may constitute a widespread mechanism of allostery. We argue that strain-particularly shear-is the most appropriate quantity for analysis of local protein deformations. This analysis can reveal mechanical and biological properties of many proteins.

The Art of Compiling Protein Binding Site Ensembles

Structure-based drug design starts with the collection, preparation, and initial analysis of protein structures. With more than 115,000 structures publically available in the Protein Data Bank (PDB), fully automated processes reliably performing these important preprocessing steps are needed. Several tools are available for these tasks, however, most of them do not address the special needs of scientists interested in protein-ligand interactions. In this paper, we summarize our research activities towards an automated processing pipeline from raw PDB data towards ready-to-use protein binding site ensembles. Starting from a single protein structure, the pipeline covers the following phases: Extracting structurally related binding sites from the PDB, aligning disconnected binding site sequences, resolving tautomeric forms and protonation, orienting hydrogens and flippable side-chains, structurally aligning the multitude of binding sites, and performing a reasonable reduction of ensemble structures. The pipeline, named SIENA, creates protein-structural ensembles for the analysis of protein flexibility, molecular design efforts like docking or de novo design within seconds. For the first time, we are able to process the whole PDB in order to create a large collection of protein binding site ensembles. SIENA is available as part of the ZBH ProteinsPlus webserver under http://proteinsplus.zbh.uni-hamburg.de.

SIENA: Efficient Compilation of Selective Protein Binding Site Ensembles

Structural flexibility of proteins has an important influence on molecular recognition and enzymatic function. In modeling, structure ensembles are therefore often applied as a valuable source of alternative protein conformations. However, their usage is often complicated by structural artifacts and inconsistent data annotation. Here, we present SIENA, a new computational approach for the automated assembly and preprocessing of protein binding site ensembles. Starting with an arbitrarily defined binding site in a single protein structure, SIENA searches for alternative conformations of the same or sequentially closely related binding sites. The method is based on an indexed database for identifying perfect k-mer matches and a recently published algorithm for the alignment of protein binding site conformations. Furthermore, SIENA provides a new algorithm for the interaction-based selection of binding site conformations which aims at covering all known ligand-binding geometries. Various experiments highlight that SIENA is able to generate comprehensive and well selected binding site ensembles improving the compatibility to both known and unconsidered ligand molecules. Starting with the whole PDB as data source, the computation time of the whole ensemble generation takes only a few seconds. SIENA is available via a Web service at www.zbh.uni-hamburg.de/siena .

Inexpensive Method for Selecting Receptor Structures for Virtual Screening

This article introduces a screening performance index (SPI) to help select from a number of experimental structures one or a few that are more likely to identify more actives among its top hits from virtual screening of a compound library. It achieved this by docking only known actives to the experimental structures without considering a large number of decoys to reduce computational costs. The SPI is calculated by using the docking energies of the actives to all the receptor structures. We evaluated the performance of the SPI by applying it to study eight protein systems: fatty acid binding protein adipocyte FABP4, serine/threonine-protein kinase BRAF, beta-1 adrenergic receptor ADRB1, TGF-beta receptor type I TGFR1, adenosylhomocysteinase SAHH, thyroid hormone receptor beta-1 THB, phospholipase A2 group IIA PA2GA, and cytochrome P450 3a4 CP3A4. We found that the SPI agreed with the results from other popular performance metrics such as Boltzmann-Enhanced Discrimination Receiver Operator Characteristics (BEDROC), Robust Initial Enhancement (RIE), Area Under Accumulation Curve (AUAC), and Enrichment Factor (EF) but is less expensive to calculate. SPI also performed better than the best docking energy, the molecular volume of the bound ligand, and the resolution of crystal structure in selecting good receptor structures for virtual screening. The implications of these findings were further discussed in the context of ensemble docking, in situations when no experimental structure for the targeted protein was available, or under circumstances when quick choices of receptor structures need to be made before quantitative indexes such as the SPI and BEDROC can be calculated.

BAX-induced apoptosis can be initiated through a conformational selection mechanism

BAX protein plays a key role in the mitochondria-mediated apoptosis. However, it remains unclear by what mechanism BAX is triggered to initiate apoptosis. Here, we reveal the mechanism using electron spin resonance (ESR) techniques. An inactive BAX monomer was found to exhibit conformational heterogeneity and exist at equilibrium in two conformations, one of which has never been reported. We show that upon apoptotic stimulus by BH3-only peptides, BAX can be induced to convert into either a ligand-bound monomer or an oligomer through a conformational selection mechanism. The kinetics of reaction is studied by means of time-resolved ESR, allowing a direct in situ observation for the transformation of BAX from the native to the bound states. In vitro mitochondrial assays provide further discrimination between the proposed BAX states, thereby revealing a population-shift allosteric mechanism in the process. BAX's apoptotic function is shown to critically depend on excursions between different structural conformations.

The expanded FindCore method for identification of a core atom set for assessment of protein structure prediction

Maximizing the scientific impact of NMR-based structure determination requires robust and statistically sound methods for assessing the precision of NMR-derived structures. In particular, a method to define a core atom set for calculating superimpositions and validating structure predictions is critical to the use of NMR-derived structures as targets in the CASP competition. FindCore (Snyder and Montelione, Proteins 2005;59:673-686) is a superimposition independent method for identifying a core atom set and partitioning that set into domains. However, as FindCore optimizes superimposition by sensitively excluding not-well-defined atoms, the FindCore core may not comprise all atoms suitable for use in certain applications of NMR structures, including the CASP assessment process. Adapting the FindCore approach to assess predicted models against experimental NMR structures in CASP10 required modification of the FindCore method. This paper describes conventions and a standard protocol to calculate an "Expanded FindCore" atom set suitable for validation and application in biological and biophysical contexts. A key application of the Expanded FindCore method is to identify a core set of atoms in the experimental NMR structure for which it makes sense to validate predicted protein structure models. We demonstrate the application of this Expanded FindCore method in characterizing well-defined regions of 18 NMR-derived CASP10 target structures. The Expanded FindCore protocol defines "expanded core atom sets" that match an expert's intuition of which parts of the structure are sufficiently well defined to use in assessing CASP model predictions. We also illustrate the impact of this analysis on the CASP GDT assessment scores.

Crystallographic phasing with NMR models: an envelope approach

X-ray crystallography and NMR are complementary tools in structural biology. However, it is often difficult to use NMR structures as search models in molecular replacement (MR) to phase crystallographic data. In this study, a new approach is reported utilizing a molecular envelope of NMR structures for MR phasing with the programFSEARCHat low resolution (about 6 Å). Several targets with both crystallographic and NMR structures available have been tested.FSEARCHwas able to find the correct translation and orientation of the search model in the crystallographic unit cell, while conventional MR procedures were unsuccessful.

An overview of tools for the validation of protein NMR structures

Biomolecular structures at atomic resolution present a valuable resource for the understanding of biology. NMR spectroscopy accounts for 11% of all structures in the PDB repository. In response to serious problems with the accuracy of some of the NMR-derived structures and in order to facilitate proper analysis of the experimental models, a number of program suites are available. We discuss nine of these tools in this review: PROCHECK-NMR, PSVS, GLM-RMSD, CING, Molprobity, Vivaldi, ResProx, NMR constraints analyzer and QMEAN. We evaluate these programs for their ability to assess the structural quality, restraints and their violations, chemical shifts, peaks and the handling of multi-model NMR ensembles. We document both the input required by the programs and output they generate. To discuss their relative merits we have applied the tools to two representative examples from the PDB: a small, globular monomeric protein (Staphylococcal nuclease from S. aureus, PDB entry 2kq3) and a small, symmetric homodimeric protein (a region of human myosin-X, PDB entry 2lw9).

Recommendations of the wwPDB NMR Validation Task Force

As methods for analysis of biomolecular structure and dynamics using nuclear magnetic resonance spectroscopy (NMR) continue to advance, the resulting 3D structures, chemical shifts, and other NMR data are broadly impacting biology, chemistry, and medicine. Structure model assessment is a critical area of NMR methods development, and is an essential component of the process of making these structures accessible and useful to the wider scientific community. For these reasons, the Worldwide Protein Data Bank (wwPDB) has convened an NMR Validation Task Force (NMR-VTF) to work with wwPDB partners in developing metrics and policies for biomolecular NMR data harvesting, structure representation, and structure quality assessment. This paper summarizes the recommendations of the NMR-VTF, and lays the groundwork for future work in developing standards and metrics for biomolecular NMR structure quality assessment.

A robust method for quantitative identification of ordered cores in an ensemble of biomolecular structures by non-linear multi-dimensional scaling using inter-atomic distance variance matrix

Superpositioning of atoms in an ensemble of biomolecules is a common task in a variety of fields in structural biology. Although several automated tools exist based on previously established methods, manual operations to define the atoms in the ordered regions are usually preferred. The task is difficult and lacks output efficiency for multi-core proteins having complicated folding topology. The new method presented here can systematically and quantitatively achieve the identification of ordered cores even for molecules containing multiple cores linked with flexible loops. In contrast to established methods, this method treats the variance of inter-atomic distances in an ensemble as information content using a non-linear (NL) function, and then subjects it to multi-dimensional scaling (MDS) to embed the row vectors in the inter-atomic distance variance matrix into a lower dimensional matrix. The plots of the identified atom groups in a one or two-dimensional map enables users to visually and intuitively infer well-ordered atoms in an ensemble, as well as to automatically identify them by the standard clustering methods. The performance of the NL-MDS method has been examined for number of structure ensembles studied by nuclear magnetic resonance, demonstrating that the method can be more suitable for structural analysis of multi-core proteins in comparison to previously established methods.

Quality assessment of protein NMR structures

Biomolecular NMR structures are now routinely used in biology, chemistry, and bioinformatics. Methods and metrics for assessing the accuracy and precision of protein NMR structures are beginning to be standardized across the biological NMR community. These include both knowledge-based assessment metrics, parameterized from the database of protein structures, and model versus data assessment metrics. On line servers are available that provide comprehensive protein structure quality assessment reports, and efforts are in progress by the world-wide Protein Data Bank (wwPDB) to develop a biomolecular NMR structure quality assessment pipeline as part of the structure deposition process. These quality assessment metrics and standards will aid NMR spectroscopists in determining more accurate structures, and increase the value and utility of these structures for the broad scientific community.

Ab initio modeling and experimental assessment of Janus Kinase 2 (JAK2) kinase-pseudokinase complex structure

The Janus Kinase 2 (JAK2) plays essential roles in transmitting signals from multiple cytokine receptors, and constitutive activation of JAK2 results in hematopoietic disorders and oncogenesis. JAK2 kinase activity is negatively regulated by its pseudokinase domain (JH2), where the gain-of-function mutation V617F that causes myeloproliferative neoplasms resides. In the absence of a crystal structure of full-length JAK2, how JH2 inhibits the kinase domain (JH1), and how V617F hyperactivates JAK2 remain elusive. We modeled the JAK2 JH1-JH2 complex structure using a novel informatics-guided protein-protein docking strategy. A detailed JAK2 JH2-mediated auto-inhibition mechanism is proposed, where JH2 traps the activation loop of JH1 in an inactive conformation and blocks the movement of kinase αC helix through critical hydrophobic contacts and extensive electrostatic interactions. These stabilizing interactions are less favorable in JAK2-V617F. Notably, several predicted binding interfacial residues in JH2 were confirmed to hyperactivate JAK2 kinase activity in site-directed mutagenesis and BaF3/EpoR cell transformation studies. Although there may exist other JH2-mediated mechanisms to control JH1, our JH1-JH2 structural model represents a verifiable working hypothesis for further experimental studies to elucidate the role of JH2 in regulating JAK2 in both normal and pathological settings.

Four complete turns of a curved 3₁₀-helix at atomic resolution: the crystal structure of the peptaibol trichovirin I-4A in a polar environment suggests a transition to α-helix for membrane function

The first crystal structure of a member of peptaibol antibiotic subfamily 4, trichovirin I-4A (14 residues), has been determined by direct methods and refined at atomic resolution. The monoclinic unit cell has two molecules in the asymmetric unit. Both molecules assume a 3₁₀ right-handed helical conformation and are significantly bent. The molecules pack loosely along the crystallographic twofold axis, forming two large tunnels between symmetry-related molecules in which no ordered solvent could be located. Carbonyl O atoms which are not involved in intramolecular hydrogen bonding participate in close van der Waals interactions with apolar groups. The necessary amphipathicity for biological activity of peptaibols is not realised in the crystal structure. Hence, a structural change of trichovirin to an α-helical conformation is proposed for membrane integration and efficient water/ion transportation across the lipid bilayer.

Application of conformational clustering in protein-ligand docking

Protein-Ligand docking is a powerful technique routinely employed in structure-based drug design. Despite many reported success stories, docking is not always able to provide an accurate and easily interpretable prediction of the structure of the bound complex formed by a small organic molecule and a pharmacologically relevant target. Cluster analysis can represent a versatile and readily available postprocessing tool to be employed in combination with protein-ligand docking to simplify the evaluation of the results and help to overcome present limitations of docking protocols.

Improved technologies now routinely provide protein NMR structures useful for molecular replacement

Molecular replacement (MR) is widely used for addressing the phase problem in X-ray crystallography. Historically, crystallographers have had limited success using NMR structures as MR search models. Here, we report a comprehensive investigation of the utility of protein NMR ensembles as MR search models, using data for 25 pairs of X-ray and NMR structures solved and refined using modern NMR methods. Starting from NMR ensembles prepared by an improved protocol, FindCore, correct MR solutions were obtained for 22 targets. Based on these solutions, automatic model rebuilding could be done successfully. Rosetta refinement of NMR structures provided MR solutions for another two proteins. We also demonstrate that such properly prepared NMR ensembles and X-ray crystal structures have similar performance when used as MR search models for homologous structures, particularly for targets with sequence identity >40%.

Objective identification of residue ranges for the superposition of protein structures

Background

The automation of objectively selecting amino acid residue ranges for structure superpositions is important for meaningful and consistent protein structure analyses. So far there is no widely-used standard for choosing these residue ranges for experimentally determined protein structures, where the manual selection of residue ranges or the use of suboptimal criteria remain commonplace.

Results

We present an automated and objective method for finding amino acid residue ranges for the superposition and analysis of protein structures, in particular for structure bundles resulting from NMR structure calculations. The method is implemented in an algorithm, CYRANGE, that yields, without protein-specific parameter adjustment, appropriate residue ranges in most commonly occurring situations, including low-precision structure bundles, multi-domain proteins, symmetric multimers, and protein complexes. Residue ranges are chosen to comprise as many residues of a protein domain that increasing their number would lead to a steep rise in the RMSD value. Residue ranges are determined by first clustering residues into domains based on the distance variance matrix, and then refining for each domain the initial choice of residues by excluding residues one by one until the relative decrease of the RMSD value becomes insignificant. A penalty for the opening of gaps favours contiguous residue ranges in order to obtain a result that is as simple as possible, but not simpler. Results are given for a set of 37 proteins and compared with those of commonly used protein structure validation packages. We also provide residue ranges for 6351 NMR structures in the Protein Data Bank.

Conclusions

The CYRANGE method is capable of automatically determining residue ranges for the superposition of protein structure bundles for a large variety of protein structures. The method correctly identifies ordered regions. Global structure superpositions based on the CYRANGE residue ranges allow a clear presentation of the structure, and unnecessary small gaps within the selected ranges are absent. In the majority of cases, the residue ranges from CYRANGE contain fewer gaps and cover considerably larger parts of the sequence than those from other methods without significantly increasing the RMSD values. CYRANGE thus provides an objective and automatic method for standardizing the choice of residue ranges for the superposition of protein structures.

X-ray vs. NMR structures as templates for computational protein design

Certain protein-design calculations involve using an experimentally determined high-resolution structure as a template to identify new sequences that can adopt the same fold. This approach has led to the successful design of many novel, well-folded, native-like proteins. Although any atomic-resolution structure can serve as a template in such calculations, most successful designs have used high-resolution crystal structures. Because there are many proteins for which crystal structures are not available, it is of interest whether nuclear magnetic resonance (NMR) templates are also appropriate. We have analyzed differences between using X-ray and NMR templates in side-chain repacking and design calculations. We assembled a database of 29 proteins for which both a high-resolution X-ray structure and an ensemble of NMR structures are available. Using these pairs, we compared the rotamericity, chi(1)-angle recovery, and native-sequence recovery of X-ray and NMR templates. We carried out design using RosettaDesign on both types of templates, and compared the energies and packing qualities of the resulting structures. Overall, the X-ray structures were better templates for use with Rosetta. However, for approximately 20% of proteins, a member of the reported NMR ensemble gave rise to designs with similar properties. Re-evaluating RosettaDesign structures with other energy functions indicated much smaller differences between the two types of templates. Ultimately, experiments are required to confirm the utility of particular X-ray and NMR templates. But our data suggest that the lack of a high-resolution X-ray structure should not preclude attempts at computational design if an NMR ensemble is available.

PDBe: Protein Data Bank in Europe

The Protein Data Bank in Europe (PDBe) (http://www.ebi.ac.uk/pdbe/) is actively working with its Worldwide Protein Data Bank partners to enhance the quality and consistency of the international archive of bio-macromolecular structure data, the Protein Data Bank (PDB). PDBe also works closely with its collaborators at the European Bioinformatics Institute and the scientific community around the world to enhance its databases and services by adding curated and actively maintained derived data to the existing structural data in the PDB. We have developed a new database infrastructure based on the remediated PDB archive data and a specially designed database for storing information on interactions between proteins and bound molecules. The group has developed new services that allow users to carry out simple textual queries or more complex 3D structure-based queries. The newly designed 'PDBeView Atlas pages' provide an overview of an individual PDB entry in a user-friendly layout and serve as a starting point to further explore the information available in the PDBe database. PDBe's active involvement with the X-ray crystallography, Nuclear Magnetic Resonance spectroscopy and cryo-Electron Microscopy communities have resulted in improved tools for structure deposition and analysis.

Alpha7 nicotinic acetylcholine receptor agonists: prediction of their binding affinity through a molecular mechanics Poisson-Boltzmann surface area approach

A group of agonists for the alpha7 neuronal nicotinic acetylcholine receptors (nAChRs) was investigated, and their free energies of binding DeltaG(bind) were calculated by applying the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) approach. This method, based on molecular dynamics simulations of fully solvated protein-ligand complexes, allowed us to estimate the contribution of both polar and nonpolar terms as well as the entropy to the overall free energy of binding. The calculated results were in a good agreement with the experimentally determined DeltaG(bind) values, thereby pointing to the MM-PBSA protocol as a valuable computational tool for the rational design of specific agents targeting the neuronal alpha7 nAChR subtypes.

Mixture models for protein structure ensembles

MOTIVATION

Protein structure ensembles provide important insight into the dynamics and function of a protein and contain information that is not captured with a single static structure. However, it is not clear a priori to what extent the variability within an ensemble is caused by internal structural changes. Additional variability results from overall translations and rotations of the molecule. And most experimental data do not provide information to relate the structures to a common reference frame. To report meaningful values of intrinsic dynamics, structural precision, conformational entropy, etc., it is therefore important to disentangle local from global conformational heterogeneity.

RESULTS

We consider the task of disentangling local from global heterogeneity as an inference problem. We use probabilistic methods to infer from the protein ensemble missing information on reference frames and stable conformational sub-states. To this end, we model a protein ensemble as a mixture of Gaussian probability distributions of either entire conformations or structural segments. We learn these models from a protein ensemble using the expectation-maximization algorithm. Our first model can be used to find multiple conformers in a structure ensemble. The second model partitions the protein chain into locally stable structural segments or core elements and less structured regions typically found in loops. Both models are simple to implement and contain only a single free parameter: the number of conformers or structural segments. Our models can be used to analyse experimental ensembles, molecular dynamics trajectories and conformational change in proteins.

AVAILABILITY

The Python source code for protein ensemble analysis is available from the authors upon request.

Characterization of the binding properties of SIRT2 inhibitors with a N-(3-phenylpropenoyl)-glycine tryptamide backbone

SIRT2 inhibitors with a N-(3-phenylpropenoyl)-glycine tryptamide backbone were studied. This backbone has been developed in our group, and it is derived from a compound originally found by virtual screening. In addition, compounds with a smaller 3-phenylpropenoic acid tryptamide backbone were also included in the study. Binding modes for the new compounds and the previously reported compounds were analyzed with molecular modelling methods. The approach, which included a combination of molecular dynamics, molecular docking and cluster analysis, showed that certain docking poses were favourable despite the conformational variation in the target protein. The N-(3-phenylpropenoyl)-glycine tryptamide backbone is also a good backbone for SIRT2 inhibitors, and the series of compounds includes several potent SIRT2 inhibitors.

Conformational analysis of alternative protein structures

MOTIVATION

Alternative structural models determined experimentally are available for an increasing number of proteins. Structural and functional studies of these proteins need to take these models into consideration as they can present considerable structural differences. The characterization of the structural differences and similarities between these models is a fundamental task in structural biology requiring appropriate methods.

RESULTS

We propose a method for characterizing sets of alternative structural models. Three types of analysis are performed: grouping according to structural similarity, visualization and detection of structural variation and comparison of subsets for identifying and locating distinct conformational states. The alpha carbon atoms are used in order to analyse the backbone conformations. Alternatively, side-chain atoms are used for detailed conformational analysis of specific sites. The method takes into account estimates of atom coordinate uncertainty. The invariant regions are used to generate optimal superpositions of these models. We present the results obtained for three proteins showing different degrees of conformational variability: relative motion of two structurally conserved subdomains, a disordered subdomain and flexibility in the functional site associated with ligand binding. The method has been applied in the analysis of the alternative models available in SCOP. Considerable structural variability can be observed for most proteins.

AVAILABILITY

The results of the analysis of the SCOP alternative models, the estimates of coordinate uncertainty as well as the source code of the implementation are available in the STRuster web site: http://struster.bioinf.mpi-inf.mpg.de.

Novel class of quinone-bearing polyamines as multi-target-directed ligands to combat Alzheimer's disease

One of the characteristics of Alzheimer's disease (AD) that hinders the discovery of effective disease-modifying therapies is the multifactorial nature of its etiopathology. To circumvent this drawback, the use of multi-target-directed ligands (MTDLs) has recently been proposed as a means of simultaneously hitting several targets involved in the development of the AD syndrome. In this paper, a new class of MTDLs based on a polyamine-quinone skeleton, whose lead (memoquin, 2) showed promising properties in preclinical investigations (Cavalli et al. Angew. Chem., Int. Ed. 2007, 46, 3689-3692), is described. 3-29 were tested in vitro against a number of isolated AD-related targets, namely, AChE and BChE, and Abeta aggregation (both AChE-mediated and self-induced). Furthermore, the ability of the compounds to counteract the oxidative stress in a human neuronal-like cellular system (SH-SY5Y cells) was assayed, in both the presence and absence of NQO1, an enzyme able to generate and maintain the reduced form of quinone.

Strategies for dealing with conformational sampling in structural calculations of flexible or kinked transmembrane peptides

Peptides corresponding to transmembrane (TM) segments from membrane proteins provide a potential route for the determination of membrane protein structure. We have determined that 2 functionally critical TM segments from the mammalian Na+/H+ exchanger display well converged structure in regions separated by break points. The flexibility of these break points results in conformational sampling in solution. A brief review of available NMR structures of helical membrane proteins demonstrates that there are a number of published structures showing similar properties. Such flexibility is likely indicative of kinks in the full-length protein. This minireview focuses on methods and protocols for NMR structure calculation and analysis of peptide structures under conditions of conformational sampling. The methods outlined allow the identification and analysis of structured peptides containing break points owing to conformational sampling and the differentiation between oligomerization and ensemble-averaged observation of multiple peptide conformations.

Clustering algorithms for identifying core atom sets and for assessing the precision of protein structure ensembles

An important open question in the field of NMR-based biomolecular structure determination is how best to characterize the precision of the resulting ensemble of structures. Typically, the RMSD, as minimized in superimposing the ensemble of structures, is the preferred measure of precision. However, the presence of poorly determined atomic coordinates and multiple "RMSD-stable domains"--locally well-defined regions that are not aligned in global superimpositions--complicate RMSD calculations. In this paper, we present a method, based on a novel, structurally defined order parameter, for identifying a set of core atoms to use in determining superimpositions for RMSD calculations. In addition we present a method for deciding whether to partition that core atom set into "RMSD-stable domains" and, if so, how to determine partitioning of the core atom set. We demonstrate our algorithm and its application in calculating statistically sound RMSD values by applying it to a set of NMR-derived structural ensembles, superimposing each RMSD-stable domain (or the entire core atom set, where appropriate) found in each protein structure under consideration. A parameter calculated by our algorithm using a novel, kurtosis-based criterion, the epsilon-value, is a measure of precision of the superimposition that complements the RMSD. In addition, we compare our algorithm with previously described algorithms for determining core atom sets. The methods presented in this paper for biomolecular structure superimposition are quite general, and have application in many areas of structural bioinformatics and structural biology.

A Comparative Study on the Application of Hierarchical−Agglomerative Clustering Approaches to Organize Outputs of Reiterated Docking Runs

Reiterated runs of standard docking protocols usually provide a collection of possible binding modes rather than pinpoint a single solution. Usually, this ensemble is then ranked by means of an energy-based scoring function. However, since many degrees of approximation have to be introduced in the computation of the binding free energy, scoring functions cannot always rank the experimental pose among the top scorers. Cluster analysis might help to overcome this limit, provided that data clusterability has been earlier assessed. In this paper, first, we present a modified version of a test earlier developed by Hopkins to assess whether or not docking outputs show the natural tendency to be grouped in clusters. Then, we report the results of a comparative study on the application of different hierarchical-agglomerative cluster rules to partition docking outputs. The rule that was able to best manage the observed data was finally applied to the whole ensemble of poses collected from several docking tools. The combination of the average linkage rule with the cutting function developed by Sutcliffe and co-workers turned out to be an approach that meets all of the criteria required for a robust clustering protocol. Furthermore, a consensus clustering allowed us to identify the pose closest to the experimental one within a statistically significant cluster, whose number was always of few units.

Structural and Functional Characterization of Transmembrane Segment IV of the NHE1 Isoform of the Na+/H+ Exchanger

The Na(+)/H(+) exchanger isoform 1 is a ubiquitously expressed integral membrane protein that regulates intracellular pH in mammals. We characterized the structural and functional aspects of the critical transmembrane (TM) segment IV. Each residue was mutated to cysteine in cysteine-less NHE1. TM IV was exquisitely sensitive to mutation with 10 of 23 mutations causing greatly reduced expression and/or activity. The Phe(161) --> Cys mutant was inhibited by treatment with the water-soluble sulfhydryl-reactive compounds [2-(trimethylammonium)ethyl]methanethiosulfonate and [2-sulfonatoethyl]methanethiosulfonate, suggesting it is a pore-lining residue. The structure of purified TM IV peptide was determined using high resolution NMR in a CD(3)OH:CDCl(3):H(2)O mixture and in Me(2)SO. In CD(3)OH: CDCl(3):H(2)O, TM IV was structured but not as a canonical alpha-helix. Residues Asp(159)-Leu(162) were a series of beta-turns; residues Leu(165)-Pro(168) showed an extended structure, and residues Ile(169)-Phe(176) were helical in character. These three structured regions rotated quite freely with respect to the others. In Me(2)SO, the structure was much less defined. Our results demonstrate that TM IV is an unusually structured transmembrane segment that is exquisitely sensitive to mutagenesis and that Phe(161) is a pore-lining residue.

Dynamics of protamine 1 binding to single DNA molecules

Protamine molecules bind to and condense DNA in the sperm of most vertebrates, packaging the sperm genome in an inactive state until it can be reactivated following fertilization. By using methods that enable the analysis of protamine binding to individual DNA molecules, we have monitored the kinetics of DNA condensation and decondensation by protamine 1 (P1) and synthetic peptides corresponding to specific segments of the bull P1 DNA binding domain. Our results show that the number of clustered arginine residues present in the DNA binding domain is the most important factor affecting the condensation and stability of the DNA-protamine complex prior to the formation of inter-protamine disulfide cross-links. The high affinity of P1 for DNA is achieved by the coordinated binding of three anchoring domains, which together in bull P1 contain 19 Arg residues. The single DNA molecule experiments show that sequences containing two or more anchoring domains have an off-rate that is at least 3 orders of magnitude slower than those containing a single domain. The use of Arg, rather than Lys residues, and the inclusion of Tyr or Phe residues in the hinge regions between anchoring domains provide additional stability to the complex.

Molecular dynamics studies of the HIV-1 TAR and its complex with argininamide

The dynamic behavior of HIV-1 TAR and its complex with argininamide is investigated by means of molecular dynamics simulations starting from NMR structures, with explicit inclusion of water and periodic boundary conditions particle mesh Ewald representation of the electrostatic energy. During simulations of free and argininamide-bound TAR, local structural patterns, as determined by NMR experiments, were reproduced. An interdomain motion was observed in the simulations of free TAR, which is absent in the case of bound TAR, leading to the conclusion that the free conformation of TAR is intrinsically more flexible than the bound conformation. In particular, in the bound conformation the TAR-argininamide interface is very well ordered, as a result of the formation of a U.A.U base triple, which imposes structural constraints on the global conformation of the molecule. Free energy analysis, which includes solvation contributions, was used to evaluate the influence of van der Waals and electrostatic terms on formation of the complex and on the conformational rearrangement from free to bound TAR.

Conformational analysis of long spacers in PROSITE patterns

To determine if variable sequences (spacers) between conserved positions in a sequence motif or pattern share a consensus structure, three-dimensional structures containing PROSITE patterns with spacers of fixed length greater than three residues were analyzed. Structural similarities of a given pattern were evaluated by computing the backbone phi, psi and side-chain chi1 dihedral order parameters. The exact bias information in analyzing the conformational variability of the patterns was taken into account by introducing a new parameter, the bias coefficient, which describes the number and distribution of residue types found at each position of a pattern in the structures. The results of the analyses show that backbone conformational heterogeneity at a given position in a sequence motif does not necessarily correlate with the residue-type variability at that position, and the long spacer region can adopt a well-defined backbone conformation, in addition to the conserved residues. Furthermore, a PROSITE pattern may be redefined to yield two or more "refined" regular expressions, each corresponding to a distinct backbone conformation. A way in which the observed structural consensus in a pattern may be employed to improve the accuracy of function prediction from sequence is suggested.

Overall rotational diffusion and internal mobility in domain II of protein G from Streptococcus determined from 15N relaxation data

The backbone dynamics and overall tumbling of protein G have been investigated using 15N relaxation. Comparison of measured R2/R1 relaxation rate ratios with known three-dimensional coordinates of the protein show that the rotational diffusion tensor is significantly asymmetric, exhibiting a prolate axial symmetry. Extensive Monte Carlo simulations have been used to estimate the uncertainty due to experimental error in the relaxation rates to be D(parallel)/D(perpendicular) = 1.68 +/- 0.08, while the dispersion in the NMR ensemble leads to a variation of D(parallel)/D(perpendicular) = 1.65 +/- 0.03. Incorporation of this tensorial description into a Lipari-Szabo type analysis of internal motion has allowed us to accurately describe the local dynamics of the molecule. This analysis differs from an earlier study where the overall rotational diffusion was described by a spherical top. In this previous analysis, exchange parameters were fitted to many of the residues in the alpha helix. This was interpreted as reflecting a small motion of the alpha helix with respect to the beta sheet. We propose that the differential relaxation properties of this helix compared to the beta sheet are due to the near-orthogonality of the NH vectors in the two structural motifs with respect to the unique axis of the diffusion tensor. Our analysis shows that when anisotropic rotational diffusion is taken into account NH vectors in these structural motifs appear to be equally rigid. This study underlines the importance of a correct description of the rotational diffusion tensor if internal motion is to be accurately investigated.

Are there non-trivial dynamic cross-correlations in proteins?

The analysis of internal motion in ensembles of flexible molecules in coordinate space requires the removal of overall motion by a least-squares fitting procedure of the Cartesian coordinates. It has been demonstrated that the choice of the atom set used for fitting influences the picture of the internal motion of BPTI. We have performed essential dynamics analyses of a 1 ns molecular dynamics trajectory of the single-stranded DNA-binding protein from the Pf3 phage using either all alpha-carbon atoms or the least mobile ones for fitting the trajectory prior to the analysis. We found that covariances of atoms separated by long distances were significantly reduced in the latter case; the overall overlap of essential spaces was still high. In the second part we present a method that does not introduce and bias caused by overall motion: principal component analysis in distance space. Non-trivial dynamic cross-correlations were preserved in distance space, which answers the question posed in the title in the affirmative. However, cross-correlations were throughout smaller than those detected by standard essential dynamics analyses.

Determinants of the substrate specificity of human cytochrome P-450 CYP2D6: design and construction of a mutant with testosterone hydroxylase activity

Cytochrome P-450 CYP2D6, human debrisoquine hydroxylase, metabolizes more than 30 prescribed drugs, the vast majority of which are small molecules containing a basic nitrogen atom. In contrast, the similar mouse protein Cyp2d-9 was first characterized as a testosterone 16alpha-hydroxylase. No common substrates have been reported for the two enzymes. Here we investigate the structural basis of this difference in substrate specificity. We have earlier used a combination of NMR data and homology modelling to generate a three-dimensional model of CYP2D6 [Modi, Paine, Sutcliffe, Lian, Primrose, Wolf, C.R. and Roberts (1996) Biochemistry 35, 4541-4550]. We have now generated a homology model of Cyp2d-9 and compared the two models to identify specific amino acid residues that we believe form the substrate-binding site in each protein and therefore influence catalytic selectivity. Although there are many similarities in active site structure, the most notable difference is a phenylalanine residue (Phe-483) in CYP2D6, which in the model is located such that the bulky phenyl ring is positioned across the channel mouth, thus limiting the size of substrate that can access the active site. In Cyp2d-9, the corresponding position is occupied by an isoleucine residue, which imposes fewer steric restraints on the size of substrate that can access the active site. To investigate whether the amino acid residue at this position does indeed influence the catalytic selectivity of these enzymes, site-directed mutagenesis was used to change Phe-483 in CYP2D6 to isoleucine and also to tryptophan. CYP2D6, Cyp2d-9 and both mutant CYP2D6 proteins were co-expressed with NADPH cytochrome P-450 reductase as a functional mono-oxygenase system in Escherichia coli and their relative catalytic activities towards bufuralol and testosterone were determined. All four proteins exhibited catalytic activity towards bufuralol but only Cyp2d-9 catalysed the formation of 16alpha-hydroxytesterone. Uniquely, the CYP2D6F483I mutant acquired the ability to metabolize testosterone to a novel product, which was identified by MS and proton NMR spectroscopy as 15alpha-hydroxytestosterone. NMR spin relaxation experiments were used to measure distances between the haem iron and protons of testosterone bound to the CYP2D6F483I mutant. These experiments demonstrate that very minor modifications to the active site structure of CYP2D6 can have a profound influence on the substrate specificity of the enzyme.

OLDERADO: on-line database of ensemble representatives and domains. On Line Database of Ensemble Representatives And DOmains

In cases where the structure of a single protein is represented by an ensemble of conformations, there is often a need to determine the common features and to choose a "representative" conformation. This occurs, for example, with structures determined by NMR spectroscopy, analysis of the trajectory from a molecular dynamics simulation, or an ensemble of structures produced by comparative modeling. We reported previously automatic methods for (1) defining the atoms with low spatial variance across an ensemble (i.e., the "core" atoms) and the domains in which these atoms lie, and (2) clustering an ensemble into conformationally related subfamilies. To extend the utility of these methods, we have developed a freely available server on the World Wide Web at http:/(/)neon.chem.le.ac.uk/olderado/. This (1) contains an automatically generated database of representative structures, core atoms, and domains determined for 449 ensembles of NMR-derived protein structures in the Protein Data Bank (PDB) in May 1997, and (2) allows the user to upload a PDB-formatted file containing the coordinates of an ensemble of structures. The server returns in real time: (1) information on the residues constituting domains: (2) the structures that constitute each conformational subfamily; and (3) an interactive java-based three-dimensional viewer to visualise the domains and clusters. Such information is useful, for example, when selecting conformations to be used in comparative modeling and when choosing parts of structures to be used in molecular replacement. Here we describe the OLDERADO server.