Determination of Conformational Equilibria in Proteins Using Residual Dipolar Couplings

In order to carry out their functions, proteins often undergo significant conformational fluctuations that enable them to interact with their partners. The accurate characterization of these motions is key in order to understand the mechanisms by which macromolecular recognition events take place. Nuclear magnetic resonance spectroscopy offers a variety of powerful methods to achieve this result. We discuss a method of using residual dipolar couplings as replica-averaged restraints in molecular dynamics simulations to determine large amplitude motions of proteins, including those involved in the conformational equilibria that are established through interconversions between different states. By applying this method to ribonuclease A, we show that it enables one to characterize the ample fluctuations in interdomain orientations expected to play an important functional role

Characterization of the Interdomain Motions in Hen Lysozyme Using Residual Dipolar Couplings as Replica-Averaged Structural Restraints in Molecular Dynamics Simulations

Hen lysozyme is an enzyme characterized by the presence of two domains whose relative motions are involved in the mechanism of binding and release of the substrates. By using residual dipolar couplings as replica-averaged structural restraints in molecular dynamics simulations, we characterize the breathing motions describing the interdomain fluctuations of this protein. We found that the ensemble of conformations that we determined spans the entire range of structures of hen lysozyme deposited in the Protein Data Bank, including both the free and bound states, suggesting that the thermal motions in the free state provide access to the structures populated upon binding. The approach that we present illustrates how the use of residual dipolar couplings as replica-averaged structural restraints in molecular dynamics simulations makes it possible to explore conformational fluctuations of a relatively large amplitude in proteins

Determination of Secondary Structure Populations in Disordered States of Proteins Using Nuclear Magnetic Resonance Chemical Shifts

One of the major open challenges in structural biology is to achieve effective descriptions of disordered states of proteins. This problem is difficult because these states are conformationally highly heterogeneous and cannot be represented as single structures, and therefore it is necessary to characterize their conformational properties in terms of probability distributions. Here we show that it is possible to obtain highly quantitative information about particularly important types of probability distributions, the populations of secondary structure elements (α-helix, β-strand, random coil, and polyproline II), by using the information provided by backbone chemical shifts. The application of this approach to mammalian prions indicates that for these proteins a key role in molecular recognition is played by disordered regions characterized by highly conserved polyproline II populations. We also determine the secondary structure populations of a range of other disordered proteins that are medically relevant, including p53, α-synuclein, and the Aβ peptide, as well as an oligomeric form of αB-crystallin. Because chemical shifts are the nuclear magnetic resonance parameters that can be measured under the widest variety of conditions, our approach can be used to obtain detailed information about secondary structure populations for a vast range of different protein states

The intrinsic stability of the human prion β-sheet region investigated by molecular dynamics

<div><p>Human prion diseases are neurodegenerative disorders associated to the misfolding of the prion protein (PrP). Common features of prion disorders are the fibrillar amyloid deposits and the formation of prefibrillar oligomeric species also suggested as the origin of cytotoxicity associated with diseases. Although the process of PrP misfolding has been extensively investigated, many crucial aspects of this process remain unclear. We have here carried out a molecular dynamics study to evaluate the intrinsic dynamics of PrP β-sheet, a region that is believed to play a crucial role in prion aggregation. Moreover, as this region mediates protein association in dimeric assemblies frequently observed in prion crystallographic investigations, we also analyzed the dynamics of these intermolecular interactions. The extensive sampling of replica exchange shows that the native antiparallel β-structure of the prion is endowed with a remarkable stability. Therefore, upon unfolding, the persistence of a structured β-region may seed molecular association and influence the subsequent phases of the aggregation process. The analysis of the four-stranded β-sheet detected in the dimeric assemblies of PrP shows a tendency of this region to form dynamical structured states. The impact on the β-sheet structure and dynamics of disease associated point mutations has also been evaluated.</p> </div

Characterization of the Conformational Equilibrium between the Two Major Substates of RNase A Using NMR Chemical Shifts

Following the recognition that NMR chemical shifts can be used for protein structure determination, rapid advances have recently been made in methods for extending this strategy for proteins and protein complexes of increasing size and complexity. A remaining major challenge is to develop approaches to exploit the information contained in the chemical shifts about conformational fluctuations in native states of proteins. In this work we show that it is possible to determine an ensemble of conformations representing the free energy surface of RNase A using chemical shifts as replica-averaged restraints in molecular dynamics simulations. Analysis of this surface indicates that chemical shifts can be used to characterize the conformational equilibrium between the two major substates of this protein

Release of High-Energy Water as an Essential Driving Force for the High-Affinity Binding of Cucurbit[<i>n</i>]urils

