Structure refinement of protein model decoys requires accurate side‐chain placement

In this study, the application of temperature‐based replica‐exchange (T‐ReX) simulations for structure refinement of decoys taken from the I‐TASSER dataset was examined. A set of eight nonredundant proteins was investigated using self‐guided Langevin dynamics (SGLD) with a generalized Born implicit solvent model to sample conformational space. For two of the protein test cases, a comparison of the SGLD/T‐ReX method with that of a hybrid explicit/implicit solvent molecular dynamics T‐ReX simulation model is provided. Additionally, the effect of side‐chain placement among the starting decoy structures, using alternative rotamer conformations taken from the SCWRL4 modeling program, was investigated. The simulation results showed that, despite having near‐native backbone conformations among the starting decoys, the determinant of their refinement is side‐chain packing to a level that satisfies a minimum threshold of native contacts to allow efficient excursions toward the downhill refinement regime on the energy landscape. By repacking using SCWRL4 and by applying the RWplus statistical potential for structure identification, the SGLD/T‐ReX simulations achieved refinement to an average of 38% increase in the number of native contacts relative to the original I‐TASSER decoy sets and a 25% reduction in values of C_α root‐mean‐square deviation. The hybrid model succeeded in obtaining a sharper funnel to low‐energy states for a modeled target than the implicit solvent SGLD model; yet, structure identification remained roughly the same. Without meeting a threshold of near‐native packing of side chains, the T‐ReX simulations degrade the accuracy of the decoys, and subsequently, refinement becomes tantamount to the protein folding problem. Proteins 2013. 2012 Published by Wiley Periodicals, Inc.     

Molecular dynamics simulation of intrinsically disordered proteins

Intrinsically disordered proteins are biomolecules that do not have a definite 3D structure; therefore, their dynamical simulation cannot start from a known list of atomistic positions, such as a Protein Data Bank file. We describe a method to start a computer simulation of these proteins. The first step of the procedure is the creation of a multi-rod configuration of the molecule, derived from its primary sequence. This structure is dynamically evolved in vacuo until its gyration radius reaches the experimental average value; at this point solvent molecules, in explicit or implicit implementation, are added to the protein and a regular molecular dynamics simulation follows. We have applied this procedure to the simulation of tau, one of the largest totally disordered proteins.     

Structure refinement of membrane proteins via molecular dynamics simulations

A refinement protocol based on physics‐based techniques established for water soluble proteins is tested for membrane protein structures. Initial structures were generated by homology modeling and sampled via molecular dynamics simulations in explicit lipid bilayer and aqueous solvent systems. Snapshots from the simulations were selected based on scoring with either knowledge‐based or implicit membrane‐based scoring functions and averaged to obtain refined models. The protocol resulted in consistent and significant refinement of the membrane protein structures similar to the performance of refinement methods for soluble proteins. Refinement success was similar between sampling in the presence of lipid bilayers and aqueous solvent but the presence of lipid bilayers may benefit the improvement of lipid‐facing residues. Scoring with knowledge‐based functions (DFIRE and RWplus) was found to be as good as scoring using implicit membrane‐based scoring functions suggesting that differences in internal packing is more important than orientations relative to the membrane during the refinement of membrane protein homology models.     

Ultrafast molecular dynamics of biofuel extraction for microalgae and bacteria milking: blocking membrane folding pathways to damaged lipid-bilayer conformations with nanomicelles

