Advances in Molecular Dynamics Simulations and Enhanced Sampling Methods for the Study of Protein Systems

Flexible RNA aptamers as inhibitors of Bacillus anthracis ribosomal protein S8: Insights from molecular dynamics simulations

Bacillus anthracis, the causative agent of anthrax, poses a substantial threat to public health and national security, and is recognized as a potential bioweapon due to its capacity to form resilient spores with enduring viability. Inhalation or ingestion of even minute quantities of aerosolized spores can lead to widespread illness and fatalities, underscoring the formidable lethality of the bacterium. With an untreated mortality rate of 100%, Bacillus anthracis is a disconcerting candidate for bioterrorism. In response to this critical scenario, we employed state-of-the-art computational tools to conceive and characterize flexible RNA aptamer therapeutics tailored for anthrax. The foundational structure of the flexible RNA aptamers was designed by removing the C2'-C3' in each nucleotide unit. Leveraging the crystal structure of Bacillus anthracis ribosomal protein S8 complexed with an RNA aptamer, we explored the structural, dynamic, and energetic aspects of the modified RNA aptamer - S8 protein complexes through extensive all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicas each), followed by drawing comparisons to the control system. Our findings demonstrate the enhanced binding competencies of the flexible RNA aptamers to the S8 protein via better shape complementarity and improved H-bond network compared to the control RNA aptamer. This research offers valuable insights into the development of RNA aptamer therapeutics targeting Bacillus anthracis, paving the way for innovative strategies to mitigate the impact of this formidable pathogen.

Human cardiac β-myosin powerstroke energetics: Thin filament, Pi displacement, and mutation effects

The powerstroke of human cardiac β-myosin is an important stage of the cross-bridge cycle that generates force for muscle contraction. However, the starting structure of this process has never been resolved, and the relative timing of the powerstroke and inorganic phosphate (Pi) release is still controversial. In this study, we generated an atomistic model of myosin on the thin filament and utilized metadynamics simulations to predict the absent starting structure of the powerstroke. We demonstrated that the displacement of Pi from the active site during the powerstroke is likely necessary, reducing the energy barrier of the conformation change. The effects of the presence of the thin filament, the hypertrophic cardiomyopathy mutation R712L, and the binding of mavacamten on the powerstroke process were also investigated.

Computational Insights into Papaveroline as an In Silico Drug Candidate for Alzheimer's Disease via Fyn Tyrosine Kinase Inhibition

Alzheimer's disease (AD) poses a significant global health challenge, necessitating the exploration of novel therapeutic strategies. Fyn Tyrosine Kinase has emerged as a key player in AD pathogenesis, making it an attractive target for drug development. This study focuses on investigating the potential of Papaveroline as a drug candidate for AD by targeting Fyn Tyrosine Kinase. The research employed high-throughput virtual screening and QSAR analysis were conducted to identify compounds with optimal drug-like properties, emphasizing adherence to ADMET parameters for further evaluation. Molecular dynamics simulations to analyze the binding interactions between Papaveroline and Staurosporine with Fyn Tyrosine Kinase over a 200-ns period. The study revealed detailed insights into the binding mechanisms and stability of the Papaveroline-Fyn complex, showcasing the compound's potential as an inhibitor of Fyn Tyrosine Kinase. Comparative analysis with natural compounds and a reference compound highlighted Papaveroline's unique characteristics and promising therapeutic implications for AD treatment. Overall, the findings underscore Papaveroline's potential as a valuable drug candidate for targeting Fyn Tyrosine Kinase in AD therapy, offering new avenues for drug discovery in neurodegenerative diseases. This study contributes to advancing our understanding of molecular interactions in AD pathogenesis and paves the way for further research and development in this critical area.

Molecular dynamics simulation and purification of chimeric L1/L2 protein from human papillomavirus type 52 expressed in Escherichia coli BL21 (DE3)

The available prophylactic vaccines for human papillomavirus (HPV) in the market are only effective against specific types of HPV, rendering them ineffective for other types of HPV infections. The objective of this research is to investigate the stability of the recombinant protein constructed, namely chimeric L1/L2 protein from HPV type 52, with improved cross-neutralization ability. The 3D model, predicted using Alphafold, Robetta, I-Tasser, and refined with Galaxy Refinement, is validated using Ramachandran plot analysis. The stability is verified through molecular dynamics simulations, considering parameters such as RMSD, RMSF, Rg, and SASA, where stable conditions are observed. The chimeric L1/L2 protein from HPV type 52 is purified using affinity chromatography, and the His-tag is cleaved using SUMO protease to obtain pure chimeric protein with the size of ~ 55 kDa. Western blot analysis confirms binding to anti-L1 HPV type 52 polyclonal antibody. The obtained vaccine candidate can be utilized as an effective prophylactic vaccine against HPV.

Molecular Dynamic Simulations for Biopolymers with Biomedical Applications

Computational modeling (CM) is a versatile scientific methodology used to examine the properties and behavior of complex systems, such as polymeric materials for biomedical bioengineering. CM has emerged as a primary tool for predicting, setting up, and interpreting experimental results. Integrating in silico and in vitro experiments accelerates scientific advancements, yielding quicker results at a reduced cost. While CM is a mature discipline, its use in biomedical engineering for biopolymer materials has only recently gained prominence. In biopolymer biomedical engineering, CM focuses on three key research areas: (A) Computer-aided design (CAD/CAM) utilizes specialized software to design and model biopolymers for various biomedical applications. This technology allows researchers to create precise three-dimensional models of biopolymers, taking into account their chemical, structural, and functional properties. These models can be used to enhance the structure of biopolymers and improve their effectiveness in specific medical applications. (B) Finite element analysis, a computational technique used to analyze and solve problems in engineering and physics. This approach divides the physical domain into small finite elements with simple geometric shapes. This computational technique enables the study and understanding of the mechanical and structural behavior of biopolymers in biomedical environments. (C) Molecular dynamics (MD) simulations involve using advanced computational techniques to study the behavior of biopolymers at the molecular and atomic levels. These simulations are fundamental for better understanding biological processes at the molecular level. Studying the wide-ranging uses of MD simulations in biopolymers involves examining the structural, functional, and evolutionary aspects of biomolecular systems over time. MD simulations solve Newton's equations of motion for all-atom systems, producing spatial trajectories for each atom. This provides valuable insights into properties such as water absorption on biopolymer surfaces and interactions with solid surfaces, which are crucial for assessing biomaterials. This review provides a comprehensive overview of the various applications of MD simulations in biopolymers. Additionally, it highlights the flexibility, robustness, and synergistic relationship between in silico and experimental techniques.

Predicting binding events in very flexible, allosteric, multi-domain proteins

