Protein conformational changes in the bacteriorhodopsin photocycle

We report a comprehensive electron crystallographic analysis of conformational changes in the photocycle of wild-type bacteriorhodopsin and in a variety of mutant proteins with kinetic defects in the photocycle. Specific intermediates that accumulate in the late stages of the photocycle of wild-type bacteriorhodopsin, the single mutants D38R, D96N, D96G, T46V, L93A and F219L, and the triple mutant D96G/F171C/F219L were trapped by freezing two-dimensional crystals in liquid ethane at varying times after illumination with a light flash. Electron diffraction patterns recorded from these crystals were used to construct projection difference Fourier maps at 3.5 A resolution to define light-driven changes in protein conformation.Our experiments demonstrate that in wild-type bacteriorhodopsin, a large protein conformational change occurs within approximately 1 ms after illumination. Analysis of structural changes in wild-type and mutant bacteriorhodopsins under conditions when either the M or the N intermediate is preferentially accumulated reveals that there are only small differences in structure between M and N intermediates trapped in the same protein. However, a considerably larger variation is observed when the same optical intermediate is trapped in different mutants. In some of the mutants, a partial conformational change is present even prior to illumination, with additional changes occurring upon illumination. Selected mutations, such as those in the D96G/F171C/F219L triple mutant, can sufficiently destabilize the wild-type structure to generate almost the full extent of the conformational change in the dark, with minimal additional light-induced changes. We conclude that the differences in structural changes observed in mutants that display long-lived M, N or O intermediates are best described as variations of one fundamental type of conformational change, rather than representing structural changes that are unique to the optical intermediate that is accumulated. Our observations thus support a simplified view of the photocycle of wild-type bacteriorhodopsin in which the structures of the initial state and the early intermediates (K, L and M1) are well approximated by one protein conformation, while the structures of the later intermediates (M2, N and O) are well approximated by the other protein conformation. We propose that in wild-type bacteriorhodopsin and in most mutants, this conformational change between the M1 and M2 states is likely to make an important contribution towards efficiently switching proton accessibility of the Schiff base from the extracellular side to the cytoplasmic side of the membrane.     

Conformational changes induced by ionic strength and pH in two bovine myelin basic proteins.

The structures of two biologically different myelin proteins, A1 from the central nervous system and P2 from the peripheral nervous system, were investigated. Both proteins were isolated from nerve tissues. Conformational changes in the homogeneous proteins were examined in aqueous solutions by means of circular dichroism measurements. The secondary structures of both proteins proved to be very stable between pH 2.5 and pH 11.7. Unlike the P2 protein, the A1 protein is stable up to pH 13 without detectable conformational changes. The stereochemistry of the polypeptide chains of both proteins is markedly different in the presence of urea. While the value of theta222 for the A1 protein changes linearly with increasing urea concentration, a sigmoidal curve was obtained for the P2 protein. The observed differences in the dichroic properties of the basic myelin proteins A1 and P2 indicate the possibility of further structure - function correlations.     

Probing Substates in Sperm Whale Myoglobin Using High-Pressure Crystallography

Pressures in the 100 MPa range are known to have an enormous number of effects on the action of proteins, but straightforward means for determining the structural basis of these effects have been lacking. Here, crystallography has been used to probe effects of pressure on sperm whale myoglobin structure. A comparison of pressure effects with those seen at low pH suggests that structural changes under pressure are interpretable as a shift in the populations of conformational substates. Furthermore, a novel high-pressure protein crystal-cooling method has been used to show low-temperature metastability, providing an alternative to room temperature, beryllium pressure cell-based techniques. The change in protein structure due to pressure is not purely compressive and involves conformational changes important to protein activity. Correlation with low-pH structures suggests observed structural changes are associated with global conformational substates. Methods developed here open up a direct avenue for exploration of the effects of pressure on proteins.     

Structures of an unliganded neurophysin and its vasopressin complex: implications for binding and allosteric mechanisms

