Search

Isotactic polypropylene (iPP) is a frequently used polymer in industrial processing and scientific research, used as a model material. This review highlights some current and forefront issues in polymer science and engineering. This review alludes to classical deformation schemes to interpret the complicated deformation in isotactic polypropylene as a semicrystalline polymer. Also, due to the complex hierarchical structures of crystallized polymers, the progress on the multiscale and multistage microstructural evolution covering each deformation stage from the initiation of plasticity until failure is mainly concerned, particularly in terms of crystalline morphology changes. Specifically, the characterization of the microstructure evolution of a single spherulite, lamella, and chain of isotactic polypropylene will be emphasized and mentioned. Here, the authors are most interested in the results obtained from stretching (strain‐induced) orientation and crystallization studies accompanied by the in situ scattering and spectroscopy techniques offering a detection window from nano to micrometer scale, which should be able to reflect the microstructural changes at different length scales of interest. During polymer processing, the melts are frequently subjected to shear or/and elongation flow fields, producing flow‐induced molecular chain orientation in the melt. The oriented molecular chains crystallize differently than those encountered under quiescent conditions. Thus, orientation‐induced crystallization has long been an object of intense interest. A vast body of research is reported about the effects of preorientation on the crystallization of various polymers. Among them, iPP is most frequently studied since its diversified structure and morphology are very sensitive to changes in processing conditions and molecular parameters. [ABSTRACT FROM AUTHOR]

Background: Sequence matching is extremely important for applications throughout biology, particularly for discovering information such as functional and evolutionary relationships, and also for discriminating between unimportant and disease mutants. At present the functions of a large fraction of genes are unknown; improvements in sequence matching will improve gene annotations. Universal amino acid substitution matrices such as Blosum62 are used to measure sequence similarities and to identify distant homologues, regardless of the structure class. However, such single matrices do not take into account important structural information evident within the different topologies of proteins and treats substitutions within all protein folds identically. Others have suggested that the use of structural information can lead to significant improvements in sequence matching but this has not yet been very effective. Here we develop novel substitution matrices that include not only general sequence information but also have a topology specific component that is unique for each CATH topology. This novel feature of using a combination of sequence and structure information for each protein topology significantly improves the sequence matching scores for the sequence pairs tested. We have used a novel multi-structure alignment method for each homology level of CATH in order to extract topological information. Results: We obtain statistically significant improved sequence matching scores for 73 % of the alpha helical test cases. On average, 61 % of the test cases showed improvements in homology detection when structure information was incorporated into the substitution matrices. On average z-scores for homology detection are improved by more than 54 % for all cases, and some individual cases have z-scores more than twice those obtained using generic matrices. Our topology specific similarity matrices also outperform other traditional similarity matrices and single matrix based structure methods. When default amino acid substitution matrix in the Psi-blast algorithm is replaced by our structure-based matrices, the structure matching is significantly improved over conventional Psi-blast. It also outperforms results obtained for the corresponding HMM profiles generated for each topology. Conclusions: We show that by incorporating topology-specific structure information in addition to sequence information into specific amino acid substitution matrices, the sequence matching scores and homology detection are significantly improved. Our topology specific similarity matrices outperform other traditional similarity matrices, single matrix based structure methods, also show improvement over conventional Psi-blast and HMM profile based methods in sequence matching. The results support the discriminatory ability of the new amino acid similarity matrices to distinguish between distant homologs and structurally dissimilar pairs. [ABSTRACT FROM AUTHOR]

Tyrosinase is a key enzyme target to design new chemical ligands against melanogenesis. In the current review, different chemical derivatives are explored which have been used as anti-melanogenic compounds. These are different chemical compounds naturally present in plants and semi-synthetic and synthetic compounds inspired by these natural products, such as kojic acid produced by several species of fungi; arbutin—a glycosylated hydroquinone extracted from the bearberry plant; vanillin—a phenolic aldehyde extracted from the vanilla bean, etc. After enzyme inhibition screening, various chemical compounds showed different therapeutic effects as tyrosinase inhibitors with different values of the inhibition constant and IC50. We show how appropriately designed scaffolds inspired by the structures of natural compounds are used to develop novel synthetic inhibitors. We review the results of numerous studies, which could lead to the development of effective anti-tyrosinase agents with increased efficiency and safety in the near future, with many applications in the food, pharmaceutical and cosmetics industries. [ABSTRACT FROM AUTHOR]

We studied the refolding free energy landscape of 26 proteins using the Go-like model. The distance between the denaturated state and the transition state, XF, was calculated using the Bell theory and the nonlinear Dudko-Hummer-Szabo theory, and its relation to the geometrical properties of the native state was considered in detail. We showed that none of the structural parameters, such as the contact order, protein length, and radius of cross section, correlate with XF for all classes of proteins. To overcome this problem, we have introduced the nematic order parameter P02, which describes the ordering of the structured elements of the native state. Due to its topologically global nature, P02 is better than other structural parameters in describing the folding free energy landscape. In particular, P02 displays a good correlation with XF extracted from the nonlinear theory for all three classes of proteins. Therefore, this parameter can be used to predict XF for any protein, if its native structure is known. [ABSTRACT FROM AUTHOR]

Drug designing is high-priced and time taking process with low success rate. To overcome this obligation, computational drug repositioning technique is being promptly used to predict the possible therapeutic effects of FDA approved drugs against multiple diseases. In this computational study, protein modeling, shape-based screening, molecular docking, pharmacogenomics, and molecular dynamic simulation approaches have been utilized to retrieve the FDA approved drugs against AD. The predicted MADD protein structure was designed by homology modeling and characterized through different computational resources. Donepezil and galantamine were implanted as standard drugs and drugs were screened out based on structural similarities. Furthermore, these drugs were evaluated and based on binding energy (Kcal/mol) profiles against MADD through PyRx tool. Moreover, pharmacogenomics analysis showed good possible associations with AD mediated genes and confirmed through detail literature survey. The best 6 drug (darifenacin, astemizole, tubocurarine, elacridar, sertindole and tariquidar) further docked and analyzed their interaction behavior through hydrogen binding. Finally, MD simulation study were carried out on these drugs and evaluated their stability behavior by generating root mean square deviation and fluctuations (RMSD/F), radius of gyration (Rg) and soluble accessible surface area (SASA) graphs. Taken together, darifenacin, astemizole, tubocurarine, elacridar, sertindole and tariquidar displayed good lead like profile as compared with standard and can be used as possible therapeutic agent in the treatment of AD after in-vitro and in-vivo assessment. [ABSTRACT FROM AUTHOR]