Knowledge of the structures formed by proteins and small ligands is of fundamental importance for understanding molecular principles of chemotherapy and for designing new and more effective drugs. Due to the still high costs and to the several limitations of experimental techniques, it is most often desirable to predict these ligand-protein complexes in silico, particularly when screening for new putative drugs from databases of millions of compounds. While virtual screening based on molecular docking is widely used for this purpose, it generally fails in mimicking binding events associated with large conformational changes in the protein, particularly when the latter involve multiple domains. In this work, we describe a new methodology aimed at generating bound-like conformations of very flexible and allosteric proteins bearing multiple binding sites. Validation was performed on the enzyme adenylate kinase (ADK), a paradigmatic example of proteins that undergo very large conformational changes upon ligand binding. By only exploiting the unbound structure and the putative binding sites of the protein, we generated a significant fraction of bound-like structures, which employed in ensemble-docking calculations allowed to find native-like poses of substrates, inhibitors, and catalytically incompetent binders. Our protocol provides a general framework for the generation of bound-like conformations of flexible proteins that are suitable to host different ligands, demonstrating high sensitivity to the fine chemical details that regulate protein's activity. We foresee applications in virtual screening for difficult targets, prediction of the impact of amino acid mutations on structure and dynamics, and protein engineering.

Thermodynamics of a hyperthermostable carboxylesterase from Anoxybacillus geothermalis D9

Hyperthermostable enzymes are highly desirable biocatalysts due to their exceptional stability at extreme temperatures. Recently, a hyperthermostable carboxylesterase EstD9 from Anoxybacillus geothermalis D9 was biochemically characterized. The enzyme exhibited remarkable stability at high temperature. In this study, we attempted to probe the conformational adaptability of EstD9 under extreme conditions via in silico approaches. Circular dichroism revealed that EstD9 generated new β-sheets at 80 °C, making the core of the hydrolase fold more stable. Interestingly, the profiles of molecular dynamics simulation showed the lowest scores of radius of gyration and solvent accessible surface area (SASA) at 80 °C. Three loops were responsible for protecting the catalytic site, which resided at the interface between the large and cap domains. To further investigate the structural adaptation in extreme conditions, the intramolecular interactions of the native structure were investigated. EstD9 revealed 18 hydrogen bond networks, 7 salt bridges, and 9 hydrophobic clusters, which is higher than the previously reported thermostable Est30. Collectively, the analysis indicates that intramolecular interactions and structural dynamics play distinct roles in preserving the overall EstD9 structure at elevated temperatures. This work is relevant to both fundamental and applied research involving protein engineering of industrial thermostable enzymes.

Identification of potential anti-mucor agents by targeting endothelial cell receptor glucose-regulated protein-78 using in silico approach

Mucormycosis is a fungal infection of the sinuses, brain and lungs that is the cause of approximately 50% mortality rate despite the available first-line therapy. Glucose-Regulated Protein 78 (GRP78) is already reported to be a novel host receptor that mediates invasion and damage of human endothelial cells by Rhizopus oryzae and Rhizopus delemar, the most common etiologic species of Mucorales. The expression of GRP78 is also regulated by the levels of iron and glucose in the blood. There are several antifungal drugs in the market but they pose a serious side effect to the vital organs of the body. Therefore, there is an immediate need to discover effective drug molecules having increased efficacy with no side effects. With the help of various computational tools, the current study was attempted to determine potential antimucor agents against GRP78. The receptor molecule GRP78 was screened against 8820 known drugs deposited in DrugBank library using high-throughput virtual screening method. Total top 10 compounds were selected based on the binding energies greater than the reference co-crystal molecule. Furthermore, molecular dynamic (MD) simulations using AMBER were performed to calculate the stability of the top-ranked compounds in the active site of GRP78. After extensive computational studies, we propose that two compounds (CID439153 and CID5289104) have inhibitory potency against mucormycosis and can serve as potential drugs that can form the basis of treating mucormycosis disease.Communicated by Ramaswamy H. Sarma.

Accurate Characterization of Binding Kinetics and Allosteric Mechanisms for the HSP90 Chaperone Inhibitors Using AI-Augmented Integrative Biophysical Studies

The binding kinetics of drugs to their targets are gradually being recognized as a crucial indicator of the efficacy of drugs in vivo, leading to the development of various computational methods for predicting the binding kinetics in recent years. However, compared with the prediction of binding affinity, the underlying structure and dynamic determinants of binding kinetics are more complicated. Efficient and accurate methods for predicting binding kinetics are still lacking. In this study, quantitative structure-kinetics relationship (QSKR) models were developed using 132 inhibitors targeting the ATP binding domain of heat shock protein 90α (HSP90α) to predict the dissociation rate constant (koff), enabling a direct assessment of the drug-target residence time. These models demonstrated good predictive performance, where hydrophobic and hydrogen bond interactions significantly influence the koff prediction. In subsequent applications, our models were used to assist in the discovery of new inhibitors for the N-terminal domain of HSP90α (N-HSP90α), demonstrating predictive capabilities on an experimental validation set with a new scaffold. In X-ray crystallography experiments, the loop-middle conformation of apo N-HSP90α was observed for the first time (previously, the loop-middle conformation had only been observed in holo-N-HSP90α structures). Interestingly, we observed different conformations of apo N-HSP90α simultaneously in an asymmetric unit, which was also observed in a holo-N-HSP90α structure, suggesting an equilibrium of conformations between different states in solution, which could be one of the determinants affecting the binding kinetics of the ligand. Different ligands can undergo conformational selection or alter the equilibrium of conformations, inducing conformational rearrangements and resulting in different effects on binding kinetics. We then used molecular dynamics simulations to describe conformational changes of apo N-HSP90α in different conformational states. In summary, the study of the binding kinetics and molecular mechanisms of N-HSP90α provides valuable information for the development of more targeted therapeutic approaches.

ELViM: Exploring Biomolecular Energy Landscapes through Multidimensional Visualization

Molecular dynamics (MD) simulations provide a powerful means of exploring the dynamic behavior of biomolecular systems at the atomic level. However, analyzing the vast data sets generated by MD simulations poses significant challenges. This article discusses the energy landscape visualization method (ELViM), a multidimensional reduction technique inspired by the energy landscape theory. ELViM transcends one-dimensional representations, offering a comprehensive analysis of the effective conformational phase space without the need for predefined reaction coordinates. We apply the ELViM to study the folding landscape of the antimicrobial peptide Polybia-MP1, showcasing its versatility in capturing complex biomolecular dynamics. Using dissimilarity matrices and a force-scheme approach, the ELViM provides intuitive visualizations, revealing structural correlations and local conformational signatures. The method is demonstrated to be adaptable, robust, and applicable to various biomolecular systems.

Explainable Supervised Machine Learning Model To Predict Solvation Gibbs Energy

Many challenges persist in developing accurate computational models for predicting solvation free energy (ΔGsol). Despite recent developments in Machine Learning (ML) methodologies that outperformed traditional quantum mechanical models, several issues remain concerning explanatory insights for broad chemical predictions with an acceptable speed-accuracy trade-off. To overcome this, we present a novel supervised ML model to predict the ΔGsol for an array of solvent-solute pairs. Using two different ensemble regressor algorithms, we made fast and accurate property predictions using open-source chemical features, encoding complex electronic, structural, and surface area descriptors for every solvent and solute. By integrating molecular properties and chemical interaction features, we have analyzed individual descriptor importance and optimized our model though explanatory information form feature groups. On aqueous and organic solvent databases, ML models revealed the predictive relevance of solutes with increasing polar surface area and decreasing polarizability, yielding better results than state-of-the-art benchmark Neural Network methods (without complex quantum mechanical or molecular dynamic simulations). Both algorithms successfully outperformed previous ΔGsol predictions methods, with a maximum absolute error of 0.22 ± 0.02 kcal mol-1, further validated in an external benchmark database and with solvent hold-out tests. With these explanatory and statistical insights, they allow a thoughtful application of this method for predicting other thermodynamic properties, stressing the relevance of ML modeling for further complex computational chemistry problems.

