Balancing Cost and Accuracy in Quantum Mechanical Simulations on Collagen Protein Models

Abstract

Collagen proteins are spread in almost every vertebrate's tissue with mechanical function. The defining feature of this fundamental family of proteins is its well-known collagen triple-helical domain. This helical domain can have different geometries, varying in helical elongation and interstrands contact, as a function of the amino acidic composition. The helical geometrical features play an important role in the interaction of the collagen protein with cell receptors, but for the vast majority of collagen compositions, these geometrical features are unknown. Quantum mechanical (QM) simulations based on the density functional theory (DFT) provide a robust approach to characterize the scenario on the collagen composition−structure relationships. In this work, we analyze the role of the adopted computational method in predicting the collagen structure for two purposes. First, we look for a cost-effective computational approach to apply to a large-scale composition− structure analysis. Second, we attempt to assess the robustness of the predictions by varying the QM methods. Therefore, we have run geometry optimization on periodic models of the collagen protein using a variety of approaches based on the most commonly used DFT functionals (PBE, HSE06, and B3LYP) with and without dispersion correction (D3 ABC ). We have coupled these methods with several different basis sets, looking for the highest accuracy/cost ratio. Furthermore, we have studied the performance of the composite HF-3c method and the semiempirical GFN1-xTB method. Our results identify a computational recipe that is potentially capable of predicting collagen structural features in line with DFT simulations, with orders of magnitude reduced computational cost, encouraging further investigations on the topic.

Nanoscale MXene Interlayer and Substrate Adhesion for Lubrication: A Density Functional Theory Study

“…However, most of the available dispersion correction methods for DFT have been developed for organic molecules. Although they have been holistically tested for periodic organic systems, such as polymers − and molecular crystals, they should be carefully applied to inorganic solid-state materials . More advanced, parameter-free methodologies such as the Møller–Plesset perturbation theory (MP2) and the random phase approximation (RPA) can capture the elusive dispersion forces in an accurate way in solid inorganic systems.…”

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree–Fock Methods Corrected with Atom-Centered Potentials

“…The biomolecule–biomolecule subset contains model systems representative of nucleotide–nucleotide interactions as well as protein fragments interacting with carbohydrates, nucleotides, drugs, water, and with other proteins. Such interactions are relevant in applications like protein folding, , protein structure refinement, , protein−ligand binding, − intercalation, nucleobase stacking, and protein hydration, to name a few. Uncorrected minimal or double-ζ basis set HF-D3 in general overestimate the interaction energies in this subset, and application of the ACPs reduces this overestimation and decreases the error spread and SDs.…”

Modeling phosphorene and $$\hbox {MoS}_{2}$$ interacting with iron: lubricating effects compared to graphene

“…We took into account London dispersion forces by adopting the D2 scheme [29], using the value for the s 6 scaling factor suggested in the original paper ( s 6 = 0.75). This type of dispersion scheme gave very good results for computing not only organic [30][31][32][33][34] but also inorganic materials properties [35,36]. Recent investigations indicate that more modern dispersion correction schemes are also appropriate for simulating the chemisorption of 2D material on metal and ceramic surfaces [37,38].…”