Follicular Fluid Proteomic Analysis to Identify Predictive Markers of Normal Embryonic Development.

Ageing populations, mass "baby-free" policies and children born to mothers at the age at which they are biologically expected to become grandmothers are growing problems in most developed societies. Therefore, any opportunity to improve the quality of infertility treatments seems important for the survival of societies. The possibility of indirectly studying the quality of developing oocytes by examining their follicular fluids (hFFs) offers new opportunities for progress in our understanding the processes of final oocyte maturation and, consequently, for predicting the quality of the resulting embryos and personalising their culture. Using mass spectrometry, we studied follicular fluids collected individually during in vitro fertilisation and compared their composition with the quality of the resulting embryos. We analysed 110 follicular fluids from 50 oocyte donors, from which we obtained 44 high-quality, 39 medium-quality, and 27 low-quality embryos. We identified 2182 proteins by Sequential Window Acquisition of all Theoretical Mass Spectra (SWATH-MS) using a TripleTOF 5600+ hybrid mass spectrometer, of which 484 were suitable for quantification. We were able to identify several proteins whose concentrations varied between the follicular fluids of different oocytes from the same patient and between patients. Among them, the most important appear to be immunoglobulin heavy constant alpha 1 (IgA1hc) and dickkopf-related protein 3. The first one is found at higher concentrations in hFFs from which oocytes develop into poor-quality embryos, the other one exhibits the opposite pattern. None of these have, so far, had any specific links to fertility disorders. In light of these findings, these proteins should be considered a primary target for research aimed at developing a diagnostic tool for oocyte quality control and pre-fertilisation screening. This is particularly important in cases where the fertilisation of each egg is not an option for ethical or other reasons, or in countries where it is prohibited by law.

Comparison of Peptidomes Extracted from Healthy Tissue and Tumor Tissue of the Parotid Glands and Saliva Samples.

Salivary gland tumors are highly variable in clinical presentation and histology. The World Health Organization (WHO) classifies 22 types of malignant and 11 types of benign tumors of the salivary glands. Diagnosis of salivary gland tumors is based on imaging (ultrasound, magnetic resonance imaging) and fine-needle aspiration biopsy, but the final diagnosis is based on histopathological examination of the removed tumor tissue. In this pilot study, we are testing a new approach to identifying peptide biomarkers in saliva that can be used to diagnose salivary gland tumors. The research material for the peptidomic studies was extracts from washings of neoplastic tissues and healthy tissues (control samples). At the same time, saliva samples from patients and healthy individuals were analyzed. The comparison of the peptidome composition of tissue extracts and saliva samples may allow the identification of potential peptide markers of salivary gland tumors in patients' saliva. The peptidome compositions extracted from 18 tumor and 18 healthy tissue samples, patients' saliva samples (11 samples), and healthy saliva samples (8 samples) were analyzed by LC-MS tandem mass spectrometry. A group of 109 peptides was identified that were present only in the tumor tissue extracts and in the patients' saliva samples. Some of the identified peptides were derived from proteins previously suggested as potential biomarkers of salivary gland tumors (ANXA1, BPIFA2, FGB, GAPDH, HSPB1, IGHG1, VIM) or tumors of other tissues or organs (SERPINA1, APOA2, CSTB, GSTP1, S100A8, S100A9, TPI1). Unfortunately, none of the identified peptides were present in all samples analyzed. This may be due to the high heterogeneity of this type of cancer. The surprising result was that extracts from tumor tissue did not contain peptides derived from salivary gland-specific proteins (STATH, SMR3B, HTN1, HTN3). These results could suggest that the developing tumor suppresses the production of proteins that are essential components of saliva.

Compatibility of Distinct Label-Free Proteomic Workflows in Absolute Quantification of Proteins Linked to the Oocyte Quality in Human Follicular Fluid.

