Next Article in Journal
Mechanical Properties of Hybrid Carbonized Plant Fibers Reinforced Bio-Based Epoxy Laminates
Next Article in Special Issue
Flows of Linear Polymer Solutions and Other Suspensions of Rod-like Particles: Anisotropic Micropolar-Fluid Theory Approach
Previous Article in Journal
Production and Surface Modification of Cellulose Bioproducts
Previous Article in Special Issue
Taguchi Optimization of Roundness and Concentricity of a Plastic Injection Molded Barrel of a Telecentric Lens
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Polymers
Volume 13
Issue 19
10.3390/polym13193434
polymers-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvent Effect on the Structure and Properties of RGD Peptide (1FUV) at Body Temperature (310 K) Using Ab Initio Molecular Dynamics

by
Khagendra Baral
1,
Puja Adhikari
1,
Bahaa Jawad
1,
Rudolf Podgornik
2,3,4,5 and
Wai-Yim Ching
1,*
1
Department of Physics and Astronomy, University of Missouri-Kansas City, Kansas City, MO 64110, USA
2
School of Physical Sciences, Kavli Institute of Theoretical Science, University of Chinese Academy of Sciences, Beijing 100049, China
3
CAS Key Laboratory of Soft Matter Physics, Institute of Physics, Chinese Academy of Sciences, Beijing 100090, China
4
Wenzhou Institute of the University of Chinese Academy of Sciences, Wenzhou 325000, China
5
Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, SI-1000 Ljubljana, Slovenia
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(19), 3434; https://doi.org/10.3390/polym13193434
Submission received: 17 September 2021 / Revised: 27 September 2021 / Accepted: 29 September 2021 / Published: 7 October 2021
(This article belongs to the Special Issue Advanced Polymer Simulation and Processing)
Browse Figures
Review Reports Versions Notes

Abstract

:
The structure and properties of the arginine-glycine-aspartate (RGD) sequence of the 1FUV peptide at 0 K and body temperature (310 K) are systematically investigated in a dry and aqueous environment using more accurate ab initio molecular dynamics and density functional theory calculations. The fundamental properties, such as electronic structure, interatomic bonding, partial charge distribution, and dielectric response function at 0 and 310 K are analyzed, comparing them in dry and solvated models. These accurate microscopic parameters determined from highly reliable quantum mechanical calculations are useful to define the range and strength of complex molecular interactions occurring between the RGD peptide and the integrin receptor. The in-depth bonding picture analyzed using a novel quantum mechanical metric, the total bond order (TBO), quantifies the role played by hydrogen bonds in the internal cohesion of the simulated structures. The TBO at 310 K decreases in the dry model but increases in the solvated model. These differences are small but extremely important in the context of conditions prevalent in the human body and relevant for health issues. Our results provide a new level of understanding of the structure and properties of the 1FUV peptide and help in advancing the study of RGD containing other peptides.

1. Introduction

Biomolecular materials containing the arginine-glycine-aspartate (RGD) sequence are always at the center of biophysics research in their application such as in the bone scaffold, synthesis, and regeneration of tissue and cartilage [1,2], in imaging as radiotracers [3,4,5,6], for cancer therapy [7,8,9], and in targeted drug delivery [10,11,12,13,14,15]. The RGD motif peptide serves as a primary integrin recognition site in extracellular matrix proteins since it has a strong binding affinity to integrins, which are heterodimeric cell surface receptors and mediate cell-extracellular matrix adhesion [16,17,18,19]. Out of 24 known integrins, one third of them bind to the RGD motif as the primary recognition sequence in their ligands, which makes RGD an attractive target in numerous drug delivery systems [20,21,22,23,24]. The increased application of integrins in drug development and their functions in the physiological processes requires a complete understanding of the structure and properties of the RGD peptides at body temperature (310K), crucial in the design of selective inhibitors. Therefore, a detailed study of the structure and properties of 1FUV in RGD peptides at 310 K is of particular interest in delineating the modifications that occur with a rise or fall of body temperature. Moreover, the short- and long-range interactions are crucial for the molecular recognition and self-assembly of biological macromolecules, and consequently determining the microscopic parameters such as partial charge and frequency-dependent dielectric functions based on more accurate quantum mechanical calculation will facilitate the fundamental understanding of electrostatic, polar, and van der Waals-London dispersion interaction in biological (macro)molecules.
The peptide-water interfacial reactions are of paramount importance in biological systems as solvation is ubiquitous in biological interaction with water as an inevitable component in blood and multicomponent body fluids. Understanding the peptide-water interaction is also essential regarding the RGD-targeted drug delivery mechanisms. It is, therefore, crucial to investigate the consequences of peptide-water interaction and aqueous solvation and its effects on internal bonding, charge distribution, and dielectric response. The solubility of 1FUV in water is generally attributed to the hydrogen bonding between water molecules and oxygen from the peptide backbone. To explore the details of peptide solvation properties, molecular dynamics (MD), molecular mechanics, and Monte Carlo methods, all based on specific force field parameters, are widely used. Nevertheless, the information on the role played by the hydrogen bond (HB) within RGD and at the RGD-water interface is in general still missing, as are the relevant ab initio MD (AIMD) studies [12,25,26,27]. Classical MD studies cannot provide a quantitative description of essentially quantum HBs (O–H and N–H) since they depend on the details of the force field parameterization. Consequently, the need for AIMD simulation in the study of dry and aqueous solvated peptide/protein over classical MD has been repeatedly recognized [28,29,30,31] and a more accurate AIMD calculation with a sufficiently large number of water molecules is highly desirable for a complete understanding of 1FUV aqueous solvation.
Although experimental techniques such as nuclear magnetic resonance (NMR), surface-enhanced Raman spectroscopy, X-ray photoelectron spectroscopy, etc., are used to investigate the structure and properties of RGD containing peptides [16,32,33,34,35], the experimental studies specifically involving 1FUV are scant. Moreover, although experimental studies can probe the consequent changes of the peptide structure due to water contact, they cannot provide critical information and details of microscopic properties and interfacial reaction mechanisms at the atomistic scale. Computer-assisted atomistic simulation thus seems the best alternative to investigate the molecular behavior and provide new insights to complement the meager experimental data [36,37], in many cases guiding the experiments as well. Nevertheless, only a few computational studies have so far been reported on the 1FUV peptide [38].
To answer the issues raised above and provide new insights into the study of the 1FUV peptide, this study is designed to simulate four models of dry and solvated 1FUV at 0 K and body temperature, namely 1FUV0, 1FUV310, 1FUVS0, and 1FUVS310, respectively. Our study provides a deeper understanding of the temperature-dependent structure, fundamental properties, peptide-water interfacial reaction mechanism, and solvation effect. To the best of our knowledge, our effort is the first example of a study dedicated to the 1FUV peptide using AIMD that ensures the reliability of the simulated structures and calculated properties. The rest of the paper is structured as follows. In the next section, a brief description of the modeling technique followed by the simulation workflow is described. The main part is the result and discussion section, in which we discuss our findings and articulate the future prospects for applying AIMD to other complex biomolecular systems.