Deep learning for advancing peptide drug development: Tools and methods in structure prediction and design

Peptides can bind challenging disease targets with high affinity and specificity, offering enormous opportunities for addressing unmet medical needs. However, peptides' unique features, including smaller size, increased structural flexibility, and limited data availability, pose additional challenges to the design process compared to proteins. This review explores the dynamic field of peptide therapeutics, leveraging deep learning to enhance structure prediction and design. Our exploration encompasses various facets of peptide research, ranging from dataset curation handling to model development. As deep learning technologies become more refined, we channel our efforts into peptide structure prediction and design, aligning with the fundamental principles of structure-activity relationships in drug development. To guide researchers in harnessing the potential of deep learning to advance peptide drug development, our insights comprehensively explore current challenges and future directions of peptide therapeutics.

Mapping secondary substrate-binding sites on the GH11 xylanase from Bacillus subtilis

Xylanases are of significant interest for biomass conversion technologies. Here, we investigated the allosteric regulation of xylan hydrolysis by the Bacillus subtilis GH11 endoxylanase. Molecular dynamics simulations (MDS) in the presence of xylobiose identified binding to the active site and two potential secondary binding sites (SBS) around surface residues Asn54 and Asn151. Arabinoxylan titration experiments with single cysteine mutants N54C and N151C labeled with the thiol-reactive fluorophore acrylodan or the ESR spin-label MTSSL validated the MDS results. Ligand binding at the SBS around Asn54 confirms previous reports, and analysis of the second SBS around N151C discovered in the present study includes residues Val98/Ala192/Ser155/His156. Understanding the regulation of xylanases contributes to efforts for industrial decarbonization and to establishing a sustainable energy matrix.

Trapping a non-cognate nucleotide upon initial binding for replication fidelity control in SARS-CoV-2 RNA dependent RNA polymerase

The RNA dependent RNA polymerase (RdRp) in SARS-CoV-2 is a highly conserved enzyme responsible for viral genome replication/transcription. To understand how the viral RdRp achieves fidelity control during such processes, here we computationally investigate the natural non-cognate vs. cognate nucleotide addition and selectivity during viral RdRp elongation. We focus on the nucleotide substrate initial binding (RdRp active site open) to the prechemical insertion (active site closed) of the RdRp. The current studies were first carried out using microsecond ensemble equilibrium all-atom molecular dynamics (MD) simulations. Due to the slow conformational changes (from open to closed) during nucleotide insertion and selection, enhanced or umbrella sampling methods have been further employed to calculate the free energy profiles of the nucleotide insertion. Our studies find notable stability of noncognate dATP and GTP upon initial binding in the active-site open state. The results indicate that while natural cognate ATP and Remdesivir drug analogue (RDV-TP) are biased toward stabilization in the closed state to facilitate insertion, the natural non-cognate dATP and GTP can be well trapped in off-path initial binding configurations and prevented from insertion so that to be further rejected. The current work thus presents the intrinsic nucleotide selectivity of SARS-CoV-2 RdRp for natural substrate fidelity control, which should be considered in antiviral drug design.

Generalized open-source workflows for atomistic molecular dynamics simulations of viral helicases

Viral helicases are promising targets for the development of antiviral therapies. Given their vital function of unwinding double-stranded nucleic acids, inhibiting them blocks the viral replication cycle. Previous studies have elucidated key structural details of these helicases, including the location of substrate binding sites, flexible domains, and the discovery of potential inhibitors. Here we present a series of new Galaxy tools and workflows for performing and analyzing molecular dynamics simulations of viral helicases. We first validate them by demonstrating recapitulation of data from previous simulations of Zika (NS3) and SARS-CoV-2 (NSP13) helicases in apo and complex with inhibitors. We further demonstrate the utility and generalizability of these Galaxy workflows by applying them to new cases, proving their usefulness as a widely accessible method for exploring antiviral activity.

Integrating biocomputational techniques for Breast cancer drug discovery via the HER-2, BCRA, VEGF and ER protein targets

Computational modelling remains an indispensable technique in drug discovery. With myriad of high computing resources, and improved modelling algorithms, there has been a high-speed in the drug development cycle with promising success rate compared to the traditional route. For example, lapatinib; a well-known anticancer drug with clinical applications was discovered with computational drug design techniques. Similarly, molecular modelling has been applied to various disease areas ranging from cancer to neurodegenerative diseases. The techniques ranges from high-throughput virtual screening, molecular mechanics with generalized Born and surface area solvation (MM/GBSA) to molecular dynamics simulation. This review focuses on the application of computational modelling tools in the identification of drug candidates for Breast cancer. First, we begin with a succinct overview of molecular modelling in the drug discovery process. Next, we take note of special efforts on the developments and applications of combining these techniques with particular emphasis on possible breast cancer therapeutic targets such as estrogen receptor (ER), human epidermal growth factor receptor 2 (HER2), vascular endothelial growth factor (VEGF), breast cancer gene 1 (BRCA1), and breast cancer gene 2 (BRCA2). Finally, we discussed the search for covalent inhibitors against these receptors using computational techniques, advances, pitfalls, possible solutions, and future perspectives.

Fn-HMGB1 Adsorption Behavior Initiates Early Immune Recognition and Subsequent Osteoinduction of Biomaterials

Implantable biomaterials are widely used in bone tissue engineering, but little is still known about how they initiate early immune recognition and the initial dynamics. Herein, the early immune recognition and subsequent osteoinduction of biphasic calcium phosphate (BCP) after implantation to the protein adsorption behavior is attributed. By liquid chromatography tandem mass spectrometry (LC-MS/MS) analysis, the biomaterial-related molecular patterns (BAMPs) formed after BCP implantation are mapped, dominated by the highly expressed extracellular matrix protein fibronectin (Fn) and the high mobility group box 1 (HMGB1). Molecular dynamics simulations show that Fn has the ability to bind more readily to the BCP surface than HMGB1. The preferential binding of Fn provides a higher adsorption energy for HMGB1. Furthermore, multiple hydrogen bonding sites between HMGB1 and Fn are demonstrated using a molecular docking approach. Ultimately, the formation of BAMPs through HMGB1 antagonist glycyrrhizic acid (GA), resulting in impaired immune recognition of myeloid differentiation factor 88 (MYD88) mediated dendritic cells (DCs) and macrophages (Mφs), as well as failed osteoinduction processes is obstructed. This study introduces a mechanism for early immune recognition of implant materials based on protein adsorption, providing perspectives for future design and application of tissue engineering materials.

Exploring the role of microbial proteins in controlling environmental pollutants based on molecular simulation