The structures of des 1-6 bovine neurophysin-II in the unliganded state and as its complex with lysine vasopressin were determined crystallographically at resolutions of 2.4 A and 2.3 A, respectively. The structure of the protein component of the vasopressin complex was, with some local differences, similar to that determined earlier of the full-length protein complexed with oxytocin, but relatively large differences, probably intrinsic to the hormones, were observed between the structures of bound oxytocin and bound vasopressin at Gln 4. The structure of the unliganded protein is the first structure of an unliganded neurophysin. Comparison with the liganded state indicated significant binding-induced conformational changes that were the largest in the loop region comprising residues 50-58 and in the 7-10 region. A subtle binding-induced tightening of the subunit interface of the dimer also was shown, consistent with a role for interface changes in neurophysin allosteric mechanism, but one that is probably not predominant. Interface changes are suggested to be communicated from the binding site through the strands of beta-sheet that connect these two regions, in part with mediation by Gly 23. Comparison of unliganded and liganded states additionally reveals that the binding site for the hormone alpha-amino group is largely preformed and accessible in the unliganded state, suggesting that it represents the initial site of hormone protein recognition. The potential molecular basis for its thermodynamic contribution to binding is discussed.     

Insights into nucleotide signal transduction in nitrogenase: structure of an iron protein with MgADP bound

Coupling the energy of nucleoside triphosphate binding and hydrolysis to conformational changes is a common mechanism for a number of proteins with disparate cellular functions, including those involved in DNA replication, protein synthesis, and cell differentiation. Unique to this class of proteins is the dimeric Fe protein component of nitrogenase in which the binding and hydrolysis of MgATP controls intermolecular electron transfer and reduction of nitrogen to ammonia. In the work presented here, the MgADP-bound (or "off") conformational state of the nitrogenase Fe protein has been captured and a 2.15 A resolution X-ray crystal structure is presented. The structure described herein reveals likely mechanisms for long-range communication from the nucleotide-binding sites for controlling the affinity of association with the MoFe protein component. Two pathways, termed switches I and II, appear to be integral to this nucleotide signal transduction mechanism. In addition, the structure provides the basis for the changes in the biophysical properties of the [4Fe-4S] cluster observed when Fe protein binds nucleotides. The structure of the MgADP-bound Fe protein provides important insights into the respective contributions of nucleotide interaction and complex formation in defining the conformational states that are the keys to nitrogenase catalysis.     

RosettaDock in CAPRI rounds 6-12

A challenge in protein-protein docking is to account for the conformational changes in the monomers that occur upon binding. The RosettaDock method, which incorporates sidechain flexibility but keeps the backbone fixed, was found in previous CAPRI rounds (4 and 5) to generate docking models with atomic accuracy, provided that conformational changes were mainly restricted to protein sidechains. In the recent rounds of CAPRI (6-12), large backbone conformational changes occur upon binding for several target complexes. To address these challenges, we explicitly introduced backbone flexibility in our modeling procedures by combining rigid-body docking with protein structure prediction techniques such as modeling variable loops and building homology models. Encouragingly, using this approach we were able to correctly predict a significant backbone conformational change of an interface loop for Target 20 (12 A rmsd between those in the unbound monomer and complex structures), but accounting for backbone flexibility in protein-protein docking is still very challenging because of the significantly larger conformational space, which must be surveyed. Motivated by these CAPRI challenges, we have made progress in reformulating RosettaDock using a "fold-tree" representation, which provides a general framework for treating a wide variety of flexible-backbone docking problems.     

Protein conformational dynamics of crambin in crystal, solution and in the trajectories of molecular dynamics simulations

Atomic displacement parameters — B factors of the eight crambin crystal structures obtained at 0.54–1.5 Å resolution and temperatures of 100–293 K have been analyzed. The comparable contributions to the B factor values are the intramolecular motions which are modeled by the harmonic vibration calculations and derived from the molecular dynamics simulation (MD) as well as rigid body changes in the position of a protein molecule as a whole. In solution for the average NMR structure of crambin the amplitudes of the backbone atomic fluctuations of the most residues of the segments with the regular backbone conformations are close to the amplitudes of the small scale harmonic vibrations. For the same residues the probability of the medium scale fluctuations fixed by the hydrogen exchange method is very low. The restricted conformational mobility of those segments is coupled with the depressed amplitudes of the fluctuation changes of the tertiary structure registered by the residue accessibility changes in an ensemble of NMR structures that forms the average NMR structure of crambin. The amplitudes of temperature fluctuations of backbone atoms and the tertiary structure raise in the segment with the irregular conformations, turn and loops. In the same segments the amplitudes of the calculated harmonic vibrations also increase, but to a lesser extent and especially in the interhelical loop with the most strong and complicated fluctuation changes of the backbone conformation. In solution for the NMR structure in this loop the conformational transitions occur between the conformational substates separated by the energy barriers, but they are not observed even in the long 100 ns trajectories from the MD simulation of crambin. These strong local fluctuation changes of the structure may play a key role in the protein functioning and modern performance improvements in the MD simulation techniques are oriented to increase the probability of protein appearance in the trajectories from the MD simulations.     

