Predicting fibrinogen adsorption to polymeric surfaces in silico: a combined method approach

Protein Adsorption on Solid Surfaces: Data Mining, Database, Molecular Surface-Derived Properties, and Semiempirical Relationships

Protein adsorption on solid surfaces is a process relevant to biological, medical, industrial, and environmental applications. Despite this wide interest and advancement in measurement techniques, the complexity of protein adsorption has frustrated its accurate prediction. To address this challenge, here, data regarding protein adsorption reported in the last four decades was collected, checked for completeness and correctness, organized, and archived in an upgraded, freely accessible Biomolecular Adsorption Database, which is equivalent to a large-scale, ad hoc, crowd-sourced multifactorial experiment. The shape and physicochemical properties of the proteins present in the database were quantified on their molecular surfaces using an in-house program (ProMS) operating as an add-on to the PyMol software. Machine learning-based analysis indicated that protein adsorption on hydrophobic and hydrophilic surfaces is modulated by different sets of operational, structural, and molecular surface-based physicochemical parameters. Separately, the adsorption data regarding four "benchmark" proteins, i.e., lysozyme, albumin, IgG, and fibrinogen, was processed by piecewise linear regression with the protein monolayer acting as breakpoint, using the linearization of the Langmuir isotherm formalism, resulting in semiempirical relationships predicting protein adsorption. These relationships, derived separately for hydrophilic and hydrophobic surfaces, described well the protein concentration on the surface as a function of the protein concentration in solution, adsorbing surface contact angle, ionic strength, pH, and temperature of the carrying fluid, and the difference between pH and the isoelectric point of the protein. When applying the semiempirical relationships derived for benchmark proteins to two other "test" proteins with known PDB structure, i.e., β-lactoglobulin and α-lactalbumin, the errors of this extrapolation were found to be in a linear relationship with the dissimilarity between the benchmark and the test proteins. The work presented here can be used for the estimation of operational parameters modulating protein adsorption for various applications such as diagnostic devices, pharmaceuticals, biomaterials, or the food industry.

High-Throughput Routes to Biomaterials Discovery

Many existing clinical treatments are limited in their ability to completely restore decreased or lost tissue and organ function, an unenviable situation only further exacerbated by a globally aging population. As a result, the demand for new medical interventions has increased substantially over the past 20 years, with the burgeoning fields of gene therapy, tissue engineering, and regenerative medicine showing promise to offer solutions for full repair or replacement of damaged or aging tissues. Success in these fields, however, inherently relies on biomaterials that are engendered with the ability to provide the necessary biological cues mimicking native extracellular matrixes that support cell fate. Accelerating the development of such "directive" biomaterials requires a shift in current design practices toward those that enable rapid synthesis and characterization of polymeric materials and the coupling of these processes with techniques that enable similarly rapid quantification and optimization of the interactions between these new material systems and target cells and tissues. This manuscript reviews recent advances in combinatorial and high-throughput (HT) technologies applied to polymeric biomaterial synthesis, fabrication, and chemical, physical, and biological screening with targeted end-point applications in the fields of gene therapy, tissue engineering, and regenerative medicine. Limitations of, and future opportunities for, the further application of these research tools and methodologies are also discussed.

Automation and data-driven design of polymer therapeutics

Polymers are uniquely suited for drug delivery and biomaterial applications due to tunable structural parameters such as length, composition, architecture, and valency. To facilitate designs, researchers may explore combinatorial libraries in a high throughput fashion to correlate structure to function. However, traditional polymerization reactions including controlled living radical polymerization (CLRP) and ring-opening polymerization (ROP) require inert reaction conditions and extensive expertise to implement. With the advent of air-tolerance and automation, several polymerization techniques are now compatible with well plates and can be carried out at the benchtop, making high throughput synthesis and high throughput screening (HTS) possible. To avoid HTS pitfalls often described as "fishing expeditions," it is crucial to employ intelligent and big data approaches to maximize experimental efficiency. This is where the disruptive technologies of machine learning (ML) and artificial intelligence (AI) will likely play a role. In fact, ML and AI are already impacting small molecule drug discovery and showing signs of emerging in drug delivery. In this review, we present state-of-the-art research in drug delivery, gene delivery, antimicrobial polymers, and bioactive polymers alongside data-driven developments in drug design and organic synthesis. From this insight, important lessons are revealed for the polymer therapeutics community including the value of a closed loop design-build-test-learn workflow. This is an exciting time as researchers will gain the ability to fully explore the polymer structural landscape and establish quantitative structure-property relationships (QSPRs) with biological significance.