Molecular simulation has been widely used to study microbial proteins' structural composition and dynamic properties, such as volatility, flexibility, and stability at the microscopic scale. Herein, this review describes the key elements of molecular docking and molecular dynamics (MD) simulations in molecular simulation; reviews the techniques combined with molecular simulation, such as crystallography, spectroscopy, molecular biology, and machine learning, to validate simulation results and bridge information gaps in the structure, microenvironmental changes, expression mechanisms, and intensity quantification; illustrates the application of molecular simulation, in characterizing the molecular mechanisms of interaction of microbial proteins with four different types of contaminants, namely heavy metals (HMs), pesticides, dyes and emerging contaminants (ECs). Finally, the review outlines the important role of molecular simulations in the study of microbial proteins for controlling environmental contamination and provides ideas for the application of molecular simulation in screening microbial proteins and incorporating targeted mutagenesis to obtain more effective contaminant control proteins.

In Silico and In Vivo Evaluation of the Maqui Berry (Aristotelia chilensis (Mol.) Stuntz) on Biochemical Parameters and Oxidative Stress Markers in a Metabolic Syndrome Model

Metabolic syndrome (MetS) is a complex disease that includes metabolic and physiological alterations in various organs such as the heart, pancreas, liver, and brain. Reports indicate that blackberry consumption, such as maqui berry, has a beneficial effect on chronic diseases such as cardiovascular disease, obesity, and diabetes. In the present study, in vivo and in silico studies have been performed to evaluate the molecular mechanisms implied to improve the metabolic parameters of MetS. Fourteen-day administration of maqui berry reduces weight gain, blood fasting glucose, total blood cholesterol, triacylglycerides, insulin resistance, and blood pressure impairment in the diet-induced MetS model in male and female rats. In addition, in the serum of male and female rats, the administration of maqui berry (MB) improved the concentration of MDA, the activity of SOD, and the formation of carbonyls in the group subjected to the diet-induced MetS model. In silico studies revealed that delphinidin and its glycosylated derivatives could be ligands of some metabolic targets such as α-glucosidase, PPAR-α, and PPAR-γ, which are related to MetS parameters. The experimental results obtained in the study suggest that even at low systemic concentrations, anthocyanin glycosides and aglycones could simultaneously act on different targets related to MetS. Therefore, these molecules could be used as coadjuvants in pharmacological interventions or as templates for designing new multitarget molecules to manage patients with MetS.

Computational Analysis of CD46 Protein Interaction with SARS-CoV-2 Structural Proteins: Elucidating a Putative Viral Entry Mechanism into Human Cells

This study examines an unexplored aspect of SARS-CoV-2 entry into host cells, which is widely understood to occur via the viral spike (S) protein's interaction with human ACE2-associated proteins. While vaccines and inhibitors targeting this mechanism are in use, they may not offer complete protection against reinfection. Hence, we investigate putative receptors and their cofactors. Specifically, we propose CD46, a human membrane cofactor protein, as a potential putative receptor and explore its role in cellular invasion, acting possibly as a cofactor with other viral structural proteins. Employing computational techniques, we created full-size 3D models of human CD46 and four key SARS-CoV-2 structural proteins-EP, MP, NP, and SP. We further developed 3D models of CD46 complexes interacting with these proteins. The primary aim is to pinpoint the likely interaction domains between CD46 and these structural proteins to facilitate the identification of molecules that can block these interactions, thus offering a foundation for novel pharmacological treatments for SARS-CoV-2 infection.

The multi-protein targeting potential of bioactive syringin in inflammatory diseases: using molecular modelling and in-silico analysis of regulatory elements

Inflammation plays a crucial role in the onset or progression of a variety of acute and chronic diseases. Non-steroidal anti-inflammatory drugs (NSAIDs) are the only available FDA-approved therapy. The therapeutic outcome of NSAIDs is still finite due to off-target effects and extreme side effects on other vital organs. Bioactive syringin has been manifested to hold anti-osteoporosis, cardiac hypertrophy, alter autophagy, anti-cancer, neuro-preventive effects, etc. However, its multi-protein targeting potential in inflammation mostly remains unexplored. In the present work, we have checked the multi-protein targeting potential of bioactive glycoside syringin in inflammatory diseases. Based on the binding score of protein-ligand complexes, glycoside syringin scored greater than -7 kcal/mol against 12 inflammatory proteins. Our molecular dynamic simulation study (200 ns) confirmed that bioactive syringin remained inside the binding cavity of inflammatory proteins (JAK1, TYK2, and COX1) in a stable conformation. Further, our co-expression analysis suggests that these genes play an essential role in multiple pathways and are regulated by multiple miRNAs. Our study demonstrates that bioactive glycoside syringin might be a multi-protein targeting potential against inflammatory diseases and could be further investigated utilizing different preclinical approaches.Communicated by Ramaswamy H. Sarma.

Probing the Nucleic Acid Flexibility to Disarm the Viral Counter-Defense Machinery: Design and Characterization of Potent p19 Inhibitors

Plant viruses are highly destructive and significant contributors to several global pandemics and epidemics in plants. A viral disease outbreak in plants can cause a scarcity of food supply and is a severe concern to humanity. The siRNA (small interfering RNA)-mediated RNA-induced silencing complex (RISC) formation is a primary defense mechanism in plants against viruses, where the RISC binds and degrades viral mRNAs. As a counter-defense, many viruses encode RNA-silencing suppressor proteins (e.g., the p19 protein from the Tombusviridae family) for viral proliferation in plants. The functional form of p19 (homodimer) binds to plant siRNA with high affinities, thereby interrupting the RISC formation and thus preventing the viral mRNA silencing in plants. By altering the RISC formation, the p19 protein helps the virus invasion in the plant and ultimately stunts host growth. In this study, we designed several modified siRNA-based molecules for p19 inhibition. The viral p19 protein is known to interact predominantly through H-bonds with 2'-OH and phosphates of the plant siRNA. We utilized this information and in silico-designed flexible substituents of siRNA, where we removed the C2'-C3' bond in each nucleotide unit. We performed all-atom explicit-solvent molecular dynamics simulations (400 ns, 3 replicates each) for control/modified siRNA─p19 complexes (8 in total) followed by energetic estimations. Strikingly, in a few modified complexes, the siRNA not only retained the double-helical structural integrity but also displayed remarkably enhanced p19 binding compared to the control siRNA; hence, we consider it important to perform biological and chemical in vitro and in vivo studies on proposed flexible nucleic acids as p19 inhibitors for crop protection.

On the Role of Temperature in the Depolymerization of PET by FAST-PETase: An Atomistic Point of View on Possible Active Site Pre-Organization and Substrate-Destabilization Effects

Enzyme FAST-PETase, recently obtained by a machine learning approach, can depolymerize poly(ethylene terephthalate) (PET), a synthetic resin employed in plastics and in clothing fibers. Therefore it represents a promising solution for the recycling of PET-based materials. In this study, a model of PET was adopted to describe the substrate, and all-atoms classical molecular dynamics (MD) simulations on apo- and substrate-bound FAST-PETase were carried out at 30 and 50 °C to provide atomistic details on the binding step of the catalytic cycle. Comparative analysis shed light on the interactions occurring between the FAST-PETase and 4PET at 50 °C, the optimal working conditions of the enzyme. Pre-organization of the enzyme active and binding sites has been highlighted, while MD simulations of FAST-PETase:4PET pointed out the occurrence of solvent-inaccessible conformations of the substrate promoted by the enzyme. Indeed, neither of these conformations was observed during MD simulations of the substrate alone in solution performed at 30, 50 and 150 °C. The analysis led us to propose that, at 50 °C, the FAST-PETase is pre-organized to bind the PET and that the interactions occurring in the binding site can promote a more reactive conformation of PET substrate, thus enhancing the catalytic activity of the enzyme.