Conformational states of the glucocorticoid receptor DNA-binding domain from molecular dynamics simulations

Molecular dynamics simulations (MD) have been performed on variant crystal and NMR-derived structures of the glucocorticoid receptor DNA-binding domain (GR DBD). A loop region five residues long, the so-called D-box, exhibits significant flexibility, and transient perturbations of the tetrahedral geometry of two structurally important Cys4 zinc finger are seen, coupled to conformational changes in the D-box. In some cases, one of the Cys ligands to zinc exchanges with water, although no global distortion of the protein structure is observed. Thus, from MD simulation, dynamics of the D-box could partly be explained by solvent effects in conjunction with structural reformation of the zinc finger.     

Solution structure of porcine pancreatic phospholipase A2.

The lipolytic enzyme phospholipase A2 (PLA2) is involved in the degradation of high-molecular weight phospholipid aggregates in vivo. The enzyme has very high catalytic activities on aggregated substrates compared with monomeric substrates, a phenomenon called interfacial activation. Crystal structures of PLA2s in the absence and presence of inhibitors are identical, from which it has been concluded that enzymatic conformational changes do not play a role in the mechanism of interfacial activation. The high-resolution NMR structure of porcine pancreatic PLA2 free in solution was determined with heteronuclear multidimensional NMR methodology using doubly labeled 13C, 15N-labeled protein. The solution structure of PLA2 shows important deviations from the crystal structure. In the NMR structure the Ala1 alpha-amino group is disordered and the hydrogen bonding network involving the N-terminus and the active site is incomplete. The disorder observed for the N-terminal region of PLA2 in the solution structure could be related to the low activity of the enzyme towards monomeric substrates. The NMR structure of PLA2 suggests, in contrast to the crystallographic work, that conformational changes do play a role in the interfacial activation of this enzyme.     

Conformational states of cytochrome P450cam revealed by trapping of synthetic molecular wires

Members of the ubiquitous cytochrome P450 family catalyze a vast range of biologically significant reactions in mammals, plants, fungi, and bacteria. Some P450s display a remarkable promiscuity in substrate recognition, while others are very specific with respect to substrate binding or regio and stereo-selective catalysis. Recent results have suggested that conformational flexibility in the substrate access channel of many P450s may play an important role in controlling these effects. Here, we report the X-ray crystal structures at 1.8A and 1.5A of cytochrome P450cam complexed with two synthetic molecular wires, D-4-Ad and D-8-Ad, consisting of a dansyl fluorophore linked to an adamantyl substrate analog via an alpha,omega-diaminoalkane chain of varying length. Both wires bind with the adamantyl moiety in similar positions at the camphor-binding site. However, each wire induces a distinct conformational response in the protein that differs from the camphor-bound structure. The changes involve significant movements of the F, G, and I helices, allowing the substrate access channel to adapt to the variable length of the probe. Wire-induced opening of the substrate channel also alters the I helix bulge and Thr252 at the active site with binding of water that has been proposed to assist in peroxy bond cleavage. The structures suggest that the coupling of substrate-induced conformational changes to active-site residues may be different in P450cam and recently described mammalian P450 structures. The wire-induced changes may be representative of the conformational intermediates that must exist transiently during substrate entry and product egress, providing a view of how substrates enter the deeply buried active site. They also support observed examples of conformational plasticity that are believed be responsible for the promiscuity of drug metabolizing P450s. Observation of such large changes in P450cam suggests that substrate channel plasticity is a general property inherent to all P450 structures.     

Close correspondence between the motions from principal component analysis of multiple HIV-1 protease structures and elastic network modes