Cell milking is a 100% renewable green energy for CO₂ by extraction of biofuels inside the cytosol of photosynthetic micro-organisms as microalgae and bacteria. The cells are exposed to a hydrophobic solvent forming holes and cracks through their membranes from which the biofuels can leak out. In protein folding, the goal would be to find pathways to the unique functional protein conformer. However, in the lipid-bilayer interaction with the solvent for milking, the objective is to block the pathways for damaged membrane conformations of low free energy with undesired nanostructures, using the solvent properties, as shown with an ab initio structural bioinformatic model. Statistical thermodynamics is used to compute the free energy (including entropy) from the molecular dynamics trajectory of the biomolecular system with many conformational changes. This model can be extended to the general problem of biomolecules folding as for proteins and nucleic acids. Using an adaptation of the Einstein diffusion law, the conformational change dynamics of the lipid bilayer depends on the two diffusion coefficients of the solvent: D₁ before the irreversible folding transition time and the much smaller D₂ thereafter. In contrast to the n-hexane and n-heptane hydrocarbons of smaller size, the residual D₂?=?4.7?×?10^?7?cm²/s of the n-decane solvent, with the highest partition coefficient among the three extractors, is the only to present a D₂ value that is significantly below the critical threshold of 10^?6?cm²/s. Therefore, the membrane would resist to long hydrocarbons and the exposed cells would remain viable for milking.     

Distance restraints from crosslinking mass spectrometry: Mining a molecular dynamics simulation database to evaluate lysine–lysine distances

Integrative structural biology attempts to model the structures of protein complexes that are challenging or intractable by classical structural methods (due to size, dynamics, or heterogeneity) by combining computational structural modeling with data from experimental methods. One such experimental method is chemical crosslinking mass spectrometry (XL‐MS), in which protein complexes are crosslinked and characterized using liquid chromatography‐mass spectrometry to pinpoint specific amino acid residues in close structural proximity. The commonly used lysine‐reactive N‐hydroxysuccinimide ester reagents disuccinimidylsuberate (DSS) and bis(sulfosuccinimidyl)suberate (BS³) have a linker arm that is 11.4 Å long when fully extended, allowing C_α (alpha carbon of protein backbone) atoms of crosslinked lysine residues to be up to ～24 Å apart. However, XL‐MS studies on proteins of known structure frequently report crosslinks that exceed this distance. Typically, a tolerance of ～3 Å is added to the theoretical maximum to account for this observation, with limited justification for the chosen value. We used the Dynameomics database, a repository of high‐quality molecular dynamics simulations of 807 proteins representative of diverse protein folds, to investigate the relationship between lysine–lysine distances in experimental starting structures and in simulation ensembles. We conclude that for DSS/BS³, a distance constraint of 26–30 Å between C_α atoms is appropriate. This analysis provides a theoretical basis for the widespread practice of adding a tolerance to the crosslinker length when comparing XL‐MS results to structures or in modeling. We also discuss the comparison of XL‐MS results to MD simulations and known structures as a means to test and validate experimental XL‐MS methods.     

Comparison of Generalised Born/Surface Area with Periodic Boundary Simulations to Study Protein Unfolding

Simulations comparing the rapid unfolding behaviour of the model protein barnase under explicit and implicit solvent systems have been undertaken in order to validate a faster implicit solvent method for studying proteins which are kinetically stable in silico. A comparison is made between all-atom explicitly solvated simulations of barnase undertaken using Particle Mesh Ewald electrostatic interactions with all-atom implicit solvent simulations undertaken using the generalised born/surface area (GBSA) method with a long non-bonded cut-off. The two explicit solvent unfolding trajectories appear to explore slightly different pathways showing the importance of having statistically valid ensembles which are not accessible from a single trajectory. The 500?K GBSA trajectory is unsuitable for exploring intermediate structures on the unfolding pathway of barnase, as the protein almost immediately jumps to a predominately random coil conformation. However, dropping the temperature to 400?K gives rise to trajectories where the protein is unable to climb out of the energy well containing the first intermediate state, in a reasonable timescale. A similar pattern to the explicit solvent unfolding trajectories is seen in 450?K GBSA runs, with the intermediate states differing between trajectories. The development of computer simulation methods suitable for application to more kinetically stable proteins will offer insight into the atomic detail of the conformational changes associated with protein unfolding diseases.     

The effect of deuterium oxide on the stability of the collagen model peptides H‐(Pro‐Pro‐Gly)10‐OH,H‐(Gly‐Pro‐4(R)Hyp)9‐OH,and Type I collagen