We present two separate label-free quantitative workflows based on different high-resolution mass spectrometers and LC setups, which are termed after the utilized instrument: Quad-Orbitrap (nano-LC) and Triple Quad-TOF (micro-LC) and their directed adaptation toward the analysis of human follicular fluid proteome. We identified about 1000 proteins in each distinct workflow using various sample preparation methods. With assistance of the Total Protein Approach, we were able to obtain absolute protein concentrations for each workflow. In a pilot study of twenty samples linked to diverse oocyte quality status from four donors, 455 and 215 proteins were quantified by the Quad-Orbitrap and Triple Quad-TOF workflows, respectively. The concentration values obtained from both workflows correlated to a significant degree. We found reasonable agreement of both workflows in protein fold changes between tested groups, resulting in unified lists of 20 and 22 proteins linked to oocyte maturity and blastocyst development, respectively. The Quad-Orbitrap workflow was best suited for an in-depth analysis without the need of extensive fractionation, especially of low abundant proteome, whereas the Triple Quad-TOF workflow allowed a more robust approach with a greater potential to increase in effectiveness with the growing number of analyzed samples after the initial effort of building a comprehensive spectral library.

AfsK-Mediated Site-Specific Phosphorylation Regulates DnaA Initiator Protein Activity in Streptomyces coelicolor.

In all organisms, chromosome replication is regulated mainly at the initiation step. Most of the knowledge about the mechanisms that regulate replication initiation in bacteria has come from studies on rod-shaped bacteria, such as Escherichia coli and Bacillus subtilisStreptomyces is a bacterial genus that is characterized by distinctive features and a complex life cycle that shares some properties with the developmental cycle of filamentous fungi. The unusual lifestyle of streptomycetes suggests that these bacteria use various mechanisms to control key cellular processes. Here, we provide the first insights into the phosphorylation of the bacterial replication initiator protein, DnaA, from Streptomyces coelicolor We suggest that phosphorylation of DnaA triggers a conformational change that increases its ATPase activity and decreases its affinity for the replication origin, thereby blocking the formation of a functional orisome. We suggest that the phosphorylation of DnaA is catalyzed by Ser/Thr kinase AfsK, which was shown to regulate the polar growth of S. coelicolor Together, our results reveal that phosphorylation of the DnaA initiator protein functions as a negative regulatory mechanism to control the initiation of chromosome replication in a manner that presumably depends on the cellular localization of the protein.IMPORTANCE This work provides insights into the phosphorylation of the DnaA initiator protein in Streptomyces coelicolor and suggests a novel bacterial regulatory mechanism for initiation of chromosome replication. Although phosphorylation of DnaA has been reported earlier, its biological role was unknown. This work shows that upon phosphorylation, the cooperative binding of the replication origin by DnaA may be disturbed. We found that AfsK kinase is responsible for phosphorylation of DnaA. Upon upregulation of AfsK, chromosome replication occurred further from the hyphal tip. Orthologs of AfsK are exclusively found in mycelial actinomycetes that are related to Streptomyces and exhibit a complex life cycle. We propose that the AfsK-mediated regulatory pathway serves as a nonessential, energy-saving mechanism in S. coelicolor.

A general method for the derivation of the functional forms of the effective energy terms in coarse-grained energy functions of polymers. III. Determination of scale-consistent backbone-local and correlation potentials in the UNRES force field and force-field calibration and validation.

The general theory of the construction of scale-consistent energy terms in the coarse-grained force fields presented in Paper I of this series has been applied to the revision of the UNRES force field for physics-based simulations of proteins. The potentials of mean force corresponding to backbone-local and backbone-correlation energy terms were calculated from the ab initio energy surfaces of terminally blocked glycine, alanine, and proline, and the respective analytical expressions, derived by using the scale-consistent formalism, were fitted to them. The parameters of all these potentials depend on single-residue types, thus reducing their number and preventing over-fitting. The UNRES force field with the revised backbone-local and backbone-correlation terms was calibrated with a set of four small proteins with basic folds: tryptophan cage variant (TRP1; α), Full Sequence Design (FSD; α + ß), villin headpiece (villin; α), and a truncated FBP-28 WW-domain variant (2MWD; ß) (the NEWCT-4P force field) and, subsequently, with an enhanced set of 9 proteins composed of TRP1, FSD, villin, 1BDC (α), 2I18 (α), 1QHK (α + ß), 2N9L (α + ß), 1E0L (ß), and 2LX7 (ß) (the NEWCT-9P force field). The NEWCT-9P force field performed better than NEWCT-4P in a blind-prediction-like test with a set of 26 proteins not used in calibration and outperformed, in a test with 76 proteins, the most advanced OPT-WTFSA-2 version of UNRES with former backbone-local and backbone-correlation terms that contained more energy terms and more optimizable parameters. The NEWCT-9P force field reproduced the bimodal distribution of backbone-virtual-bond angles in the simulated structures, as observed in experimental protein structures.

