Deswapping bovine odorant binding protein

The X-ray structure of bovine Odorant Binding Protein (bOBP) revealed its association as a domain swapped dimer. bOBP, devoid of any cysteines, contrasts with other mammalian OBPs, which are monomeric and possess at least one disulfide bridge. We have produced a mutant of bOBP in which a glycine residue was inserted after position 121. This mutation yielded a monomeric bOBP-121Gly+ in which domain swapping has been reverted. Here, we have subsequently introduced two mutations, Trp64Cys and His155Cys, in view to stabilize the putative monomer with a disulfide bridge. We have determined the crystal structure of this triple mutant at 1.65 A resolution. The mutant protein is monomeric, stabilized by a disulfide bridge between Trp64Cys and His155Cys, with a backbone superimposable to that of native bOBP, with the exception of the hinge and of the 10 residues at the C-terminus. bOBP triple mutant binds 1-amino-anthracene, 1-octen-3-ol (bOBP co-purified ligand) and other ligands with microM Kd values comparable to those of the swapped dimer.     

Dissociation and unfolding of bovine odorant binding protein at acidic pH

Dissociation of bovine odorant binding protein (bOBP) dimers to monomers at pH 2.5 has been confirmed through size exclusion chromatography experiments. Moreover, structural and binding properties of the acidic monomer and neutral dimer have been compared using a combination of experimental (circular dichroism and fluorescence) and computational (molecular dynamics) techniques. The secondary and tertiary structures of bOBP are largely maintained at acidic pH, but molecular dynamics simulations suggest the loop regions (N-terminal residues, Omega-loop and C-terminal segments) are more relaxed and Phe36 and Tyr83 residues are involved in the regulation of the binding cavity entrance. The formation of a molten globule state at acidic pH, suggested by the strong enhancement of fluorescence of 8-anilino-1-naphtalenesulphonic acid (ANS), is not confirmed by any significant change in the near UV circular dichroism spectrum. Functionality measurements, deduced from the interaction of bOBP with 1-amino-anthracene (AMA), show that the binding capacity of the protein at acidic pH is preserved, though slightly looser than at neutral pH. Unfolding of acidic bOBP, induced by guanidinium chloride (GdnHCl), was investigated by means of CD spectroscopy, steady state fluorescence, fluorescence anisotropy and light scattering. The stability of the acidic monomer is lower than that of the neutral dimer, owing to the loss of the swapping interactions, but renaturation is completely reversible. Finally, in contrast with the neutral dimer, at low denaturant concentration some aggregation of the acidic monomer, which vanishes before the unfolding transition, has been observed.     

Constrained modeling of spin-labeled major coat protein mutants from M13 bacteriophage in a phospholipid bilayer

The family of three-dimensional molecular structures of the major coat protein from the M13 bacteriophage, which was determined in detergent micelles by NMR methods, has been analyzed by constrained geometry optimization in a phospholipid environment. A single-layer solvation shell of dioleoyl phosphatidylcholine lipids was built around the protein, after replacing single residues by cysteines with a covalently attached maleimide spin label. Both the residues substituted and the phospholipid were chosen for comparison with site-directed spin labeling EPR measurements of distance and local mobility made previously on membranous assemblies of the M13 coat protein purified from viable mutants. The main criteria for identifying promising candidate structures, out of the 300 single-residue mutant models generated for the membranous state, were 1) lack of steric conflicts with the phospholipid bilayer, 2) good match of the positions of spin-labeled residues along the membrane normal with EPR measurements, and 3) a good match between the sequence profiles of local rotational freedom and a structural restriction parameter for the spin-labeled residues obtained from the model. A single subclass of structure has been identified that best satisfies these criteria simultaneously. The model presented here is useful for the interpretation of future experimental data on membranous M13 coat protein systems. It is also a good starting point for full-scale molecular dynamics simulations and for the design of further site-specific spectroscopic experiments.     

Investigation of the dynamics of the viral immediate-early protein 1 in different conformations and oligomerization states