2. Materials and Methods

2.1. Modeling 1FUV Peptide

The RGD peptides that specifically bind to αvβ3 and αvβ5 integrins have medical significance in designing the inhibitors of tumor and retinal angiogenesis [39,40], in osteoporosis [41], and in targeted drugs for tumor vasculature [42]. Considering the importance of αvβ3 and αvβ5 integrins [16] and as 1FUV binds only to these integrins, we choose to focus on 1FUV in this work.
The initial structure of 1FUV was taken from the Research Collaboratory for Structural Bioinformatics protein data bank (PDB) generated from the NMR measurements [43]. The stored PDB file contains about 19 such models with the same composition, the number of atoms, molecular weight, and the chemical formula, but slightly different molecular volumes. In this study, we chose the first model out of 19, which has 135 atoms comprised of 11 amino acids. The initial atomic positions of 1FUV were enclosed in a sufficiently large supercell of size 37.12 Å × 34.41Å × 28.95 Å to avoid any unintended interaction due to periodic boundary conditions. This structure was then fully relaxed by applying a density functional theory (DFT)-based Vienna ab initio simulation package (VASP) [44,45] known for its efficiency in structure optimization and energy minimization to obtain 1FUV0, the structure at 0 K. We used a projector-augmented wave method [46,47] with Perdew-Burke-Ernzerhof potential for the exchange correlation functional within the generalized gradient approximation [48]. We employed a relatively high cutoff energy of 500 eV and set the electronic convergence criterion at 10−5 eV. The force convergence criterion for ionic relaxation was set at 10−3 eV/Å. A single K-point calculation was used as the simulated models are sufficiently large supercells.
The fully relaxed structure at 0 K was then heated slowly to 310 K, equivalent to human body temperature, using NVT ensemble within 4 ps. A Nose thermostat was used to control the temperature of the heat bath [49,50]. A time step of 0.5 fs was used for the ionic motion during the simulation to ensure accuracy for integration of the equation of motion, especially for the light atoms such as hydrogen. The peptide model at this elevated temperature was then equilibrated for 5 ps to obtain the equilibrium structure at 310 K, viz. 1FUV310. Out of 5 ps run, the fluctuations of temperature and pressure were minimized by the initial 4 ps run, and the averaged output structure is taken from the last 1 ps run, thus used as a production run. The solvated model at 0 K (1FUVS0) was constructed based on 1FUV0 by surrounding the peptide structure with 80 water molecules. The 1FUVS0 was then again fully relaxed by using VASP. Similar to the dry model, this solvated structure was then heated to 310 K to obtain the 1FUVS310 model using NVT ensemble within 4 ps. The same protocol as used in the dry model was adopted to obtain the equilibrated structure of the solvated model at 310 K.
To reveal the structural stability of simulated models, we present the root-mean-square deviation (RMSD) of atomic positions of the peptide structure during equilibration at 310 K in Figure 1a. The time-RMSD profile of global structure shows a smooth increase followed by a leveling off at ~ 1.90 Å in the dry model and ~ 1.15 Å in the solvated model, respectively, indicating that the equilibrium has been reached. The constant RMSD shown in the inset for the last 1 ps run, used as a production run, ensures the sufficient equilibration of simulated structures. Figure 1b shows the velocity autocorrelation function (VACF), including exponential fitting, during the equilibration of dry and solvated structures at 310 K which confirms the loss of initial velocity from the previous configuration and attest to sufficient equilibration of the systems. It also implies that the amplitude of VACF is smaller in the solvated model and dies out quicker than that in the dry model. In addition, the exponential fitting shows the solvated model reaches equilibrium faster than the dry model.

2.2. Calculation of Properties

The AIMD simulated and VASP relaxed structures discussed above are then used as input for the electronic structures and dielectric response calculations using the in-house developed package of the orthogonalized linear combination of atomic orbitals (OLCAO) methodology [51]. The combination of these two DFT methods, VASP and OLCAO, is found to be robust and highly effective in studying the electronic structure and bonding properties, especially for complex biomolecules [52,53,54,55]. In recent years, we have further demonstrated that the combination of these two methods is extremely effective in dealing with new and emerging problems for complex protein structures related to SARS-CoV-2 [56,57,58,59,60].
The OLCAO is an all-electron method using local density approximation of DFT in which the Gaussian-type orbitals are employed for the expansion of the atomic basis set. There are many advantages of using OLCAO for electronic structure calculation, such as the flexibility of the basis sets choice, lower computational cost, and ease of bonding and charge transfer analysis using the Mulliken scheme [51]. The implementation of localized atomic orbitals in the basis expansion enables us to quantify the charge transfer and interatomic bonding via effective charge (Q*) on each atom and bond order (BO) value between pairs of atoms (ραβ) using the Mulliken scheme [61,62] as
Q α * = ∑ i ∑ n , o c c ∑ j . β C i α * n C j β n S i α , j β
ρ α β = ∑ n , o c c ∑ i , j C i α * n C j β n S i α , j β
where Siα,jβ is the overlap matrix between the basis Bloch sums of the orbital index (i, j) and atomic specification (α, β). N is the band index, i, j are the orbital quantum numbers, and Cjβ is the eigenvector coefficient.
From the calculated value of Q* the charge transfer between the ions due to atomic interaction can be quantified in terms of the partial charge (PC). It is defined as the deviation of Q* from the charge of neutral atom Q0 and is given by ΔQ = Q0 − Q*. The BO value calculated from the above equation gives the direct quantitative measure of bond strength between a pair of bonded atoms. The sum of all BO values in the system is the total bond order (TBO), which is a quantum mechanically calculated parameter. It is a single quantum mechanical metric helpful in assessing the internal cohesion and strength in a material [63]. The use of TBO to characterize the internal strength of a material and correlate it with the calculated physical properties is a novel and highly appealing concept.
OLCAO is equally suitable to calculate the imaginary part, ε2 (ℏω), of the frequency-dependent complex dielectric function within the random phase approximation of inter-band optical transition theory according to the following equation:
ε 2 ( ℏ ω ) = e 2 π m ω 2 ∫ B Z d k 3 ∑ n . l | 〈 ψ n ( k → , r → ) | − i ℏ ∇ → | ψ l ( k → , r → ) 〉 | 2 × f l ( k → ) [ 1 − f n ( k → ) ] δ [ E n ( k → ) − E l ( k → ) − ℏ ω ]

3. Results

3.1. Analysis of 1FUV Structures