The collagen triple helix has a larger accessible surface area per molecular mass than globular proteins, and therefore potentially more water interaction sites. The effect of deuterium oxide on the stability of collagen model peptides and Type I collagen molecules was analyzed by circular dichroism and differential scanning calorimetry. The transition temperatures (T_m) of the protonated peptide (Pro‐Pro‐Gly)₁₀ were 25.4 and 28.7°C in H₂O and D₂O, respectively. The increase of the T_m of (Pro‐Pro‐Gly)₁₀ measured calorimetrically at 1.0°C min^?1 in a low pH solution from the protonated to the deuterated solvent was 5.1°C. The increases of the T_m for (Gly‐Pro‐4(R)Hyp)₉ and pepsin‐extracted Type I collagen were measured as 4.2 and 2.2°C, respectively. These results indicated that the increase in the T_m in the presence of D₂O is comparable to that of globular proteins, and much less than reported previously for collagen model peptides [Gough and Bhatnagar, J Biomol Struct Dyn 1999, 17, 481–491]. These experimental results suggest that the interaction of water molecules with collagen is similar to the interaction of water with globular proteins, when the ratio of collagen to water is very small and collagen is monomerically dispersed in the solvent. © 2009 Wiley Periodicals, Inc. Biopolymers 93: 93–101, 2010. This article was originally published online as an accepted preprint. The “Published Online” date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com     

Atomistic mechanisms of huntingtin N‐terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q_N) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N‐terminal 17‐amino‐acid sequence (htt^NT) positioned just before the Q_N region is important for the binding of huntingtin to membranes. Through all‐atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt^NTQ_N fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps—approach, reorganization, anchoring, and insertion—that are very diverse at the atomic level. On the membrane, the htt^NT peptide forms a stable α‐helix essentially parallel to the membrane with its nonpolar side‐chains—mainly Leu‐4, Leu‐7, Phe‐11 and Leu‐14—positioned in the hydrophobic core of the membrane. Salt‐bridges involving Glu‐5, Glu‐12, Lys‐6, and Lys‐15, as well as hydrogen bonds involving Thr‐3 and Ser‐13 with the phospholipids also stabilize the structure and orientation of the htt^NT peptide. These observations do not significantly change upon adding the Q_N region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt^NT region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q_N region could also play a key role for the oligomerization of htt^NTQ_N on phospholipid membranes. Proteins 2014; 82:1409–1427. © 2014 Wiley Periodicals, Inc.     

Contributions of water transfer energy to protein‐ligand association and dissociation barriers: Watermap analysis of a series of p38α MAP kinase inhibitors

In our previous work, we proposed that desolvation and resolvation of the binding sites of proteins can serve as the slowest steps during ligand association and dissociation, respectively, and tested this hypothesis on two protein‐ligand systems with known binding kinetics behavior. In the present work, we test this hypothesis on another kinetically‐determined protein‐ligand system—that of p38α and eight Type II BIRB 796 inhibitor analogs. The k_on values among the inhibitor analogs are narrowly distributed (10⁴ ≤ k_on ≤ 10⁵ M^?1 s^?1), suggesting a common rate‐determining step, whereas the k_off values are widely distributed (10^?1 ≤ k_off ≤ 10^?6 s^?1), suggesting a spectrum of rate‐determining steps. We calculated the solvation properties of the DFG‐out protein conformation using an explicit solvent molecular dynamics simulation and thermodynamic analysis method implemented in WaterMap to predict the enthalpic and entropic costs of water transfer to and from bulk solvent incurred upon association and dissociation of each inhibitor. The results suggest that the rate‐determining step for association consists of the transfer of a common set of enthalpically favorable solvating water molecules from the binding site to bulk solvent. The rate‐determining step for inhibitor dissociation consists of the transfer of water from bulk solvent to specific binding site positions that are unfavorably solvated in the apo protein, and evacuated during ligand association. Different sets of unfavorable solvation are evacuated by each ligand, and the observed dissociation barriers are qualitatively consistent with the calculated solvation free energies of those sets.     

Exploring membrane and cytoplasm proteomic responses of Alkalimonas amylolytica N10 to different external pHs with combination strategy of de novo peptide sequencing