The viral immediate-early protein 1 (IE1) is crucial for efficient replication of cytomegalovirus (CMV). A recent crystal structure of the IE1 protein from rhesus CMV revealed that the protein exhibits a novel fold and crystallizes in two slightly different dimeric arrangements. Molecular dynamics simulations and energetic analyses performed in this study show that both dimers are stable and allowed us to identify a common set of five residues that appear particularly important for dimer formation. These residues are distributed over the entire dimer interface and do not form a typical hot spot for protein interactions. In addition, the dimer interface of IE1 proved to include a high portion of hydrophilic interactions pointing toward the transient nature of dimer formation. Characterization of monomeric and dimeric IE1 revealed three sequentially discontinuous dynamic domains that exhibit correlated motion within the domain and are simultaneously anti-correlated to the adjacent domains. The hinge motions observed between the dynamic domains increase the shape complementarity to the coiled–coil region of tripartite motif proteins, suggesting that the detected dynamics of IE1 might be physiologically important by enabling a better interaction with its cellular target molecules.     

Pump-probe molecular dynamics as a tool for studying protein motion and long range coupling

A new method for analyzing the dynamics of proteins is developed and tested. The method, pump-probe molecular dynamics, excites selected atoms or residues with a set of oscillating forces, and the transmission of the impulse to other parts of the protein is probed using Fourier transform of the atomic motions. From this analysis, a coupling profile can be determined which quantifies the degree of interaction between pump and probe residues. Various physical properties of the method such as reciprocity and speed of transmission are examined to establish the soundness of the method. The coupling strength can be used to address questions such as the degree of interaction between different residues at the level of dynamics, and identify propagation of influence of one part of the protein on another via "pathways" through the protein. The method is illustrated by analysis of coupling between different secondary structure elements in the allosteric protein calmodulin, and by analysis of pathways of residue-residue interaction in the PDZ domain protein previously elucidated by genomics and mutational studies.     

Molecular dynamics as a tool to detect protein foldability. A mutant of domain B1 of protein G with non-native secondary structure propensities

The usefulness of molecular dynamics to assess the structural integrity of mutants containing several mutations has been investigated. Our goal was to determine whether molecular dynamics would be able to discriminate mutants of a protein having a close-to-wild-type fold, from those that are not folded under the same conditions. We used as a model the B1 domain of protein G in which we replaced the unique central alpha-helix by the sequence of the second beta-hairpin, which has a strong intrinsic propensity to form this secondary structure in solution. In the resulting protein, one-third of the secondary structure has been replaced by a non-native one. Models of the mutants were built based on the three-dimensional structure of the wild-type GB1 domain. During 2 ns of molecular dynamics simulations on these models, mutants containing up to 10 mutations in the helix retained the native fold, while another mutant with an additional mutation unfolded. This result is in agreement with our circular dichroism and NMR experiments, which indicated that the former mutants fold into a structure similar to the wild-type, as opposed to the latter mutant which is partly unfolded. Additionally, a mutant containing six mutations scattered through the surface of the domain, and which is unfolded, was also detected by the simulation. This study suggests that molecular dynamics calculations could be performed on molecular models of mutants of a protein to evaluate their foldability, prior to a mutagenesis experiment.     

Intramolecular dynamics of low molecular weight protein tyrosine phosphatase in monomer-dimer equilibrium studied by NMR: a model for changes in dynamics upon target binding

Low molecular weight protein tyrosine phosphatase (LMW-PTP) dimerizes in the phosphate-bound state in solution with a dissociation constant of K(d)=1.5(+/-0.1)mM and an off-rate on the order of 10(4)s(-1). 1H and 15N NMR chemical shifts identify the dimer interface, which is in excellent agreement with that observed in the crystal structure of the dimeric S19A mutant. Two tyrosine residues of each molecule interact with the active site of the other molecule, implying that the dimer may be taken as a model for a complex between LMW-PTP and a target protein. 15N relaxation rates for the monomeric and dimeric states were extrapolated from relaxation data acquired at four different protein concentrations. Relaxation data of satisfactory precision were extracted for the monomer, enabling model-free analyses of backbone fluctuations on pico- to nanosecond time scales. The dimer relaxation data are of lower quality due to extrapolation errors and the possible presence of higher-order oligomers at higher concentrations. A qualitative comparison of order parameters in the monomeric and apparent dimeric states shows that loops forming the dimer interface become rigidified upon dimerization. Qualitative information on monomer-dimer exchange and intramolecular conformational exchange was obtained from the concentration dependence of auto- and cross-correlated relaxation rates. The loop containing the catalytically important Asp129 fluctuates between different conformations in both the monomeric and dimeric (target bound) states. The exchange rate compares rather well with that of the catalyzed reaction step, supporting existing hypotheses that catalysis and enzyme dynamics may be coupled. The side-chain of Trp49, which is important for substrate specificity, exhibits conformational dynamics in the monomer that are largely quenched upon formation of the dimer, suggesting that binding is associated with the selection of a single side-chain conformer.     

A diverse view of protein dynamics from NMR studies of HIV-1 protease flaps

Internal motion in proteins fulfills a multitude of roles in biological processes. NMR spectroscopy has been applied to elucidate protein dynamics at the atomic level, albeit at a low resolution, and is often complemented by molecular dynamics simulation. However, it is critical to justify the consistency between simulation results and conclusions often drawn from multiple experiments in which uncertainties arising from assumed motional models may not be explicitly evaluated. To understand the role of the flaps of HIV-1 protease dimer in substrate recognition and protease function, many molecular dynamics simulations have been performed. The simulations have resulted in various proposed models of the flap dynamics, some of which are more consistent than others with our working model previously derived from experiments. However, using the working model to discriminate among the simulation results is not straightforward because the working model was derived from a combination of NMR experiments and crystal structure data. In this study, we use the NMR chemical shifts and relaxation data of the protease "monomer" rather than structural data to narrow down the possible conformations of the flaps of the "dimer". For the first time, we show that the tips of the flaps in the unliganded protease dimer interact with each other in solution. Accordingly, we discuss the consistency of the simulations with the model derived from all experimental data.     

Under Pressure That Splits a Family in Two. The Case of Lipocalin Family

The lipocalin family is typically composed of small proteins characterized by a range of different molecular recognition properties. Odorant binding proteins (OBPs) are a class of proteins of this family devoted to the transport of small hydrophobic molecules in the nasal mucosa of vertebrates. Among OBPs, bovine OBP (bOBP) is of great interest for its peculiar structural organization, characterized by a domain swapping of its two monomeric subunits. The effect of pressure on unfolding and refolding of native dimeric bOBP and of an engineered monomeric form has been investigated by theoretical and experimental studies under pressure. A coherent model explains the pressure-induced protein structural changes: i) the substrate-bound protein stays in its native configuration up to 330 MPa, where it loses its substrate; ii) the substrate-free protein dissociates into monomers at 200 MPa; and iii) the monomeric substrate-free form unfolds at 120 MPa. Molecular dynamics simulations showed that the pressure-induced tertiary structural changes that accompany the quaternary structural changes are mainly localized at the interface between the monomers. Interestingly, pressure-induced unfolding is reversible, but dimerization and substrate binding can no longer occur. The volume of the unfolding kinetic transition state of the monomer has been found to be similar to that of the folded state. This suggests that its refolding requires relatively large structural and/or hydrational changes, explaining thus the relatively low stability of the monomeric form of this class of proteins.     

Isolation, characterization, and antimicrobial properties of bovine oligosaccharide-binding protein. A microbicidal granule protein of eosinophils and neutrophils

Peptidoglycan recognition proteins (PGRPs) constitute a recently characterized family of pattern-recognition molecules that are conserved from insects to humans and are implicated in mammalian innate immunity. Here we report the isolation, characterization, cDNA cloning, and antimicrobial activities of a bovine PGRP ortholog termed bovine oligosaccharide-binding protein (bOBP). Milligram quantities of bOBP were purified from peripheral leukocytes, thus allowing for the characterization of the disulfide array and for determining the in vitro antimicrobial activities of the native protein. Of the tissues analyzed, bOBP mRNA was detected only in bone marrow where the protein is synthesized as a 190 amino acid precursor. The mature 169 amino acid protein is stored in the cytoplasmic granules of neutrophils and eosinophils but is absent from lymphocytes, monocytes, and platelets. bOBP was microbicidal for Gram-positive and Gram-negative bacteria and yeast at low micromolar concentrations. The finding that bOBP was microbicidal for organisms in which peptidoglycan is absent (Cryptococcus neoformans) or buried (Salmonella typhimurium) indicates that previous conclusions about the specificity of peptidoglycan recognition proteins must be reevaluated and suggests that other envelope components may mediate the antimicrobial action of PGRP family members.     

SAXS studies of human protein HC (alpha1-microglobulin)

Protein HC is a low molecular weight heterogeneous glycoprotein widely distributed in human body fluids and belonging to the lipocalin superfamily. The monomer contains a single (183 amino acid residues long) peptide chain with 3 cysteine residues (2 of which form a disulfide bridge) and is glycosylated. The molecular mass of the glycosylated protein is about 27 kDa. Native gel electrophoresis results revealed partial oligomerisation of protein HC, which therefore was analysed by gel filtration. Two forms (monomer and dimer) of the protein HC were isolated. The SAXS data were recorded on an X33 camera using synchrotron radiation (lambda=0.15 nm) at X33 beamline at the DORIS storage ring of DESY (Hamburg, Germany). Solution scattering results permitted determination of the structural parameters of both forms of the protein studied. The monomer of protein HC is characterised by a radius of gyration R(G)= 2.20 nm and D(max)=6.3 nm and the dimer by R(G)=2.99 nm and D(max)=9.5 nm.     

Loopless Rop: structure and dynamics of an engineered homotetrameric variant of the repressor of primer protein

The repressor of primer (Rop) protein has become a steady source of surprises concerning the relationship between the sequences and the structures of several of its mutants and variants. Here we add another piece to the puzzle of Rop by showing that an engineered deletion mutant of the protein (corresponding to a deletion of residues 30-34 of the wild-type protein and designed to restore the heptad periodicity at the turn region) results in a complete reorganization of the bundle which is converted from a homodimer to a homotetramer. In contrast (and as previously shown), a two-residue insertion, which also restores the heptad periodicity, is essentially identical with wild-type Rop. The new deletion mutant structure is a canonical, left-handed, all-antiparallel bundle with a completely different hydrophobic core and distinct surface properties. The structure agrees and qualitatively explains the results from functional, thermodynamic, and kinetic studies which indicated that this deletion mutant is a biologically inactive hyperstable homotetramer. Additional insight into the stability and dynamics of the mutant structure has been obtained from extensive molecular dynamics simulations in explicit water and with full treatment of electrostatics.     

The Quaternary Structure of the Recombinant Bovine Odorant-Binding Protein Is Modulated by Chemical Denaturants

A large group of odorant-binding proteins (OBPs) has attracted great scientific interest as promising building blocks in constructing optical biosensors for dangerous substances, such as toxic and explosive molecules. Native tissue-extracted bovine OBP (bOBP) has a unique dimer folding pattern that involves crossing the α-helical domain in each monomer over the other monomer’s β-barrel. In contrast, recombinant bOBP maintaining the high level of stability inherent to native tissue bOBP is produced in a stable native-like state with a decreased tendency for dimerization and is a mixture of monomers and dimers in a buffered solution. This work is focused on the study of the quaternary structure and the folding-unfolding processes of the recombinant bOBP in the absence and in the presence of guanidine hydrochloride (GdnHCl). Our results show that the recombinant bOBP native dimer is only formed at elevated GdnHCl concentrations (1.5 M). This process requires re-organizing the protein structure by progressing through the formation of an intermediate state. The bOBP dimerization process appears to be irreversible and it occurs before the protein unfolds. Though the observed structural changes for recombinant bOBP at pre-denaturing GdnHCl concentrations show a local character and the overall protein structure is maintained, such changes should be considered where the protein is used as a sensitive element in a biosensor system.     

Rab32 is an A-kinase anchoring protein and participates in mitochondrial dynamics

A-kinase anchoring proteins (AKAPs) tether the cAMP-dependent protein kinase (PKA) and other signaling enzymes to distinct subcellular organelles. Using the yeast two-hybrid approach, we demonstrate that Rab32, a member of the Ras superfamily of small molecular weight G-proteins, interacts directly with the type II regulatory subunit of PKA. Cellular and biochemical studies confirm that Rab32 functions as an AKAP inside cells. Anchoring determinants for PKA have been mapped to sites within the conserved alpha5 helix that is common to all Rab family members. Subcellular fractionation and immunofluorescent approaches indicate that Rab32 and a proportion of the cellular PKA pool are associated with mitochondria. Transient transfection of a GTP binding-deficient mutant of Rab32 promotes aberrant accumulation of mitochondria at the microtubule organizing center. Further analysis of this mutant indicates that disruption of the microtubule cytoskeleton results in aberrantly elongated mitochondria. This implicates Rab32 as a participant in synchronization of mitochondrial fission. Thus, Rab32 is a dual function protein that participates in both mitochondrial anchoring of PKA and mitochondrial dynamics.     

Exploration of disorder in protein structures by X-ray restrained molecular dynamics

Conformational disorder in crystal structures of ribonuclease-A and crambin is studied by including two independent structures in least-squares optimizations against X-ray data. The optimizations are carried out by X-ray restrained molecular dynamics (simulated annealing refinement) and by conventional least-squares optimization. Starting from two identical structures, the optimizations against X-ray data lead to significant deviations between the two, with rms backbone displacements of 0.45 A for refinement of ribonuclease at 1.53 A resolution, and 0.31 A for crambin at 0.945 A. More than 15 independent X-ray restrained molecular dynamics runs have been carried out for ribonuclease, and the displacements between the resulting structures are highly reproducible for most atoms. These include residues with two or more conformations with significant dihedral angle differences and alternative hydrogen bonding, as well as groups of residues that undergo displacements that are suggestive of rigid-body librations. The crystallographic R-values obtained are approximately 13%, as compared to 15.3% for a comparable refinement with a single structure. Least-squares optimization without an intervening restrained molecular dynamics stage is sufficient to reproduce most of the observed displacements. Similar results are obtained for crambin, where the higher resolution of the X-ray data allows for refinement of unconstrained individual anisotropic temperature factors. These are shown to be correlated with the displacements in the two-structure refinements.     

Identifying the protein folding nucleus using molecular dynamics

Molecular dynamics simulations of folding in an off-lattice protein model reveal a nucleation scenario, in which a few well-defined contacts are formed with high probability in the transition state ensemble of conformations. Their appearance determines folding cooperativity and drives the model protein into its folded conformation. Amino acid residues participating in those contacts may serve as "accelerator pedals" used by molecular evolution to control protein folding rate.     

Partition of protein solvation into group contributions from molecular dynamics simulations

Linear response theory coupled to molecular dynamics simulations with an explicit solvent representation is used to derive fractional contributions of amino acid residues to the solvation of proteins. The new fractional methods developed here are compared with standard approaches based on empirical 1D and 3D statistical potentials, as well as with estimates obtained from the analysis of classical molecular interaction potentials. The new fractional methods, which have a clear physical basis and explicitly account for the effects due to protein structure and flexibility, provide an accurate picture of the contribution to solvation of different regions of the protein.     

The transmembrane domain of the oncogenic mutant ErbB-2 receptor: a structure obtained from site-specific infrared dichroism and molecular dynamics

ErbB-2 is a member of the family of epidermal growth factor receptors, which shows an oncogenic mutation in the rat gene neu, Val664Glu in the transmembrane domain that causes permanent dimerisation and subsequently leads to uncontrollable cell division and tumour formation. We have obtained the alpha-helical structure of the mutant transmembrane domain dimer experimentally with site-specific infrared dichroism (SSID) based on six transmembrane peptides with 13C18O carbonyl group-labelled residues. The derived orientational data indicate a local helix tilt ranging from 28(+/-6) degrees to 22(+/-4) degrees. Altogether using orientational constraints from SSID and experimental alpha-helical constraints while performing a systematic conformational search including molecular dynamics simulation in a lipid bilayer, we have obtained a unique experimentally defined atomic structure. The resulting structure consists of a right handed alpha-helical bundle with the residues Ile659, Val663, Leu667, Ile671, Val674 and Leu679 in the dimerisation interface. The right-handed bundle is in contrast to the left-handed structures obtained in previous modelling efforts. In order to facilitate tight helical packing, the spacious Glu664 residues do not interact directly but with water molecules that enter the bilayer.