The large number of available HIV-1 protease structures provides a remarkable sampling of conformations of the different conformational states, which can be viewed as direct structural information about the dynamics of the HIV-1 protease. After structure matching, we apply principal component analysis (PCA) to obtain the important apparent motions for both bound and unbound structures. There are significant similarities between the first few key motions and the first few low-frequency normal modes calculated from a static representative structure with an elastic network model (ENM), strongly suggesting that the variations among the observed structures and the corresponding conformational changes are facilitated by the low-frequency, global motions intrinsic to the structure. Similarities are also found when the approach is applied to an NMR ensemble, as well as to molecular dynamics (MD) trajectories. Thus, a sufficiently large number of experimental structures can directly provide important information about protein dynamics, but ENM can also provide similar sampling of conformations.     

Calculating three-dimensional changes in protein structure due to amino-acid substitutions: the variable region of immunoglobulins

A procedure (coupled perturbation procedure, CPP) is introduced as a specific method for calculating the detailed three-dimensional structure of a protein molecule which has a number of amino-acid substitutions relative to some previously determined "parent" protein structure. The accuracy of the procedure is tested by calculating the conformation of a region of the human immunoglobulin fragment Fab Kol based on the analogous region of the human immunoglobulin fragment Fab New. Both structures have previously been determined crystallographically. The calculated model is accurate to the extent that both of the sequence differences in the region are modeled correctly and that conformational changes in a number of nearby residues are correctly identified. CPP is shown to give better results than other commonly used modeling procedures when applied to the same problem.     

A Correspondence Between Solution-State Dynamics of an Individual Protein and the Sequence and Conformational Diversity of its Family

Conformational ensembles are increasingly recognized as a useful representation to describe fundamental relationships between protein structure, dynamics and function. Here we present an ensemble of ubiquitin in solution that is created by sampling conformational space without experimental information using “Backrub” motions inspired by alternative conformations observed in sub-Angstrom resolution crystal structures. Backrub-generated structures are then selected to produce an ensemble that optimizes agreement with nuclear magnetic resonance (NMR) Residual Dipolar Couplings (RDCs). Using this ensemble, we probe two proposed relationships between properties of protein ensembles: (i) a link between native-state dynamics and the conformational heterogeneity observed in crystal structures, and (ii) a relation between dynamics of an individual protein and the conformational variability explored by its natural family. We show that the Backrub motional mechanism can simultaneously explore protein native-state dynamics measured by RDCs, encompass the conformational variability present in ubiquitin complex structures and facilitate sampling of conformational and sequence variability matching those occurring in the ubiquitin protein family. Our results thus support an overall relation between protein dynamics and conformational changes enabling sequence changes in evolution. More practically, the presented method can be applied to improve protein design predictions by accounting for intrinsic native-state dynamics.     

Molecular crowding inhibits intramolecular breathing motions in proteins

In aqueous solution some proteins undergo large-scale movements of secondary structures, subunits or domains, referred to as protein “breathing”, that define a native-state ensemble of structures. These fluctuations are sensitive to the nature and concentration of solutes and other proteins and are thereby expected to be different in the crowded interior of a cell than in dilute solution. Here we use a combination of wide angle X-ray scattering (WAXS) and computational modeling to derive a quantitative measure of the spatial scale of conformational fluctuations in a protein solution. Concentration-dependent changes in the observed scattering intensities are consistent with a model of structural fluctuations in which secondary structures undergo rigid-body motions relative to one another. This motion increases with decreasing protein concentration or increasing temperature. Analysis of a set of five structurally and functionally diverse proteins reveals a diversity of kinetic behaviors. Proteins with multiple disulfide bonds exhibit little or no increase in breathing in dilute solutions. The spatial extent of structural fluctuations appears highly dependent on both protein structure and concentration and is universally suppressed at very high protein concentrations.     

Solution NMR studies of colicin E1 C-terminal thermolytic peptide. Structural comparison with colicin A and the effects of pH changes

The aqueous solution structure of the C-terminal thermolytic peptide of colicin E1 has been investigated using both one- and two-dimensional NMR techniques. The NMR data are consistent with a fold for the peptide very similar to that reported for the colicin A C-terminal peptide in the crystalline state, although some differences have been noted. The one-dimensional NMR spectrum of the peptide has been used to follow changes in both the structure and dynamics of the peptide on changing pH. The in vitro functionally competent form of the peptide (present in solution only below pH 6) does not differ in structure significantly from the higher pH form. However, small local conformational changes are observed together with an increase in mobility in some of the more hydrophilic regions. This suggests that the effect of lower pH is to change the ease with which the major conformational changes during insertion into a membrane can occur.     