The RGD peptide, 1FUV, studied in this work contains 11 amino acid sequence ACDCRGDCFCG, namely alanine (Ala, A), cysteine (Cys, C), asparagine (Asp, D), arginine (Arg, R), glycine (Gly, G), and phenylalanine (Phe, F). As this amino acid sequence contains four C, this peptide is also known as RGD-4C, in which the cysteine position is at 2, 4, 8, and 10th place of the sequence, named as Cys1, Cys2, Cys3, and Cys4, respectively. Each Cys amino acid contains the S atom and the arrangement of two Cys, and hence the S–S bond plays an important role in the configuration and dipole moment of this peptide [64]. Although the presence of four Cys groups in RGD-4C allows three possible combinations of S-S bonds, only two natural configurations of this peptide exist based on the S–S bonding arrangement [16]. The peptide with Cys1–Cys4 and Cys2–Cys3 disulfide bonding arrangement is known as RGD-A isomer, while the one with Cys1–Cys3, Cys2–Cys4 bonding arrangement is known as RGD-B isomer [16]. As RGD-A is a far stronger binder to integrin αvβ3 than RGD-B [16], it is reasonable to explore details of the structure and properties of RGD-A.
The relaxed structures for RGD-A peptide (1FUV) at 0 K and 310 K in the dry and aqueous environment are shown in Figure 2. 1FUV0 shows only a minor change in its structure due to VASP relaxation compared with the original NMR structure [16] taken from the PDB. However, there are noticeable changes in its conformation after AIMD simulation at 310 K. It shows a rotation of amino acid groups and changes in the bond length (BL) and BO value, which will be discussed in detail later. In the case of the solvated model, there again appear notable changes in the peptide structure. More water molecules interact with the peptide and there is a rotation of amino acid groups at 310 K as compared with 0 K. In both dry and solvated models, the S–S bond containing amino acids come closer at 310 K as compared with those at 0 K. However, the average S–S BL increases and hence its strength decreases at 310 K. It is noticed that the configurations of HBs, O–H, change at 310 K as compared with 0 K, resulting in the overall decrease in its strength in the dry model. However, in the solvated model, a large contribution of O–H bonding strength comes from peptide-water interaction due to which its value increases at 310 K. These differences are, although small, important on issues related to human health.

3.2. Electronic Structure and Interatomic Bonding

Electronic structures are the fundamental properties of a material that help to understand many other physical properties. The calculated total density of states (TDOS) for the simulated models is shown in Figure 3. The top of the valence band is set at 0 eV. The overall features of TDOS at 0 K and 310 K in dry and solvated models look similar as they contain the same amino acid groups; however, there is a decrease in the highest occupied molecular orbital (HOMO)-lowest unoccupied molecular orbital (LUMO) gap with the increase in temperature. The position of the energy difference between HOMO and LUMO plays the main role in various chemical reactions. The number of peaks decreases in solvated models as compared with the dry model, but the overall intensity increases. The HOMO-LUMO gap for 1FUV0 is ~ 3.40 eV, while its value in 1FUV310 decreases to ~2.10 eV due to a temperature rise, indicating that the molecule is chemically more reactive at a higher temperature. The HOMO-LUMO gap decreases largely in solvated models as compared with dry models. In 1FUVS0, this gap is ~0.3 eV, but it disappears in 1FUVS310.
The nature of bonding and its quantitative assessment is important in the electronic structure study of the material. While most of the MD studies interpret the bonding analysis based exclusively on the geometric distance configuration, OLCAO can quantify the strength of each bonded pair by providing the BO value. We calculate the TBO values in all 4 models and resolve them into partial pair-resolved BO (PBO) components. This is an unorthodox approach that generalizes the bonding between individual bonded pairs simply and effectively to reveal the hidden details which cannot be extricated by other methods, either experimental or computational.
In general, the BO value scales inversely with BL; however, the actual BO value of a bonded pair depends on its bonding environment and the atoms surrounding it. The calculated values of TBO and PBO are listed in Table 1, showing that at 310 K the TBO value decreases in the dry model, resulting from the elongation of some bonds. The TBO value of the solvated model increases at 310 K, opposite to the case of the dry model, which is due to the formation of more bonds between the peptide backbone and water at this elevated temperature. Figure 4 shows the comparison of PBO values for each bonded pair in simulated models at 0 K and 310 K in the form of a vertical bar diagram. The scale breaks are used along the vertical axis to delineate the minor differences between PBO values and show their magnified image. Figure 4 shows that there is a minor change in the PBO values due to the rise of temperature, except in the case of the O–H pair in the solvated model. The PBO for the O–H pair increases sharply in the solvated model at 310 K as compared with the model at 0 K, which arises due to the interaction of more water molecules with the peptide, and an increase in the HBs population at 310 K. It should be emphasized again that the small changes in these numbers at different temperature can be significant in the context of the human body.
We provide further information on the strength of individual bonded pairs presenting a complex distribution of BO versus BL in Figure 5. It clearly shows a wide variety of bonding pairs ranging from a strong covalent bond originating from N–H, N–C, C–H, C–C, C–O, and O–H pairs to weaker covalent bonds at longer BL. HBs are extremely weak but not entirely negligible because of their sheer numbers. The first group of strong covalent bonds arises from N–H, C–H, and O–H pairs with BL less than 1.2 Å. The next group of covalent BO pairs arises from C–C, C–O, and N–C bonds, whose BL lies in the range 1.2–1.7 Å. The C–C bonds show two clusterings, one having a high BO greater than 0.5 e from C–C double bonds, whereas the other group has a BO value less than 0.4 e from single C–C bonds. The BO values for N–C pairs also exhibit a group separation, having BO values greater than 0.4 e and less than 0.4 e.
The HBs are distributed in the range of 1.5–2.5 Å, with some O–H pairs having relatively higher BO values up to 0.1 e. In solvated models, the covalent O–H bonds within the water molecule have BO values higher than 0.2 e, while those bonds originating from the next neighbor or weak HBs beyond 2.5 Å have almost zero BO values. The bonding distributions due to C–S and S–S pairs have BO values comparable or slightly higher than those of strong HBs. The S–S bonding configurations are important in peptide structure, and they show significant changes in BO values due to the rise of temperature and solvation. Especially in the solvated model at 310 K, the BO of S–S bonds has a large difference. Wang et al. studied the influence of disulfide bond on properties of RGD peptide using MD and found that this bond causes the restriction on peptide mobility and affects its dipole moment [64]. The bonded pairs H–H and H–S with BL larger than 2.5 Å have BO values close to zero.

3.3. Hydrogen Bonding Analysis