Qualitative and Quantitative Analysis of Proteome and Peptidome of Human Follicular Fluid Using Multiple Samples from Single Donor with LC-MS and SWATH Methodology.

Human follicular fluid (hFF) is a natural environment of oocyte maturation, and some components of hFF could be used to judge oocyte capability for fertilization and further development. In our pilot small-scale study three samples from four donors (12 samples in total) were analyzed to determine which hFF proteins/peptides could be used to differentiate individual oocytes and which are patient-specific. Ultrafiltration was used to fractionate hFF to high-molecular-weight (HMW) proteome (>10 kDa) and low-molecular-weight (LMW) peptidome (<10 kDa) fractions. HMW and LMW compositions were analyzed using LC-MS in SWATH data acquisition and processing methodology. In total we were able to identify 158 proteins, from which 59 were never reported before as hFF components. 55 (45 not reported before) proteins were found by analyzing LMW fraction, 67 (14 not reported before) were found by analyzing HMW fraction, and 36 were identified in both fractions of hFF. We were able to perform quantitative analysis for 72 proteins from HMW fraction of hFF. We found that concentrations of 11 proteins varied substantially among hFF samples from single donors, and those proteins are promising targets to identify biomarkers useful in oocyte quality assessment.

Maximum Likelihood Calibration of the UNRES Force Field for Simulation of Protein Structure and Dynamics.

By using the maximum likelihood method for force-field calibration recently developed in our laboratory, which is aimed at achieving the agreement between the simulated conformational ensembles of selected training proteins and the corresponding ensembles determined experimentally at various temperatures, the physics-based coarse-grained UNRES force field for simulations of protein structure and dynamics was optimized with seven small training proteins exhibiting a variety of secondary and tertiary structures. Four runs of optimization, in which the number of optimized force-field parameters was gradually increased, were carried out, and the resulting force fields were subsequently tested with a set of 22 α-, 12 ß-, and 12 α + ß-proteins not used in optimization. The variant in which energy-term weights, local, and correlation potentials, side-chain radii, and anisotropies were optimized turned out to be the most transferable and outperformed all previous versions of UNRES on the test set.

Lessons from application of the UNRES force field to predictions of structures of CASP10 targets.

The performance of the physics-based protocol, whose main component is the United Residue (UNRES) physics-based coarse-grained force field, developed in our laboratory for the prediction of protein structure from amino acid sequence, is illustrated. Candidate models are selected, based on probabilities of the conformational families determined by multiplexed replica-exchange simulations, from the 10th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP10). For target T0663, classified as a new fold, which consists of two domains homologous to those of known proteins, UNRES predicted the correct symmetry of packing, in which the domains are rotated with respect to each other by 180° in the experimental structure. By contrast, models obtained by knowledge-based methods, in which each domain is modeled very accurately but not rotated, resulted in incorrect packing. Two UNRES models of this target were featured by the assessors. Correct domain packing was also predicted by UNRES for the homologous target T0644, which has a similar structure to that of T0663, except that the two domains are not rotated. Predictions for two other targets, T0668 and T0684_D2, are among the best ones by global distance test score. These results suggest that our physics-based method has substantial predictive power. In particular, it has the ability to predict domain-domain orientations, which is a significant advance in the state of the art.

A Maximum-Likelihood Approach to Force-Field Calibration.