Ligand-induced conformational changes: improved predictions of ligand binding conformations and affinities

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.     

Modeling protein conformational changes by iterative fitting of distance constraints using reoriented normal modes

Recently we have developed a normal-modes-based algorithm that predicts the direction of protein conformational changes given the initial state crystal structure together with a small number of pairwise distance constraints for the end state. Here we significantly extend this method to accurately model both the direction and amplitude of protein conformational changes. The new protocol implements a multisteps search in the conformational space that is driven by iteratively minimizing the error of fitting the given distance constraints and simultaneously enforcing the restraint of low elastic energy. At each step, an incremental structural displacement is computed as a linear combination of the lowest 10 normal modes derived from an elastic network model, whose eigenvectors are reorientated to correct for the distortions caused by the structural displacements in the previous steps. We test this method on a list of 16 pairs of protein structures for which relatively large conformational changes are observed (root mean square deviation >3 angstroms), using up to 10 pairwise distance constraints selected by a fluctuation analysis of the initial state structures. This method has achieved a near-optimal performance in almost all cases, and in many cases the final structural models lie within root mean square deviation of 1 approximately 2 angstroms from the native end state structures.     

Structural insights into the molecular mechanism of H‐NOX activation

Nitric oxide (NO) signaling in mammals controls important processes such as smooth muscle relaxation and neurotransmission by the activation of soluble guanylate cyclase (sGC). NO binding to the heme domain of sGC leads to dissociation of the iron–histidine (Fe–His) bond, which is required for enzyme activity. The heme domain of sGC belongs to a larger class of proteins called H‐NOX (Heme‐Nitric oxide/OXygen) binding domains. Previous crystallographic studies on H‐NOX domains demonstrate a correlation between heme bending and protein conformation. It was unclear, however, whether these structural changes were important for signal transduction. Subsequent NMR solution structures of H‐NOX proteins show a conformational change upon disconnection of the heme and proximal helix, similar to those observed in the crystallographic studies. The atomic details of these conformational changes, however, are lacking in the NMR structures especially at the heme pocket. Here, a high‐resolution crystal structure of an H‐NOX mutant mimicking a broken Fe–His bond is reported. This mutant exhibits specific changes in heme conformation and major N‐terminal displacements relative to the wild‐type H‐NOX protein. Fe–His ligation is ubiquitous in all H‐NOX domains, and therefore, the heme and protein conformational changes observed in this study are likely to occur throughout the H‐NOX family when NO binding leads to rupture of the Fe–His bond.     

Progress in protein-protein docking: atomic resolution predictions in the CAPRI experiment using RosettaDock with an improved treatment of side-chain flexibility

RosettaDock uses real-space Monte Carlo minimization (MCM) on both rigid-body and side-chain degrees of freedom to identify the lowest free energy docked arrangement of 2 protein structures. An improved version of the method that uses gradient-based minimization for off-rotamer side-chain optimization and includes information from unbound structures was used to create predictions for Rounds 4 and 5 of CAPRI. First, large numbers of independent MCM trajectories were carried out and the lowest free energy docked configurations identified. Second, new trajectories were started from these lowest energy structures to thoroughly sample the surrounding conformation space, and the lowest energy configurations were submitted as predictions. For all cases in which there were no significant backbone conformational changes, a small number of very low-energy configurations were identified in the first, global search and subsequently found to be close to the center of the basin of attraction in the free energy landscape in the second, local search. Following the release of the experimental coordinates, it was found that the centers of these free energy minima were remarkably close to the native structures in not only the rigid-body orientation but also the detailed conformations of the side-chains. Out of 8 targets, the lowest energy models had interface root-mean-square deviations (RMSDs) less than 1.1 A from the correct structures for 6 targets, and interface RMSDs less than 0.4 A for 3 targets. The predictions were top submissions to CAPRI for Targets 11, 12, 14, 15, and 19. The close correspondence of the lowest free energy structures found in our searches to the experimental structures suggests that our free energy function is a reasonable representation of the physical chemistry, and that the real space search with full side-chain flexibility to some extent solves the protein-protein docking problem in the absence of significant backbone conformational changes. On the other hand, the approach fails when there are significant backbone conformational changes as the steric complementarity of the 2 proteins cannot be modeled without incorporating backbone flexibility, and this is the major goal of our current work.