Search

Protein filaments are dynamic entities that respond to external stimuli by slightly or substantially modifying the internal binding geometries between successive protomers. This results in overall changes in the filament architecture, which are difficult to model due to the helical character of the system. Here, we describe how distortions in RecA nucleofilaments and their consequences on the filament-DNA and bound DNA-DNA interactions at different stages of the homologous recombination process can be modeled using the PTools/Heligeom software and subsequent molecular dynamics simulation with NAMD. Modeling methods dealing with helical macromolecular objects typically rely on symmetric assemblies and take advantage of known symmetry descriptors. Other methods dealing with single objects, such as MMTK or VMD, do not integrate the specificities of regular assemblies. By basing the model building on binding geometries at the protomer-protomer level, PTools/Heligeom frees the building process from a priori knowledge of the system topology and enables irregular architectures and symmetry disruption to be accounted for. Graphical abstract: [Image: see text] Model of ATP hydrolysis-induced distortions in the recombinant nucleoprotein, obtained by combining RecA-DNA and two RecA-RecA binding geometries.

Reproduction

To understand molecular recognition, one needs to characterize the conformational changes that take place when a free (unbound) molecule forms a complex with another protein. An analysis of 281 protein components is presented, delineating such changes in terms of accessible surface area, secondary structure, crystallographic B factors and disorder-to-order transitions., Protein interactions are essential in all biological processes. The changes brought about in the structure when a free component forms a complex with another molecule need to be characterized for a proper understanding of molecular recognition as well as for the successful implementation of docking algorithms. Here, unbound (U) and bound (B) forms of protein structures from the Protein–Protein Interaction Affinity Database are compared in order to enumerate the changes that occur at the interface atoms/residues in terms of the solvent-accessible surface area (ASA), secondary structure, temperature factors (B factors) and disorder-to-order transitions. It is found that the interface atoms optimize contacts with the atoms in the partner protein, which leads to an increase in their ASA in the bound interface in the majority (69%) of the proteins when compared with the unbound interface, and this is independent of the root-mean-square deviation between the U and B forms. Changes in secondary structure during the transition indicate a likely extension of helices and strands at the expense of turns and coils. A reduction in flexibility during complex formation is reflected in the decrease in B factors of the interface residues on going from the U form to the B form. There is, however, no distinction in flexibility between the interface and the surface in the monomeric structure, thereby highlighting the potential problem of using B factors for the prediction of binding sites in the unbound form for docking another protein. 16% of the proteins have missing (disordered) residues in the U form which are observed (ordered) in the B form, mostly with an irregular conformation; the data set also shows differences in the composition of interface and non-interface residues in the disordered polypeptide segments as well as differences in their surface burial.

The buried surface area (BSA), which measures the size of the interface in a protein-protein complex may differ from the accessible surface area (ASA) lost upon association (which we call DSA), if conformation changes take place. To evaluate the DSA, we measure the ASA of the interface atoms in the bound and unbound states of the components of 144 protein-protein complexes taken from the Protein-Protein Interaction Affinity Database of Kastritis et al. (2011). We observe differences exceeding 20%, and a systematic bias in the distribution. On average, the ASA calculated in the bound state of the components is 3.3% greater than in their unbound state, and the BSA, 7% greater than the DSA. The bias is observed even in complexes where the conformation changes are small. An examination of the bound and unbound structures points to a possible origin: local movements optimize contacts with the other component at the cost of internal contacts, and presumably also the binding free energy.

We compare the changes in side chain conformations that accompany the formation of protein–protein complexes, in residues forming either the interface or the remainder of the solvent-accessible surface of the proteins in the Docking Benchmark 3.0. We find that the interface residues undergo significantly more changes than other surface residues, and these changes are more likely to convert them from a high-energy torsion angle state to a lower-energy one than the reverse. Moreover, in both the unbound proteins and the complexes, the interface residues are more frequently found to be in a high-energy torsion angle state than the noninterface residues. As these differences exist before the binding step, they may be relevant to specificity and help in identifying binding sites for docking predictions. Proteins 2010. © 2010 Wiley-Liss, Inc.