Mushrooms are new potential sources of valuable medicines, long neglected because of difficulties experienced in their cultivation. There is a large variety of medicinal mushrooms which possess significant therapeutic properties and are used as medications for various diseases because they contain several novel highly bioactive components. Medicinal mushrooms can be identified based on their morphology, size, mass, and the color of the stalk, cap and spore, and attachment to the stalk. Medicinal mushrooms possess a variety of important biological activities and are used as antioxidants, hepatoprotectors, anticancer, antidiabetic, anti-inflammatory, antiaging, antiviral, antiparasitic, and antimicrobial agents, among others. This review provides a basic overview of the chemical scaffolds present in mushrooms and their therapeutic implications in the human body. [ABSTRACT FROM AUTHOR]

Charcot–Marie–Tooth disease (CMT) is a heritable neurodegenerative disease of peripheral nervous system diseases in which more than 100 genes and their mutations are associated. Two consanguineous families Dera Ghazi Khan (PAK-CMT1-DG KHAN) and Layyah (PAK-CMT2-LAYYAH) with multiple CMT-affected subjects were enrolled from Punjab province in Pakistan. Basic epidemiological data were collected for the subjects. Nerve conduction study (NCS) and electromyography (EMG) were performed for the patients. Whole-exome sequencing (WES) followed by Sanger sequencing was applied to report the genetic basic of CMT. The NCS findings revealed that sensory and motor nerve conduction velocities for both families were <38 m/s. EMG presented denervation, neuropathic motor unit potential, and reduced interference pattern of peripheral nerves. WES identified that a novel nonsense mutation (c. 226 G>T) in GADP1 gene and a previously known missense mutation in MFN2 gene (c. 334 G>A) cause CMT4A (Charcot–Marie–Tooth disease type 4A) in the PAK-CMT1-DG KHAN family and CMT2A (Charcot–Marie–Tooth disease type 2A) in the PAK-CMT2-LAYYAH family, respectively. Mutations followed Mendelian pattern with autosomal recessive mode of inheritance. Multiple sequence alignment by Clustal Omega indicated that mutation-containing domain in both genes is highly conserved, and in situ analysis revealed that both mutations are likely to be pathogenic. We reported that a novel nonsense mutation and a previously known missense mutation in GAPD1 gene and MFN2 gene, respectively, cause CMT in consanguineous Pakistani families. [ABSTRACT FROM AUTHOR]

Drug designing is high-priced and time taking process with low success rate. To overcome this obligation, computational drug repositioning technique is being promptly used to predict the possible therapeutic effects of FDA approved drugs against multiple diseases. In this computational study, protein modeling, shape-based screening, molecular docking, pharmacogenomics, and molecular dynamic simulation approaches have been utilized to retrieve the FDA approved drugs against AD. The predicted MADD protein structure was designed by homology modeling and characterized through different computational resources. Donepezil and galantamine were implanted as standard drugs and drugs were screened out based on structural similarities. Furthermore, these drugs were evaluated and based on binding energy (Kcal/mol) profiles against MADD through PyRx tool. Moreover, pharmacogenomics analysis showed good possible associations with AD mediated genes and confirmed through detail literature survey. The best 6 drug (darifenacin, astemizole, tubocurarine, elacridar, sertindole and tariquidar) further docked and analyzed their interaction behavior through hydrogen binding. Finally, MD simulation study were carried out on these drugs and evaluated their stability behavior by generating root mean square deviation and fluctuations (RMSD/F), radius of gyration (Rg) and soluble accessible surface area (SASA) graphs. Taken together, darifenacin, astemizole, tubocurarine, elacridar, sertindole and tariquidar displayed good lead like profile as compared with standard and can be used as possible therapeutic agent in the treatment of AD after in-vitro and in-vivo assessment. [ABSTRACT FROM AUTHOR]

Ewing's sarcoma is a rare type of cancer that forms in bones and soft tissues in the body, affecting mostly children and young adults. Current treatments for ES are limited to chemotherapy and/or radiation, followed by surgery. Recently, microRNAs have shown favourable results as latent diagnostic and prognostic biomarkers in various cancers. Furthermore, microRNAs have shown to be a good therapeutic agent due to their involvement in the dysregulation of various molecular pathways linked to tumour progression, invasion, angiogenesis, and metastasis. In this review, comprehensive data mining was employed to explore various microRNAs that might have therapeutic potential as target molecules in the treatment of ES. [ABSTRACT FROM AUTHOR]

In this study, a series of 1,3,4-thiadiazole-thiazolidinone derivatives (7a–j) were synthesized as carbonic anhydrase inhibitors. The in vitro carbonic anhydrase inhibitory activity of these synthesized compounds was determined and the results revealed that compound 7i (IC50 = 0.402 ± 0.017 μM) is more potent than the standard reference acetazolamide (IC50 = 0.998 ± 0.046 μM). Moreover, kinetic analysis was also performed to see the mode of binding with the target enzyme, and it was shown that compound 7i bound with the target in a competitive manner. All compounds were screened for antioxidant activity against DPPH and anticancer activity against A549 (human lung cancer cells), and the MTT assay results indicated that compounds 7a and 7d had more inhibitory growth. Moreover, compound 7i exhibited a good conformational state with hydrogen bond interactions within the receptor binding pocket according to molecular docking studies. Therefore, the novel compound 7i might be employed as a promising pharmacophore for potent and selective carbonic anhydrase inhibitors. [ABSTRACT FROM AUTHOR]

Dosage-sensitive sex reversal, adrenal hypoplasia critical region, on chromosome X, gene 1 (DAX1) is an orphan nuclear receptor encoded by the NR0B1 gene. The functional study showed that DAX1 is a physiologically significant target for EWS/FLI1-mediated oncogenesis, particularly Ewing Sarcoma (ES). In this study, a three-dimensional DAX1 structure was modeled by employing a homology modeling approach. Furthermore, the network analysis of genes involved in Ewing Sarcoma was also carried out to evaluate the association of DAX1 and other genes with ES. Moreover, a molecular docking study was carried out to check the binding profile of screened flavonoid compounds against DAX1. Therefore, 132 flavonoids were docked in the predicted active binding pocket of DAX1. Moreover, the pharmacogenomics analysis was performed for the top ten docked compounds to evaluate the ES-related gene clusters. As a result, the five best flavonoid-docked complexes were selected and further evaluated by Molecular Dynamics (MD) simulation studies at 100 ns. The MD simulation trajectories were evaluated by generating RMSD, hydrogen bond plot analysis, and interaction energy graphs. Our results demonstrate that flavonoids showed interactive profiles in the active region of DAX1 and can be used as potential therapeutic agents against DAX1-mediated augmentation of ES after in-vitro and in-vivo evaluations. [ABSTRACT FROM AUTHOR]