Consensus QSPR modelling for the prediction of cellular response and fibrinogen adsorption to the surface of polymeric biomaterials

In the current study, we have developed predictive quantitative structure-activity relationship (QSAR) models for cellular response (foetal rate lung fibroblast proliferation) and protein adsorption (fibrinogen adsorption (FA)) on the surface of tyrosine-derived biodegradable polymers designed for tissue engineering purpose using a dataset of 66 and 40 biodegradable polymers, respectively, employing two-dimensional molecular descriptors. Best four individual models have been selected for each of the endpoints. These models are developed using partial least squares regression with a unique combination of six and four descriptors for cellular response and protein adsorption, respectively. The generated models were strictly validated using internal and external metrics to determine the predictive ability and robustness of proposed models. Subsequently, the validated individual models for each response endpoints were used for the generation of 'intelligent' consensus models ( http://teqip.jdvu.ac.in/QSAR_Tools/DTCLab/ ) to improve the quality of predictions for the external data set. These models may help in prediction of virtual polymer libraries for rational design/optimization for properties relevant to biomedical applications prior to their synthesis.

Quantitative design rules for protein-resistant surface coatings using machine learning

Preventing biological contamination (biofouling) is key to successful development of novel surface and nanoparticle-based technologies in the manufacturing industry and biomedicine. Protein adsorption is a crucial mediator of the interactions at the bio - nano -materials interface but is not well understood. Although general, empirical rules have been developed to guide the design of protein-resistant surface coatings, they are still largely qualitative. Herein we demonstrate that this knowledge gap can be addressed by using machine learning approaches to extract quantitative relationships between the material surface chemistry and the protein adsorption characteristics. We illustrate how robust linear and non-linear models can be constructed to accurately predict the percentage of protein adsorbed onto these surfaces using lysozyme or fibrinogen as prototype common contaminants. Our computational models could recapitulate the adsorption of proteins on functionalised surfaces in a test set with an r2 of 0.82 and standard error of prediction of 13%. Using the same data set that enabled the development of the Whitesides rules, we discovered an extension to the original rules. We describe a workflow that can be applied to large, consistently obtained data sets covering a broad range of surface functional groups and protein types.

Feasibility of attenuated total reflection-fourier transform infrared (ATR-FTIR) chemical imaging and partial least squares regression (PLSR) to predict protein adhesion on polymeric surfaces

PLSR with ATR-FTIR chemical imaging predicts protein adhesion on polymeric surfaces well (R² = 0.99, RMSECV = 0.16).

Developing a Suitable Model for Water Uptake for Biodegradable Polymers Using Small Training Sets

Prediction of the dynamic properties of water uptake across polymer libraries can accelerate polymer selection for a specific application. We first built semiempirical models using Artificial Neural Networks and all water uptake data, as individual input. These models give very good correlations (R² > 0.78 for test set) but very low accuracy on cross-validation sets (less than 19% of experimental points within experimental error). Instead, using consolidated parameters like equilibrium water uptake a good model is obtained (R² = 0.78 for test set), with accurate predictions for 50% of tested polymers. The semiempirical model was applied to the 56-polymer library of L-tyrosine-derived polyarylates, identifying groups of polymers that are likely to satisfy design criteria for water uptake. This research demonstrates that a surrogate modeling effort can reduce the number of polymers that must be synthesized and characterized to identify an appropriate polymer that meets certain performance criteria.

Stepping into the omics era: Opportunities and challenges for biomaterials science and engineering

The research paradigm in biomaterials science and engineering is evolving from using low-throughput and iterative experimental designs towards high-throughput experimental designs for materials optimization and the evaluation of materials properties. Computational science plays an important role in this transition. With the emergence of the omics approach in the biomaterials field, referred to as materiomics, high-throughput approaches hold the promise of tackling the complexity of materials and understanding correlations between material properties and their effects on complex biological systems. The intrinsic complexity of biological systems is an important factor that is often oversimplified when characterizing biological responses to materials and establishing property-activity relationships. Indeed, in vitro tests designed to predict in vivo performance of a given biomaterial are largely lacking as we are not able to capture the biological complexity of whole tissues in an in vitro model. In this opinion paper, we explain how we reached our opinion that converging genomics and materiomics into a new field would enable a significant acceleration of the development of new and improved medical devices. The use of computational modeling to correlate high-throughput gene expression profiling with high throughput combinatorial material design strategies would add power to the analysis of biological effects induced by material properties. We believe that this extra layer of complexity on top of high-throughput material experimentation is necessary to tackle the biological complexity and further advance the biomaterials field.

STATEMENT OF SIGNIFICANCE

In this opinion paper, we postulate that converging genomics and materiomics into a new field would enable a significant acceleration of the development of new and improved medical devices. The use of computational modeling to correlate high-throughput gene expression profiling with high throughput combinatorial material design strategies would add power to the analysis of biological effects induced by material properties. We believe that this extra layer of complexity on top of high-throughput material experimentation is necessary to tackle the biological complexity and further advance the biomaterials field.

Analysis of high-throughput screening reveals the effect of surface topographies on cellular morphology

Surface topographies of materials considerably impact cellular behavior as they have been shown to affect cell growth, provide cell guidance, and even induce cell differentiation. Consequently, for successful application in tissue engineering, the contact interface of biomaterials needs to be optimized to induce the required cell behavior. However, a rational design of biomaterial surfaces is severely hampered because knowledge is lacking on the underlying biological mechanisms. Therefore, we previously developed a high-throughput screening device (TopoChip) that measures cell responses to large libraries of parameterized topographical material surfaces. Here, we introduce a computational analysis of high-throughput materiome data to capture the relationship between the surface topographies of materials and cellular morphology. We apply robust statistical techniques to find surface topographies that best promote a certain specified cellular response. By augmenting surface screening with data-driven modeling, we determine which properties of the surface topographies influence the morphological properties of the cells. With this information, we build models that predict the cellular response to surface topographies that have not yet been measured. We analyze cellular morphology on 2176 surfaces, and find that the surface topography significantly affects various cellular properties, including the roundness and size of the nucleus, as well as the perimeter and orientation of the cells. Our learned models capture and accurately predict these relationships and reveal a spectrum of topographies that induce various levels of cellular morphologies. Taken together, this novel approach of high-throughput screening of materials and subsequent analysis opens up possibilities for a rational design of biomaterial surfaces.

A modular versatile chip carrier for high-throughput screening of cell-biomaterial interactions

The field of biomaterials research is witnessing a steady rise in high-throughput screening approaches, comprising arrays of materials of different physico-chemical composition in a chip format. Even though the cell arrays provide many benefits in terms of throughput, they also bring new challenges. One of them is the establishment of robust homogeneous cell seeding techniques and strong control over cell culture, especially for long time periods. To meet these demands, seeding cells with low variation per tester area is required, in addition to robust cell culture parameters. In this study, we describe the development of a modular chip carrier which represents an important step in standardizing cell seeding and cell culture conditions in array formats. Our carrier allows flexible and controlled cell seeding and subsequent cell culture using dynamic perfusion. To demonstrate the application of our device, we successfully cultured and evaluated C2C12 premyoblast cell viability under dynamic conditions for a period of 5 days using an automated pipeline for image acquisition and analysis. In addition, using computational fluid dynamics, lactate and BMP-2 as model molecules, we estimated that there is good exchange of nutrients and metabolites with the flowing medium, whereas no cross-talk between adjacent TestUnits should be expected. Moreover, the shear stresses to the cells can be tailored uniformly over the entire chip area. Based on these findings, we believe our chip carrier may be a versatile tool for high-throughput cell experiments in biomaterials sciences.

Computational modeling of in vitro biological responses on polymethacrylate surfaces

The objective of this research was to examine the capabilities of QSPR (Quantitative Structure Property Relationship) modeling to predict specific biological responses (fibrinogen adsorption, cell attachment and cell proliferation index) on thin films of different polymethacrylates. Using 33 commercially available monomers it is theoretically possible to construct a library of over 40,000 distinct polymer compositions. A subset of these polymers were synthesized and solvent cast surfaces were prepared in 96 well plates for the measurement of fibrinogen adsorption. NIH 3T3 cell attachment and proliferation index were measured on spin coated thin films of these polymers. Based on the experimental results of these polymers, separate models were built for homo-, co-, and terpolymers in the library with good correlation between experiment and predicted values. The ability to predict biological responses by simple QSPR models for large numbers of polymers has important implications in designing biomaterials for specific biological or medical applications.

Logical Analysis of Data in Structure-Activity Investigation of Polymeric Gene Delivery

To date semi-empirical or surrogate modeling has demonstrated great success in the prediction of the biologically relevant properties of polymeric materials. For the first time, a correlation between the chemical structures of poly(β-amino esters) and their efficiency in transfecting DNA was established using the novel technique of logical analysis of data (LAD). Linear combination and explicit representation models were introduced and compared in the framework of the present study. The most successful regression model yielded satisfactory agreement between the predicted and experimentally measured values of transfection efficiency (Pearson correlation coefficient, 0.77; mean absolute error, 3.83). It was shown that detailed analysis of the rules provided by the LAD algorithm offered practical utility to a polymer chemist in the design of new biomaterials.

Neural network analysis identifies scaffold properties necessary for in vitro chondrogenesis in elastin-like polypeptide biopolymer scaffolds

The successful design of biomaterial scaffolds for articular cartilage tissue engineering requires an understanding of the impact of combinations of material formulation parameters on diverse and competing functional outcomes of biomaterial performance. This study sought to explore the use of a type of unsupervised artificial network, a self-organizing map, to identify relationships between scaffold formulation parameters (crosslink density, molecular weight, and concentration) and 11 such outcomes (including mechanical properties, matrix accumulation, metabolite usage and production, and histological appearance) for scaffolds formed from crosslinked elastin-like polypeptide (ELP) hydrogels. The artificial neural network recognized patterns in functional outcomes and provided a set of relationships between ELP formulation parameters and measured outcomes. Mapping resulted in the best mean separation amongst neurons for mechanical properties and pointed to crosslink density as the strongest predictor of most outcomes, followed by ELP concentration. The map also grouped formulations together that simultaneously resulted in the highest values for matrix production, greatest changes in metabolite consumption or production, and highest histological scores, indicating that the network was able to recognize patterns amongst diverse measurement outcomes. These results demonstrated the utility of artificial neural network tools for recognizing relationships in systems with competing parameters, toward the goal of optimizing and accelerating the design of biomaterial scaffolds for articular cartilage tissue engineering.

The BAD project: data mining, database and prediction of protein adsorption on surfaces

Protein adsorption at solid-liquid interfaces is critical to many applications, including biomaterials, protein microarrays and lab-on-a-chip devices. Despite this general interest, and a large amount of research in the last half a century, protein adsorption cannot be predicted with an engineering level, design-orientated accuracy. Here we describe a Biomolecular Adsorption Database (BAD), freely available online, which archives the published protein adsorption data. Piecewise linear regression with breakpoint applied to the data in the BAD suggests that the input variables to protein adsorption, i.e., protein concentration in solution; protein descriptors derived from primary structure (number of residues, global protein hydrophobicity and range of amino acid hydrophobicity, isoelectric point); surface descriptors (contact angle); and fluid environment descriptors (pH, ionic strength), correlate well with the output variable-the protein concentration on the surface. Furthermore, neural network analysis revealed that the size of the BAD makes it sufficiently representative, with a neural network-based predictive error of 5% or less. Interestingly, a consistently better fit is obtained if the BAD is divided in two separate sub-sets representing protein adsorption on hydrophilic and hydrophobic surfaces, respectively. Based on these findings, selected entries from the BAD have been used to construct neural network-based estimation routines, which predict the amount of adsorbed protein, the thickness of the adsorbed layer and the surface tension of the protein-covered surface. While the BAD is of general interest, the prediction of the thickness and the surface tension of the protein-covered layers are of particular relevance to the design of microfluidics devices.

Network models in drug discovery and regenerative medicine

Network motifs and modelling paradigms are attracting increasing attention as modelling tools in drug design and development, and in regenerative medicine. There is a gradual but inexorable convergence between these hitherto disparate disciplines. This review summarizes some very recent work in these areas, leading to an understanding of the complementary roles networks play and factors driving this convergence: network paradigms can be excellent ways of modelling and understanding drug molecules and their action, an understanding of the robustness and vulnerabilities of biological targets may improve the efficacy of drug design and discovery, drug design has an increasingly large role to play in directing stem cell properties, stem cell regulatory networks can be modelled in useful ways using network models at a reasonable level of scale, and the network tools of drug design are also very useful for the design of biomaterials used in regenerative medicine.

A new approach to the rationale discovery of polymeric biomaterials

This paper attempts to illustrate both the need for new approaches to biomaterials discovery as well as the significant promise inherent in the use of combinatorial and computational design strategies. The key observation of this Leading Opinion Paper is that the biomaterials community has been slow to embrace advanced biomaterials discovery tools such as combinatorial methods, high-throughput experimentation, and computational modeling in spite of the significant promise shown by these discovery tools in materials science, medicinal chemistry and the pharmaceutical industry. It seems that the complexity of living cells and their interactions with biomaterials has been a conceptual as well as a practical barrier to the use of advanced discovery tools in biomaterials science. However, with the continued increase in computer power, the goal of predicting the biological response of cells in contact with biomaterials surfaces is within reach. Once combinatorial synthesis, high-throughput experimentation, and computational modeling are integrated into the biomaterials discovery process, a significant acceleration is possible in the pace of development of improved medical implants, tissue regeneration scaffolds, and gene/drug delivery systems.

Glass transition temperature prediction of polymers through the mass-per-flexible-bond principle

A semi-empirical method based on the mass-per-flexible-bond (M/f) principle was used to quantitatively explain the large range of glass transition temperatures (T(g)) observed in a library of 132 L-tyrosine derived homo, co- and terpolymers containing different functional groups. Polymer class specific behavior was observed in T(g) vs. M/f plots, and explained in terms of different densities, steric hindrances and intermolecular interactions of chemically distinct polymers. The method was found to be useful in the prediction of polymer T(g). The predictive accuracy was found to range from 6.4 to 3.7 K, depending on polymer class. This level of accuracy compares favorably with (more complicated) methods used in the literature. The proposed method can also be used for structure prediction of polymers to match a target T(g) value, by keeping the thermal behavior of a terpolymer constant while independently choosing its chemistry. Both applications of the method are likely to have broad applications in polymer and (bio)material science.

Prediction of Fibrinogen Adsorption for Biodegradable Polymers: Integration of Molecular Dynamics and Surrogate Modeling

This work is a part of a series of publications devoted to the development of surrogate (semi-empirical) models for the prediction of fibrinogen adsorption onto polymer surfaces. Since fibrinogen is one of the key proteins involved in platelet activation and the formation of thrombosis, the modeling of fibrinogen adsorption on the surface of blood contacting medical devices is of high theoretical and practical significance. We report here, for the first time, on the incorporation three-dimensional structures of polymers obtained from atomistic simulations into conventional mesoscopic-scale calculations. Low energy conformations derived from Molecular Dynamics simulations for 45 representatives of a combinatorial library of polyarylates were used in an improved modeling procedure (referred to as "3D surrogate model") instead of simplistic two-dimensional representations of polymer structures, which were used in several previous models (collectively referred to as "2D surrogate models"). In the framework of this 3D model we created 12 model sets of polymers to account for their chirality, conformational diversity and the structural influence of a solvent. For each polymer set, three-dimensional molecular descriptors were generated and then ranked with respect to the experimental fibrinogen adsorption data by means of a Monte Carlo Decision Tree. The most significant descriptors identified by Decision Tree and the experimental dataset were utilized to predict fibrinogen adsorption using an Artificial Neural Network (ANN). The best prediction achieved by the 3D surrogate model demonstrated a noticeable improvement in the predictive quality as compared to the previously used 2D model (as evidenced by the increase in the average Pearson correlation coefficient from 0.67+/-0.13 to 0.54+/-0.12). The predictive quality of the 3D surrogate model compares favorably with the best results previously reported for extended 2D model that combines an ANN with Partial Least Squares (PLS) regression and principal component (PC) analysis. The significance of the newly developed 3D model is that it allows high accuracy prediction of fibrinogen adsorption without the need for experimentally derived descriptors and it has better predictive quality than the original 2D surrogate model due to utilization of realistic polymer representations.

Mini-review: combinatorial approaches for the design of novel coating systems

Combinatorial and high throughput experimental methods are being applied to the design and development of novel polymers and coatings used in a number of application areas. Methods have been developed for polymer synthesis and screening and for the development of polymer thin film and coating libraries and the screening of these libraries for key properties such as surface energy and modulus. Combinatorial and high throughput methods enable the efficient exploration of a large number of compositional variables over a wide range. In the development of coatings for use in the marine environment, the key challenge is in the development of screening methods that can predict good performance. A number of assays are under development that will permit the rapid screening of the interaction of coatings with representative marine organisms.