Hydrogen bonding is especially important in a biomolecular system since it provides the key information to understand many intriguing phenomena, making the detailed analysis of HBs in the solvated models crucial. Unfortunately, in most MD studies, the analysis of HBs is carried out based exclusively on structural data, such as the BL and the location of HBs, failing to provide a quantitative assessment of HB strength and their local bonding arrangements. On the contrary, we bring forth quantitative information and deeper analysis of HBs using ab initio calculation in terms of BO-BL plot, as shown in Figure 5. It shows there exists some strong HB with BO value ~ 0.1 e, while most of the HBs have BO less than 0.1 e. At the body temperature, the overall strength of HB from the O–H bonding pair slightly decreases in the dry model while it increases in the solvated model.
In the case of solvated models, the change in the HB pattern is more complicated due to the presence of H2O molecules which increase the population and strength of O–H bonds at 310 K. The more HBs imply the stronger binding force of biomolecule with solvent. There is one HB originating from the N–H pair within the peptide backbone whose BO value at 310 K increases in the dry model and decreases in the solvated model because of the rotation and changes in the configuration of the peptide structure. The results on HB analysis in solvated models indicate a significant difference in the human body as compared with the vicinal environment of a dry peptide at zero temperature.

3.4. Partial Charge Distribution

The PC distribution of biomolecules is important for determining the long-range electrostatic and polar interaction involving the molecule. It entails a deep qualitative understanding of the structure and reactivity of a molecule. Rather than just making a rudimentary judgment by dividing PC into positive, negative, and neutral as adopted in many simplified calculations, we provide a quantitative assessment of PC of every atom in the simulated models. As the PC gives the quantity of charge transfer in atomic interaction, a positive value implies a loss, and a negative value means the gain of electronic charge.
The distribution of calculated PC for each atom in dry and solvated models is shown in Figure 6 and Figure 7, respectively, which provide a wealth of information that corroborates the electronic structure results discussed above. The solid symbols represent PC of ions at 0 K, while the hollow symbols represent PC at the body temperature. O (P) and H (P) denote O and H atoms from the peptide backbone while O (W) and H (W) denote those from the water molecule. In both dry and solvated models at 310 K, small changes appear in the PC of ions compared with their values in corresponding models at 0 K. Figure 6 and Figure 7 show that H and S lose charges, N gains charge, while C is both losing or gaining charges. In the dry model, O gains charge, while in the solvated model (Figure 7), few O atoms from the peptide backbone lose charge, and the remaining gain charge. The calculated PC values for each atom are scattered in a certain range rather than clustering and overlapping to each other, which implies the PC value is not the same even for the same atom in the same model, as their values depend on the local bonding environment. The PC of O and H in water and peptide backbone is not the same, implying that their local bonding characteristics are completely different (Figure 7). The dispersed distribution of calculated PC for each atom implies that the quantum mechanical calculation is essential to deal with such type of problem rather than using MD studies in which the PC values are predefined and kept constant in the simulation where the local environment changes.

3.5. Dielectric Response

The dielectric response of biomolecules in an aqueous solution has been studied for decades [65,66,67] and remains a subject of active research. The static dielectric constant of peptide/protein plays a crucial role in the calculation of the electrostatic field obtained by solving the Poisson–Boltzmann equation [68]. However, it is interesting to note that the dielectric properties of the 1FUV are seldom reported. We present a detailed analysis of the imaginary part, ε2 (ℏω), of the frequency-dependent dielectric function, based on the ab initio calculation under random phase approximation such as those used in inorganic crystals and glasses [69,70,71,72].
Figure 8a shows the spectrum for ε2 (ℏω) in 1FUV0 calculated using OLCAO, displaying a sharp peak at ~5.70 eV and one broad peak at ~14.30 eV. A similar spectrum is observed in the dry model at 310 K; however, the broad peak position shifts slightly to the lower energy side. In solvated models, in addition to the two peaks observed in dry models, there appears another sharp peak at ~0.1 eV, but with slightly less intensity (Figure 8b). One of the striking features observed is that at the body temperature, the absorption end of ε2 (ℏω) shifts to lower energy which is consistent with the decrease in the HOMO-LUMO gap at 310 K. We want to make it clear that our aim in this work is to compare the full dielectric spectra at 0 K and body temperature in case of a dry and solvated environment and the actual calculation of the dielectric constant is out of the scope of this work.

4. Discussions and Summary

In this section, we discuss the results presented above and summarize the main findings to answer the queries raised in the introduction. More accurate AIMD and DFT calculations adopted in this work provide us a wealth of quantitative information on the fundamental properties of 1FUV, which are not yet available or cannot be attained by an experimental protocol. Our findings provide a comparative analysis of the structure and properties of 1FUV and a better understanding of its solvation properties at 0 K and body temperature, 310 K. The key message is that there are discernible changes in the peptide structure and properties at 310 K, which could point to the potential implications for human health issues if such calculations can eventually be applied to large biomolecular systems.
Our results show noticeable changes in the conformation of amino acids at body temperature than at 0 K. The major structural differences observed in dry and solvated peptides at body temperature help in designing the specificity of integrins. The disulfide bridge connectivity greatly determines the orientation of the RGD motif and peptide structure which may dictate selectivity towards different integrins. The HOMO-LUMO gap of 1FUV decreases due to the rise of temperature and solvation effect. Such information is quite helpful to understand the chemical reactions involving biomolecules.
The bonding analysis is a significant result that helps to understand the interatomic interaction within 1FUV and other biomolecules. The TBO value decreases in the dry model while it increases in the solvated model at 310 K. The dry model at 0 K is more cohesive than that at 310 K, characterized by a higher TBO value, while it is just the opposite in the solvated models. The BO-BL plot analysis provides a detailed picture of interatomic bonding within the peptide structure and its explicit solvent environment. The solvated model shows a considerable peptide-water interactions occur at 310 K, increasing the HBs population than at 0 K. The detailed analysis of HBs and quantification of their strength is crucial in understanding the biological interactions inside the human body and for the design and delivery of RGD peptide-targeted drugs. The partial charge analysis provides another set of significant parameters to define the reactivity of biomolecules and their interactions with the integrins. It shows that the PC value strongly depends on the temperature as well as the solvated environment. Furthermore, the calculated PC values scatter in a certain range, even for the same atom in the same model. Quantum mechanical calculation is, therefore, absolutely necessary to capture all these crucial differences and to provide more accurate results inaccessible to classical MD studies.
Another important result is the calculation of the dielectric response function and its relation to the dielectric constant. Although this topic has attracted researchers from various backgrounds for a long time, unified and reliable ab initio calculations are still out of reach. As an initial assessment, we analyze the imaginary part of the dielectric function of the 1FUV peptide and study the differences in the dielectric spectra due to the rise in temperature and the solvation effect. Such optical spectra are helpful to estimate the long-range van der Waals–London interaction in biomolecules, as shown in the study of carbon nanotubes [73]. Overall, our results show apparent differences in bonding, PC, and dielectric function between the 1FUV at 0 K and the body temperature in a dry and solvated environment. These results can be quite helpful to develop and refine the force field parameters used in MD simulation. Unfortunately, as there is a lack of theoretical studies in 1FUV and the focus of the existing experimental study [16] is different, we are unable to compare our results with the previous results. However, we anticipate that our results will inspire future experimental and/or theoretical works to investigate further the 1FUV and related peptides, including the effect of salts and different numbers of water molecules, which significantly influence peptide structure and properties.
In the past, we have successfully demonstrated the applicability of AIMD in the study of oxide glassy materials and fluoride salts [74,75,76,77]. Our accomplishment in the present work is applying AIMD to a relatively small peptide structure that looks quite promising. It ensures the feasibility and encourages a further study of other possibly more complex biomolecules by using AIMD, of course assuming sufficient computing resources to be available.