A well‐known key enzyme in melanogenesis and hyperpigmentation is tyrosinase. The present study introduces a novel series of thiophenyl‐pyrazolylthiazole‐coumarin hybrids (6a–6h) as tyrosinase inhibitors. The in‐vitro tyrosinase inhibition results indicated that all compounds have strong tyrosinase inhibitory activity, particularly compound 6g (IC50 = 0.043 ± 0.006 μM), was identified as the most active compound compared to the positive control (kojic acid, IC50 = 18.521 ± 1.162 μM). Lineweaver‐Burk plots were employed to analyze the kinetic mechanism, and compound 6g formed an enzyme‐inhibitor complex by inhibiting tyrosinase non‐competitively. Furthermore, all compounds demonstrated excellent antioxidant activity against DPPH. MTT assay was used to screen the cytotoxicity of all compounds on B16F10 melanoma cells, and they had no toxic effect on the cells. The binding affinity of compounds with tyrosinase was also investigated using molecular docking, and the ligands displayed good binding energy values. These molecules could be a promising lead for skin pigmentation and associated diseases as nontoxic pharmacological scaffolds. [ABSTRACT FROM AUTHOR]

The COVID-19 outbreak was a major event that greatly impacted the economy and the health systems around the world. Understanding the behavior of the virus and being able to perform long-term and short-term future predictions of the daily new cases is a working field for machine learning methods and mathematical models. This paper compares Verhulst's, Gompertz´s, and SIR models from the point of view of their efficiency to describe the behavior of COVID-19 in Spain. These mathematical models are used to predict the future of the pandemic by first solving the corresponding inverse problems to identify the model parameters in each wave separately, using as observed data the daily cases in the past. The posterior distributions of the model parameters are then inferred via the Metropolis–Hastings algorithm, comparing the robustness of each prediction model and making different representations to visualize the results obtained concerning the posterior distribution of the model parameters and their predictions. The knowledge acquired is used to perform predictions about the evolution of both the daily number of infected cases and the total number of cases during each wave. As a main conclusion, predictive models are incomplete without a corresponding uncertainty analysis of the corresponding inverse problem. The invariance of the output (posterior prediction) with respect to the forward predictive model that is used shows that the methodology shown in this paper can be used to adopt decisions in real practice (public health). [ABSTRACT FROM AUTHOR]

Background: Protein secondary structure prediction (SSP) has been an area of intense research interest. Despite advances in recent methods conducted on large datasets, the estimated upper limit accuracy is yet to be reached. Since the predictions of SSP methods are applied as input to higher-level structure prediction pipelines, even small errors may have large perturbations in final models. Previous works relied on cross validation as an estimate of classifier accuracy. However, training on large numbers of protein chains compromises the classifier ability to generalize to new sequences. This prompts a novel approach to training and an investigation into the possible structural factors that lead to poor predictions. Here, a small group of 55 proteins termed the compact model is selected from the CB513 dataset using a heuristics-based approach. In a prior work, all sequences were represented as probability matrices of residues adopting each of Helix, Sheet and Coil states, based on energy calculations using the C-Alpha, C-Beta, Side-chain (CABS) algorithm. The functional relationship between the conformational energies computed with CABS force-field and residue states is approximated using a classifier termed the Fully Complex-valued Relaxation Network (FCRN). The FCRN is trained with the compact model proteins. Results: The performance of the compact model is compared with traditional cross-validated accuracies and blind-tested on a dataset of G Switch proteins, obtaining accuracies of ∼81 %. The model demonstrates better results when compared to several techniques in the literature. A comparative case study of the worst performing chain identifies hydrogen bond contacts that lead to Coil↔Sheet misclassifications. Overall, mispredicted Coil residues have a higher propensity to participate in backbone hydrogen bonding than correctly predicted Coils. Conclusions: The implications of these findings are: (i) the choice of training proteins is important in preserving the generalization of a classifier to predict new sequences accurately and (ii) SSP techniques sensitive in distinguishing between backbone hydrogen bonding and side-chain or water-mediated hydrogen bonding might be needed in the reduction of Coil↔Sheet misclassifications. [ABSTRACT FROM AUTHOR]

Ewing sarcoma (ES) is a highly malignant carcinoma prevalent in children and most frequent in the second decade of life. It mostly occurs due to t(11;22) (q24;q12) translocation. This translocation encodes the oncogenic fusion protein EWS/FLI (Friend leukemia integration 1 transcription factor), which acts as an aberrant transcription factor to deregulate target genes essential for cancer. Traditionally, flavonoids from plants have been investigated against viral and cancerous diseases and have shown some promising results to combat these disorders. In the current study, representative flavonoid compounds from various subclasses are selected and used to disrupt the RNA-binding motif of EWS, which is required for EWS/FLI fusion. By blocking the RNA-binding motif of EWS, it might be possible to combat ES. Therefore, molecular docking experiments validated the binding interaction patterns and structural behaviors of screened flavonoid compounds within the active region of the Ewing sarcoma protein (EWS). Furthermore, pharmacogenomics analysis was used to investigate potential drug interactions with Ewing sarcoma-associated genes. Finally, molecular dynamics simulations were used to investigate the stability of the best selected docked complexes. Taken together, daidzein, kaempferol, and genistein exhibited a result comparable to ifosfamide in the proposed in silico study and can be further analyzed as possible candidate compounds in biological in vitro studies against ES. [ABSTRACT FROM AUTHOR]

MicroRNAs are hidden players in complex psychophysical phenomena such as depression and anxiety related disorders though the activation and deactivation of multiple proteins in signaling cascades. Depression is classified as a mood disorder and described as feelings of sadness, loss, or anger that interfere with a person's everyday activities. In this review, we have focused on exploration of the significant role of miRNAs in depression by affecting associated target proteins (cellular and synaptic) and their signaling pathways which can be controlled by the attachment of miRNAs at transcriptional and translational levels. Moreover, miRNAs have potential role as biomarkers and may help to cure depression through involvement and interactions with multiple pharmacological and physiological therapies. Taken together, miRNAs might be considered as promising novel therapy targets themselves and may interfere with currently available antidepressant treatments. [ABSTRACT FROM AUTHOR]

To develop new alkaline phosphatase inhibitors (ALP), a series of pyrazolo-oxothiazolidine derivatives were synthesized and biologically assessed, and the results showed that all of the synthesized compounds significantly inhibited ALP. Specifically, compound 7g displayed the strongest inhibitory activity (IC50 = 0.045 ± 0.004 μM), which is 116-fold more active than monopotassium phosphate (IC50 = 5.242 ± 0.472 μM) as a standard reference. The most potent compound among the series (7g) was checked for its mode of binding with the enzyme and shown as non-competitively binding with the target enzyme. The antioxidant activity of these compounds was examined to investigate the radical scavenging effect. Moreover, the MTT assay method was performed to evaluate their toxic effects on the viability of MG-63 human osteosarcoma cells, and all compounds have no toxic effect on the cells at 4 μM. Computational research was also conducted to examine the binding affinity of the ligands with alkaline phosphatase, and the results revealed that all compounds showed good binding energy values within the active site of the target. Therefore, these novel pyrazolo-oxothiazolidine derivatives might be employed as promising pharmacophores for potent and selective alkaline phosphatase inhibitors. [ABSTRACT FROM AUTHOR]

Noise is a basic ingredient in data, since observed data are always contaminated by unwanted deviations, i.e., noise, which, in the case of overdetermined systems (with more data than model parameters), cause the corresponding linear system of equations to have an imperfect solution. In addition, in the case of highly underdetermined parameterization, noise can be absorbed by the model, generating spurious solutions. This is a very undesirable situation that might lead to incorrect conclusions. We presented mathematical formalism based on the inverse problem theory combined with artificial intelligence methodologies to perform an enhanced sampling of noisy biomedical data to improve the finding of meaningful solutions. Random sampling methods fail for high-dimensional biomedical problems. Sampling methods such as smart model parameterizations, forward surrogates, and parallel computing are better suited for such problems. We applied these methods to several important biomedical problems, such as phenotype prediction and a problem related to predicting the effects of protein mutations, i.e., if a given single residue mutation is neutral or deleterious, causing a disease. We also applied these methods to de novo drug discovery and drug repositioning (repurposing) through the enhanced exploration of huge chemical space. The purpose of these novel methods that address the problem of noise and uncertainty in biomedical data is to find new therapeutic solutions, perform drug repurposing, and accelerate and optimize drug discovery, thus reestablishing homeostasis. Finding the right target, the right compound, and the right patient are the three bottlenecks to running successful clinical trials from the correct analysis of preclinical models. Artificial intelligence can provide a solution to these problems, considering that the character of the data restricts the quality of the prediction, as in any modeling procedure in data analysis. The use of simple and plain methodologies is crucial to tackling these important and challenging problems, particularly drug repositioning/repurposing in rare diseases. [ABSTRACT FROM AUTHOR]

Highly designable structures can be distinguished based on certain geometric graphical features of the interactions, confirming the fact that the topology of a protein structure and its residue-residue interaction network are important determinants of its designability. The most designable structures and least designable structures obtained for sets of proteins having the same number of residues are compared. It is shown that the most designable structures predicted by the graph features of the contact diagrams are more densely packed, whereas the poorly designable structures are more open structures or structures that are loosely packed. Interestingly enough, it can also be seen that the highly designable identified are also common structural motifs found in nature. [ABSTRACT FROM AUTHOR]

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer’s and Parkinson’s diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aβ40 and Aβ42 and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics. [ABSTRACT FROM AUTHOR]

Hematopoietic stem cells (HSCs) reside in a specialized microenvironment in a peculiar anatomic location which regulates the maintenance of stem cells and controls its functions. Recent scientific progress in experimental technologies have enabled the specific detection of epigenetic factors responsible for the maintenance and quiescence of the hematopoietic niche, which has improved our knowledge of regulatory mechanisms. The aberrant role of RNA-binding proteins and their impact on the disruption of stem cell biology have been reported by a number of recent studies. Despite recent modernization in hematopoietic microenvironment research avenues, our comprehension of the signaling mechanisms and interactive pathways responsible for integration of the hematopoietic niche is still limited. In the past few decades, zebrafish usage with regards to exploratory studies of the hematopoietic niche has expanded our knowledge for deeper understanding of novel cellular interactions. This review provides an update on the functional roles of different genetic and epigenetic factors and molecular signaling events at different sections of the hematopoietic microenvironment. The explorations of different molecular approaches and interventions of latest web-based tools being used are also outlined. This will help us to get more mechanistic insights and develop therapeutic options for the malignancies. [ABSTRACT FROM AUTHOR]

Big data in health care is a fast-growing field and a new paradigm that is transforming case-based studies to large-scale, data-driven research. As big data is dependent on the advancement of new data standards, technology, and relevant research, the future development of big data applications holds foreseeable promise in the modern day health care revolution. Enormously large, rapidly growing collections of biomedical omics-data (genomics, proteomics, transcriptomics, metabolomics, glycomics, etc.) and clinical data create major challenges and opportunities for their analysis and interpretation and open new computational gateways to address these issues. The design of new robust algorithms that are most suitable to properly analyze this big data by taking into account individual variability in genes has enabled the creation of precision (personalized) medicine. We reviewed and highlighted the significance of big data analytics for personalized medicine and health care by focusing mostly on machine learning perspectives on personalized medicine, genomic data models with respect to personalized medicine, the application of data mining algorithms for personalized medicine as well as the challenges we are facing right now in big data analytics. [ABSTRACT FROM AUTHOR]

A series of nine novel 1,2,4-triazole based compounds were synthesized through a multistep reaction pathway and their structures were scrutinized by using spectral methods such as FTIR, LC-MS, 1H NMR, and 13C NMR. The synthesized derivatives were screened for inhibitory activity against the mushroom tyrosinase and we found that all the synthesized compounds demonstrated decent inhibitory activity against tyrosinase. However, among the series of compounds, N-(4- fluorophenyl)-2-(5-(2-fluorophenyl)-4-(4-fluorophenyl)-4H-1,2,4-triazol-3-ylthio) acetamide exhibited more prominent activity when accompanied with the standard drug kojic acid. Furthermore, the molecular docking studies identified the interaction profile of all synthesized derivatives at the active site of tyrosinase. Based on these results, N-(4-fluorophenyl)-2-(5-(2-fluorophenyl)-4-(4-fluorophenyl)- 4H-1,2,4-triazol-3-ylthio) acetamide could be used as a novel scaffold to design some new drugs against melanogenesis. [ABSTRACT FROM AUTHOR]

The prediction of the dynamics of the COVID-19 outbreak and the corresponding needs of the health care system (COVID-19 patients' admissions, the number of critically ill patients, need for intensive care units, etc.) is based on the combination of a limited growth model (Verhulst model) and a short-term predictive model that allows predictions to be made for the following day. In both cases, the uncertainty analysis of the prediction is performed, i.e., the set of equivalent models that adjust the historical data with the same accuracy. This set of models provides the posterior distribution of the parameters of the predictive model that adjusts the historical series. It can be extrapolated to the same analyzed time series (e.g., the number of infected individuals per day) or to another time series of interest to which it is correlated and used, e.g., to predict the number of patients admitted to urgent care units, the number of critically ill patients, or the total number of admissions, which are directly related to health needs. These models can be regionalized, that is, the predictions can be made at the local level if data are disaggregated. We show that the Verhulst and the Gompertz models provide similar results and can be also used to monitor and predict new outbreaks. However, the Verhulst model seems to be easier to interpret and to use. [ABSTRACT FROM AUTHOR]

Proteins are tiny players involved in the activation and deactivation of multiple signaling cascades through interactions in cells. The TNFR1 and MADD interact with each other and mediate downstream protein signaling pathways which cause neuronal cell death and Alzheimer's disease. In the current study, a molecular docking approach was employed to explore the interactive behavior of TNFR1 and MADD proteins and their role in the activation of downstream signaling pathways. The computational sequential and structural conformational results revealed that Asp400, Arg58, Arg59 were common residues of TNFR1 and MADD which are involved in the activation of downstream signaling pathways. Aspartic acid in negatively charged residues is involved in the biosynthesis of protein. However, arginine is a positively charged residue with the potential to interact with oppositely charged amino acids. Furthermore, our molecular dynamic simulation results also ensured the stability of the backbone of TNFR1 and MADD death domains (DDs) in binding interactions. This DDs interaction mediates some conformational changes in TNFR1 which leads to the activation of mediators proteins in the cellular signaling pathways. Taken together, a better understanding of TNFR1 and MADD receptors and their activated signaling cascade may help treat Alzheimer's disease. The death domains of TNFR1 and MADD could be used as a novel pharmacological target for the treatment of Alzheimer's disease by inhibiting the MAPK pathway. [ABSTRACT FROM AUTHOR]

We discuss the use of the regularized linear discriminant analysis (LDA) as a model reduction technique combined with particle swarm optimization (PSO) in protein tertiary structure prediction, followed by structure refinement based on singular value decomposition (SVD) and PSO. The algorithm presented in this paper corresponds to the category of template-based modeling. The algorithm performs a preselection of protein templates before constructing a lower dimensional subspace via a regularized LDA. The protein coordinates in the reduced spaced are sampled using a highly explorative optimization algorithm, regressive–regressive PSO (RR-PSO). The obtained structure is then projected onto a reduced space via singular value decomposition and further optimized via RR-PSO to carry out a structure refinement. The final structures are similar to those predicted by best structure prediction tools, such as Rossetta and Zhang servers. The main advantage of our methodology is that alleviates the ill-posed character of protein structure prediction problems related to high dimensional optimization. It is also capable of sampling a wide range of conformational space due to the application of a regularized linear discriminant analysis, which allows us to expand the differences over a reduced basis set. [ABSTRACT FROM AUTHOR]

Accurate prediction of protein stability changes resulting from amino acid substitutions is of utmost importance in medicine to better understand which mutations are deleterious, leading to diseases, and which are neutral. Since conducting wet lab experiments to get a better understanding of protein mutations is costly and time consuming, and because of huge number of possible mutations the need of computational methods that could accurately predict effects of amino acid mutations is of greatest importance. In this research, we present a robust methodology to predict the energy changes of a proteins upon mutations. The proposed prediction scheme is based on two step algorithm that is a Holdout Random Sampler followed by a neural network model for regression. The Holdout Random Sampler is utilized to analysis the energy change, the corresponding uncertainty, and to obtain a set of admissible energy changes, expressed as a cumulative distribution function. These values are further utilized to train a simple neural network model that can predict the energy changes. Results were blindly tested (validated) against experimental energy changes, giving Pearson correlation coeffcients of 0.66 for Single Point Mutations and 0.77 for Multiple Point Mutations. These results confirm the successfulness of our method, since it outperforms majority of previous studies in this field. [ABSTRACT FROM AUTHOR]

• An efficient one-pot synthesis of dihydropyrimidinones was performed with cysteine as a bio-organic catalyst. • These compounds (4a – 4j) displayed significant in vitro urease inhibitory activity. • MTT assay was used to evaluate the cytotoxicity of all compounds in HT1080 human fibroblast cell lines. • The kinetic mechanism was investigated by Lineweaver-Burk plots. • In silico molecular docking studies supported the experimental results. The aim of the present research was to develop a one-pot, three-component method for the synthesis of some 3,4-dihydro pyrimidine-2-(1 H)-one (DHPM) derivatives using aldehydes, urea, and ethyl acetoacetate in the presence of cysteine as a green, bio-organic catalyst. This procedure has several key advantages, including a facile and convenient method, green and low-cost catalyst, high yield, short reaction time, and avoidance of the use of environmentally harmful solvents. It is important to note that the synthesized compounds (4a – j) were evaluated for in vitro urease inhibitory activity, and all were found to be more effective than the standard thiourea (IC 50 = 4.745 ± 0.054 µM), with urease inhibitory potentials between IC 50 = 0.058 ± 0.006 µM to 0.297 ± 0.066 µM. Lineweaver-Burk plots were used to study the kinetics of the most potent compound (4e), and the results showed that the compound non-competitively inhibited urease by forming an enzyme-inhibitor complex. In addition, HT1080 human fibroblast cell lines were used to test the toxicity of the synthesized compounds on cell viability using the MTT assay method, and all compounds except 4g, 4i, and 4j exhibited no significant toxic effects on cells even at high concentrations. Moreover, molecular docking was used to determine the binding interactions between the molecules (ligands) and the active site of the urease enzyme, and the ligands indicated acceptable binding energy values. [Display omitted] [ABSTRACT FROM AUTHOR]

[Display omitted] • A series of 1,3,4-thiadiazol-1,3,4-oxadiazol-acetamide derivatives (7a – 7j) were synthesized as elastase inhibitors. • The in-vitro elastase inhibitory activity of all synthesized compounds was determined successfully. • The kinetic mechanism was investigated by Lineweaver-Burk plots. • MTT assay was used to evaluate the cytotoxicity of all compounds on B16F10 melanoma cell lines. • Molecular docking was applied to investigate the binding affinity of compounds with elastase and the ligands showed good binding energy values. Considering the biological significance of 1,3,4-thiadiazole/oxadiazole heterocyclic scaffolds, a novel series of 1,3,4-thiadiazole-1,3,4-oxadiazole-acetamide derivatives (7a – j) was designed and synthesized using molecular hybridization. The inhibitory effects of the target compounds on elastase were evaluated, and all of these molecules were found to be potent inhibitors compared to the standard reference oleanolic acid. Compound 7f exhibited the excellent inhibitory activity (IC 50 = 0.06 ± 0.02 μM), which is 214-fold more active than oleanolic acid (IC 50 = 12.84 ± 0.45 μM). Kinetic analysis was also performed on the most potent compound (7f) to determine the mode of binding with the target enzyme, and it was discovered that 7f inhibits the enzyme in a competitive manner. Furthermore, the MTT assay method was used to assess their toxicity on the viability of B16F10 melanoma cell lines, and all compounds did not display any toxic effect on the cells even at high concentrations. The molecular docking studies of all compounds also justified with their good docking score and among them, compound 7f had a good conformational state with hydrogen bond interactions within the receptor binding pocket, which is consistent with the experimental inhibition studies. [ABSTRACT FROM AUTHOR]

We discuss applicability of principal component analysis (PCA) for protein tertiary structure prediction from amino acid sequence. The algorithm presented in this paper belongs to the category of protein refinement models and involves establishing a low-dimensional space where the sampling (and optimization) is carried out via particle swarm optimizer (PSO). The reduced space is found via PCA performed for a set of low-energy protein models previously found using different optimization techniques. A high frequency term is added into this expansion by projecting the best decoy into the PCA basis set and calculating the residual model. This term is aimed at providing high frequency details in the energy optimization. The goal of this research is to analyze how the dimensionality reduction affects the prediction capability of the PSO procedure. For that purpose, different proteins from the Critical Assessment of Techniques for Protein Structure Prediction experiments were modeled. In all the cases, both the energy of the best decoy and the distance to the native structure have decreased. Our analysis also shows how the predicted backbone structure of native conformation and of alternative low energy states varies with respect to the PCA dimensionality. Generally speaking, the reconstruction can be successfully achieved with 10 principal components and the high frequency term. We also provide a computational analysis of protein energy landscape for the inverse problem of reconstructing structure from the reduced number of principal components, showing that the dimensionality reduction alleviates the ill-posed character of this high-dimensional energy optimization problem. The procedure explained in this paper is very fast and allows testing different PCA expansions. Our results show that PSO improves the energy of the best decoy used in the PCA when the adequate number of PCA terms is considered. [ABSTRACT FROM AUTHOR]

Background: Over the past 50,000 years, shifts in human-environmental or human-human interactions shaped genetic differences within and among human populations, including variants under positive selection. Shaped by environmental factors, such variants influence the genetics of modern health, disease, and treatment outcome. Because evolutionary processes tend to act on gene regulation, we test whether regulatory variants are under positive selection. We introduce a new approach to enhance detection of genetic markers undergoing positive selection, using conditional entropy to capture recent local selection signals. Results We use conditional logistic regression to compare our Adjusted Haplotype Conditional Entropy (H|H) measure of positive selection to existing positive selection measures. H|H and existing measures were applied to published regulatory variants acting in cis (cis-eQTLs), with conditional logistic regression testing whether regulatory variants undergo stronger positive selection than the surrounding gene. These cis-eQTLs were drawn from six independent studies of genotype and RNA expression. The conditional logistic regression shows that, overall, H|H is substantially more powerful than existing positive-selection methods in identifying cis-eQTLs against other Single Nucleotide Polymorphisms (SNPs) in the same genes. When broken down by Gene Ontology, H|H predictions are particularly strong in some biological process categories, where regulatory variants are under strong positive selection compared to the bulk of the gene, distinct from those GO categories under overall positive selection. . However, cis-eQTLs in a second group of genes lack positive selection signatures detectable by H|H, consistent with ancient short haplotypes compared to the surrounding gene (for example, in innate immunity GO:0042742); under such other modes of selection, H|H would not be expected to be a strong predictor.. These conditional logistic regression models are adjusted for Minor allele frequency(MAF); otherwise, ascertainment bias is a huge factor in all eQTL data sets. Relationships between Gene Ontology categories, positive selection and eQTL specificity were replicated with H|H in a single larger data set. Our measure, Adjusted Haplotype Conditional Entropy (H|H), was essential in generating all of the results above because it: 1) is a stronger overall predictor for eQTLs than comparable existing approaches, and 2) shows low sequential auto-correlation, overcoming problems with convergence of these conditional regression statistical models. Conclusions: Our new method, H|H, provides a consistently more robust signal associated with cis-eQTLs compared to existing methods. We interpret this to indicate that some cis-eQTLs are under positive selection compared to their surrounding genes. Conditional entropy indicative of a selective sweep is an especially strong predictor of eQTLs for genes in several biological processes of medical interest. Where conditional entropy is a weak or negative predictor of eQTLs, such as innate immune genes, this would be consistent with balancing selection acting on such eQTLs over long time periods. Different measures of selection may be needed for variant prioritization under other modes of evolutionary selection. [ABSTRACT FROM AUTHOR]

Background: Traditional cancer treatments have centered on cytotoxic drugs and general purpose chemotherapy that may not be tailored to treat specific cancers. Identification of molecular markers that are related to different types of cancers might lead to discovery of drugs that are patient and disease specific. This study aims to use microarray gene expression cancer data to identify biomarkers that are indicative of different types of cancers. Our aim is to provide a multi-class cancer classifier that can simultaneously differentiate between cancers and identify type-specific biomarkers, through the application of the Binary Coded Genetic Algorithm (BCGA) and a neural network based Extreme Learning Machine (ELM) algorithm. Results: BCGA and ELM are combined and used to select a subset of genes that are present in the Global Cancer Mapping (GCM) data set. This set of candidate genes contains over 52 biomarkers that are related to multiple cancers, according to the literature. They include APOA1, VEGFC, YWHAZ, B2M, EIF2S1, CCR9 and many other genes that have been associated with the hallmarks of cancer. BCGA-ELM is tested on several cancer data sets and the results are compared to other classification methods. BCGA-ELM compares or exceeds other algorithms in terms of accuracy. We were also able to show that over 50% of genes selected by BCGA-ELM on GCM data are cancer related biomarkers. Conclusions: We were able to simultaneously differentiate between 14 different types of cancers, using only 92 genes, to achieve a multi-class classification accuracy of 95.4% which is between 21.6% and 38% higher than other results in the literature for multi-class cancer classification. Our findings suggest that computational algorithms such as BCGA-ELM can facilitate biomarker-driven integrated cancer research that can lead to a detailed understanding of the complexities of cancer. [ABSTRACT FROM AUTHOR]

Vascular smooth muscle α2C-adrenoceptors (α2C-ARs) mediate vasoconstriction of small blood vessels, especially arterioles. Studies of endogenous receptors in human arteriolar smooth muscle cells (referred to as microVSM) and transiently transfected receptors in heterologous HEK293 cells show that the α2C-ARs are perinuclear receptors that translocate to the cell surface under cellular stress and elicit a biological response. Recent studies in microVSM unraveled a crucial role of Rap1A-Rho-ROCK-F-actin pathways in receptor translocation, and identified protein-protein interaction of α2C-ARs with the actin binding protein filamin-2 as an essential step in the process. To better understand the molecular nature and specificity of this interaction, in this study, we constructed comparative models of human α2C-AR and human filamin-2 proteins. Finally, we performed in silico protein-protein docking to provide a structural platform for the investigation of human α2C-AR and filamin-2 interactions. We found that electrostatic interactions seem to play a key role in this complex formation which manifests in interactions between the C-terminal arginines of α2C-ARs (particularly R454 and R456) and negatively charged residues from filamin-2 region between residues 1979 and 2206. Phylogenetic and sequence analysis showed that these interactions have evolved in warm-blooded animals. [ABSTRACT FROM AUTHOR]

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications. [ABSTRACT FROM AUTHOR]

HIV-2 infection is frequently neglected in HIV/AIDS campaigns. However, a special emphasis must be given to HIV-2 as an untreated infection that also leads to AIDS and death, and for which the efficacy of most available drugs is limited against HIV-2. HIV envelope glycoproteins mediate binding to the receptor CD4 and co-receptors at the surface of the target cell, enabling fusion with the cell membrane and viral entry. Here, we developed and optimized a computer-assisted drug design approach of an important HIV-2 glycoprotein that allows us to explore and gain further insights at the molecular level into protein structures and interactions crucial for the inhibition of HIV-2 cell entry. The 3D structure of a key HIV-2ROD gp125 region was generated by a homology modeling campaign. To disclose the importance of the main structural features and compare them with experimental results, 3D-models of six mutants were also generated. These mutations revealed the selective impact on the behavior of the protein. Furthermore, molecular dynamics simulations were performed to optimize the models, and the dynamic behavior was tackled to account for structure flexibility and interactions network formation. Structurally, the mutations studied lead to a loss of aromatic features, which is very important for the establishment of π-π interactions and could induce a structural preference by a specific coreceptor. These new insights into the structure-function relationship of HIV-2 gp125 V3 and surrounding regions will help in the design of better models and the design of new small molecules capable to inhibit the attachment and binding of HIV with host cells. [ABSTRACT FROM AUTHOR]

SARS-CoV-2 exploits angiotensin-converting enzyme 2 (ACE2) as a receptor to invade cells. It has been reported that the UK and South African strains may have higher transmission capabilities, eventually in part due to amino acid substitutions on the SARS-CoV-2 Spike protein. The pathogenicity seems modified but is still under investigation. Here we used the experimental structure of the Spike RBD domain co-crystallized with part of the ACE2 receptor, several in silico methods and numerous experimental data reported recently to analyze the possible impacts of three amino acid replacements (Spike K417N, E484K, N501Y) with regard to ACE2 binding. We found that the N501Y replacement in this region of the interface (present in both the UK and South African strains) should be favorable for the interaction with ACE2, while the K417N and E484K substitutions (South African strain) would seem neutral or even unfavorable. It is unclear if the N501Y substitution in the South African strain could counterbalance the K417N and E484K Spike replacements with regard to ACE2 binding. Our finding suggests that the UK strain should have higher affinity toward ACE2 and therefore likely increased transmissibility and possibly pathogenicity. If indeed the South African strain has a high transmission level, this could be due to the N501Y replacement and/or to substitutions in regions located outside the direct Spike–ACE2 interface but not so much to the K417N and E484K replacements. Yet, it should be noted that amino acid changes at Spike position 484 can lead to viral escape from neutralizing antibodies. Further, these amino acid substitutions do not seem to induce major structural changes in this region of the Spike protein. This structure–function study allows us to rationalize some observations made for the UK strain but raises questions for the South African strain. [ABSTRACT FROM AUTHOR]

Glucocorticoid (GC) resistance complicates the treatment of ~10–20% of children with nephrotic syndrome (NS), yet the molecular basis for resistance remains unclear. We used RNAseq analysis and in silico algorithm-based approaches on peripheral blood leukocytes from 12 children both at initial NS presentation and after ~7 weeks of GC therapy to identify a 12-gene panel able to differentiate steroid resistant NS (SRNS) from steroid-sensitive NS (SSNS). Among this panel, subsequent validation and analyses of one biologically relevant candidate, sulfatase 2 (SULF2), in up to a total of 66 children, revealed that both SULF2 leukocyte expression and plasma arylsulfatase activity Post/Pre therapy ratios were greater in SSNS vs. SRNS. However, neither plasma SULF2 endosulfatase activity (measured by VEGF binding activity) nor plasma VEGF levels, distinguished SSNS from SRNS, despite VEGF's reported role as a downstream mediator of SULF2's effects in glomeruli. Experimental studies of NS-related injury in both rat glomeruli and cultured podocytes also revealed decreased SULF2 expression, which were partially reversible by GC treatment of podocytes. These findings together suggest that SULF2 levels and activity are associated with GC resistance in NS, and that SULF2 may play a protective role in NS via the modulation of downstream mediators distinct from VEGF. [ABSTRACT FROM AUTHOR]

