Mutant bovine odorant-binding protein: Temperature affects the protein stability and dynamics as revealed by infrared spectroscopy and molecular dynamics simulations

Bovine odorant-binding protein (bOBP), a member of the lipocalin family, presents the so-called 3D "domain-swapped" protein structure. In fact, in solution, it appears as a dimer in which each monomer is composed by the classical lipocalin fold, with a central beta-barrel followed by a stretch of residues and the alpha-helix domain protruding out of the barrel and crossing the dimer interface. Recently, a deswapped mutant form of bOBP was obtained, in which a Gly residue was inserted after position 121 and the two residues in position 64 and 156 were replaced by Cys residues for restoring the disulfide bridge common to the lipocalin family. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the effect of temperature on the structural stability and conformational dynamics of the mutant bOBP. The spectroscopic and molecular simulation data pointed out that the hydrophobic regions of the protein matrix appear to be an important factor for the protein stability and integrity. In addition, it was also found that the mutant bOBP is significantly stabilized by the binding of the ligand, which may have an impact on the biological function of bOBP. The obtained results will allow for a better use of this protein as probe for the design of advanced protein-based biosensors for the detection of compounds used in the fabrication of explosive powders.     

Probing protein structure and dynamics by second-derivative ultraviolet absorption analysis of cation-{pi} interactions

We describe an alternate approach for studying protein structure using the detection of ultraviolet (UV) absorbance peak shifts of aromatic amino acid side chains induced by the presence of salts. The method is based on the hypothesis that salt cations (Li+, Na+, and Cs+) of varying sizes can differentially diffuse through protein matrices and interact with benzyl, phenyl, and indole groups through cation-pi interactions. We have investigated the potential of this method to probe protein dynamics by measuring high resolution second-derivative UV spectra as a function of salt concentration for eight proteins of varying physical and chemical properties and the N-acetylated C-ethyl esterified amino acids to represent totally exposed side chains. We show that small shifts in the wavelength maxima for Phe, Tyr, and Trp in the presence of high salt concentrations can be reliably measured and that the magnitude and direction of the peak shifts are influenced by several factors, including protein size, charge, and the local environment and solvent accessibility of the aromatic groups. Evaluating the empirical UV spectral data in light of known protein structural information shows that probing cation-pi interactions in proteins reveals unique information about the influence of structure on aromatic side chain spectroscopic behavior.     

Dark-interval relaxation kinetics (DIRK) of absorbance changes as a quantitative probe of steady-state electron transfer

We introduce a new, non-invasive technique to measure linear electron transfer in intact leaves under steady-state illumination. Dark-interval relaxation kinetic or ‘DIRK’ analysis is based on measurements of the initial rates of relaxation of steady-state absorbance signals upon a rapid light-dark transition. We show that estimates of electron flux by DIRK analysis of absorbance signals, reflecting redox changes in the photosynthetic electron transfer chain, can yield quantitative information about photosynthetic flux when the light-dependent partitioning of electrons among redox components of the electron transfer chain are considered. This concept is modeled in computer simulations and then demonstrated in vivo with tobacco plants under non-photorespiratory conditions resulting in linear relationships between DIRK analysis and gross carbon assimilation (A_G). Estimation based on DIRK analysis of the number of electrons transferred through the photosynthetic apparatus for each CO₂ fixed was within 20% of the theoretical value. Possible errors and future improvements are discussed. We conclude that the DIRK method represents a useful tool to address issues such as plant stress and photosynthetic regulation. This revised version was published online in June 2006 with corrections to the Cover Date.     

Comparison of the dynamics of the membrane-bound form of fd coat protein in micelles and in bilayers by solution and solid-state nitrogen-15 nuclear magnetic resonance spectroscopy