A new approach to the calibration of the force fields is proposed, in which the force-field parameters are obtained by maximum-likelihood fitting of the calculated conformational ensembles to the experimental ensembles of training system(s). The maximum-likelihood function is composed of logarithms of the Boltzmann probabilities of the experimental conformations, calculated with the current energy function. Because the theoretical distribution is given in the form of the simulated conformations only, the contributions from all of the simulated conformations, with Gaussian weights in the distances from a given experimental conformation, are added to give the contribution to the target function from this conformation. In contrast to earlier methods for force-field calibration, the approach does not suffer from the arbitrariness of dividing the decoy set into native-like and non-native structures; however, if such a division is made instead of using Gaussian weights, application of the maximum-likelihood method results in the well-known energy-gap maximization. The computational procedure consists of cycles of decoy generation and maximum-likelihood-function optimization, which are iterated until convergence is reached. The method was tested with Gaussian distributions and then applied to the physics-based coarse-grained UNRES force field for proteins. The NMR structures of the tryptophan cage, a small α-helical protein, determined at three temperatures (T = 280, 305, and 313 K) by Halabis et al. ( J. Phys. Chem. B 2012 , 116 , 6898 - 6907 ), were used. Multiplexed replica-exchange molecular dynamics was used to generate the decoys. The iterative procedure exhibited steady convergence. Three variants of optimization were tried: optimization of the energy-term weights alone and use of the experimental ensemble of the folded protein only at T = 280 K (run 1); optimization of the energy-term weights and use of experimental ensembles at all three temperatures (run 2); and optimization of the energy-term weights and the coefficients of the torsional and multibody energy terms and use of experimental ensembles at all three temperatures (run 3). The force fields were subsequently tested with a set of 14 α-helical and two α + ß proteins. Optimization run 1 resulted in better agreement with the experimental ensemble at T = 280 K compared with optimization run 2 and in comparable performance on the test set but poorer agreement of the calculated folding temperature with the experimental folding temperature. Optimization run 3 resulted in the best fit of the calculated ensembles to the experimental ones for the tryptophan cage but in much poorer performance on the training set, suggesting that use of a small α-helical protein for extensive force-field calibration resulted in overfitting of the data for this protein at the expense of transferability. The optimized force field resulting from run 2 was found to fold 13 of the 14 tested α-helical proteins and one small α + ß protein with the correct topologies; the average structures of 10 of them were predicted with accuracies of about 5 Å C(α) root-mean-square deviation or better. Test simulations with an additional set of 12 α-helical proteins demonstrated that this force field performed better on α-helical proteins than the previous parametrizations of UNRES. The proposed approach is applicable to any problem of maximum-likelihood parameter estimation when the contributions to the maximum-likelihood function cannot be evaluated at the experimental points and the dimension of the configurational space is too high to construct histograms of the experimental distributions.

Mean-field interactions between nucleic-acid-base dipoles can drive the formation of a double helix.

A proposed coarse-grained model of nucleic acids demonstrates that average interactions between base dipoles, together with chain connectivity and excluded-volume interactions, are sufficient to form double-helical structures of DNA and RNA molecules. Additionally, local interactions determine helix handedness and direction of strand packing. This result, and earlier research on reduced protein models, suggests that mean-field multipole-multipole interactions are the principal factors responsible for the formation of regular structure of biomolecules.

Substitutions at position 263 within the von Willebrand factor type A domain determine the functionality of complement C2 protein.

The complement system is one of the first defense lines protecting from invading pathogens. However, it may turn offensive to the body's own cells and tissues when deregulated by the presence of rare genetic variants that impair physiological regulation and/or provoke abnormal activity of key enzymatic components. Factor B and complement C2 are examples of paralogs engaged in the alternative and classical/lectin complement pathway, respectively. Pathogenic mutations in the von Willebrand factor A domain (vWA) of FB have been known for years. Despite substantial homology between two proteins and the demonstration that certain substitutions in FB translated to C2 result in analogous phenotype, there was a limited number of reports on pathogenic C2 variants in patients. Recently, we studied a cohort of patients suffering from rare kidney diseases and confirmed the existence of two gain-of-function and three loss-of-function mutations within the C2 gene sequences coding for the vWA domain (amino acids 254-452) or nearly located unstructured region (243-253) of C2 protein. Herein, we report the functional consequences of amino acid substitution of glutamine at position 263. The p.Q263G variant resulted in the gain-of-function phenotype, similarly to a homologous mutation p.D279G in FB. Conversely, the p.Q263P variant found in a patient with C3 glomerulopathy resulted in the loss of C2 function. Our results confirm that the N-terminal part of the vWA domain is a hot spot crucial for the complement C2 function.