Novel computational and drug design strategies for inhibition of human papillomavirus-associated cervical cancer and DNA polymerase theta receptor by Apigenin derivatives

The present study deals with the advanced in-silico analyses of several Apigenin derivatives to explore human papillomavirus-associated cervical cancer and DNA polymerase theta inhibitor properties by molecular docking, molecular dynamics, QSAR, drug-likeness, PCA, a dynamic cross-correlation matrix and quantum calculation properties. The initial literature study revealed the potent antimicrobial and anticancer properties of Apigenin, prompting the selection of its potential derivatives to investigate their abilities as inhibitors of human papillomavirus-associated cervical cancer and DNA polymerase theta. In silico molecular docking was employed to streamline the findings, revealing promising energy-binding interactions between all Apigenin derivatives and the targeted proteins. Notably, Apigenin 4'-O-Rhamnoside and Apigenin-4'-Alpha-L-Rhamnoside demonstrated higher potency against the HPV45 oncoprotein E7 (PDB ID 2EWL), while Apigenin and Apigenin 5-O-Beta-D-Glucopyranoside exhibited significant binding energy against the L1 protein in humans. Similarly, a binding affinity range of - 7.5 kcal/mol to - 8.8 kcal/mol was achieved against DNA polymerase theta, indicating the potential of Apigenin derivatives to inhibit this enzyme (PDB ID 8E23). This finding was further validated through molecular dynamic simulation for 100 ns, analyzing parameters such as RMSD, RMSF, SASA, H-bond, and RoG profiles. The results demonstrated the stability of the selected compounds during the simulation. After passing the stability testing, the compounds underwent screening for ADMET, pharmacokinetics, and drug-likeness properties, fulfilling all the necessary criteria. QSAR, PCA, dynamic cross-correlation matrix, and quantum calculations were conducted, yielding satisfactory outcomes. Since this study utilized in silico computational approaches and obtained outstanding results, further validation is crucial. Therefore, additional wet-lab experiments should be conducted under in vivo and in vitro conditions to confirm the findings.

Perspectives on Ligand/Protein Binding Kinetics Simulations: Force Fields, Machine Learning, Sampling, and User-Friendliness

Computational techniques applied to drug discovery have gained considerable popularity for their ability to filter potentially active drugs from inactive ones, reducing the time scale and costs of preclinical investigations. The main focus of these studies has historically been the search for compounds endowed with high affinity for a specific molecular target to ensure the formation of stable and long-lasting complexes. Recent evidence has also correlated the in vivo drug efficacy with its binding kinetics, thus opening new fascinating scenarios for ligand/protein binding kinetic simulations in drug discovery. The present article examines the state of the art in the field, providing a brief summary of the most popular and advanced ligand/protein binding kinetics techniques and evaluating their current limitations and the potential solutions to reach more accurate kinetic models. Particular emphasis is put on the need for a paradigm change in the present methodologies toward ligand and protein parametrization, the force field problem, characterization of the transition states, the sampling issue, and algorithms' performance, user-friendliness, and data openness.

Protein-Inspired Control over Synthetic Polymer Folding for Structured Functional Nanoparticles in Water

The folding of proteins into functional nanoparticles with defined 3D structures has inspired chemists to create simple synthetic systems mimicking protein properties. The folding of polymers into nanoparticles in water proceeds via different strategies, resulting in the global compaction of the polymer chain. Herein, we review the different methods available to control the conformation of synthetic polymers and collapse/fold them into structured, functional nanoparticles, such as hydrophobic collapse, supramolecular self-assembly, and covalent cross-linking. A comparison is made between the design principles of protein folding to synthetic polymer folding and the formation of structured nanocompartments in water, highlighting similarities and differences in design and function. We also focus on the importance of structure for functional stability and diverse applications in complex media and cellular environments.

The effect of mutation on neurotoxicity reduction of new chimeric reteplase, a computational study

Background and purpose

Excitotoxicity in nerve cells is a type of neurotoxicity in which excessive stimulation of receptors (such as N-methyl-d-aspartate glutamate receptors (NMDAR)) leads to the influx of high-level calcium ions into cells and finally cell damage or death. This complication can occur after taking some of the plasminogen activators like tissue plasminogen activator and reteplase. The interaction of the kringle2 domain in such plasminogen activator with the amino-terminal domain (ATD) of the NR1 subunit of NMDAR finally leads to excitotoxicity. In this study, we assessed the interaction of two new chimeric reteplase, mutated in the kringle2 domain, with ATD and compared the interaction of wild-type reteplase with ATD, computationally.

Experimental approach

Homology modeling, protein docking, molecular dynamic simulation, and molecular dynamics trajectory analysis were used for the assessment of this interaction.

Findings/Results

The results of the free energy analysis between reteplase and ATD (wild reteplase: -2127.516 ± 0.0, M1-chr: -1761.510 ± 0.0, M2-chr: -521.908 ± 0.0) showed lower interaction of this chimeric reteplase with ATD compared to the wild type.

Conclusion and implications

The decreased interaction between two chimeric reteplase and ATD of NR1 subunit in NMDAR which leads to lower neurotoxicity related to these drugs, can be the start of a way to conduct more tests and if the results confirm this feature, they can be considered potential drugs in acute ischemic stroke treatment.

Computational Modeling of DNA 3D Structures: From Dynamics and Mechanics to Folding

DNA carries the genetic information required for the synthesis of RNA and proteins and plays an important role in many processes of biological development. Understanding the three-dimensional (3D) structures and dynamics of DNA is crucial for understanding their biological functions and guiding the development of novel materials. In this review, we discuss the recent advancements in computer methods for studying DNA 3D structures. This includes molecular dynamics simulations to analyze DNA dynamics, flexibility, and ion binding. We also explore various coarse-grained models used for DNA structure prediction or folding, along with fragment assembly methods for constructing DNA 3D structures. Furthermore, we also discuss the advantages and disadvantages of these methods and highlight their differences.

Past, Present, and Future Perspectives on Computer-Aided Drug Design Methodologies

The application of computational approaches in drug discovery has been consolidated in the last decades. These families of techniques are usually grouped under the common name of "computer-aided drug design" (CADD), and they now constitute one of the pillars in the pharmaceutical discovery pipelines in many academic and industrial environments. Their implementation has been demonstrated to tremendously improve the speed of the early discovery steps, allowing for the proficient and rational choice of proper compounds for a desired therapeutic need among the extreme vastness of the drug-like chemical space. Moreover, the application of CADD approaches allows the rationalization of biochemical and interactive processes of pharmaceutical interest at the molecular level. Because of this, computational tools are now extensively used also in the field of rational 3D design and optimization of chemical entities starting from the structural information of the targets, which can be experimentally resolved or can also be obtained with other computer-based techniques. In this work, we revised the state-of-the-art computer-aided drug design methods, focusing on their application in different scenarios of pharmaceutical and biological interest, not only highlighting their great potential and their benefits, but also discussing their actual limitations and eventual weaknesses. This work can be considered a brief overview of computational methods for drug discovery.

Accelerating therapeutic protein design with computational approaches toward the clinical stage