Identification of differentially proteomic responses to external pHs would pave an access for understanding of survival mechanisms of bacteria living at extreme pH environment. We cultured Alkalimonas amylolytica N10 (N10), a novel alkaliphilic bacterium found in Lake Chahannor, in media with three different pHs and extracted the correspondent membrane and cytoplasm proteins for proteomic analysis through 2‐DE. The differential 2‐DE spots corresponding to the altered pHs were delivered to MALDI TOF/TOF MS for protein identification. Since the genomic data of strain N10 was unavailable, we encountered a problem at low rate of protein identification with 18.1%. We employed, therefore, a combined strategy of de novo sequencing to analyze MS/MS signals generated from MALDI TOF/TOF MS. A significantly improved rate of protein identification was thus achieved at over than 70.0%. Furthermore, we extensively investigated the expression of these pH‐dependent N10 genes using Western blot and real‐time PCR. The conclusions drawn from immunoblot and mRNA measurements were mostly in agreement with the proteomic observations. We conducted the bioinformatic analysis to all the pH‐dependent N10 proteins and found that some membrane proteins participated in iron transport were differentially expressed as external pH elevated and most of differential proteins with increased or bell‐shape mode of pH‐dependence were involved in bioenergetic process and metabolism of carbohydrates, fatty acid, amino acids, and nucleotides. Our data thus provide a functional profile of the pH‐responsive proteins in alkaliphiles, leading to elucidation of alkaliphilic‐adaptive mechanism.     

Characterization of the overall and local dynamics of a protein with intermediate rotational anisotropy: Differentiating between conformational exchange and anisotropic diffusion in the B3 domain of protein G

Because the overall tumbling provides a major contribution to protein spectral densities measured in solution, the choice of a proper model for this motion is critical for accurate analysis of protein dynamics. Here we study the overall and backbone dynamics of the B3 domain of protein G using ¹⁵N relaxation measurements and show that the picture of local motions is markedly dependent on the model of overall tumbling. The main difference is in the interpretation of the elevated R ₂ values in the -helix: the isotropic model results in conformational exchange throughout the entire helix, whereas no exchange is predicted by anisotropic models that place the longitudinal axis of diffusion tensor almost parallel to the helix axis. Due to small size (fast tumbling) of the protein, the T ₁ values have low sensitivity to NH bond orientation. The diffusion tensor derived from orientation dependence of R ₂/R ₁ is anisotropic (D _par/D _perp=1.4), with a small rhombic component. In order to distinguish the correct picture of motion, we apply model-independent methods that are sensitive to conformational exchange and do not require knowledge of protein structure or assumptions about its dynamics. A comparison of the CSA/dipolar cross-correlation rate constants with ¹⁵N relaxation rates and the estimation of R _ex terms from relaxation data at 9.4 and 14.1 T indicate no conformational exchange in the helix, in support of the anisotropic models. The experimentally derived diffusion tensor is in excellent agreement with theoretical predictions from hydrodynamic calculations; a detailed comparison with various hydrodynamic models revealed optimal parameters for hydrodynamic calculations.     

Water scaffolding in collagen: Implications on protein dynamics as revealed by solid‐state NMR

Solid‐state NMR studies of collagen samples of various origins confirm that the amplitude of collagen backbone and sidechain motions increases significantly on increasing the water content. This conclusion is supported by the changes observed in three different NMR observables: (i) the linewidth dependence on the ¹H decoupling frequency; (ii) ¹³C CSA changes for the peptide carbonyl groups, and (iii) dephasing rates of ¹H‐¹³C dipolar couplings. In particular, a nearly threefold increase in motional amplitudes of the backbone librations about C‐C^α or N‐C^α bonds was found on increasing the added water content up to 47 wt%D₂O. On the basis of the frequencies of NMR observables involved, the timescale of the protein motions dependent on the added water content is estimated to be of the order of microseconds. This estimate agrees with that from wideline T₂ ¹H NMR measurements. Also, our wideline ¹H NMR measurements revealed that the timescale of the microsecond motions in proteins reduces significantly on increasing the added water content, i.e., an ～15‐fold increase in protein motional frequencies is observed on increasing the added water content to 45 wt% D₂O. The observed changes in collagen dynamics is attributed to the increase in water translational diffusion on increasing the amount of added water, which leads to more frequent “bound water/free water” exchange on the protein surface, accompanied by the breakage and formation of new hydrogen bonds with polar functionalities of protein. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 246–256, 2014.     