Molecular dynamics simulations and isothermal titration calorimetry (ITC) experiments with neutral guests illustrate that the release of high-energy water from the cavity of cucurbit­[<i>n</i>]­uril (CB<i>n</i>) macrocycles is a major determinant for guest binding in aqueous solutions. The energy of the individual encapsulated water molecules decreases with increasing cavity size, because larger cavities allow for the formation of more stable H-bonded networks. Conversely, the total energy of internal water increases with the cavity size because the absolute number of water molecules increases. For CB7, which has emerged as an ultrahigh affinity binder, these counteracting effects result in a maximum energy gain through a complete removal of water molecules from the cavity. A new design criterion for aqueous synthetic receptors has therefore emerged, which is the optimization of the size of cavities and binding pockets with respect to the energy and number of residing water molecules

Molecular determinants of inactivation of the resuscitation promoting factor B from <i>Mycobacterium tuberculosis</i>

<div><p>Inactivation of revival of <i>Mycobacterium tuberculosis</i> from dormancy is one of the main goals of the WHO Global Plan to stop tuberculosis (TB) 2011–2015, given the huge reservoir of latently infected individuals. This process requires a group of secreted proteins, denoted as resuscitation-promoting factors (Rpfs). Of these, RpfB is the sole member indispensable for resuscitation <i>in vivo</i>. The first class of inhibitors of RpfB was identified among 2-nitrophenylthiocyanates. However, their inactivation mechanism is hitherto not known. To gain insight into the inactivation mechanism of one of the most promising RpfB inhibitors, 4-benzoyl-2-nitrophenyl thiocyanate, NPT7, we have performed replica exchange molecular dynamics (REMD) simulations, starting from the crystal structure of RpfB catalytic domain, derived in this study. We validated our results by resuscitation experiments of <i>M</i>. <i>tuberculosis</i> cultures. The atomic resolution crystal structure of RpfB catalytic domain identified the potential of the enzyme catalytic cleft to bind benzene rings. REMD simulations, 48 replicas, identified the key interactions for the binding of NPT7 to RpfB catalytic site. Of these, an important role is played by the thiocyanate group of NPT7. Consistently, we prove that the substitution of this group implies a complete loss of RpfB inactivation. Our results provide valuable information for modifications of NPT7 structure to enhance its binding affinity to RpfB, with the final aim of developing second-generation inhibitors of therapeutic interest in TB eradication strategy.</p> </div

Structure and dynamics of the multi-domain resuscitation promoting factor RpfB from <i>Mycobacterium tuberculosis</i>

<p>RpfB is multidomain protein that is crucial for <i>Mycobacterium tuberculosis</i> resuscitation from dormancy. This protein cleaves cell wall peptidoglycan, an essential bacterial cell wall polymer formed by glycan chains of β-(1-4)-linked-N-acetylglucosamine (GlcNAc) and N-acetylmuramic acid (MurNAc) cross-linked by short peptide stems. RpfB is structurally complex being composed of five distinct domains, namely a catalytic, a G5 and three DUF348 domains. Here, we have undertaken a combined experimental and computation structural investigations on the entire protein to gain insights into its structure–function relationships. CD spectroscopy and light scattering experiments have provided insights into the protein fold stability and into its oligomeric state. Using the available structure information, we modeled the entire protein structure, which includes the two DUF348 domains whose structure is experimentally unknown, and we analyzed the dynamic behavior of RpfB using molecular dynamics simulations. Present results highlight an intricate mutual influence of the dynamics of the different protein domains. These data provide interesting clues on the functional role of non-catalytic domains of RpfB and on the mechanism of peptidoglycan degradation necessary to resuscitation of <i>M. tuberculosis</i>.</p

Structures of the Excited States of Phospholamban and Shifts in Their Populations upon Phosphorylation

Phospholamban is an integral membrane protein that controls the calcium balance in cardiac muscle cells. As the function and regulation of this protein require the active involvement of low populated states in equilibrium with the native state, it is of great interest to acquire structural information about them. In this work, we calculate the conformations and populations of the ground state and the three main excited states of phospholamban by incorporating nuclear magnetic resonance residual dipolar couplings as replica-averaged structural restraints in molecular dynamics simulations. We then provide a description of the manner in which phosphorylation at Ser16 modulates the activity of the protein by increasing the sizes of the populations of its excited states. These results demonstrate that the approach that we describe provides a detailed characterization of the different states of phospholamban that determine the function and regulation of this membrane protein. We anticipate that the knowledge of conformational ensembles enable the design of new dominant negative mutants of phospholamban by modulating the relative populations of its conformational substates

Structural rearrangement of RpfCc monitored by NMR spectroscopy.

<p>(A, B) Overlay of the ¹H-¹⁵N HSQC spectra acquired at pH 7 (red), pH 5 (blue) and pH 3 (green). (C) Chemical shift changes (ppm) plotted versus the primary sequence. (D) Chemical shift perturbation mapping onto the RpfCc NMR structure.</p