We present the analysis of the defective genetic pathways of the Late-Onset Alzheimer's Disease (LOAD) compared to the Mild Cognitive Impairment (MCI) and Healthy Controls (HC) using different sampling methodologies. These algorithms sample the uncertainty space that is intrinsic to any kind of highly underdetermined phenotype prediction problem, by looking for the minimum-scale signatures (header genes) corresponding to different random holdouts. The biological pathways can be identified performing posterior analysis of these signatures established via cross-validation holdouts and plugging the set of most frequently sampled genes into different ontological platforms. That way, the effect of helper genes, whose presence might be due to the high degree of under determinacy of these experiments and data noise, is reduced. Our results suggest that common pathways for Alzheimer's disease and MCI are mainly related to viral mRNA translation, influenza viral RNA transcription and replication, gene expression, mitochondrial translation, and metabolism, with these results being highly consistent regardless of the comparative methods. The cross-validated predictive accuracies achieved for the LOAD and MCI discriminations were 84% and 81.5%, respectively. The difference between LOAD and MCI could not be clearly established (74% accuracy). The most discriminatory genes of the LOAD-MCI discrimination are associated with proteasome mediated degradation and G-protein signaling. Based on these findings we have also performed drug repositioning using Dr. Insight package, proposing the following different typologies of drugs: isoquinoline alkaloids, antitumor antibiotics, phosphoinositide 3-kinase PI3K, autophagy inhibitors, antagonists of the muscarinic acetylcholine receptor and histone deacetylase inhibitors. We believe that the potential clinical relevance of these findings should be further investigated and confirmed with other independent studies. [ABSTRACT FROM AUTHOR]

We present the analysis of defective pathways in multiple myeloma (MM) using two recently developed sampling algorithms of the biological pathways: The Fisher's ratio sampler, and the holdout sampler. We performed the retrospective analyses of different gene expression datasets concerning different aspects of the disease, such as the existing difference between bone marrow stromal cells in MM and healthy controls (HC), the gene expression profiling of CD34+ cells in MM and HC, the difference between hyperdiploid and non-hyperdiploid myelomas, and the prediction of the chromosome 13 deletion, to provide a deeper insight into the molecular mechanisms involved in the disease. Our analysis has shown the importance of different altered pathways related to glycosylation, infectious disease, immune system response, different aspects of metabolism, DNA repair, protein recycling and regulation of the transcription of genes involved in the differentiation of myeloid cells. The main difference in genetic pathways between hyperdiploid and non-hyperdiploid myelomas are related to infectious disease, immune system response and protein recycling. Our work provides new insights on the genetic pathways involved in this complex disease and proposes novel targets for future therapies. [ABSTRACT FROM AUTHOR]

Entropy should directly reflect the extent of disorder in proteins. By clustering structurally related proteins and studying the multiple-sequence-alignment of the sequences of these clusters, we were able to link between sequence, structure, and disorder information. We introduced several parameters as measures of fluctuations at a given MSA site and used these as representative of the sequence and structure entropy at that site. In general, we found a tendency for negative correlations between disorder and structure, and significant positive correlations between disorder and the fluctuations in the system. We also found evidence for residue-type conservation for those residues proximate to potentially disordered sites. Mutation at the disorder site itself appear to be allowed. In addition, we found positive correlation for disorder and accessible surface area, validating that disordered residues occur in exposed regions of proteins. Finally, we also found that fluctuations in the dihedral angles at the original mutated residue and disorder are positively correlated while dihedral angle fluctuations in spatially proximal residues are negatively correlated with disorder. Our results seem to indicate permissible variability in the disordered site, but greater rigidity in the parts of the protein with which the disordered site interacts. This is another indication that disordered residues are involved in protein function. [ABSTRACT FROM AUTHOR]

We discuss the relationship between the problem of protein tertiary structure prediction from the amino acid sequence and the uncertainty analysis. The algorithm presented in this paper belongs to the category of decoy-based modeling, where different known protein models are used to establish a low dimensional space via principal component analysis. The low dimensional space is utilized to perform an energy optimization via a family of very explorative particle swarm optimizers to find the global minimum. The aim of this procedure is to get a representative sample of the nonlinear equivalent region, that is, protein models that have their energy lower than a certain energy bound. The posterior analysis of this family provides very valuable information about the backbone structure of the native conformation and its possible alternate states. This methodology has the advantage of being simple and fast and can help refine the tertiary protein structure. We comprehensively illustrate the performance of our algorithm on one protein from the CASP-9 protein structure prediction experiment. We also provide a theoretical analysis of the energy landscape found in the tertiary structure protein inverse problem, explaining why model reduction techniques (principal component analysis in this case) serve to alleviate the ill-posed character of this high dimensional optimization problem. In addition, we expand the computational benchmark with a summary of other CASP-9 proteins in the Appendix. [ABSTRACT FROM AUTHOR]

Ranking protein structure models is an elusive problem in bioinformatics. These models are evaluated on both the degree of similarity to the native structure and the folding pathway. Here, we simulated the use of the coarse-grained UNited RESidue (UNRES) force field as a tool to choose the best protein structure models for a given protein sequence among a pool of candidate models, using server data from the CASP11 experiment. Because the original UNRES was optimized for Molecular Dynamics simulations, we reoptimized UNRES using a deep feed-forward neural network, and we show that introducing additional descriptive features can produce better results. Overall, we found that the reoptimized UNRES performs better in selecting the best structures and tracking protein unwinding from its native state. We also found a relatively poor correlation between UNRES values and the model's Template Modeling Score (TMS). This is remedied by reoptimization. We discuss some cases where our reoptimization procedure is useful. [ABSTRACT FROM AUTHOR]

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes. [ABSTRACT FROM AUTHOR]

Protein-peptide interactions play essential roles in many cellular processes and their structural characterization is the major focus of current experimental and theoretical research. Two decades ago, it was proposed to employ the steered molecular dynamics (SMD) to assess the strength of protein-peptide interactions. The idea behind using SMD simulations is that the mechanical stability can be used as a promising and an efficient alternative to computationally highly demanding estimation of binding affinity. However, mechanical stability defined as a peak in force-extension profile depends on the choice of the pulling direction. Here we propose an uncommon choice of the pulling direction along resultant dipole moment (RDM) vector, which has not been explored in SMD simulations so far. Using explicit solvent all-atom MD simulations, we apply SMD technique to probe mechanical resistance of ligand-receptor system pulled along two different vectors. A novel pulling direction—when ligand unbinds along the RDM vector—results in stronger forces compared to commonly used ligand unbinding along center of masses vector. Our observation that RDM is one of the factors influencing the mechanical stability of protein-peptide complex can be used to improve the ranking of binding affinities by using mechanical stability as an effective scoring function. [ABSTRACT FROM AUTHOR]