Vagal nerve stimulation improves mitochondrial dynamics via an M3 receptor/CaMKKβ/AMPK pathway in isoproterenol‐induced myocardial ischaemia

Mitochondrial dynamics—fission and fusion—are associated with ischaemic heart disease (IHD). This study explored the protective effect of vagal nerve stimulation (VNS) against isoproterenol (ISO)‐induced myocardial ischaemia in a rat model and tested whether VNS plays a role in preventing disorders of mitochondrial dynamics and function. Isoproterenol not only caused cardiac injury but also increased the expression of mitochondrial fission proteins [dynamin‐related peptide1 (Drp1) and mitochondrial fission protein1 (Fis‐1)) and decreased the expression of fusion proteins (optic atrophy‐1 (OPA1) and mitofusins1/2 (Mfn1/2)], thereby disrupting mitochondrial dynamics and leading to increase in mitochondrial fragments. Interestingly, VNS restored mitochondrial dynamics through regulation of Drp1, Fis‐1, OPA1 and Mfn1/2; enhanced ATP content and mitochondrial membrane potential; reduced mitochondrial permeability transition pore (MPTP) opening; and improved mitochondrial ultrastructure and size. Furthermore, VNS reduced the size of the myocardial infarction and ameliorated cardiomyocyte apoptosis and cardiac dysfunction induced by ISO. Moreover, VNS activated AMP‐activated protein kinase (AMPK), which was accompanied by phosphorylation of Ca²⁺/calmodulin‐dependent protein kinase kinase β (CaMKKβ) during myocardial ischaemia. Treatment with subtype‐3 of muscarinic acetylcholine receptor (M₃R) antagonist 4‐diphenylacetoxy‐N‐methylpiperidine methiodide or AMPK inhibitor Compound C abolished the protective effects of VNS on mitochondrial dynamics and function, suggesting that M₃R/CaMKKβ/AMPK signalling are involved in mediating beneficial effects of VNS. This study demonstrates that VNS modulates mitochondrial dynamics and improves mitochondrial function, possibly through the M₃R/CaMKKβ/AMPK pathway, to attenuate ISO‐induced cardiac damage in rats. Targeting mitochondrial dynamics may provide a novel therapeutic strategy in IHD.     

Artificial proteins as allosteric modulators of PDZ3 and SH3 in two‐domain constructs: A computational characterization of novel chimeric proteins

Artificial multidomain proteins with enhanced structural and functional properties can be utilized in a broad spectrum of applications. The design of chimeric fusion proteins utilizing protein domains or one‐domain miniproteins as building blocks is an important advancement for the creation of new biomolecules for biotechnology and medical applications. However, computational studies to describe in detail the dynamics and geometry properties of two‐domain constructs made from structurally and functionally different proteins are lacking. Here, we tested an in silico design strategy using all‐atom explicit solvent molecular dynamics simulations. The well‐characterized PDZ3 and SH3 domains of human zonula occludens (ZO‐1) (3TSZ), along with 5 artificial domains and 2 types of molecular linkers, were selected to construct chimeric two‐domain molecules. The influence of the artificial domains on the structure and dynamics of the PDZ3 and SH3 domains was determined using a range of analyses. We conclude that the artificial domains can function as allosteric modulators of the PDZ3 and SH3 domains. Proteins 2016; 84:1358–1374. © 2016 Wiley Periodicals, Inc.     

FF12MC: A revised AMBER forcefield and new protein simulation protocol