Author Contributions

W.-Y.C., K.B. and P.A. conceived the project. K.B. and W.-Y.C. performed the calculations and made figures. K.B. and W.-Y.C. drafted the paper with inputs from P.A., R.P. and B.J. All authors participated in the discussion and interpretation of the results. All authors edited and proofread the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This project is funded by the National Science Foundation of USA: RAPID DMR/CMMT-2028803.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data that supports the results in this study are available from the corresponding author upon reasonable request.

Acknowledgments

This research used the resources of the National Energy Research Scientific Computing Center supported by DOE under Contract No. DE-AC03-76SF00098 and the Research Computing Support Services (RCSS) of the University of Missouri System. R.P. would like to acknowledge the support from the 1000-Talents Program of the Chinese Foreign Experts Bureau, the School of Physics, the University of the Chinese Academy of Sciences, Beijing, and the Institute of Physics of the Chinese Academy of Sciences, Beijing.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alipour, M.; Baneshi, M.; Hosseinkhani, S.; Mahmoudi, R.; Jabari Arabzadeh, A.; Akrami, M.; Mehrzad, J.; Bardania, H. Recent progress in biomedical applications of RGD-based ligand: From precise cancer theranostics to biomaterial engineering: A systematic review. J. Biomed. Mater. Res. A 2020, 108, 839–850. [Google Scholar] [CrossRef] [PubMed]
  2. Jeschke, B.; Meyer, J.; Jonczyk, A.; Kessler, H.; Adamietz, P.; Meenen, N.M.; Kantlehner, M.; Goepfert, C.; Nies, B. RGD-peptides for tissue engineering of articular cartilage. Biomaterials 2002, 23, 3455–3463. [Google Scholar] [CrossRef]
  3. Dijkgraaf, I.; Beer, A.J.; Wester, H.-J. Application of RGD-containing peptides as imaging probes for alphavbeta3 expression. Front. Biosci. 2009, 14, 887–899. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Ko, H.Y.; Choi, K.-J.; Lee, C.H.; Kim, S. A multimodal nanoparticle-based cancer imaging probe simultaneously targeting nucleolin, integrin αvβ3 and tenascin-C proteins. Biomaterials 2011, 32, 1130–1138. [Google Scholar] [CrossRef] [PubMed]
  5. Liu, S. Radiolabeled multimeric cyclic RGD peptides as integrin αvβ3 targeted radiotracers for tumor imaging. Mol. Pharm. 2006, 3, 472–487. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, S. Radiolabeled cyclic RGD peptide bioconjugates as radiotracers targeting multiple integrins. Biocon. Chem. 2015, 26, 1413–1438. [Google Scholar] [CrossRef] [Green Version]
  7. Danhier, F.; Le Breton, A.; Préat, V. RGD-based strategies to target Alpha(v) Beta (3) integrin in cancer therapy and diagnosis. Mol. Pharm. 2012, 9, 2961–2973. [Google Scholar] [CrossRef]
  8. Yin, H.-Q.; Bi, F.-L.; Gan, F. Rapid synthesis of cyclic RGD conjugated gold nanoclusters for targeting and fluorescence imaging of Melanoma A375 Cells. Biocon. Chem. 2015, 26, 243–249. [Google Scholar] [CrossRef]
  9. Fu, S.; Xu, X.; Ma, Y.; Zhang, S.; Zhang, S. RGD peptide-based non-viral gene delivery vectors targeting integrin αvβ3 for cancer therapy. J. Drug Target. 2019, 27, 1–11. [Google Scholar] [CrossRef]
  10. Sun, Y.; Kang, C.; Liu, F.; Zhou, Y.; Luo, L.; Qiao, H. RGD peptide-based target drug delivery of doxorubicin nanomedicine. Drug Dev. Res. 2017, 78, 283–291. [Google Scholar] [CrossRef]
  11. Nikfar, Z.; Shariatinia, Z. Tripeptide arginyl-glycyl-aspartic acid (RGD) for delivery of cyclophosphamide anticancer drug: A computational approach. Int. J. Nano Dimen. 2020, 11, 312–336. [Google Scholar]
  12. Nikfar, Z.; Shariatinia, Z. The RGD tripeptide anticancer drug carrier: DFT computations and molecular dynamics simulations. J. Mol. Liq. 2019, 281, 565–583. [Google Scholar] [CrossRef]
  13. Jin, H.; Varner, J. Integrins: Roles in cancer development and as treatment targets. Br. J. Cancer 2004, 90, 561–565. [Google Scholar] [CrossRef]
  14. Sun, Y.-X.; Zeng, X.; Meng, Q.-F.; Zhang, X.-Z.; Cheng, S.-X.; Zhuo, R.-X. The influence of RGD addition on the gene transfer characteristics of disulfide-containing polyethyleneimine/DNA complexes. Biomaterials 2008, 29, 4356–4365. [Google Scholar] [CrossRef] [PubMed]
  15. Temming, K.; Schiffelers, R.; Molema, G.; Kok, R.J. RGD-based strategies for selective delivery of therapeutics and imaging agents to the tumour vasculature. Drug Resist. Updates 2005, 8, 381–402. [Google Scholar] [CrossRef]
  16. Assa-Munt, N.; Jia, X.; Laakkonen, P.; Ruoslahti, E. Solution structures and integrin binding activities of an RGD peptide with two isomers. Biochemistry 2001, 40, 2373–2378. [Google Scholar] [CrossRef]
  17. Ruoslahti, E.; Pierschbacher, M.D. New perspectives in cell adhesion: RGD and integrins. Science 1987, 238, 491–497. [Google Scholar] [CrossRef]
  18. Paradise, R.K.; Lauffenburger, D.A.; Van Vliet, K.J. Acidic extracellular pH promotes activation of integrin αvβ3. PLoS ONE 2011, 6, e15746. [Google Scholar] [CrossRef] [Green Version]
  19. Craig, D.; Gao, M.; Schulten, K.; Vogel, V. Structural insights into how the MIDAS ion stabilizes integrin binding to an RGD peptide under force. Structure 2004, 12, 2049–2058. [Google Scholar] [CrossRef] [Green Version]
  20. Ruoslahti, E. Integrins. J. Clin. Investig. 1991, 87, 1–5. [Google Scholar] [CrossRef]
  21. Barczyk, M.; Carracedo, S.; Gullberg, D. Integrins. Cell Tiss. Res. 2010, 339, 269–280. [Google Scholar] [CrossRef] [Green Version]
  22. Takada, Y.; Ye, X.; Simon, S. The integrins. Genome Biol. 2007, 8, 1–9. [Google Scholar] [CrossRef] [Green Version]
  23. Li, J.; Wang, Y.; Liang, R.; An, X.; Wang, K.; Shen, G.; Tu, Y.; Zhu, J.; Tao, J. Recent advances in targeted nanoparticles drug delivery to melanoma. Nanomed. Nanotechnol. Biol. Med. 2015, 11, 769–794. [Google Scholar] [CrossRef]
  24. Marinelli, L.; Lavecchia, A.; Gottschalk, K.-E.; Novellino, E.; Kessler, H. Docking studies on αvβ3 integrin ligands: Pharma-cophore refinement and implications for drug design. J. Med. Chem. 2003, 46, 4393–4404. [Google Scholar] [CrossRef]
  25. Li, N.; Qiu, S.; Fang, Y.; Wu, J.; Li, Q. Comparison of Linear vs. Cyclic RGD Pentapeptide Interactions with Integrin αvβ3 by molecular dynamics simulations. Biology 2021, 10, 688. [Google Scholar] [CrossRef]
  26. Damjanovic, J.; Miao, J.; Huang, H.; Lin, Y.-S. Elucidating solution structures of cyclic peptides using molecular dynamics simulations. Chem. Rev. 2021, 121, 2292–2324. [Google Scholar] [CrossRef] [PubMed]
  27. Muir, J.; Costa, D.; Idriss, H. DFT computational study of the RGD peptide interaction with the rutile TiO2 (110) surface. Surf. Sci. 2014, 624, 8–14. [Google Scholar] [CrossRef]
  28. Xu, M.; Zhu, T.; Zhang, J.Z.H. A force balanced fragmentation method for ab initio molecular dynamic simulation of protein. Front. Chem. 2018, 6, 189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Isborn, C.M.; Mar, B.D.; Curchod, B.F.E.; Tavernelli, I.; Martínez, T.J. The charge transfer problem in density functional theory calculations of aqueously solvated molecules. J. Phys. Chem. B 2013, 117, 12189–12201. [Google Scholar] [CrossRef]
  30. Ufimtsev, I.S.; Luehr, N.; Martinez, T.J. Charge transfer and polarization in solvated proteins from ab initio molecular dynamics. J. Phys. Chem. Let. 2011, 2, 1789–1793. [Google Scholar] [CrossRef]
  31. Wei, D.; Guo, H.; Salahub, D. Conformational dynamics of an alanine dipeptide analog: An ab initio molecular dynamics study. Phys. Rev. E 2001, 64, 011907. [Google Scholar] [CrossRef]
  32. Jois, S.D.; Tambunan, U.S.; Chakrabarti, S.; Siahaan, T.J. Solution structure of a cyclic RGD peptide that inhibits platelet ag-gregation. J. Biomol. Struct. Dyn. 1996, 14, 1–11. [Google Scholar] [CrossRef]
  33. Portela, A.; Yano, T.-A.; Santschi, C.; Martin, O.J.F.; Tabata, H.; Hara, M. Highly sensitive SERS analysis of the cyclic Arg-Gly-Asp peptide ligands of cells using nanogap antennas. J. Biophotonics 2016, 10, 294–302. [Google Scholar] [CrossRef]
  34. Stevens, J.S.; de Luca, A.C.; Pelendritis, M.; Terenghi, G.; Downes, S.; Schroeder, S.L.M. Quantitative analysis of complex amino acids and RGD peptides by X-ray photoelectron spectroscopy (XPS). Surf. Interf. Anal. 2013, 45, 1238–1246. [Google Scholar] [CrossRef] [Green Version]
  35. Zhang, R.; Huang, G.; Chen, L.; Wu, W. Studies on the conformations and Hydrogen-bonding interactions of RGD tri-peptide in aqueous solutions by molecular dynamics simulations and 2D-NOESY spectroscopy. J. Solut. Chem. 2015, 44, 1281–1291. [Google Scholar] [CrossRef]
  36. Karplus, M.; McCammon, J.A. Molecular dynamics simulations of biomolecules. Nat. Struct. Biol. 2002, 9, 646–652. [Google Scholar] [CrossRef] [PubMed]
  37. Wakefield, A.E.; Wuest, W.M.; Voelz, V.A. Molecular simulation of conformational pre-organization in cyclic RGD peptides. J. Chem. Inf. Model. 2015, 55, 806–813. [Google Scholar] [CrossRef]
  38. Adhikari, P.; Wen, A.M.; French, R.H.; Parsegian, V.A.; Steinmetz, N.F.; Podgornik, R.; Ching, W.-Y. Electronic structure, dielectric response and surface charge distribution of RGD (1FUV) peptide. Sci. Rep. 2014, 4, 5605. [Google Scholar] [CrossRef] [PubMed]
  39. Brooks, P.C.; Strömblad, S.; Klemke, R.; Visscher, D.; Sarkar, F.H.; Cheresh, D.A. Antiintegrin alpha v beta 3 blocks human breast cancer growth and angiogenesis in human skin. J. Clin. Investig. 1995, 96, 1815–1822. [Google Scholar] [CrossRef] [Green Version]
  40. Hammes, H.-P.; Brownlee, M.; Jonczyk, A.; Sutter, A.; Preissner, K.T. Subcutaneous injection of a cyclic peptide antagonist of vitronectin receptor–type integrins inhibits retinal neovascularization. Nat. Med. 1996, 2, 529–533. [Google Scholar] [CrossRef]
  41. Yamamoto, M.; Fisher, J.E.; Gentile, M.; Seedor, J.G.; Leu, C.-T.; Rodan, S.B.; Rodan, G.A. The integrin ligand echistatin prevents bone loss in ovariectomized mice and rats. Endocrinology 1998, 139, 1411–1419. [Google Scholar] [CrossRef]
  42. Arap, W.; Pasqualini, R.; Ruoslahti, E. Cancer treatment by targeted drug delivery to tumor vasculature in a mouse model. Sciences 1998, 279, 377–380. [Google Scholar] [CrossRef] [Green Version]
  43. Solution Structure of an RGD Peptide Isomer-A. Available online: https://www.rcsb.org/structure/1FUV (accessed on 7 June 2021).
  44. Kresse, G.; Furthmüller, J. Efficient iterative schemes for ab initio total-energy calculations using a plane-wave basis set. Phys. Rev. B 1996, 54, 169–186. [Google Scholar] [CrossRef]
  45. Kresse, G.; Furthmüller, J. Vienna Ab-Initio Simulation Package (VASP); University of Vienna: Vienna, Austria, 2001. [Google Scholar]
  46. Blöchl, P.E. Projector augmented-wave method. Phys. Rev. B 1994, 50, 17953. [Google Scholar] [CrossRef] [Green Version]
  47. Kresse, G.; Joubert, D. From ultrasoft pseudopotentials to the projector augmented-wave method. Phys. Rev. B 1999, 59, 1758. [Google Scholar] [CrossRef]
  48. Perdew, J.P.; Burke, K.; Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 1996, 77, 3865. [Google Scholar] [CrossRef] [Green Version]
  49. Nosé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511–519. [Google Scholar] [CrossRef] [Green Version]
  50. Nosé, S. A molecular dynamics method for simulations in the canonical ensemble. Mol. Phys. 1984, 52, 255–268. [Google Scholar] [CrossRef]
  51. Ching, W.-Y.; Rulis, P. Electronic Structure Methods for Complex Materials: The Orthogonalized Linear Combination of Atomic Orbitals; Oxford University Press: Oxford, UK, 2012. [Google Scholar]
  52. Eifler, J.; Podgornik, R.; Steinmetz, N.F.; French, R.H.; Parsegian, V.A.; Ching, W. Charge distribution and hydrogen bonding of a collagen α2-chain in vacuum, hydrated, neutral, and charged structural models. Int. J. Quantum Chem. 2016, 116, 681–691. [Google Scholar] [CrossRef]
  53. Poudel, L.; Steinmetz, N.F.; French, R.; Parsegian, V.A.; Podgornik, R.; Ching, W.-Y. Implication of the solvent effect, metal ions and topology in the electronic structure and hydrogen bonding of human telomeric G-quadruplex DNA. Phys. Chem. Chem. Phys. 2016, 18, 21573–21585. [Google Scholar] [CrossRef] [PubMed]
  54. Poudel, L.; Tamerler, C.; Misra, A.; Ching, W.-Y. Atomic-scale quantification of interfacial binding between peptides and in-organic crystals: The case of calcium carbonate binding peptide on aragonite. J. Phys. Chem. C 2017, 121, 28354–28363. [Google Scholar] [CrossRef]
  55. Jawad, B.; Poudel, L.; Podgornik, R.; Steinmetz, N.F.; Ching, W.-Y. Molecular mechanism and binding free energy of doxoru-bicin intercalation in DNA. Phys. Chem. Chem. Phys. 2019, 21, 3877–3893. [Google Scholar] [CrossRef]
  56. Adhikari, P.; Ching, W.-Y. Amino acid interacting network in the receptor-binding domain of SARS-CoV-2 spike protein. RSC Adv. 2020, 10, 39831–39841. [Google Scholar] [CrossRef]
  57. Adhikari, P.; Li, N.; Shin, M.; Steinmetz, N.F.; Twarock, R.; Podgornik, R.; Ching, W.-Y. Intra-and intermolecular atomic-scale interactions in the receptor binding domain of SARS-CoV-2 spike protein: Implication for ACE2 receptor binding. Phys. Chem. Chem. Phys. 2020, 22, 18272–18283. [Google Scholar] [CrossRef]
  58. Ching, W.-Y.; Adhikari, P.; Jawad, B.; Podgornik, R. Ultra-large-scale ab initio quantum chemical computation of bio-molecular systems: The case of spike protein of SARS-CoV-2 virus. Comput. Struct. Biotech. J. 2021, 19, 1288–1301. [Google Scholar] [CrossRef] [PubMed]
  59. Jawad, B.; Adhikari, P.; Podgornik, R.; Ching, W.-Y. Key interacting residues between RBD of SARS-CoV-2 and ACE2 receptor: Combination of molecular dynamic simulation and density functional calculation. J. Chem. Inf. Model. 2021, 61, 4425–4441. [Google Scholar] [CrossRef] [PubMed]
  60. Jawad, B.; Poudel, L.; Podgornik, R.; Ching, W.-Y. Thermodynamic dissection of the intercalation binding process of doxoru-bicin to dsDNA with implications of ionic and solvent effects. J. Phys. Chem. B 2020, 124, 7803–7818. [Google Scholar] [CrossRef] [PubMed]
  61. Mulliken, R.S. Electronic population analysis on LCAO–MO molecular wave functions. I. J. Chem. Phys. 1955, 23, 1833–1840. [Google Scholar] [CrossRef] [Green Version]
  62. Mulliken, R.S. Electronic population analysis on LCAO–MO molecular wave functions. II. overlap populations, bond orders, and covalent bond energies. J. Chem. Phys. 1955, 23, 1841–1846. [Google Scholar] [CrossRef]
  63. Dharmawardhana, C.C.; Misra, A.; Ching, W.-Y. Quantum mechanical metric for internal cohesion in cement crystals. Sci. Rep. 2014, 4, 7332. [Google Scholar] [CrossRef]
  64. Wang, Y.; Goh, S.-Y.; Kuczera, K. Molecular dynamics study of disulfide bond influence on properties of an RGD peptide. J. Pept. Res. 1999, 53, 188–200. [Google Scholar] [CrossRef]
  65. Löffler, G.; Schreiber, H.; Steinhauser, O. Calculation of the dielectric properties of a protein and its solvent: Theory and a case study. J. Mol. Biol. 1997, 270, 520–534. [Google Scholar] [CrossRef] [Green Version]
  66. Schutz, C.N.; Warshel, A. What are the dielectric?constants? of proteins and how to validate electrostatic models? Proteins Struct. Funct. Bioinform. 2001, 44, 400–417. [Google Scholar] [CrossRef]
  67. King, G.; Lee, F.S.; Warshel, A. Microscopic simulations of macroscopic dielectric constants of solvated proteins. J. Chem. Phys. 1991, 95, 4366–4377. [Google Scholar] [CrossRef]
  68. Gilson, M.K.; Honig, B.H. Calculation of electrostatic potentials in an enzyme active site. Nature 1987, 330, 84–86. [Google Scholar] [CrossRef] [PubMed]
  69. Aryal, S.; Rulis, P.; Ouyang, L.; Ching, W.Y. Structure and properties of the low-density phase ι-Al2O3 from first principles. Phys. Rev. B 2011, 84, 174123. [Google Scholar] [CrossRef] [Green Version]
  70. Baral, K.; Adhikari, P.; Ching, W.Y. Ab initio modeling of the electronic structures and physical properties of a-Si1− xGexO2 Glass (x = 0 to 1). J. Am. Ceram. Soc. 2016, 99, 3677–3684. [Google Scholar] [CrossRef]
  71. Hasan, S.; Baral, K.; Li, N.; Ching, W.-Y. Structural and physical properties of 99 complex bulk chalcogenides crystals using first-principles calculations. Sci. Rep. 2021, 11, 1–18. [Google Scholar] [CrossRef]
  72. Hussain, A.; Mehmood, S.; Rasool, M.; Aryal, S.; Rulis, P.; Ching, W. Electronic structure, mechanical, and optical properties of CaO Al2O3 system: A first principles approach. Ind. J. Phys. 2016, 90, 917–929. [Google Scholar] [CrossRef]
  73. Rajter, R.F.; French, R.; Ching, W.Y.; Carter, W.C.; Chiang, Y.M. Calculating van der Waals-London dispersion spectra and Hamaker coefficients of carbon nanotubes in water from ab initio optical properties. J. Appl. Phys. 2007, 101, 054303. [Google Scholar] [CrossRef]
  74. Baral, K.; Li, A.; Ching, W.-Y. Ab initio modeling of structure and properties of single and mixed alkali silicate glasses. J. Phys. Chem. A 2017, 121, 7697–7708. [Google Scholar] [CrossRef] [PubMed]
  75. Baral, K.; Li, A.; Ching, W.-Y. Ab initio molecular dynamics simulation of Na-doped aluminosilicate glasses and glass-water interaction. AIP Adv. 2019, 9, 075218. [Google Scholar] [CrossRef] [Green Version]
  76. Baral, K.; Li, A.; Ching, W.-Y. Ab initio study of hydrolysis effects in single and ion-exchanged alkali aluminosilicate glasses. J. Phys. Chem. B 2020, 124, 8418–8433. [Google Scholar] [CrossRef] [PubMed]
  77. Baral, K.; San, S.; Sakidja, R.; Couet, A.; Sridharan, K.; Ching, W.-Y. Temperature-dependent properties of molten Li2BeF4 Salt using Ab initio molecular dynamics. ACS Omega 2021, 6, 19822–19835. [Google Scholar] [CrossRef] [PubMed]