Therapeutic protein, represented by antibodies, is of increasing interest in human medicine. However, clinical translation of therapeutic protein is still largely hindered by different aspects of developability, including affinity and selectivity, stability and aggregation prevention, solubility and viscosity reduction, and deimmunization. Conventional optimization of the developability with widely used methods, like display technologies and library screening approaches, is a time and cost-intensive endeavor, and the efficiency in finding suitable solutions is still not enough to meet clinical needs. In recent years, the accelerated advancement of computational methodologies has ushered in a transformative era in the field of therapeutic protein design. Owing to their remarkable capabilities in feature extraction and modeling, the integration of cutting-edge computational strategies with conventional techniques presents a promising avenue to accelerate the progression of therapeutic protein design and optimization toward clinical implementation. Here, we compared the differences between therapeutic protein and small molecules in developability and provided an overview of the computational approaches applicable to the design or optimization of therapeutic protein in several developability issues.

Fast, Accurate, and System-Specific Variable-Resolution Modeling of Proteins

In recent years, a few multiple-resolution modeling strategies have been proposed, in which functionally relevant parts of a biomolecule are described with atomistic resolution, with the remainder of the system being concurrently treated using a coarse-grained model. In most cases, the parametrization of the latter requires lengthy reference all-atom simulations and/or the usage of off-shelf coarse-grained force fields, whose interactions have to be refined to fit the specific system under examination. Here, we overcome these limitations through a novel multiresolution modeling scheme for proteins, dubbed coarse-grained anisotropic network model for variable resolution simulations, or CANVAS. This scheme enables a user-defined modulation of the resolution level throughout the system structure; a fast parametrization of the potential without the necessity of reference simulations; and the straightforward usage of the model on the most commonly used molecular dynamics platforms. The method is presented and validated with two case studies, the enzyme adenylate kinase and the therapeutic antibody pembrolizumab, by comparing the results obtained with the CANVAS model against fully atomistic simulations. The modeling software, implemented in Python, is made freely available for the community on a collaborative github repository.

Comparative Binding Ability of Human Monoclonal Antibodies against Omicron Variants of SARS-CoV-2: An In Silico Investigation

Mutation(s) in the spike protein is the major characteristic trait of newly emerged SARS-CoV-2 variants such as Alpha (B.1.1.7), Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2), and Delta-plus. Omicron (B.1.1.529) is the latest addition and it has been characterized by high transmissibility and the ability to escape host immunity. Recently developed vaccines and repurposed drugs exert limited action on Omicron strains and hence new therapeutics are immediately needed. Herein, we have explored the efficiency of twelve therapeutic monoclonal antibodies (mAbs) targeting the RBD region of the spike glycoprotein against all the Omicron variants bearing a mutation in spike protein through molecular docking and molecular dynamics simulation. Our in silico evidence reveals that adintivimab, beludivimab, and regadanivimab are the most potent mAbs to form strong biophysical interactions and neutralize most of the Omicron variants. Considering the efficacy of mAbs, we incorporated CDRH3 of beludavimab within the framework of adintrevimab, which displayed a more intense binding affinity towards all of the Omicron variants viz. BA.1, BA.2, BA.2.12.1, BA.4, and BA.5. Furthermore, the cDNA of chimeric mAb was cloned in silico within pET30ax for recombinant production. In conclusion, the present study represents the candidature of human mAbs (beludavimab and adintrevimab) and the therapeutic potential of designed chimeric mAb for treating Omicron-infected patients.

Modelling eNvironment for Isoforms (MoNvIso): A general platform to predict structural determinants of protein isoforms in genetic diseases

The seamless integration of human disease-related mutation data into protein structures is an essential component of any attempt to correctly assess the impact of the mutation. The key step preliminary to any structural modelling is the identification of the isoforms onto which mutations should be mapped due to there being several functionally different protein isoforms from the same gene. To handle large sets of data coming from omics techniques, this challenging task needs to be automatized. Here we present the MoNvIso (Modelling eNvironment for Isoforms) code, which identifies the most useful isoform for computational modelling, balancing the coverage of mutations of interest and the availability of templates to build a structural model of both the wild-type isoform and the related variants.

Across atoms to crossing continents: Application of similarity measures to biological location data

Biological processes involve movements across all measurable scales. Similarity measures can be applied to compare and analyze these movements but differ in how differences in movement are aggregated across space and time. The present study reviews frequently-used similarity measures, such as the Hausdorff distance, Fréchet distance, Dynamic Time Warping, and Longest Common Subsequence, jointly with several measures less used in biological applications (Wasserstein distance, weak Fréchet distance, and Kullback-Leibler divergence), and provides computational tools for each of them that may be used in computational biology. We illustrate the use of the selected similarity measures in diagnosing differences within two extremely contrasting sets of biological data, which, remarkably, may both be relevant for magnetic field perception by migratory birds. Specifically, we assess and discuss cryptochrome protein conformational dynamics and extreme migratory trajectories of songbirds between Alaska and Africa. We highlight how similarity measures contrast regarding computational complexity and discuss those which can be useful in noise elimination or, conversely, are sensitive to spatiotemporal scales.

Accuracy and precision of binding free energy prediction for a tacrine related lead inhibitor of acetylcholinesterase with an arsenal of supercomputerized molecular modelling methods: a comparative study

Nowadays, advanced computational chemistry methods offer various strategies for revealing prospective hit structures in drug development essentially through accurate binding free energy predictions. After the era of molecular docking and quantitative structure-activity relationships, much interest has been lately oriented to perturbed molecular dynamic approaches like replica exchange with solute tempering and free energy perturbation (REST/FEP) and the potential of the mean force with adaptive biasing and accelerated weight histograms (PMF/AWH). Both of these receptor-based techniques can exploit exascale CPU&GPU supercomputers to achieve high throughput performance. In this fundamental study, we have compared the predictive power of a panel of supercomputerized molecular modelling methods to distinguish the major binding modes and the corresponding binding free energies of a promising tacrine related potential antialzheimerics in human acetylcholinesterase. The binding free energies were estimated using flexible molecular docking, molecular mechanics/generalized Born surface area/Poisson-Boltzmann surface area (MM/GBSA/PBSA), transmutation REST/FEP with 12 x 5 ns/λ windows, annihilation FEP with 20 x 5 ns/λ steps, PMF with weight histogram analysis method (WHAM) and 40 x 5 ns samples, and PMF/AWH with 10 x 100 ns replicas. Confrontation of the classical approaches such as canonical molecular dynamics and molecular docking with alchemical calculations and steered molecular dynamics enabled us to show how large errors in ΔG predictions can be expected if these in silico methods are employed in the elucidation of a common case of enzyme inhibition.Communicated by Ramaswamy H. Sarma.

In Vitro, In Silico and Network Pharmacology Mechanistic Approach to Investigate the α-Glucosidase Inhibitors Identified by Q-ToF-LCMS from Phaleria macrocarpa Fruit Subcritical CO2 Extract