Modeling the Structure, Dynamics, and Transformations of Proteins with the UNRES Force Field.

The physics-based united-residue (UNRES) model of proteins ( www.unres.pl ) has been designed to carry out large-scale simulations of protein folding. The force field has been derived and parameterized based on the principles of statistical-mechanics, which makes it independent of structural databases and applicable to treat nonstandard situations such as, proteins that contain D-amino-acid residues. Powered by Langevin dynamics and its replica-exchange extensions, UNRES has found a variety of applications, including ab initio and database-assisted protein-structure prediction, simulating protein-folding pathways, exploring protein free-energy landscapes, and solving biological problems. This chapter provides a summary of UNRES and a guide for potential users regarding the application of the UNRES package in a variety of research tasks.

Computational techniques for efficient conformational sampling of proteins.

In this review, we summarize the computational methods for sampling the conformational space of biomacromolecules. We discuss the methods applicable to find only lowest energy conformations (global minimization of the potential-energy function) and to generate canonical ensembles (canonical Monte Carlo method and canonical molecular dynamics method and their extensions). Special attention is devoted to the use of coarse-grained models that enable simulations to be enhanced by several orders of magnitude.

Gain-of-Function Mutations R249C and S250C in Complement C2 Protein Increase C3 Deposition in the Presence of C-Reactive Protein.

The impairment of the alternative complement pathway contributes to rare kidney diseases such as atypical hemolytic uremic syndrome (aHUS) and C3 glomerulopathy (C3G). We recently described an aHUS patient carrying an exceptional gain-of-function (GoF) mutation (S250C) in the classical complement pathway component C2 leading to the formation of hyperactive classical convertases. We now report the identification of the same mutation and another C2 GoF mutation R249C in two other patients with a glomerulopathy of uncertain etiology. Both mutations stabilize the classical C3 convertases by a similar mechanism. The presence of R249C and S250C variants in serum increases complement-dependent cytotoxicity (CDC) in antibody-sensitized human cells and elevates deposition of C3 on ELISA plates coated with C-reactive protein (CRP), as well as on the surface of glomerular endothelial cells. Our data justify the inclusion of classical pathway genes in the genetic analysis of patients suspected of complement-driven renal disorders. Also, we point out CRP as a potential antibody-independent trigger capable of driving excessive complement activation in carriers of the GoF mutations in complement C2.

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. III. Dynamics of long-range hydrophobic interactions.

A 20-residue peptide, IG(42-61), derived from the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptoccocus was studied using circular dichroism, nuclear magnetic resonance (NMR) spectroscopy at various temperatures and by differential scanning calorimetry (DSC). Unlike other related peptides studied so far, this peptide displays two heat capacity peaks in DSC measurements (at a scanning rate of 1.5 deg/min at a peptide concentration of 0.07 mM), which suggests a three-state folding/unfolding process. The results from DSC and NMR measurements suggest the formation of a dynamic network of hydrophobic interactions stabilizing the structure, which resembles a beta-hairpin shape over a wide range of temperatures (283-313 K). Our results show that IG (42-61) possesses a well-organized three-dimensional structure stabilized by long-range hydrophobic interactions (Tyr50 ... Phe57 and Trp48 ... Val59) at T = 283 K and (Trp48 ... Val59) at 305 and 313 K. The mechanism of beta-hairpin folding and unfolding, as well as the influence of peptide length on its conformational properties, are also discussed.

Evidence, from simulations, of a single state with residual native structure at the thermal denaturation midpoint of a small globular protein.