Specialized to simulate proteins in molecular dynamics (MD) simulations with explicit solvation, FF12MC is a combination of a new protein simulation protocol employing uniformly reduced atomic masses by tenfold and a revised AMBER forcefield FF99 with (i) shortened C? H bonds, (ii) removal of torsions involving a nonperipheral sp³ atom, and (iii) reduced 1–4 interaction scaling factors of torsions ? and ψ. This article reports that in multiple, distinct, independent, unrestricted, unbiased, isobaric–isothermal, and classical MD simulations FF12MC can (i) simulate the experimentally observed flipping between left‐ and right‐handed configurations for C14–C38 of BPTI in solution, (ii) autonomously fold chignolin, CLN025, and Trp‐cage with folding times that agree with the experimental values, (iii) simulate subsequent unfolding and refolding of these miniproteins, and (iv) achieve a robust Z score of 1.33 for refining protein models TMR01, TMR04, and TMR07. By comparison, the latest general‐purpose AMBER forcefield FF14SB locks the C14–C38 bond to the right‐handed configuration in solution under the same protein simulation conditions. Statistical survival analysis shows that FF12MC folds chignolin and CLN025 in isobaric–isothermal MD simulations 2–4 times faster than FF14SB under the same protein simulation conditions. These results suggest that FF12MC may be used for protein simulations to study kinetics and thermodynamics of miniprotein folding as well as protein structure and dynamics. Proteins 2016; 84:1490–1516. © 2016 The Authors Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Kinetic and mechanistic differences in the interactions between caldendrin and calmodulin with AKAP79 suggest different roles in synaptic function

The kinetic and mechanistic details of the interaction between caldendrin, calmodulin and the B‐domain of AKAP79 were determined using a biosensor‐based approach. Caldendrin was found to compete with calmodulin for binding at AKAP79, indicating overlapping binding sites. Although the AKAP79 affinities were similar for caldendrin (K_D = 20 n m ) and calmodulin (K_D = 30 n m ), their interaction characteristics were different. The calmodulin interaction was well described by a reversible one‐step model, but was only detected in the presence of Ca²⁺. Caldendrin interacted with a higher level of complexity, deduced to be an induced fit mechanism with a slow relaxation back to the initial encounter complex. It interacted with AKAP79 also in the absence of Ca²⁺, but with different kinetic rate constants. The data are consistent with a similar initial Ca²⁺‐dependent binding step for the two proteins. For caldendrin, a second Ca²⁺‐independent rearrangement step follows, resulting in a stable complex. The study shows the importance of establishing the mechanism and kinetics of protein–protein interactions and that minor differences in the interaction of two homologous proteins can have major implications in their functional characteristics. These results are important for the further elucidation of the roles of caldendrin and calmodulin in synaptic function. Copyright © 2012 John Wiley & Sons, Ltd.     

Photosynthesis and calcification in the calcifying algae Halimeda discoidea studied with microsensors

With microsensors, we measured the steady‐state microprofiles of O₂, pH and Ca²⁺ on the topside of young segments of Halimeda discoidea, as well as the surface dynamics upon light–dark shifts. The effect of several inhibitors was studied. The steady‐state measurements showed that under high light intensity, calcium and protons were taken up, while O₂ was produced. In the dark, O₂ was consumed, the pH decreased to below seawater level and Ca²⁺ uptake was reduced to 50%. At low light intensity (12 mmol photons m^‐2 s^‐1), Ca²⁺ efflux was observed. Upon light–dark shifts, a complicated pattern of both the pH and calcium surface dynamics was observed. Illumination caused an initial pH decrease, followed by a gradual pH increase: this indicated that the surface pH of H. discoidea is determined by more than one light‐induced process. When photosynthesis was inhibited by dichlorophenyl dimethyl urea (DCMU), a strong acidification was observed upon illumination. The nature and physiological function of this putative pump is not known. The calcium dynamics followed all pH dynamics closely, both in the presence and absence of DCMU. The Ca‐channel blockers verapamil and nifedipine had no effect on the Ca²⁺ dynamics and steady‐state profiles. Thus, in H. discoidea, calcification is not regulated by the alga, but is a consequence of pH increase during photosynthesis. Acetazolamide had no effect on photosynthesis, whereas ethoxyzolamide inhibited photosynthesis at higher light intensities. Therefore, all carbonic anhydrase activity is intracellular. Carbonic anhydrase is required to alleviate the CO₂ limitation. Calcification cannot supply sufficient protons and CO₂ to sustain photosynthesis.     

Dynamic Profiling of Protein Mole Synthesis Rates during C2C12 Myoblast Differentiation

Mole (MSR) and fractional (FSR) synthesis rates of proteins during C2C12 myoblast differentiation are investigated. Myoblast cultures supplemented with D₂O during 0–24 h or 72–96 h of differentiation are analyzed by LC‐MS/MS to calculate protein FSR and MSR after samples are spiked with yeast alcohol dehydrogenase (ADH1). Profiling of 153 proteins detected 70 significant (p ≤ 0.05, FDR ≤ 1%) differences in abundance between cell states. Early differentiation is enriched by clusters of ribosomal and heat shock proteins, whereas later differentiation is associated with actin filament binding. The median (first–third quartile) FSR (%/h) during early differentiation 4.1 (2.7–5.3) is approximately twofold greater than later differentiation 1.7 (1.0–2.2), equating to MSR of 0.64 (0.38–1.2) and 0.28 (0.1–0.5) fmol h^?1 µg^?1 total protein, respectively. MSR corresponds more closely with abundance data and highlights proteins associated with glycolytic processes and intermediate filament protein binding that are not evident among FSR data. Similarly, MSR during early differentiation accounts for 78% of the variation in protein abundance during later differentiation, whereas FSR accounts for 4%. Conclusively, the interpretation of protein synthesis data differs when reported in mole or fractional terms, which has consequences when studying the allocation of cellular resources.     

Protein solution structure calculations in solution: Solvated molecular dynamics refinement of calbindin D9k

The three-dimensional solution structures of proteins determinedwith NMR-derived constraints are almost always calculated in vacuo. Thesolution structure of (Ca²⁺)_{_2}-calbindinD_9k has been redetermined by new restrained molecular dynamics(MD) calculations that include Ca²⁺ ions and explicit solventmolecules. Four parallel sets of MD refinements were run to provide accuratecomparisons of structures produced in vacuo, in vacuo withCa²⁺ ions, and with two different protocols in a solvent bathwith Ca²⁺ ions. The structural ensembles were analyzed interms of structural definition, molecular energies, packing density,solvent-accessible surface, hydrogen bonds, and the coordination of calciumions in the two binding loops. Refinement including Ca²⁺ ionsand explicit solvent results in significant improvements in the precisionand accuracy of the structure, particularly in the binding loops. Theseresults are consistent with results previously obtained in free MDsimulations of proteins in solution and show that the rMD refinedNMR-derived solution structures of proteins, especially metalloproteins, canbe significantly improved by these strategies.     

Quantitative analysis of protein interaction network dynamics in yeast

Many cellular functions are mediated by protein–protein interaction networks, which are environment dependent. However, systematic measurement of interactions in diverse environments is required to better understand the relative importance of different mechanisms underlying network dynamics. To investigate environment‐dependent protein complex dynamics, we used a DNA‐barcode‐based multiplexed protein interaction assay in Saccharomyces cerevisiae to measure in vivo abundance of 1,379 binary protein complexes under 14 environments. Many binary complexes (55%) were environment dependent, especially those involving transmembrane transporters. We observed many concerted changes around highly connected proteins, and overall network dynamics suggested that “concerted” protein‐centered changes are prevalent. Under a diauxic shift in carbon source from glucose to ethanol, a mass‐action‐based model using relative mRNA levels explained an estimated 47% of the observed variance in binary complex abundance and predicted the direction of concerted binary complex changes with 88% accuracy. Thus, we provide a resource of yeast protein interaction measurements across diverse environments and illustrate the value of this resource in revealing mechanisms of network dynamics.