The fruit of Phaleria macrocarpa have been traditionally used as an antidiabetic remedy in Malaysia and neighbouring countries. Despite its potential for diabetes treatment, no scientific study has ever been conducted to predict the inhibitor interaction of the protein α-glucosidase identified in an extract prepared with a non-conventional extraction technique. Hence, the major aim of this research was to evaluate the in vitro antioxidant, the α-glucosidase inhibitors, and the molecular dynamic simulations of the α-glucosidase inhibitors identified by Quadrupole Time-of-Flight Liquid Chromatography Mass Spectrometry (Q-ToF-LCMS) analysis. Initially, dry fruit were processed using non-conventional and conventional extraction methods to obtain subcritical carbon dioxide extracts (SCE-1 and SCE-2) and heating under reflux extract (HRE), respectively. Subsequently, all extracts were evaluated for their in vitro antioxidative and α-glucosidase inhibitory potentials. Subsequently, the most bioactive extract (SCE-2) was subjected to Q-ToF-LCMS analysis to confirm the presence of α-glucosidase inhibitors, which were then analysed through molecular dynamic simulations and network pharmacology approaches to confirm their possible mechanism of action. The highest inhibitory effects of the 2,2-diphenyl-1-picrylhydrazyl (DPPH) radical and α-glucosidase on SCE-2 was found as 75.36 ± 0.82% and 81.79 ± 0.82%, respectively, compared to the SCE-1 and HRE samples. The Q-ToF-LCMS analysis tentatively identified 14 potent α-glucosidase inhibitors. Finally, five identified compounds, viz., lupenone, swertianolin, m-coumaric acid, pantothenic acid, and 8-C-glucopyranosyleriodictylol displayed significant stability, compactness, stronger protein-ligand interaction up to 100 ns further confirming their potential as α-glucosidase inhibitors. Consequently, it was concluded that the SCE-2 possesses a strong α-glucosidase inhibitory effect due to the presence of these compounds. The findings of this study might prove useful to develop these compounds as alternative safe α-glucosidase inhibitors to manage diabetes more effectively.

Stochastic Lag Time Parameterization for Markov State Models of Protein Dynamics

Markov state models (MSMs) play a key role in studying protein conformational dynamics. A sliding count window with a fixed lag time is widely used to sample sub-trajectories for transition counting and MSM construction. However, sub-trajectories sampled with a fixed lag time may not perform well under different selections of lag time, which requires strong prior practice and leads to less robust estimation. To alleviate it, we propose a novel stochastic method from a Poisson process to generate perturbative lag time for sub-trajectory sampling and utilize it to construct a Markov chain. Comprehensive evaluations on the double-well system, WW domain, BPTI, and RBD-ACE2 complex of SARS-CoV-2 reveal that our algorithm significantly increases the robustness and power of a constructed MSM without disturbing the Markovian properties. Furthermore, the superiority of our algorithm is amplified for slow dynamic modes in complex biological processes.

Structural study of the uPA-nafamostat complex reveals a covalent inhibitory mechanism of nafamostat

Nafamostat mesylate (NM) is a synthetic compound that inhibits various serine proteases produced during the coagulation cascade and inflammation. Previous studies showed that NM was a highly safe drug for the treatment of different cancers, but the precise functions and mechanisms of NM are not clear. In this study, we determined a series of crystal structures of NM and its hydrolysates in complex with a serine protease (urokinase-type plasminogen activator [uPA]). These structures reveal that NM was cleaved by uPA and that a hydrolyzed product (4-guanidinobenzoic acid [GBA]) remained covalently linked to Ser195 of uPA, and the other hydrolyzed product (6-amidino-2-naphthol [6A2N]) released from uPA. Strikingly, in the inactive uPA (uPA-S195A):NM structure, the 6A2N side of intact NM binds to the specific pocket of uPA. Molecular dynamics simulations and end-point binding free-energy calculations show that the conf1 of NM (6A2N as P1 group) in the uPA-S195A:NM complex may be more stable than conf2 of NM (GBA as P1 group). Moreover, in the structure of uPA:NM complex, the imidazole group of His57 flips further away from Ser195 and disrupts the stable canonical catalytic triad conformation. These results not only reveal the inhibitory mechanism of NM as an efficient serine protease inhibitor but also might provide the structural basis for the further development of serine protease inhibitors.

Conformational Stability and Denaturation Processes of Proteins Investigated by Electrophoresis under Extreme Conditions

The functional structure of proteins results from marginally stable folded conformations. Reversible unfolding, irreversible denaturation, and deterioration can be caused by chemical and physical agents due to changes in the physicochemical conditions of pH, ionic strength, temperature, pressure, and electric field or due to the presence of a cosolvent that perturbs the delicate balance between stabilizing and destabilizing interactions and eventually induces chemical modifications. For most proteins, denaturation is a complex process involving transient intermediates in several reversible and eventually irreversible steps. Knowledge of protein stability and denaturation processes is mandatory for the development of enzymes as industrial catalysts, biopharmaceuticals, analytical and medical bioreagents, and safe industrial food. Electrophoresis techniques operating under extreme conditions are convenient tools for analyzing unfolding transitions, trapping transient intermediates, and gaining insight into the mechanisms of denaturation processes. Moreover, quantitative analysis of electrophoretic mobility transition curves allows the estimation of the conformational stability of proteins. These approaches include polyacrylamide gel electrophoresis and capillary zone electrophoresis under cold, heat, and hydrostatic pressure and in the presence of non-ionic denaturing agents or stabilizers such as polyols and heavy water. Lastly, after exposure to extremes of physical conditions, electrophoresis under standard conditions provides information on irreversible processes, slow conformational drifts, and slow renaturation processes. The impressive developments of enzyme technology with multiple applications in fine chemistry, biopharmaceutics, and nanomedicine prompted us to revisit the potentialities of these electrophoretic approaches. This feature review is illustrated with published and unpublished results obtained by the authors on cholinesterases and paraoxonase, two physiologically and toxicologically important enzymes.

Solvent paramagnetic relaxation enhancement as a versatile method for studying structure and dynamics of biomolecular systems

Solvent paramagnetic relaxation enhancement (sPRE) is a versatile nuclear magnetic resonance (NMR)-based method that allows characterization of the structure and dynamics of biomolecular systems through providing quantitative experimental information on solvent accessibility of NMR-active nuclei. Addition of soluble paramagnetic probes to the solution of a biomolecule leads to paramagnetic relaxation enhancement in a concentration-dependent manner. Here we review recent progress in the sPRE-based characterization of structural and dynamic properties of biomolecules and their complexes, and aim to deliver a comprehensive illustration of a growing number of applications of the method to various biological systems. We discuss the physical principles of sPRE measurements and provide an overview of available co-solute paramagnetic probes. We then explore how sPRE, in combination with complementary biophysical techniques, can further advance biomolecular structure determination, identification of interaction surfaces within protein complexes, and probing of conformational changes and low-population transient states, as well as deliver insights into weak, nonspecific, and transient interactions between proteins and co-solutes. In addition, we present examples of how the incorporation of solvent paramagnetic probes can improve the sensitivity of NMR experiments and discuss the prospects of applying sPRE to NMR metabolomics, drug discovery, and the study of intrinsically disordered proteins.

Protein Design: From the Aspect of Water Solubility and Stability