The folding of the B-domain of staphylococcal protein A has been studied by coarse-grained canonical and multiplexed replica-exchange molecular dynamics simulations with the UNRES force field in a broad range of temperatures (270 K < or = T < or = 350 K). In canonical simulations, the folding was found to occur either directly to the native state or through kinetic traps, mainly the topological mirror image of the native three-helix bundle. The latter folding scenario was observed more frequently at low temperatures. With increase of temperature, the frequency of the transitions between the folded and misfolded/unfolded states increased and the folded state became more diffuse with conformations exhibiting increased root-mean-square deviations from the experimental structure (from about 4 A at T = 300 K to 8.7 A at T = 325 K). An analysis of the equilibrium conformational ensemble determined from multiplexed replica exchange simulations at the folding-transition temperature (T(f) = 325 K) showed that the conformational ensemble at this temperature is a collection of conformations with residual secondary structures, which possess native or near-native clusters of nonpolar residues in place, and not a 50-50% mixture of fully folded and fully unfolded conformations. These findings contradict the quasi-chemical picture of two- or multistate protein folding, which assumes an equilibrium between the folded, unfolded, and intermediate states, with equilibrium shifting with temperature but with the native conformations remaining essentially unchanged. Our results also suggest that long-range hydrophobic contacts are the essential factor to keep the structure of a protein thermally stable.

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. I. Importance of hydrophobic interactions in stabilization of beta-hairpin structure.

We previously studied a 16-amino acid-residue fragment of the C-terminal beta-hairpin of the B3 domain (residues 46-61), [IG(46-61)] of the immunoglobulin binding protein G from Streptoccocus, and found that hydrophobic interactions and the turn region play an important role in stabilizing the structure. Based on these results, we carried out systematic structural studies of peptides derived from the sequence of IG (46-61) by systematically shortening the peptide by one residue at a time from both the C- and the N-terminus. To determine the structure and stability of two resulting 12- and 14-amino acid-residue peptides, IG(48-59) and IG(47-60), respectively, we carried out circular dichroism, NMR, and calorimetric studies of these peptides in pure water. Our results show that IG(48-59) possesses organized three-dimensional structure stabilized by hydrophobic interactions (Tyr50-Phe57 and Trp48-Val59) at T = 283 and 305 K. At T = 313 K, the structure breaks down because of increased chain entropy, but the turn region is preserved in the same position observed for the structure of the whole protein. The breakdown of structure occurs near the melting temperature of this peptide (T(m) = 310 K) measured by differential scanning calorimetry (DSC). The melting temperature of IG(47-60) determined by DSC is T(m) = 330 K and its structure is similar to that of the native beta-hairpin at all (lower) temperatures examined (283-313 K). Both of these truncated sequences are conserved in all known amino acid sequences of the B domains of the immunoglobulin binding protein G from bacteria. Thus, this study contributes to an understanding of the mechanism of folding of this whole family of proteins, and provides information about the mechanism of formation and stabilization of a beta-hairpin structural element.

Mechanism of formation of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus. II. Interplay of local backbone conformational dynamics and long-range hydrophobic interactions in hairpin formation.

Two peptides, corresponding to the turn region of the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G from Streptococcus, consisting of residues 51-56 [IG(51-56)] and 50-57 [IG(50-57)], respectively, were studied by circular dichroism and NMR spectroscopy at various temperatures and by differential scanning calorimetry. Our results show that the part of the sequence corresponding to the beta-turn in the native structure (DDATKT) of the B3 domain forms bent conformations similar to those observed in the native protein. The formation of a turn is observed for both peptides in a broad range of temperatures (T = 283-323 K), which confirms the conclusion drawn from our previous studies of longer sequences from the C-terminal beta-hairpin of the B3 domain of the immunoglobulin binding protein G (16, 14, and 12 residues), that the DDATKT sequence forms a nucleation site for formation of the beta-hairpin structure of peptides corresponding to the C-terminal part of all the B domains of the immunoglobulin binding protein G. We also show and discuss the role of long-range hydrophobic interactions as well as local conformational properties of polypeptide chains in the mechanism of formation of the beta-hairpin structure.