International audience; Describing biological macromolecular energetics from computer simulations can pose major challenges, and often necessitates enhanced conformational sampling. We describe the calculation of conformational free-energy profiles along carefully chosen collective coordinates: “consensus” normal modes, developed recently as robust alternatives to conventional normal modes. In an application to the HIV-1 protease, we obtain efficient sampling of significant flap opening movements governing inhibitor binding from relatively short simulations, in close correspondence with experimental results.

International audience; While proteins and nucleic acids are the fundamental components of an organism, Biology itself is based on the interactions they make with each other. Analyzing macromolecular interactions typically requires handling systems involving from two to hundreds of polypeptide chains. After a brief overview of the modeling challenges faced in computational structural biology, this chapter presents concepts and tools aiming at improving our understanding of the link between the static structures of macromolecular complexes and their biophysical/biological properties. We discuss geometrical approaches suited to atomic-resolution complexes and to large protein assemblies; for each, we also present examples of their successful application in quantifying and interpreting biological data. This methodology includes state-of-the-art geometric analyses of surface area, volume, curvature, and topological properties (isolated components, cavities, voids, cycles) related to Voronoï constructions in the context of structure analysis. On the applied side, we present novel insights into real biological problems gained thanks to these modeling tools.

Unfixed chicken erythrocyte chromatin fibers in very low salt have been imaged with a scanning force microscope operating in the tapping mode in air at ambient humidity. These images reveal a three-dimensional organization of the fibers. The planar "zig-zag" conformation is rare, and extended "beads-on-a-string" fibers are seen only in chromatin depleted of histones H1 and H5. Glutaraldehyde fixation reveals very similar structures. Fibers fixed in 10 mM salt appear somewhat more compacted. These results, when compared with modeling studies, suggest that chromatin fibers may exist as irregular three-dimensional arrays of nucleosomes even at low ionic strength.

International audience; The CAPRI (Critical Assessment of Predicted Interactions) and CASP (Critical Assessment of protein Structure Prediction) experiments have demonstrated the power of community-wide tests of methodology in assessing the current state of the art and spurring progress in the very challenging areas of protein docking and structure prediction. We sought to bring the power of community-wide experiments to bear on a very challenging protein design problem that provides a complementary but equally fundamental test of current understanding of protein-binding thermodynamics. We have generated a number of designed protein-protein interfaces with very favorable computed binding energies but which do not appear to be formed in experiments, suggesting that there may be important physical chemistry missing in the energy calculations. A total of 28 research groups took up the challenge of determining what is missing: we provided structures of 87 designed complexes and 120 naturally occurring complexes and asked participants to identify energetic contributions and/or structural features that distinguish between the two sets. The community found that electrostatics and solvation terms partially distinguish the designs from the natural complexes, largely due to the nonpolar character of the designed interactions. Beyond this polarity difference, the community found that the designed binding surfaces were, on average, structurally less embedded in the designed monomers, suggesting that backbone conformational rigidity at the designed surface is important for realization of the designed function. These results can be used to improve computational design strategies, but there is still much to be learned; for example, one designed complex, which does form in experiments, was classified by all metrics as a nonbinder.

Background Cathepsin B (catB) is a promising target for anti-cancer drug design due to its implication in several steps of tumorigenesis. catB activity and inhibition are pH-dependent, making it difficult to identify efficient inhibitor candidates for clinical trials. In addition it is known that heparin binding stabilizes the enzyme in alkaline conditions. However, the molecular mechanism of stabilization is not well understood, indicating the need for more detailed structural and dynamic studies in order to clarify the influence of pH and heparin binding on catB stability. Results Our pKa calculations of catB titratable residues revealed distinct protonation states under different pH conditions for six key residues, of which four lie in the crucial interdomain interface. This implies changes in the overall charge distribution at the catB surface, as revealed by calculation of the electrostatic potential. We identified two basic surface regions as possible heparin binding sites, which were confirmed by docking calculations. Molecular dynamics (MD) of both apo catB and catB-heparin complexes were performed using protonation states for catB residues corresponding to the relevant acidic or alkaline conditions. The MD of apo catB at pH 5.5 was very stable, and presented the highest number and occupancy of hydrogen bonds within the inter-domain interface. In contrast, under alkaline conditions the enzyme's overall flexibility was increased: interactions between active site residues were lost, helical content decreased, and domain separation was observed as well as high-amplitude motions of the occluding loop – a main target of drug design studies. Essential dynamics analysis revealed that heparin binding modulates large amplitude motions promoting rearrangement of contacts between catB domains, thus favoring the maintenance of helical content as well as active site stability. Conclusions The results of our study contribute to unraveling the molecular events involved in catB inactivation in alkaline pH, highlighting the fact that protonation changes of few residues can alter the overall dynamics of an enzyme. Moreover, we propose an allosteric role for heparin in the regulation of catB stability in such a manner that the restriction of enzyme flexibility would allow the establishment of stronger contacts and thus the maintenance of overall structure.

International audience; Small G-proteins of the superfamily Ras function as molecular switches, interacting with different cellular partners according to their activation state. G-protein activation involves the dissociation of bound GDP and its replacement by GTP, in an exchange reaction that is accelerated and regulated in the cell by guanine-nucleotide exchange factors (GEFs). Large conformational changes accompany the exchange reaction, and our understanding of the mechanism is correspondingly incomplete. However, much knowledge has been derived from structural studies of blocked or inactive mutant GEFs, which presumably closely represent intermediates in the exchange reaction and yet which are by design incompetent for carrying out the nucleotide exchange reaction. In this study we have used comparative modelling to recreate an exchange-competent form of a late, pre-GDP-ejection intermediate species in Arf1, a well-characterized small G-protein. We extensively characterized three distinct models of this intermediate using molecular dynamics simulations, allowing us to address ambiguities related to the mutant structural studies. We observed in particular the unfavorable nature of Mg 2z associated forms of the complex and the establishment of closer Arf1-GEF contacts in its absence. The results of this study shed light on GEF-mediated activation of this small G protein and on predicting the fate of the Mg ion at a critical point in the exchange reaction. The structural models themselves furnish additional targets for interfacial inhibitor design, a promising direction for exploring potentially druggable targets with high biological specificity.

International audience; Protein flexibility is essential for enzymatic function, ligand binding, and protein–protein or protein–nucleic acid interactions. Normal mode analysis has increasingly been shown to be well suited for studying such flexibility, as it can be used to identify favorable structural deformations that correspond to functional motions. However, normal modes are strictly relevant to a single structure, reflecting a particular minimum on a complex energy surface, and are thus susceptible to artifacts. We describe a new theoretical framework for determining ‘‘consensus’’ normal modes from a set of related structures, such as those issuing from a short molecular dynamics simulation. This approach is more robust than standard normal mode analysis, and provides higher collectivity and symmetry properties. In an application to HIV-1 protease, the low-frequency consensus modes describe biologically relevant motions including flap openingand closing that can be used in interpreting structural changes accompanying the bindingof widely differing inhibitors.

International audience; Recent studies of the mechanisms involved in the regulation of gene expression in eukaryotic organisms depict a highly complex process requiring a coordinated rearrangement of numerous molecules to mediate DNA accessibility. Silencing in Saccharomyces cerevisiae involves the Sir family of proteins. Sir3p, originally described as repressing key areas of the yeast genome through interactions with the tails of histones H3 and H4, appears to have additional roles in that process, including involvement with a DNA binding component. Our in vitro studies focused on the characterization of Sir3p-nucleic acid interactions and their biological functions in Sir3p-mediated silencing using binding assays, EM imaging, and theoretical modeling. Our results suggest that the initial Sir3p recruitment is partially DNA-driven, highly cooperative, and dependent on nucleosomal features other than histone tails. The initial step appears to be rapidly followed by the spreading of silencing using linker DNA as a track.

International audience; In mammals, the olfactory epithelium secretes odorant-binding proteins (OBPs), which are lipocalins found freely dissolved in the mucus layer protecting the olfactory neurons. OBPs may act as passive transporters of predominantly hydrophobic odorant molecules across the aqueous mucus layer, or they may play a more active role in which the olfactory neuronal receptor recognizes the OBP-ligand complex. To better understand the molecular events accompanying the initial steps in the olfaction process, we have performed molecular dynamics studies of rat and pig OBPs with the odorant molecule thymol. These calculations provide an atomic level description of conformational changes and pathway intermediates that remain difficult to study directly. A series of eight independent molecular dynamics trajectories of rat OBP permitted the observation of a consensus pathway for ligand unbinding and the calculation of the potential of mean force (PMF) along this path. Titration microcalorimetry confirmed the specific binding of thymol to this protein with a strong hydrophobic component. In both rat and pig OBPs we observed lipocalin strand pair opening in the presence of ligand, consistent with potential roles of these proteins in olfactive receptor recognition.

The proteins Arno and Gea2 of the Sec7 family can promote GDP–GTP exchange on Arf1, a small GTP-binding protein, which coordinates coated vesicle formation for protein transport within the cell. Crystal structures of the essential Sec7 domain (Sec7d) of Gea2 in the free and Arf1-bound forms suggest that conformational dynamics of the Sec7d as well as those of the G-protein play a role in nucleotide exchange. Starting from a set of complementary crystal structures, we compared the collective movements of unbound Gea2 and Arno Sec7 domains, Arf1-GDP, and the Arf1–Gea2(Sec7d) nucleotide-free complex using normal modes analyses. In all unbound Sec7d analyses, significant low-energy movements were found to lead to closure of the hydrophobic groove towards the form seen in the Arf1–Gea2(Sec7d) complex, suggesting that groove closure is a general feature of the Sec7 family. Low-energy movements in Arf1-GDP implicate critical switch 1 and 2 residues which are coupled to modifications in the myristoylated N-terminal-helix binding site at the other end of the “interswitch” beta hairpin. It is suggested that Sec7d groove closure upon docking of the two molecules may permit extraction of switch 1 from Arf1-GDP and prepare the complex for movement of the interswitch, which is central to the membrane-linked exchange activity. Large-scale collective movements in the Arf1–Sec7d complex appear to participate in the insertion of the Sec7d Glu finger into the GDP binding site to promote actual nucleotide release.

We discuss the efflux of entrapped marker material from liposomes or cells through pores in the membrane, being monitored by the time course of a certain signal F (e.g., fluorescence emission). This is expressed in terms of an appropriate normalized function of time, the so-called efflux function E(t). Under conditions frequently encountered in practice the measured E(t) can be easily related to the forward rate of pore formation if the liposomes/cells are monodisperse in size. In the basic case of a time-independent rate law it turns out that E(t) must be single exponential. Deviations from such a simple functional behavior might be due to a fairly broad distribution of liposome/cell sizes and/or a more complicated pore formation mechanism. A relevant evaluation of original data is demonstrated making use of experimental results obtained with small unilamellar lipid vesicles where pores are induced by the antibiotic peptide alamethicin. This includes the application of a general method to eliminate the effect of a given liposome/cell size distribution.

We demonstrate the utility of normal mode analysis in correctly predicting the binding modes of inhibitors in the active sites of matrix metalloproteinases (MMPs). We show the accuracy in predicting the positions of MMP-3 inhibitors is strongly dependent on which structure is used as the target, especially when it has been energy minimized. This dependency can be overcome by using intermediate structures generated along one of the normal modes previously calculated for a given target. These results may be of prime importance for further in silico drug discovery.