Polymers 13 03434 g001 550
Figure 1. Calculated (a) RMSD and (b) VACF including exponential fitting, of simulated dry and solvated 1FUV peptide at 310 K.
Figure 1. Calculated (a) RMSD and (b) VACF including exponential fitting, of simulated dry and solvated 1FUV peptide at 310 K.
Polymers 13 03434 g001
Polymers 13 03434 g002 550
Figure 2. Structures of dry (a,b) and solvated (c,d) 1FUV peptide at 0 K and 310 K, respectively. (N = Blue, C = Brown, O = Red, H = Black, and S = Yellow).
Figure 2. Structures of dry (a,b) and solvated (c,d) 1FUV peptide at 0 K and 310 K, respectively. (N = Blue, C = Brown, O = Red, H = Black, and S = Yellow).
Polymers 13 03434 g002
Polymers 13 03434 g003 550
Figure 3. Calculated TDOS for dry and solvated 1FUV peptide at 0 K and 310 K.
Figure 3. Calculated TDOS for dry and solvated 1FUV peptide at 0 K and 310 K.
Polymers 13 03434 g003
Polymers 13 03434 g004 550
Figure 4. Comparison of PBO values for each bonded pair in dry and solvated 1FUV at 0 K and 310 K.
Figure 4. Comparison of PBO values for each bonded pair in dry and solvated 1FUV at 0 K and 310 K.
Polymers 13 03434 g004
Polymers 13 03434 g005 550
Figure 5. BO versus BL plot for dry and solvated 1FUV at 0 K and 310 K.
Figure 5. BO versus BL plot for dry and solvated 1FUV at 0 K and 310 K.
Polymers 13 03434 g005
Polymers 13 03434 g006 550
Figure 6. Calculated partial charge for each atom in dry 1FUV at 0 K and 310 K. (The subscripts 0 and 310 denote the temperature 0 K and 310 K, respectively).
Figure 6. Calculated partial charge for each atom in dry 1FUV at 0 K and 310 K. (The subscripts 0 and 310 denote the temperature 0 K and 310 K, respectively).
Polymers 13 03434 g006
Polymers 13 03434 g007 550
Figure 7. Calculated partial charge for each atom in solvated models of 1FUV at 0 K and 310 K. (The subscripts 0 and 310 denote the temperature 0 K and 310 K, respectively).
Figure 7. Calculated partial charge for each atom in solvated models of 1FUV at 0 K and 310 K. (The subscripts 0 and 310 denote the temperature 0 K and 310 K, respectively).
Polymers 13 03434 g007
Polymers 13 03434 g008 550
Figure 8. Calculated imaginary part of the complex dielectric function for (a) dry and (b) solvated 1FUV at 0 K and 310 K.
Figure 8. Calculated imaginary part of the complex dielectric function for (a) dry and (b) solvated 1FUV at 0 K and 310 K.
Polymers 13 03434 g008
Table
Table 1. Calculated TBO and PBO values for each bonded pair in simulated models.
Table 1. Calculated TBO and PBO values for each bonded pair in simulated models.
ModelsTBOPBO
C‒CC‒HC‒OC‒SH‒HH‒SH‒ON‒CN‒HN‒OS‒S
1FUV054.27211.99116.6807.6690.6370.2320.0370.5299.9796.1810.0010.336
1FUV31053.88711.94416.7057.5810.6170.1980.0270.5309.9266.0340.0000.323
1FUVS093.1669.99015.9479.2340.4000.2510.02542.2078.8595.9050.0030.346
1FUVS31095.66310.47916.1199.1250.4340.2920.00543.7249.1116.1190.0000.246
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Baral, K.; Adhikari, P.; Jawad, B.; Podgornik, R.; Ching, W.-Y. Solvent Effect on the Structure and Properties of RGD Peptide (1FUV) at Body Temperature (310 K) Using Ab Initio Molecular Dynamics. Polymers 2021, 13, 3434. https://doi.org/10.3390/polym13193434

AMA Style

Baral K, Adhikari P, Jawad B, Podgornik R, Ching W-Y. Solvent Effect on the Structure and Properties of RGD Peptide (1FUV) at Body Temperature (310 K) Using Ab Initio Molecular Dynamics. Polymers. 2021; 13(19):3434. https://doi.org/10.3390/polym13193434

Chicago/Turabian Style

Baral, Khagendra, Puja Adhikari, Bahaa Jawad, Rudolf Podgornik, and Wai-Yim Ching. 2021. "Solvent Effect on the Structure and Properties of RGD Peptide (1FUV) at Body Temperature (310 K) Using Ab Initio Molecular Dynamics" Polymers 13, no. 19: 3434. https://doi.org/10.3390/polym13193434

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop