Validation of Molecular Dynamics Simulations for Prediction of Three-Dimensional Structures of Small Proteins

A computational odyssey: uncovering classical β-lactamase inhibitors in dry fruits

In the antibacterial arsenal, β-lactams have held a prominent position, but increasing resistance due to unauthorized use and genetic factors requires new strategies. Combining β-lactamase inhibitors with broad-spectrum β-lactams proves effective in combating this resistance. ESBL producers demand new inhibitors, leading to the exploration of plant-derived secondary metabolites for potent β-lactam antibiotics or alternative inhibitors. Using virtual screening, molecular docking, ADMET analysis, and molecular dynamic simulation, this study actively analyzed the inhibitory activity of figs, cashews, walnuts, and peanuts against SHV-1, NDM-1, KPC-2, and OXA-48 β-lactamases. Using AutoDock Vina, the docking affinities of various compounds for target enzymes were initially screened, revealing 12 bioactive compounds with higher affinities for the target enzymes compared to Avibactam and Tazobactam. Top-scoring metabolites, including Oleanolic acid, Protocatechuic acid, and Tannin, were subjected to MD simulation studies to further analyze the stability of the docked complexes using WebGro. The simulation coordinates, in terms of RMSD, RMSF, SASA, Rg, and hydrogen bonds formed, showed that these phytocompounds are stable enough to retain in the active sites at various orientations. The PCA and FEL analysis also showed the stability of the dynamic motion of Cα residues of phytochemical-bound enzymes. The pharmacokinetic analysis of the top phytochemicals was performed to analyze their bioavailability and toxicity. This study provides new insights into the therapeutic potential of phytochemicals of selected dry fruits and contributes to future experimental studies to identify βL inhibitors from plants.Communicated by Ramaswamy H. Sarma.

A new horizon in the phosphorylated sites of AGA: the structural impact of C163S mutation in aspartylglucosaminuria through molecular dynamics simulation

Aspartylglucosaminuria (AGU) is a lysosomal storage disorder caused by insufficient aspartylglucosaminidase (AGA) activity leading to chronic neurodegeneration. We utilized the PhosphoSitePlus tool to identify the AGA protein's phosphorylation sites. The phosphorylation was induced on the specific residue of the three-dimensional AGA protein, and the structural changes upon phosphorylation were studied via molecular dynamics simulation. Furthermore, the structural behaviour of C163S mutation and C163S mutation with adjacent phosphorylation was investigated. We have examined the structural impact of phosphorylated forms and C163S mutation in AGA. Molecular dynamics simulations (200 ns) exposed patterns of deviation, fluctuation, and change in compactness of Y178 phosphorylated AGA protein (Y178-p), T215 phosphorylated AGA protein (T215-p), T324 phosphorylated AGA protein (T324-p), C163S mutant AGA protein (C163S), and C163S mutation with Y178 phosphorylated AGA protein (C163S-Y178-p). Y178-p, T215-p, and C163S mutation demonstrated an increase in intramolecular hydrogen bonds, leading to greater compactness of the AGA forms. Principle component analysis (PCA) and Gibbs free energy of the phosphorylated/C163S mutation structures exhibit transition in motion/orientation than Wild type (WT). T215-p may be more dominant among these than the other studied phosphorylated forms. It might contribute to hydrolyzing L-asparagine functioning as an asparaginase, thereby regulating neurotransmitter activity. This study revealed structural insights into the phosphorylation of Y178, T215, and T324 in AGA protein. Additionally, it exposed the structural changes of the C163S mutation and C163S-Y178-p of AGA protein. This research will shed light on a better understanding of AGA's phosphorylated mechanism.Communicated by Ramaswamy H. Sarma.

Identifying molecular targets of Aspiletrein-derived steroidal saponins in lung cancer using network pharmacology and molecular docking-based assessments

Lung cancer is one of the leading cancers and causes of cancer-related deaths worldwide. Due to its high prevalence and mortality rate, its clinical management remains a significant challenge. Previously, the in vitro anticancer activity of Aspiletrein A, a steroid and a saponin from Aspidistra letreae, against non-small cell lung cancer (NSCLC) cells was reported. However, the anticancer molecular mechanism of other Aspiletreins from A. letreae remains unknown. Using in silico network pharmacology approaches, the targets of Aspiletreins were predicted using the Swiss Target Prediction database. In addition, key mediators in NSCLC were obtained from the Genetic databases. The compound-target interacting networks were constructed using the STRING database and Cytoscape, uncovering potential targets, including STAT3, VEGFA, HSP90AA1, FGF2, and IL2. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analysis demonstrated that several pathways were highly relevant to cancer pathogenesis. Additionally, molecular docking and molecular dynamic analyses revealed the interaction between key identified targets and Aspiletreins, including hydrogen bonding and Van der Waals interaction. This study provides potential targets of Aspiletreins in NSCLC, and its approach of integrating network pharmacology, bioinformatics, and molecular docking is a powerful tool for investigating the mechanism of new drug targets on a specific disease.

A comparative study of 5- fluorouracil, doxorubicin, methotrexate, paclitaxel for their inhibition ability for Mpro of nCoV: Molecular docking and molecular dynamics simulations

A new corona virus (nCoV) is aetiological agent responsible for the viral pneumonia epidemic. Three is no specific therapeutic medicines available for the treatment of this condition and also effective treatment choices are few. In this work author tried to investigate some repurposing drug such as 5- fluorouracil, doxorubicin, methotrexate and paclitaxel against the main protease (Mpro) of nCoV by the computational model. Molecular docking was performed to screen out the best compound and doxorubicin was found to have minimum binding energy −121.89 kcal/mol. To further study, MD simulations were performed at 300 K and the result successfully corroborate the energy obtained by molecular docking. Temperature dependent MD simulation of the best molecule that is doxorubicin obtained from docking result was performed to check the variation in structural changes in Mpro of nCoV at 290 K, 310 K, 320 K and 325 K. It is sound that doxorubicin binds effectively with Mpro of nCoV at 290 K. Further ADME properties of the 5- fluorouracil, doxorubicin, methotrexate and paclitaxel were also evaluated to understand the bioavailability.

Docking cyclic peptides formed by a disulfide bond through a hierarchical strategy

MOTIVATION

Cyclization is a common strategy to enhance the therapeutic potential of peptides. Many cyclic peptide drugs have been approved for clinical use, in which the disulfide-driven cyclic peptide is one of the most prevalent categories. Molecular docking is a powerful computational method to predict the binding modes of molecules. For protein-cyclic peptide docking, a big challenge is considering the flexibility of peptides with conformers constrained by cyclization.

RESULTS

Integrating our efficient peptide 3D conformation sampling algorithm MODPEP2.0 and knowledge-based scoring function ITScorePP, we have proposed an extended version of our hierarchical peptide docking algorithm, named HPEPDOCK2.0, to predict the binding modes of the peptide cyclized through a disulfide against a protein. Our HPEPDOCK2.0 approach was extensively evaluated on diverse test sets and compared with the state-of-the-art cyclic peptide docking program AutoDock CrankPep (ADCP). On a benchmark dataset of 18 cyclic peptide-protein complexes, HPEPDOCK2.0 obtained a native contact fraction of above 0.5 for 61% of the cases when the top prediction was considered, compared with 39% for ADCP. On a larger test set of 25 cyclic peptide-protein complexes, HPEPDOCK2.0 yielded a success rate of 44% for the top prediction, compared with 20% for ADCP. In addition, HPEPDOCK2.0 was also validated on two other test sets of 10 and 11 complexes with apo and predicted receptor structures, respectively. HPEPDOCK2.0 is computationally efficient and the average running time for docking a cyclic peptide is about 34 min on a single CPU core, compared with 496 min for ADCP. HPEPDOCK2.0 will facilitate the study of the interaction between cyclic peptides and proteins and the development of therapeutic cyclic peptide drugs.

AVAILABILITY AND IMPLEMENTATION

http://huanglab.phys.hust.edu.cn/hpepdock/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A Computational Approach Applied to the Study of Potential Allosteric Inhibitors Protease NS2B/NS3 from Dengue Virus

Dengue virus (DENV) is a danger to more than 400 million people in the world, and there is no specific treatment. Thus, there is an urgent need to develop an effective method to combat this pathology. NS2B/NS3 protease is an important biological target due it being necessary for viral replication and the fact that it promotes the spread of the infection. Thus, this study aimed to design DENV NS2B/NS3pro allosteric inhibitors from a matrix compound. The search was conducted using the Swiss Similarity tool. The compounds were subjected to molecular docking calculations, molecular dynamics simulations (MD) and free energy calculations. The molecular docking results showed that two compounds, ZINC000001680989 and ZINC000001679427, were promising and performed important hydrogen interactions with the Asn152, Leu149 and Ala164 residues, showing the same interactions obtained in the literature. In the MD, the results indicated that five residues, Lys74, Leu76, Asn152, Leu149 and Ala166, contribute to the stability of the ligand at the allosteric site for all of the simulated systems. Hydrophobic, electrostatic and van der Waals interactions had significant effects on binding affinity. Physicochemical properties, lipophilicity, water solubility, pharmacokinetics, druglikeness and medicinal chemistry were evaluated for four compounds that were more promising, showed negative indices for the potential penetration of the Blood Brain Barrier and expressed high human intestinal absorption, indicating a low risk of central nervous system depression or drowsiness as the the side effects. The compound ZINC000006694490 exhibited an alert with a plausible level of toxicity for the purine base chemical moiety, indicating hepatotoxicity and chromosome damage in vivo in mouse, rat and human organisms. All of the compounds selected in this study showed a synthetic accessibility (SA) score lower than 4, suggesting the ease of new syntheses. The results corroborate with other studies in the literature, and the computational approach used here can contribute to the discovery of new and potent anti-dengue agents.

Efficient 3D conformer generation of cyclic peptides formed by a disulfide bond

Cyclic peptides formed by disulfide bonds have been one large group of common drug candidates in drug development. Structural information of a peptide is essential to understand its interaction with its target. However, due to the high flexibility of peptides, it is difficult to sample the near-native conformations of a peptide. Here, we have developed an extended version of our MODPEP approach, named MODPEP2.0, to fast generate the conformations of cyclic peptides formed by a disulfide bond. MODPEP2.0 builds the three-dimensional (3D) structures of a cyclic peptide from scratch by assembling amino acids one by one onto the cyclic fragment based on the constructed rotamer and cyclic backbone libraries. Being tested on a data set of 193 diverse cyclic peptides, MODPEP2.0 obtained a considerable advantage in both accuracy and computational efficiency, compared with other sampling algorithms including PEP-FOLD, ETKDG, and modified ETKDG (mETKDG). MODPEP2.0 achieved a high sampling accuracy with an average C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\alpha$$\end{document}α RMSD of 2.20 Å and 1.66 Å when 10 and 100 conformations were considered, respectively, compared with 3.41 Å and 2.62 Å for PEP-FOLD, 3.44 Å and 3.16 Å for ETKDG, 3.09 Å and 2.72 Å for mETKDG. MODPEP2.0 also reproduced experimental peptide structures for 81.35% of the test cases when an ensemble of 100 conformations were considered, compared with 54.95%, 37.50% and 50.00% for PEP-FOLD, ETKDG, and mETKDG. MODPEP2.0 is computationally efficient and can generate 100 peptide conformations in one second. MODPEP2.0 will be useful in sampling cyclic peptide structures and modeling related protein-peptide interactions, facilitating the development of cyclic peptide drugs.

Decoding Conformational Imprint of Convoluted Molecular Interactions Between Prenylflavonoids and Aggregated Amyloid-Beta42 Peptide Causing Alzheimer's Disease

Protein misfolding occurs due to the loss of native protein structure and adopts an abnormal structure, wherein the misfolded proteins accumulate and form aggregates, which result in the formation of amyloid fibrils that are associated with neurodegenerative diseases. Amyloid beta (Aβ42) aggregation or amyloidosis is contemplated as a unique hallmark characteristic of Alzheimer’s disease (AD). Due to aberrant accrual and aggregation of Aβ42 in extracellular space, the formation of senile plaques is found in AD patients. These senile plaques occur usually in the cognitive and memory region of the brain, enfeebles neurodegeneration, hinders the signaling between synapse, and disrupts neuronal functioning. In recent years, herbal compounds are identified and characterized for their potential as Aβ42 inhibitors. Thus, understanding their structure and molecular mechanics can provide an incredible finding in AD therapeutics. To describe the structure-based molecular studies in the rational designing of drugs against amyloid fibrils, we examined various herbal compounds that belong to prenylflavonoids. The present study characterizes the trends we identified at molecular docking studies and dynamics simulation where we observed stronger binding orientation of bavachalcone, bavachin, and neobavaisoflavone with the amyloid-beta (Aβ42) fibril structure. Hence, we could postulate that these herbal compounds could be potential inhibitors of Aβ42 fibrils; these anti-aggregation agents need to be considered in treating AD.

Factual insights of the allosteric inhibition mechanism of SARS-CoV-2 main protease by quercetin: an in silico analysis

SARS-CoV-2 main protease (Mpro) cleaves the viral polypeptide 1a and 1ab in a site-specific ((L-Q|(S, A, G)) manner and produce functional enzymes for mediating viral replication. Numerous studies have reported synthetic competitive inhibitors against this target enzyme but increase in substrate concentration often reduces the effectiveness of such inhibitors. Allosteric inhibition by natural compound can provide safe and effective treatment by alleviating this limitation. Present study deals with in silico allosteric inhibition analysis of quercetin, against SARS-CoV-2-Mpro. Molecular docking of quercetin with Mpro revealed consistent binding of quercetin at a site other than active site in multiple runs, with the highest binding energy of - 8.31 kcal/mol, forming 6 H-bonds with residues Gln127, Cys128, Lys137, Asp289 and Glu290. Molecular dynamic simulation of 50 ns revealed high stability of Mpro-quercetin complex with RMSD values ranging from 0.1 to 0.25 nm. Moreover, native-Mpro and Mpro-quercetin complex conformations extracted at different time points from simulation trajectories were subjected to active site-specific docking with modelled substrate peptide (AVLQSGFR) by ZDOCK server. Results displayed site-specific cleavage of peptide when docked with native-Mpro. While substrate peptide remained intact when docked with Mpro-quercetin complex, also the binding energy of peptide reduced from 785 to 86 from 1 to 50 ns as quercetin induced alterations in the active site cavity reducing its affinity for the substrate. Further, no interactions were noticed between peptide and active site residues of Mpro-quercetin complex conformations at 40 and 50 ns. Hence, quercetin displayed effective allosteric inhibition potential against SARS-CoV-2 Mpro, and can be developed into an efficient treatment for COVID-19.

Molecular Docking and Molecular Dynamics Aided Virtual Search of OliveNet™ Directory for Secoiridoids to Combat SARS-CoV-2 Infection and Associated Hyperinflammatory Responses

Molecular docking and molecular dynamics aided virtual search of OliveNet™ directory identified potential secoiridoids that combat SARS-CoV-2 entry, replication, and associated hyperinflammatory responses. OliveNet™ is an active directory of phytochemicals obtained from different parts of the olive tree, Olea europaea (Oleaceae). Olive oil, olive fruits containing phenolics, known for their health benefits, are indispensable in the Mediterranean and Arabian diets. Secoiridoids is the largest group of olive phenols and is exclusive to the olive fruits. Functional food like olive fruits could help prevent and alleviate viral disease at an affordable cost. A systematized virtual search of 932 conformers of 78 secoiridoids utilizing Autodock Vina, followed by precision docking using Idock and Smina indicated that Nüzhenide oleoside (NZO), Oleuropein dimer (OED), and Dihydro oleuropein (DHO) blocked the SARS-CoV-2 spike (S) protein-ACE-2 interface; Demethyloleuropein (DMO), Neo-nüzhenide (NNZ), and Nüzhenide (NZE) blocked the SARS-CoV-2 main protease (Mpro). Molecular dynamics (MD) simulation of the NZO-S-protein-ACE-2 complex by Desmond revealed stability during 50 ns. RMSD of the NZO-S-protein-ACE-2 complex converged at 2.1 Å after 20 ns. During MD, the interaction fractions confirmed multiple interactions of NZO with Lys417, a crucial residue for inhibition of S protein. MD of DMO-Mpro complex proved its stability as the RMSD converged at 1.6 Å. Analysis of interactions during MD confirmed the interaction of Cys145 of Mpro with DMO and, thus, its inhibition. The docking predicted IC50 of NZO and DMO was 11.58 and 6.44 μM, respectively. Molecular docking and dynamics of inhibition of the S protein and Mpro by NZO and DMO correlated well. Docking of the six-hit secoiridoids to IL1R, IL6R, and TNFR1, the receptors of inflammatory cytokines IL1β, IL6, and TNFα, revealed the anti-inflammatory potential except for DHO. Due to intricate structures, the secoiridoids violated Lipinski's rule of five. However, the drug scores of secoiridoids supported their use as drugs. The ADMET predictions implied that the secoiridoids are non-toxic and pose low oral absorption. Secoiridoids need further optimization and are a suitable lead for the discovery of anti-SARS-CoV-2 therapeutics. For the moment, olive secoiridoids presents an accessible mode of prevention and therapy of SARS-CoV-2 infection.

Application of a Distance-Dependent Sigmoidal Dielectric Constant to the REMC/SAAP3D Simulations of Chignolin, Trp-Cage, and the G10q Mutant

The replica-exchange Monte Carlo method based on the single amino acid potential (SAAP) force field, i.e., REMC/SAAP3D, was recently developed by our group for the molecular simulation of short peptides. In this study, the method has been improved by applying a distance-dependent dielectric (DDD) constant and extended to the peptides containing D-amino acid (AA) residues. For chignolin (10 AAs), a sigmoidal DDD model reasonably allocated the native-like β-hairpin structure with all-atom root mean square deviation (RMSD) = 2.0 Å as a global energy minimum. The optimal DDD condition was subsequently applied for Trp-cage (20 AAs) and its G10q mutant. The native-like α-rich folded structures with main-chain RMSD = 3.7 and 3.8 Å were obtained as global energy minima for Trp-cage and G10q, respectively. The results suggested that the REMC/SAAP3D method with the sigmoidal DDD model is useful for structural prediction for the short peptides comprised of up to 20 AAs. In addition, the relative contributions of SAAP to the total energy (%SAAP) were evaluated by energetic component analysis. The ratios of %SAAP were about 40 and 20% for chignolin and Trp-cage (or G10q), respectively. It was proposed that SAAP is more important for the secondary structure formation than for the assembly to a higher-order folded structure, in which the attractive van der Waals interaction may play a more important role.

Influence of the conformations of αA-crystallin peptides on the isomerization rates of aspartic acid residues

The isomerization rate of aspartic acid (Asp) residue is known to be affected by the three-dimensional structures of peptides and proteins. Although the isomerized Asp residues were experimentally observed, structural features which affect the isomerization cannot be elucidated sufficiently because of protein denaturation and aggregation. In this study, molecular dynamics (MD) simulations were conducted on three αA-crystallin peptides (T6, T10, and T18), each containing a single Asp residue with different isomerization rate (T18 > T6 > T10) to clarify the structural factors of Asp isomerization tendency. For MD trajectories, distances between side-chain carboxyl carbon of Asp and main-chain amide nitrogen of (n + 1) residue (C_γ-N distances), root mean square fluctuations (RMSFs), and polar surface areas for main-chain amide nitrogen of (n + 1) residues (PSA_N) were calculated, because these structural features are considered to relate to the formations of cyclic imide intermediates. RMSFs and PSA_N are indexes of peptide backbone flexibilities and solvent exposure of the amide nitrogen, respectively. The average C_γ-N distances of T10 was longer than those of the other two peptides. In addition, the peptide containing Asp residue with a higher isomerization rate showed higher flexibility of the peptide backbone around the Asp residue. PSA_N for amide nitrogen in T18 were much larger than those of other two peptides. The computational results suggest that Asp-residue isomerization rates are affected by these factors.

Comparative modelling studies of fruit bromelain using molecular dynamics simulation

Fruit bromelain is a cysteine protease accumulated in pineapple fruits. This proteolytic enzyme has received high demand for industrial and therapeutic applications. In this study, fruit bromelain sequences QIM61759, QIM61760 and QIM61761 were retrieved from the National Center for Biotechnology Information (NCBI) Genbank Database. The tertiary structure of fruit bromelain QIM61759, QIM61760 and QIM61761 was generated by using MODELLER. The result revealed that the local stereochemical quality of the generated models was improved by using multiple templates during modelling process. Moreover, by comparing with the available papain model, structural analysis provides an insight on how pro-peptide functions as a scaffold in fruit bromelain folding and contributing to inactivation of mature protein. The structural analysis also disclosed the similarities and differences between these models. Lastly, thermal stability of fruit bromelain was studied. Molecular dynamics simulation of fruit bromelain structures at several selected temperatures demonstrated how fruit bromelain responds to elevation of temperature.

The First 3D Model of the Full-Length KIT Cytoplasmic Domain Reveals a New Look for an Old Receptor

Receptor tyrosine kinases (RTKs) are key regulators of normal cellular processes and have a critical role in the development and progression of many diseases. RTK ligand-induced stimulation leads to activation of the cytoplasmic kinase domain that controls the intracellular signalling. Although the kinase domain of RTKs has been extensively studied using X-ray analysis, the kinase insert domain (KID) and the C-terminal are partially or fully missing in all reported structures. We communicate the first structural model of the full-length RTK KIT cytoplasmic domain, a crucial target for cancer therapy. This model was achieved by integration of ab initio KID and C-terminal probe models into an X-ray structure, and by their further exploration through molecular dynamics (MD) simulation. An extended (2-µs) MD simulation of the proper model provided insight into the structure and conformational dynamics of the full-length cytoplasmic domain of KIT, which can be exploited in the description of the KIT transduction processes.

Three dimensional structures of putative, primitive proteins to investigate the origin of homochirality

Primitive proteins are likely to have been constructed from non-enzymatically generated amino acids, due to the weak enzymatic activities of primitive biomolecules such as ribozymes. On the other hand, almost all present proteins are constructed only from L-amino acids. Therefore, there must have been a mechanism early in the origins of life that selected for one of the optical isomers of amino acids. In this study, we used molecular dynamics simulations to predict the three-dimensional structures of the putative primitive proteins constructed only from glycine, alanine, aspartic acid, and valine ([GADV]-peptides). The [GADV]-peptides were generated computationally at random from L-amino acids (L-[GADV]-peptides) and from both L- and D-amino acids (DL-[GADV]-peptides). The results indicate that the tendency of secondary structure formation for L-[GADV]-peptides was larger than that for DL-[GADV]-peptides, and L-[GADV]-peptides were more rigid than DL-[GADV]-peptides. These results suggest that the proteins with rigid structure motifs were more prone to have been generated in a primordial soup that included only L-amino acids than a the soup including racemic amino acids. The tendency of the rigid structure motif formation may have played a role in selecting for the homochirality that dominates life on Earth today.

Impact of N-(2-aminoethyl) Glycine Unit on Watson-Crick Base Pairs

We aim to understand the structure and stability of the backbone tailored Watson-Crick base pairs, Guanine-Cytosine (GC), Adenine-Thymine (AT) and Adenine-Uracil (AU) by incorporating N-(2-aminoethyl) glycine units (linked by amide bonds) at the purine and pyrimidine sites of the nucleobases. Density functional theory (DFT) is employed in which B3LYP/6-311++G∗ ∗ level of theory has been used to optimize all the structures. The peptide attached base pairs are compared with the natural deoxyribose nucleic acid (DNA)/ribonucleic acid (RNA) base pairs and the calculations are carried out in both the gas and solution phases. The structural propensities of the optimized base pairs are analyzed using base pair geometries, hydrogen bond distances and stabilization energies and, compared with the standard reference data. The structural parameters were found to correlate well with the available data. The addition of peptide chain at the back bone of the DNA/RNA base pairs results only with a minimal distortion and hence does not alter the structural configuration of the base pairs. Also enhanced stability of the base pairs is spotted while adding peptidic chain at the purine site rather than the pyrimidine site of the nucleobases. The stability of the complexes is further interpreted by considering the hydrogen bonded N–H stretching frequencies of the respective base pairs. The discrimination in the interaction energies observed in both gas and solution phases are resulted due to the existence of distinct lowest unoccupied molecular orbitals (LUMO) in the solution phase. The reactivity of the base pairs is also analyzed through the in-depth examinations on the highest occupied molecular orbital (HOMO)-LUMO orbitals.

Influences of conformations of peptides on stereoinversions and/or isomerizations of aspartic acid residues

Recently, non-enzymatic stereoinversions of aspartic acid (Asp) residues in proteins and peptides have been reported. Here, we performed replica exchange molecular dynamics (REMD) simulations of model peptides (exon 6, 26A-1, and 26A-2) extracted from elastin to investigate their structural features, thereby revealing the factor that influences stereoinversions. For REMD trajectories, we calculated distances between carboxyl carbon in Asp and amide nitrogen in the (n + 1) residue (CN distances). Because bond formation between carbon and nitrogen is indispensable to the formation of a succinimide intermediate the distance between them seems to play an important role in stereoinversion. Moreover, we calculated polar surface areas (PSAs) for the trajectories, finding that CN distances and PSA were different for each peptide, with the longest CN distance and smallest PSA observed for exon 6 peptide, where stereoinversion of Asp is the slowest. Although the average CN distance was shorter for exon 26A-1 peptide than for exon 26A-2 peptide, the number of conformations with CN distances <3.0 Å was greater for exon 26A-2 peptide than for exon 26A-1 peptide. Furthermore, PSA for amide nitrogen of the (n + 1) residue was larger for exon 26A-2 peptide than for exon 26A-1 peptide. These results indicated that the flexibility of Asp and (n + 1) residues and hydrophilicity of peptides, especially in the (n + 1) residue, play important roles in the stereoinversion of Asp. This article is part of a Special Issue entitled: D-Amino acids: biology in the mirror, edited by Dr. Loredano Pollegioni, Dr. Jean-Pierre Mothet and Dr. Molla Gianluca.