Solid-state and solution 15N nuclear magnetic resonance experiments on uniformly and specifically 15N labeled coat protein in phospholipid bilayers and in detergent micelles are used to describe the dynamics of the membrane-bound form of the protein. The residues in the N- and C-terminal portions of the coat protein in both phospholipid bilayers and in detergent micelles are mobile, while those in the hydrophobic midsection are immobile. There is evidence for a gradient of mobility in the C-terminal region of the coat protein in micelles; at 25 degrees C only the last two residues are mobile on the 10(9)-Hz timescale, while the last six to eight residues appear to be mobile on slower timescales and highly mobile at higher temperatures. Since all of the C-terminal residues are immobile in the virus particles, the mobility of these residues in the membrane-bound form of the protein may be important for the formation of protein-DNA interactions in the assembly process.     

Proportion of solvent-exposed amino acids in a protein and rate of protein evolution

Translational selection, including gene expression, protein abundance, and codon usage bias, has been suggested as the single dominant determinant of protein evolutionary rate in yeast. Here, we show that protein structure is also an important determinant. Buried residues, which are responsible for maintaining protein structure or are located on a stable interaction surface between 2 subunits, are usually under stronger evolutionary constraints than solvent-exposed residues. Our partial correlation analysis shows that, when whole proteins are included, the variance of evolutionary rate explained by the proportion of solvent-exposed residues (P(exposed)) can reach two-thirds of that explained by translational selection, indicating that P(exposed) is the most important determinant of protein evolutionary rate next only to translational selection. Our result suggests that proteins with many residues under selective constraint (e.g., maintaining structure or intermolecular interaction) tend to evolve slowly, supporting the "fitness (functional) density" hypothesis.     

Characterization of threonine side chain dynamics in an antifreeze protein using natural abundance 13C NMR spectroscopy

The dynamics of threonine side chains of the Tenebrio molitor antifreeze protein (TmAFP) were investigated using natural abundance (13)C NMR. In TmAFP, the array of threonine residues on one face of the protein is responsible for conferring its ability to bind crystalline ice and inhibit its growth. Heteronuclear longitudinal and transverse relaxation rates and the [(1)H]-(13)C NOE were determined in this study. The C alpha H relaxation measurements were compared to the previously measured (15)N backbone parameters and these are found to be in agreement. For the analysis of the threonine side chain motions, the model of restricted rotational diffusion about the chi(1) dihedral angle was employed [London and Avitabile (1978) J. Am. Chem. Soc., 100, 7159-7165]. We demonstrate that the motion experienced by the ice binding threonine side chains is highly restricted, with an approximate upper limit of less than +/-25 degrees.     

Molecular dynamics studies of a DNA-binding protein: 2. An evaluation of implicit and explicit solvent models for the molecular dynamics simulation of the Escherichia coli trp repressor.

Although aqueous simulations with periodic boundary conditions more accurately describe protein dynamics than in vacuo simulations, these are computationally intensive for most proteins. Trp repressor dynamic simulations with a small water shell surrounding the starting model yield protein trajectories that are markedly improved over gas phase, yet computationally efficient. Explicit water in molecular dynamics simulations maintains surface exposure of protein hydrophilic atoms and burial of hydrophobic atoms by opposing the otherwise asymmetric protein-protein forces. This properly orients protein surface side chains, reduces protein fluctuations, and lowers the overall root mean square deviation from the crystal structure. For simulations with crystallographic waters only, a linear or sigmoidal distance-dependent dielectric yields a much better trajectory than does a constant dielectric model. As more water is added to the starting model, the differences between using distance-dependent and constant dielectric models becomes smaller, although the linear distance-dependent dielectric yields an average structure closer to the crystal structure than does a constant dielectric model. Multiplicative constants greater than one, for the linear distance-dependent dielectric simulations, produced trajectories that are progressively worse in describing trp repressor dynamics. Simulations of bovine pancreatic trypsin were used to ensure that the trp repressor results were not protein dependent and to explore the effect of the nonbonded cutoff on the distance-dependent and constant dielectric simulation models. The nonbonded cutoff markedly affected the constant but not distance-dependent dielectric bovine pancreatic trypsin inhibitor simulations. As with trp repressor, the distance-dependent dielectric model with a shell of water surrounding the protein produced a trajectory in better agreement with the crystal structure than a constant dielectric model, and the physical properties of the trajectory average structure, both with and without a nonbonded cutoff, were comparable.     