Water solubility and structural stability are key merits for proteins defined by the primary sequence and 3D-conformation. Their manipulation represents important aspects of the protein design field that relies on the accurate placement of amino acids and molecular interactions, guided by underlying physiochemical principles. Emulated designer proteins with well-defined properties both fuel the knowledge-base for more precise computational design models and are used in various biomedical and nanotechnological applications. The continuous developments in protein science, increasing computing power, new algorithms, and characterization techniques provide sophisticated toolkits for solubility design beyond guess work. In this review, we summarize recent advances in the protein design field with respect to water solubility and structural stability. After introducing fundamental design rules, we discuss the transmembrane protein solubilization and de novo transmembrane protein design. Traditional strategies to enhance protein solubility and structural stability are introduced. The designs of stable protein complexes and high-order assemblies are covered. Computational methodologies behind these endeavors, including structure prediction programs, machine learning algorithms, and specialty software dedicated to the evaluation of protein solubility and aggregation, are discussed. The findings and opportunities for Cryo-EM are presented. This review provides an overview of significant progress and prospects in accurate protein design for solubility and stability.

Iterative Landmark-Based Umbrella Sampling (ILBUS) Protocol for Sampling of Conformational Space of Biomolecules

Computer simulations of biomolecules such as molecular dynamics often suffer from insufficient sampling. Due to limited computational resources, insufficient sampling prevents obtaining proper equilibrium distributions of observed properties. To deal with this problem, we proposed a simulation protocol for efficient resampling of collected off-equilibrium trajectories. These trajectories are utilized for the initial mapping of the conformational space, which is later properly resampled by the introduced Iterative Landmark-Based Umbrella Sampling (ILBUS) method. Reconstruction of static equilibrium properties is achieved by the multistate Bennett acceptance ratio (MBAR) method, which enables efficient use of simulated data. The ILBUS protocol is geometry-based and does not demand any additional collective variable or a dimensional-reduction technique. The only requirement is a set of suitably spaced reference conformations, which serve as landmarks in the mapped conformational space. Additionally, the ILBUS protocol encompasses an iterative process that optimizes the force constant used in the umbrella sampling simulation. Such tuning is an inherent feature of the protocol and does not need to be performed by the user in advance. Furthermore, even the simulations with suboptimal force constants can be used in estimates by MBAR. We demonstrate the feasibility and the performance of this approach in the study of the conformational landscape of the alanine dipeptide, met-enkephalin, and adenylate kinase.

Deep learning-based molecular dynamics simulation for structure-based drug design against SARS-CoV-2

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus type 2 (SARS-CoV-2), has led to a global pandemic. Deep learning (DL) technology and molecular dynamics (MD) simulation are two mainstream computational approaches to investigate the geometric, chemical and structural features of protein and guide the relevant drug design. Despite a large amount of research papers focusing on drug design for SARS-COV-2 using DL architectures, it remains unclear how the binding energy of the protein-protein/ligand complex dynamically evolves which is also vital for drug development. In addition, traditional deep neural networks usually have obvious deficiencies in predicting the interaction sites as protein conformation changes. In this review, we introduce the latest progresses of the DL and DL-based MD simulation approaches in structure-based drug design (SBDD) for SARS-CoV-2 which could address the problems of protein structure and binding prediction, drug virtual screening, molecular docking and complex evolution. Furthermore, the current challenges and future directions of DL-based MD simulation for SBDD are also discussed.

Fundamental Redesign of the TIGER2hs Kernel to Address Severe Parameter Sensitivity

Replica exchange molecular dynamics simulations are one of the most popular approaches to enhance conformational sampling of molecular systems. Applications range from protein folding to protein-protein or other host-guest interactions, as well as binding free energy calculations. While these methods are computationally expensive, highly accurate results can be obtained. We recently developed TIGER2hs, an improved version of the temperature intervals with global exchange of replicas (TIGER2) algorithm. This method combines the replica-based enhanced sampling in an explicit solvent with a hybrid solvent energy evaluation. During the exchange attempts, bulk water is replaced by an implicit solvent model, allowing sampling with significantly less replicas than parallel tempering (REMD). This enables accurate enhanced sampling calculations with only a fraction of computational resources compared to REMD. Our latest results highlight several issues with sampling imbalance and parameter sensitivity within the original TIGER2 exchange algorithms that affect the overall state populations. A high sensitivity on replica number and maximum temperature is eliminated by changing to a pairwise exchange kernel (PE) without additional sorting. Simulations are controlled by adjusting the average temperature change per exchange ⟨ΔT/χ⟩ to below 30 K to mimic a controlled temperature mixing of replicas similar to REMD. Thus, this parameter provides an applicable property for selecting combinations of replica number and maximum temperature to adjust simulations for best accuracy, with flexible resource investment. This increases the robustness of the method and ensures results in excellent agreement with REMD, as demonstrated for three different peptides.

Polymorphisms in common antihypertensive targets: Pharmacogenomic implications for the treatment of cardiovascular disease

The idea of personalized medicine came to fruition with sequencing the human genome; however, aside from a few cases, the genetic revolution has yet to materialize. Cardiovascular diseases are the leading cause of death globally, and hypertension is a common prelude to nearly all cardiovascular diseases. Thus, hypertension is an ideal candidate disease to apply tenants of personalized medicine to lessen cardiovascular disease. Herein is a survey that visually depicts the polymorphisms in the top eight antihypertensive targets. Although there are numerous genome-wide association studies regarding cardiovascular disease, few studies look at the effects of receptor polymorphisms on drug treatment. With 17,000+ polymorphisms in the combined target proteins examined, it is expected that some of the clinical variability in the treatment of hypertension is due to polymorphisms in the drug targets. Recent advances in techniques and technology, such as high throughput examination of single mutations, structure prediction, computational power for modeling, and CRISPR models of point mutations, allow for a relatively rapid and comprehensive examination of the effects of known and future polymorphisms on drug affinity and effects. As hypertension is easy to measure and has a plethora of clinically viable ligands, hypertension makes an excellent disease to study pharmacogenomics in the lab and the clinic. If the promises of personalized medicine are to materialize, a concerted effort to examine the effects polymorphisms have on drugs is required. A clinician with the knowledge of a patient's genotype can then prescribe drugs that are optimal for treating that specific patient.

Structural and Biophysical Insights into SPINK1 Bound to Human Cationic Trypsin

(1) The serine protease inhibitor Kazal type 1 (SPINK1) inhibits trypsin activity in zymogen granules of pancreatic acinar cells. Several mutations in the SPINK1 gene are associated with acute recurrent pancreatitis (ARP) and chronic pancreatitis (CP). The most common variant is SPINK1 p.N34S. Although this mutation was identified two decades ago, the mechanism of action has remained elusive. (2) SPINK1 and human cationic trypsin (TRY1) were expressed in E. coli, and inhibitory activities were determined. Crystals of SPINK1-TRY1 complexes were grown by using the hanging-drop method, and phases were solved by molecular replacement. (3) Both SPINK1 variants show similar inhibitory behavior toward TRY1. The crystal structures are almost identical, with minor differences in the mutated loop. Both complexes show an unexpected rotamer conformation of the His63 residue in TRY1, which is a member of the catalytic triad. (4) The SPINK1 p.N34S mutation does not affect the inhibitory behavior or the overall structure of the protein. Therefore, the pathophysiological mechanism of action of the p.N34S variant cannot be explained mechanistically or structurally at the protein level. The observed histidine conformation is part of a mechanism for SPINK1 that can explain the exceptional proteolytic stability of this inhibitor.