A proteome-scale analysis of vertebrate protein amino acid occurrence: Thermoadaptation shows a correlation with protein solvation but less so with dynamics

Despite differences in behaviors and living conditions, vertebrate organisms share the great majority of proteins, often with subtle differences in amino acid sequence. Here, we present a simple way to analyze the difference in amino acid occurrence by comparing highly homologous proteins on a subproteome level between several vertebrate model organisms. Specifically, we use this method to identify a pattern of amino acid conservation as well as a shift in amino acid occurrence between homeotherms (warm-blooded species) and poikilotherms (cold-blooded species). Importantly, this general analysis and a specific example further establish a broad correlation, if not likely connection between the thermal adaptation of protein sequences and two of their physical features: on average a change in their protein dynamics and, even more strongly, in their solvation. For poikilotherms, such as frog and fish, the lower body temperature is expected to increase the protein–protein interaction due to a decrease in protein internal dynamics. In order to counteract the tendency for enhanced binding caused by low temperatures, poikilotherms enhance the solvation of their proteins by favoring polar amino acids. This feature appears to dominate over possible changes in dynamics for some proteins. The results suggest that a general trend for amino acid choice is part of the mechanism for thermoadaptation of vertebrate organisms at the molecular level.     

Relative stability of protein structures determined by X-ray crystallography or NMR spectroscopy: a molecular dynamics simulation study

The relative stability of protein structures determined by either X-ray crystallography or nuclear magnetic resonance (NMR) spectroscopy has been investigated by using molecular dynamics simulation techniques. Published structures of 34 proteins containing between 50 and 100 residues have been evaluated. The proteins selected represent a mixture of secondary structure types including all alpha, all beta, and alpha/beta. The proteins selected do not contain cysteine-cysteine bridges. In addition, any crystallographic waters, metal ions, cofactors, or bound ligands were removed before the systems were simulated. The stability of the structures was evaluated by simulating, under identical conditions, each of the proteins for at least 5 ns in explicit solvent. It is found that not only do NMR-derived structures have, on average, higher internal strain than structures determined by X-ray crystallography but that a significant proportion of the structures are unstable and rapidly diverge in simulations.     

Optimization and anti-inflammatory evaluation of methyl gallate derivatives as a myeloid differentiation protein 2 inhibitor

Myeloid differentiation protein 2 (MD2) is a co-receptor of toll-like receptor 4 (TLR4) responsible for the recognition of lipopolysaccharide (LPS) and mediates a series of TLR4-dependent inflammatory responses in inflammatory lung diseases including acute lung injury (ALI). Targeting MD2 thus may provide a therapeutic strategy against these lung diseases. In this study, we identified a novel compound 4k with the potent anti-inflammatory activity among 39 methyl gallate derivatives (MGDs). MGD 4k exhibited a high binding affinity to MD2, which in turn prevented the formation of the LPS/MD2/TLR4 complex. In addition, MGD 4k significantly reversed the upregulation of LPS-induced inflammatory mediators such as tumor necrosis factor-α, interleukin-6, intracellular adhesion molecule-1, vascular cell adhesion molecule-1, and monocyte chemoattractant protein-1 in vitro and in vivo. Mechanistically, MGD 4k performed anti-inflammatory function by inactivating JNK, ERK and p38 signaling pathways. Taken together, our study identified MGD 4k as a novel potential therapeutic agent for ALI through inhibiting MD2, inflammatory responses, and major inflammation-associated signaling pathways.     

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.     

The Ca2+ response of a smart forisome protein is dependent on polymerization

Forisomes are giant self‐assembling mechanoproteins that undergo reversible structural changes in response to Ca²⁺ and various other stimuli. Artificial forisomes assembled from the monomer MtSEO‐F1 can be used as smart biomaterials, but the molecular basis of their functionality is not understood. To determine the role of protein polymerization in forisome activity, we tested the Ca²⁺ association of MtSEO‐F1 dimers (the basic polymerization unit) by circular dichroism spectroscopy and microscale thermophoresis. We found that soluble MtSEO‐F1 dimers neither associate with Ca²⁺ nor undergo structural changes. However, polarization modulation infrared reflection absorption spectroscopy revealed that aggregated MtSEO‐F1 dimers and fully‐assembled forisomes associate with Ca²⁺, allowing the hydration of poorly‐hydrated protein areas. A change in the signal profile of complete forisomes indicated that Ca²⁺ interacts with negatively‐charged regions in the protein complexes that only become available during aggregation. We conclude that aggregation is required to establish the Ca²⁺ response of forisome polymers.     

Amide proton exchange measurements as a probe of the stability and dynamics of the N-terminal domain of the ribosomal protein L9: comparison with the intact protein.

Amide H/D exchange rates have been measured for the N-terminal domain of the ribosomal protein L9, residues 1-56. The rates were measured at pD 3.91, 5.03, and 5.37. At pD 5.37, 18 amides exchange slowly enough to give reliable rate measurements. At pD 3.91, seven additional residues could be followed. The exchange is shown to occur by the EX2 mechanism for all conditions studied. The rates for the N-terminal domain are very similar to those previously measured for the corresponding region in the full-length protein (Lillemoen J et al., 1997, J Mol Biol 268:482-493). In particular, the rates for the residues that we have shown to exchange via global unfolding in the N-terminal domain agree within the experimental error with the rates measured by Hoffman and coworkers, suggesting that the structure of the domain is preserved in isolation and that the stability of the isolated domain is comparable to the stability of this domain in intact L9.     

Role of the electrostatic energy between the photo-products and ionic groups on the photoinduced electron transfer rates from aromatic amino acids to the excited flavin in five single-point substitution isoforms of the charged amino acid residue-13 in the FMN-binding protein

Effect of the charge (negative, positive or neutral) of amino acid residue-13 on the photoinduced electron transfer (ET) from Trp32, Tyr35 and Trp106 to the excited isoalloxazine was evaluated in the flavin mononucleotide-binding protein from Desulfovibrio vulgaris isolate Miyazaki F (DvFBP). The protein structures of the wild type and the four isoforms where glutamic acid-13 is replaced with lysine (E13K), arginine (E13R), threonine (E13T) and glutamine (E13Q) in aqueous solution were obtained by molecular dynamics simulation. The distances between the amino acid residue-13 and isoalloxazine (Iso), and between the amino acid residue-13 and the ET donors were longer than 1 nm. The ET rates were evaluated with the Kakitani and Mataga model (KM theory) from their ultrafast fluorescence dynamics by means of a non-linear least squares method. Electrostatic (ES) energies between the photo-products and other ionic groups in the proteins markedly varied among ET donors and among the DvFBP isoforms, while the other physical quantities related to the ET rates, the solvent reorganisation and ES energies between the Iso anion and the donor cations did not vary much between the proteins and donors. A plot of the logarithmic ET rates versus either the total free energy gaps or the net ES energies between the photo-products and the other ionic groups both displayed a parabolic function and so the net ES energies are an important influential factor upon the ET rate, in addition to the donor–acceptor distance.     

NMR solution structure and backbone dynamics of domain III of the E protein of tick-borne Langat flavivirus suggests a potential site for molecular recognition

Flaviviruses cause many human diseases, including dengue fever, yellow fever, West Nile viral encephalitis, and hemorrhagic fevers, and are transmitted to their vertebrate hosts by infected mosquitoes and ticks. Domain III of the envelope protein (E-D3) is considered to be the primary viral determinant involved in the virus-host-cell receptor interaction, and thus represents an excellent target for antiviral drug development. Langat (LGT) virus is a naturally attenuated BSL-2 TBE virus and is a model for the pathogenic BSL-3 and BSL-4 viruses in the serogroup. We have determined the solution structure of LGT-E-D3 using heteronuclear NMR spectroscopy. The backbone dynamics of LGT-E-D3 have been investigated using 15N relaxation measurements. A detailed analysis of the solution structure and dynamics of LGT-E-D3 suggests potential residues that could form a surface for molecular recognition, and thereby represent a target site for antiviral therapeutics design.     

Exploring the effect of N308D mutation on protein tyrosine phosphatase-2 cause gain-of-function activity by a molecular dynamics study

One of the most common protein tyrosine phosphatase-2 (SHP2) mutations in Noonan syndrome is the N308D mutation, and it increases the activity of the protein. However, the molecular basis of the activation of N308D mutation on SHP2 conformations is poorly understood. Here, molecular dynamic simulations were performed on SHP2 and SHP2-N308D to explore the effect of N308D mutation on SHP2 cause gain of function activity, respectively. The principal component analysis, dynamic cross-correlation map, secondary structure analysis, residue interaction networks, and solvent accessible surface area analysis suggested that the N308D mutation distorted the residues interactions network between the allosteric site (residue Gly244-Gly246) and C-SH2 domain, including the hydrogen bond formation and the binding energy. Meanwhile, the activity of catalytic site (residue Gly503-Val505) located in the Q-loop in mutant increased due to this region's high fluctuations. Therefore, the substrate had more chances to access to the catalytic activity site of the precision time protocol domain of SHP2-N308D, which was easy to be exposed. In addition, we had speculated that the Lys244 located in the allosteric site was the key residue which lead to the protein conformation changes. Consequently, overall calculations presented in this study ultimately provide a useful understanding of the increased activity of SHP2 caused by the N308D mutation.     

Use of phospholipid transfer protein as a probe to study the lipid dynamics and alkaline phosphatase activity in the brush border membrane of human term placenta

Incubation of placental brush border membrane (BBM) along with sonicated vesicles of exogenous lipids (egg yolk PC) in the presence of phospholipid-transfer protein (PL-TP) showed a decrease in the alkaline phosphatase activity due to the change in the membrane micro-environment, such as fluidity. Effect of substrate concentration was tested by Lineweaver-Burk plot, which showed decreased V(max) and K(M). The effect of temperature was probed by the Arrhenius plot, which showed no change in transition temperature, but a decline in the energy of activation both below and above the transition temperature. The protein-catalyzed transfer of phospholipid from the donor unilamellar vesicles resulted in a substantial increase in the BBM phospholipid and a net decrease in cholesterol/phospholipid molar ratio. The change in membrane fluidity was assessed by translational as well as rotational diffusion of membrane extrinsic fluorescent probes, pyrene and diphenyl-hexatriene. An increased lateral mobility was recorded by the increased pyrene excimer formation. A decrease in fluorescent polarization of diphenyl-hexatriene was observed, which led to the decrease in fluorescence anisotropy and order parameter, and therefore, an increase in membrane fluidity (rotational diffusion). Mean anisotropy parameter was also decreased in the presence of PL-TP. Thus, the placental BBM alkaline phosphatase activity showed a distinct lipid dependence which may have important physiological consequences.     

Eliminating the six N-terminal amino acids of the caspase 3 large subunit improved production of a biologically active IL2-Caspase3 chimeric protein

Designing a chimeric protein and developing a procedure for its stable production as a biologically active protein, are key steps in its potential application to clinical trails. IL2-Caspase3 chimeric protein designed to target activated T lymphocytes was found to be a promising molecule for targeted treatment, however was found to be difficult to produce as a biological active molecule. Thus, we designed a new version of the molecule, IL2-Caspase3s, in which six amino acids (aa 29-34) from the N-terminus of the large subunit of caspase 3 were excluded. Repeated expressions, productions, and partial purifications of the IL2-Caspase3s yielded reproducible batches with consistent results. We found that IL2-Caspase3s causes cell death in a specific, dose-, and time-dependent manner. Cell death due to IL2-Caspase3s is caused by apoptosis. This improved and biologically stable IL2-Caspase3s chimeric protein may be developed in the future for clinical trails as a promising therapy for several pathologies involving activated T-cells. Moreover, this truncated caspase 3 sequence, lacking the N-terminal six amino acids of its large subunit, may be used in other caspase 3-based chimeric proteins targeted against various human diseases, using the appropriate targeting moiety.