Statistical Measures to Quantify Similarity between Molecular Dynamics Simulation Trajectories

Design a novel of Brucellosis preventive vaccine based on IgV_CTLA-4 and multiple epitopes via immunoinformatics approach

Brucellosis is a zoonotic disease caused by Brucella, which is difficult to be eliminated by conventional drugs. Therefore, a novel multi-epitope vaccine (MEV) was designed to prevent human Brucella infection. Based on the method of "reverse vaccinology", cytotoxic T lymphocytes epitopes (CTLEs), helper T lymphocytes epitopes (HTLEs), linear B-cell epitopes (LBEs) and conformational B-cell epitopes (CBEs) of four Brucella proteins (VirB9, VirB10, Omp 19 and Omp 25) were obtained. In order to keep the correct protein folding, the multiple epitopes was constructed by connecting epitopes through linkers. In view of the significant connection between human leukocyte antigen CTLA-4 and B7 molecules found on antigen presenting cells (APCs), a new vaccine (V_C4MEV) for preventing brucellosis was created by combining CTLA-4 immunoglobulin variable region (IgV_CTLA-4) with MEV protein. Immunoinformatics analysis showed that V_C4MEV had a good secondary and tertiary structure. Additionally, molecular docking and molecular dynamics simulation (MD) revealed a robust binding affinity between IgV_ CTLA-4 and the B7 molecule. Notably, the vaccine V_C4MEV was demonstrated favorable immunogenicity and antigenicity in both in vitro and in vivo experiments. V_C4MEV had the potential to activate defensive cells and immune responses, offering a hopeful approach for developing vaccines against Brucella in the upcoming years.

An in-silico receptor-pharmacophore based multistep molecular docking and simulation study to evaluate the inhibitory potentials against NS1 of DENV-2

DENV-2 strain is the most fatal and infectious of the five dengue virus serotypes. The non-structural protein NS1 encoded by its genome is the most significant protein required for viral pathogenesis and replication inside the host body. Thus, targeting the NS1 protein and designing an inhibitor to limit its stability and secretion is a propitious attempt in our fight against dengue. Four novel inhibitors are designed to target the conserved cysteine residues (C55, C313, C316, and C329) and glycosylation sites (N130 and N207) of the NS1 protein in an attempt to halt the spread of the dengue infection in the host body altogether. Numerous computer-aided drug designing techniques including molecular docking, molecular dynamics simulation, virtual screening, principal component analysis, and dynamic cross-correlation matrix were employed to determine the structural and functional activity of the NS1-inhibitor complexes. From our analysis, it was evident that the extent of structural and atomic level fluctuations of the ligand-bound protein exhibited a declining trend in contrast to unbound protein which was prominently noticeable through the RMSD, RMSF, Rg, and SASA graphs. The ADMET analysis of the four ligands revealed a promising pharmacokinetics and pharmacodynamic profile, along with good bioavailability and toxicity properties. The proposed drugs when bound to the targeted cavities resulted in stable conformations in comparison to their unbound state, implying they have good affinity promising effective drug action. Thus, they can be tested in vitro and used as potential anti-dengue drugs.Communicated by Ramaswamy H. Sarma.

Contriving a novel of CHB therapeutic vaccine based on IgV_CTLA-4 and L protein via immunoinformatics approach

Chronic infection induced by immune tolerance to hepatitis B virus (HBV) is one of the most common causes of hepatic cirrhosis and hepatoma. Fortunately, the application of therapeutic vaccine can not only reverse HBV-tolerance, but also serve a potentially effective therapeutic strategy for treating chronic hepatitis B (CHB). However, the clinical effect of the currently developed CHB therapeutic vaccine is not optimistic due to the weak immunogenicity. Given that the human leukocyte antigen CTLA-4 owns strong binding ability to the surface B7 molecules (CD80 and CD86) of antigen presenting cell (APCs), the immunoglobulin variable region of CTLA-4 (IgV_CTLA-4) was fused with the L protein of HBV to contrive a novel therapeutic vaccine (V_C4HBL) for CHB in this study. We found that the addition of IgV_CTLA-4 did not interfere with the formation of L protein T cell and B cell epitopes after analysis by means of immunoinformatics approaches. Meanwhile, we also found that the IgV_CTLA-4 had strong binding force to B7 molecules through molecular docking and molecular dynamics (MD) simulation. Notably, our vaccine V_C4HBL showed good immunogenicity and antigenicity by in vitro and in vivo experiments. Therefore, the V_C4HBL is promising to again effectively activate the cellular and humoral immunity of CHB patients, and provides a potentially effective therapeutic strategy for the treatment of CHB in the future.Communicated by Ramaswamy H. Sarma.

Machine Learning Integrating Protein Structure, Sequence, and Dynamics to Predict the Enzyme Activity of Bovine Enterokinase Variants

Despite recent advances in computational protein science, the dynamic behavior of proteins, which directly governs their biological activity, cannot be gleaned from sequence information alone. To overcome this challenge, we propose a framework that integrates the peptide sequence, protein structure, and protein dynamics descriptors into machine learning algorithms to enhance their predictive capabilities and achieve improved prediction of the protein variant function. The resulting machine learning pipeline integrates traditional sequence and structure information with molecular dynamics simulation data to predict the effects of multiple point mutations on the fold improvement of the activity of bovine enterokinase variants. This study highlights how the combination of structural and dynamic data can provide predictive insights into protein functionality and address protein engineering challenges in industrial contexts.

A. E. Albadri A, Patel H, Messaoudi S, Kadri A, M. Al-Hazmy S, Aouadi K. Synthesis, optical properties, DNA, β-cyclodextrin interaction, hydrogen isotope sensor and computational study of new enantiopure isoxazolidine derivative (ISoXD)

A novel isoxazolidine derivative (ISoXD) dye was successfully synthesized and comprehensively characterized. In this study, we conducted a thorough examination of its various properties, including optical characteristics, interactions with DNA and β-cyclodextrin (β-CD), molecular docking, molecular dynamic simulation, and density functional theory (DFT) calculations. Our investigation encompassed a systematic analysis of the absorption and emission spectra of ISoXD in diverse solvents. The observed variations in the spectroscopic data were attributed to the specific solvent's capacity to engage in hydrogen bonding interactions. Remarkably, the most pronounced intensities were observed in glycol, which can establish many hydrogen bonds with ISoXD. Furthermore, our study revealed a significant distinction in the fluorescence behavior of ISoXD when subjected to different solvents, particularly between CHCl3 and CDCl3. Moreover, we explored the fluorescence intensity of the ISoXD complex in the presence of various metals, both in ethanol and water. The ISoXD complex exhibited a substantial increase of fluorescence upon interaction with different metal ions. The utilization of DFT calculations allowed us to propose an intramolecular charge transfer (ICT) mechanism as a plausible explanation for this quenching phenomenon. The interaction of ISoXD with DNA and β-CD was studied using absorption spectra. The binding constant (K) and the standard Gibbs free energy change (ΔGo) for the interaction between DNA and β-CD with ISoXD were determined. In docking study, ISoXD exhibited significant docking scores (-6.511) and MM-GBSA binding free energies (-66.27 kcal/mol) within the PARP-1 binding cavity. Its binding pattern closely resembles to the co-crystal ligand veliparib, and during a 100ns MD simulation, ISoXD displayed strong stability and formed robust hydrogen bonds with key amino acids. These findings suggest ISoXD's potential as a PARP-1 inhibitor for further investigation in therapeutic development.

Elucidation of Molecular Interactions Between Drug-Polymer in Amorphous Solid Dispersion by a Computational Approach Using Molecular Dynamics Simulations

Introduction

Amorphous drug dispersion is frequently used to enhance the solubility and dissolution of poorly water-soluble drugs, thereby improving their oral bioavailability. The dispersion of these drugs into polymer matrix can inhibit their recrystallization. The inter-molecular interactions between drug and polymer plays a role in the improvement of the dissolution rate, solubility, and physical stability of drug.

Aim

This study aims to investigate the formation and interactions of ritonavir (RTV)/poloxamer (PLX) amorphous formulation using a computational approach via molecular dynamics (MD) simulations, which mimicked solvent evaporation and melt-quenching method.

Methods

TheRoot Mean Square Deviation (RMSD) value, Root Mean Square Fluctuation (RMSF), Radial Distribution Function (RDF), Radius of Gyration (Rg), Solvent Accessible Surface Area (SASA), and hydrogen bond interactions were analyzed to determine interaction mechanisms between RTV and PLX in amorphous solid dispersion.

Results

The pi-alkyl bonds between RTV and PLX were formed after simulations of solvent evaporation, while the hydrogen bond interactions of RTV-PLX was observed during melt method simulations. These results indicate the successful formulation of amorphous solid dispersion (ASD) from RTV and PLX. The RMSD values obtained from the solvent evaporation, melt-cooling-A, melt-cooling-B, and melt-cooling-C methods were 3.33 Å, 1.97 Å, 1.30 Å, and 1.29 Å, respectively, while the average RMSF values were 2.65 Å, 1.04 Å, 1.05 Å, and 1.07 Å, respectively. This indicates that the suppression of translational motion of RTV from the melt method can be stronger than solvent evaporation caused by the intermolecular interactions of RTV-PLX.

Conclusion

MD simulations helped in understanding the formation and interaction mechanisms of ASD formulations that were difficult to detect by experimental approaches.

Inhibition potential of natural flavonoids against selected omicron (B.1.19) mutations in the spike receptor binding domain of SARS-CoV-2: a molecular modeling approach

The omicron (B.1.19) variant of contagious severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is considered a variant of concern (VOC) due to its increased transmissibility and highly infectious nature. The spike receptor-binding domain (RBD) is a hotspot of mutations and is regarded as a prominent target for screening drug candidates owing to its crucial role in viral entry and immune evasion. To date, no effective therapy or antivirals have been reported; therefore, there is an urgent need for rapid screening of antivirals. An extensive molecular modelling study has been performed with the primary goal to assess the inhibition potential of natural flavonoids as inhibitors against RBD from a manually curated library. Out of 40 natural flavonoids, five natural flavonoids, namely tomentin A (-8.7 kcal/mol), tomentin C (-8.6 kcal/mol), hyperoside (-8.4 kcal/mol), catechin gallate (-8.3 kcal/mol), and corylifol A (-8.2 kcal/mol), have been considered as the top-ranked compounds based on their binding affinity and molecular interaction profiling. The state-of-the-art molecular dynamics (MD) simulations of these top-ranked compounds in complex with RBD exhibited stable dynamics and structural compactness patterns on 200 nanoseconds. Additionally, complexes of these molecules demonstrated favorable free binding energies and affirmed the docking and simulation results. Moreover, the post-simulation validation of these interacted flavonoids using principal component analysis (PCA) revealed stable interaction patterns with RBD. The integrated results suggest that tomentin A, tomentin C, hyperoside, catechin gallate, and corylifol A might be effective against the emerging variants of SARS-CoV-2 and should be further evaluated using in-vitro and in-vivo experiments.Communicated by Ramaswamy H. Sarma.

Molecular modelling of SdiA protein by selected flavonoid and terpenes compounds to attenuate virulence in Klebsiella pneumoniae

Klebsiella pneumoniae is one of the perturbing multidrug resistant (MDR) and ESKAPE pathogens contributing to the mounting morbidity, mortality and extended rate of hospitalization. Its virulence, often regulated by quorum sensing (QS) reinforces the need to explore alternative and prospective antivirulence agents, relatively from plants secondary metabolites. Computer aided drug discovery using molecular modelling techniques offers advantage to investigate prospective drugs to combat MDR pathogens. Thus, this study employed virtual screening of selected terpenes and flavonoids from medicinal plants to interrupt the QS associated SdiA protein in K. pneumoniae to attenuate its virulence. 4LFU was used as a template to model the structure of SdiA. ProCheck, Verify3D, Ramachandran plot scores, and ProSA-Web all attested to the model's good quality. Since SdiA protein in K. pneumoniae leads to the expression of virulence, 31 prospective bioactive compounds were docked for antagonistic potential. The stability of the protein-ligand complex, atomic motions and inter-atomic interactions were further investigated through molecular dynamics simulations (MDS) at 100 ns production runs. The binding free energy was estimated using the molecular mechanics/poisson-boltzmann surface area (MM/PB-SA). Furthermore, the drug-likeness properties of the studied compounds were validated. Docking studies showed phytol possesses the highest binding affinity (-9.205 kcal/mol) while glycitein had -9.752 kcal/mol highest docking score. The MDS of the protein in complex with the best-docked compounds revealed phytol with the highest binding energy of -44.2625 kcal/mol, a low root-mean-square deviation (RMSD) value of 1.54 Å and root-mean-square fluctuation (RMSF) score of 1.78 Å. Analysis of the drug-likeness properties prediction and bioavailability of these compounds revealed their conformed activity to lipinski's rules with bioavailability scores of 0.55 F. The studied terpenes and flavonoids compounds effectively thwart SdiA protein, therefore regulate inter- or intra cellular communication and associated in virulence Enterobacteriaceae, serving as prospective antivirulence drugs.Communicated by Ramaswamy H. Sarma.

Exploration of CviR-mediated quorum sensing inhibitors from Cladosporium spp. against Chromobacterium violaceum through computational studies

An opportunistic human pathogenic bacterium, Chromobacterium violaceum resists the potency of most antibiotics by exploiting the quorum sensing system within their community to control virulence factor expression. Therefore, blocking the quorum sensing mechanism could help to treat several infectious caused by this organism. The quorum sensing receptor (CviR) of C. violaceum was used as a model target in the current investigation to identify potentially novel quorum sensing inhibitors from Cladosporium spp. through in silico computational approaches. The molecular docking results confirmed the anti-quorum sensing potential of bioactive compounds from Cladosporium spp. through binding to CviR with varying docking scores between - 5.2 and - 9.5 kcal/mol. Relative to the positive control [Azithromycin (- 7.4 kcal/mol)], the top six metabolites of Cladosporium spp. had higher docking scores and were generally greater than - 8.5 kcal/mol. The thermodynamic stability and binding affinity refinement of top-ranked CviR inhibitors were further studied through a 160 ns molecular dynamic (MD) simulation. The Post-MD simulation analysis confirmed the top-ranked compounds' affinity, stability, and biomolecular interactions with CviR at 50 ns, 100 ns, and 160 ns with Coniochaetone K of the Cladosporium spp. having the highest binding free energy (- 30.87 kcal/mol) and best interactions (two consistent hydrogen bond contact) following the 160 ns simulation. The predicted pharmacokinetics properties of top selected compounds point to their drug likeliness, potentiating their chance as a possible drug candidate. Overall, the top-ranked compounds from Cladosporium spp., especially Coniochaetone K, could be identified as potential C. violaceum CviR inhibitors. The development of these compounds as broad-spectrum antibacterial medicines is thus possible in the future following the completion of further preclinical and clinical research.

Insilico generation of novel ligands for the inhibition of SARS-CoV-2 main protease (3CLpro) using deep learning

The recent emergence of novel severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) causing the coronavirus disease (COVID-19) has become a global public health crisis, and a crucial need exists for rapid identification and development of novel therapeutic interventions. In this study, a recurrent neural network (RNN) is trained and optimized to produce novel ligands that could serve as potential inhibitors to the SARS-CoV-2 viral protease: 3 chymotrypsin-like protease (3CL^pro). Structure-based virtual screening was performed through molecular docking, ADMET profiling, and predictions of various molecular properties were done to evaluate the toxicity and drug-likeness of the generated novel ligands. The properties of the generated ligands were also compared with current drugs under various phases of clinical trials to assess the efficacy of the novel ligands. Twenty novel ligands were selected that exhibited good drug-likeness properties, with most ligands conforming to Lipinski's rule of 5, high binding affinity (highest binding affinity: -9.4 kcal/mol), and promising ADMET profile. Additionally, the generated ligands complexed with 3CL^pro were found to be stable based on the results of molecular dynamics simulation studies conducted over a 100 ns period. Overall, the findings offer a promising avenue for the rapid identification and development of effective therapeutic interventions to treat COVID-19.

Novel Thiosemicarbazone Quantum Dots in the Treatment of Alzheimer's Disease Combining In Silico Models Using Fingerprints and Physicochemical Descriptors

Searching for thiosemicarbazone derivatives with the potential to inhibit acetylcholinesterase for the treatment of Alzheimer's disease (AD) is an important current goal. The QSAR_KPLS, QSAR_ANN, and QSAR_SVR models were constructed using binary fingerprints and physicochemical (PC) descriptors of 129 thiosemicarbazone compounds screened from a database of 3791 derivatives. The R ² and Q ² values for the QSAR_KPLS, QSAR_ANN, and QSAR_SVR models are greater than 0.925 and 0.713 using dendritic fingerprint (DF) and PC descriptors, respectively. The in vitro pIC₅₀ activities of four new design-oriented compounds N1, N2, N3, and N4, from the QSAR_KPLS model using DFs, are consistent with the experimental results and those from the QSAR_ANN and QSAR_SVR models. The designed compounds N1, N2, N3, and N4 do not violate Lipinski-5 and Veber rules using the ADME and BoiLED-Egg methods. The binding energy, kcal mol^-1, of the novel compounds to the 1ACJ-PDB protein receptor of the AChE enzyme was also obtained by molecular docking and dynamics simulations consistent with those predicted from the QSAR_ANN and QSAR_SVR models. New compounds N1, N2, N3, and N4 were synthesized, and the experimental in vitro pIC₅₀ activity was determined in agreement with those obtained from in silico models. The newly synthesized thiosemicarbazones N1, N2, N3, and N4 can inhibit 1ACJ-PDB, which is predicted to be able to cross the barrier. The DFT B3LYP/def-SV(P)-ECP quantization calculation method was used to calculate E _HOMO and E _LUMO to account for the activities of compounds N1, N2, N3, and N4. The quantum calculation results explained are consistent with those obtained in in silico models. The successful results here may contribute to the search for new drugs for the treatment of AD.

Docking and Molecular Dynamics Identify Leads against 5 Alpha Reductase 2 for Benign Prostate Hyperplasia Treatment

Steroid 5 alpha-reductase 2 (5αR-2) is a membrane-embedded protein that together with other isoforms plays a key role in the metabolism of steroids. This enzyme catalyzes the reduction of testosterone to the more potent ligand, dihydrotestosterone (DHT) in the prostate. Androgens, testosterone, and DHT play important roles in prostate growth, development, and function. At the same time, both testosterone and DHT have been implicated in the pathogenesis of benign prostate hyperplasia (BPH). Inhibition of the DHT formation, therefore, provides a therapeutic strategy that offers the possibility of preventing, delaying, or treating BPH. Currently, two steroidal drugs that inhibit 5αR-2, dutasteride and finasteride, have been approved for clinical use. These two come at a high cost and also portray undesirable sexual side effects which necessitate the need to find new chemotherapeutic alternatives for the disease. Based on the aforementioned, finasteride and dutasteride were subjected to scaffold hopping, fragment-based de novo design, molecular docking, and molecular dynamics simulations employing databases like ChEMBL, DrugBank, PubChem, ChemSpider, and Zinc15 in the identification of potential hits targeting 5αR-2. Altogether, ten novel compounds targeting 5αR-2 were identified with binding energies lower or comparable to finasteride and dutasteride, the main inhibitors for this target. Molecular docking and molecular dynamics simulations studies identify amino acid residues Glu57, Phe219, Phe223, and Leu224 to be critical for ligand binding and complex stability. The physicochemical and pharmacological profiling suggests the potential of the hit compounds to be drug-like and orally active. Similarly, the quality parameter assessments revealed the hits possess LELP greater than 3 implying their promise as lead-like molecules. The compounds A5, A9, and A10 were, respectively, predicted to treat prostate disorders with Pa (0.188, 0.361, and 0.270) and Pi (0.176, 0.050, and 0.093), while A8 and A9 were found to be associated with BPH treatment with Pa (0.09 and 0.127) and Pi (0.077 and 0.033), respectively. Structural similarity searches via DrugBank identified the drugs faropenem, acemetacin, estradiol valerate, and yohimbine to be useful for BPH treatment suggesting the de novo designed ligands as potential chemotherapeutic agents for treating this disease.

Homology Modeling, de Novo Design of Ligands, and Molecular Docking Identify Potential Inhibitors of Leishmania donovani 24-Sterol Methyltransferase

The therapeutic challenges pertaining to leishmaniasis due to reported chemoresistance and toxicity necessitate the need to explore novel pathways to identify plausible inhibitory molecules. Leishmania donovani 24-sterol methyltransferase (LdSMT) is vital for the synthesis of ergosterols, the main constituents of Leishmania cellular membranes. So far, mammals have not been shown to possess SMT or ergosterols, making the pathway a prime candidate for drug discovery. The structural model of LdSMT was elucidated using homology modeling to identify potential novel 24-SMT inhibitors via virtual screening, scaffold hopping, and de-novo fragment-based design. Altogether, six potential novel inhibitors were identified with binding energies ranging from −7.0 to −8.4 kcal/mol with e-LEA3D using 22,26-azasterol and S1–S4 obtained from scaffold hopping via the ChEMBL, DrugBank, PubChem, ChemSpider, and ZINC15 databases. These ligands showed comparable binding energy to 22,26-azasterol (−7.6 kcal/mol), the main inhibitor of LdSMT. Moreover, all the compounds had plausible ligand efficiency-dependent lipophilicity (LELP) scores above 3. The binding mechanism identified Tyr92 to be critical for binding, and this was corroborated via molecular dynamics simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) calculations. The ligand A1 was predicted to possess antileishmanial properties with a probability of activity (Pa) of 0.362 and a probability of inactivity (Pi) of 0.066, while A5 and A6 possessed dermatological properties with Pa values of 0.205 and 0.249 and Pi values of 0.162 and 0.120, respectively. Structural similarity search via DrugBank identified vabicaserin, daledalin, zanapezil, imipramine, and cefradine with antileishmanial properties suggesting that the de-novo compounds could be explored as potential antileishmanial agents.

Functional Dynamics of Substrate Recognition in TEM Beta-Lactamase

The beta-lactamase enzyme provides effective resistance to beta-lactam antibiotics due to substrate recognition controlled by point mutations. Recently, extended-spectrum and inhibitor-resistant mutants have become a global health problem. Here, the functional dynamics that control substrate recognition in TEM beta-lactamase are investigated using all-atom molecular dynamics simulations. Comparisons are made between wild-type TEM-1 and TEM-2 and the extended-spectrum mutants TEM-10 and TEM-52, both in apo form and in complex with four different antibiotics (ampicillin, amoxicillin, cefotaxime and ceftazidime). Dynamic allostery is predicted based on a quasi-harmonic normal mode analysis using a perturbation scan. An allosteric mechanism known to inhibit enzymatic function in TEM beta-lactamase is identified, along with other allosteric binding targets. Mechanisms for substrate recognition are elucidated using multivariate comparative analysis of molecular dynamics trajectories to identify changes in dynamics resulting from point mutations and ligand binding, and the conserved dynamics, which are functionally important, are extracted as well. The results suggest that the H10-H11 loop (residues 214-221) is a secondary anchor for larger extended spectrum ligands, while the H9-H10 loop (residues 194-202) is distal from the active site and stabilizes the protein against structural changes. These secondary non-catalytically-active loops offer attractive targets for novel noncompetitive inhibitors of TEM beta-lactamase.

JEDi: java essential dynamics inspector - a molecular trajectory analysis toolkit

Background

Principal component analysis (PCA) is commonly applied to the atomic trajectories of biopolymers to extract essential dynamics that describe biologically relevant motions. Although application of PCA is straightforward, specialized software to facilitate workflows and analysis of molecular dynamics simulation data to fully harness the power of PCA is lacking. The Java Essential Dynamics inspector (JEDi) software is a major upgrade from the previous JED software.

Results

Employing multi-threading, JEDi features a user-friendly interface to control rapid workflows for interrogating conformational motions of biopolymers at various spatial resolutions and within subregions, including multiple chain proteins. JEDi has options for Cartesian-based coordinates (cPCA) and internal distance pair coordinates (dpPCA) to construct covariance (Q), correlation (R), and partial correlation (P) matrices. Shrinkage and outlier thresholding are implemented for the accurate estimation of covariance. The effect of rare events is quantified using outlier and inlier filters. Applying sparsity thresholds in statistical models identifies latent correlated motions. Within a hierarchical approach, small-scale atomic motion is first calculated with a separate local cPCA calculation per residue to obtain eigenresidues. Then PCA on the eigenresidues yields rapid and accurate description of large-scale motions. Local cPCA on all residue pairs creates a map of all residue-residue dynamical couplings. Additionally, kernel PCA is implemented. JEDi output gives high quality PNG images by default, with options for text files that include aligned coordinates, several metrics that quantify mobility, PCA modes with their eigenvalues, and displacement vector projections onto the top principal modes. JEDi provides PyMol scripts together with PDB files to visualize individual cPCA modes and the essential dynamics occurring within user-selected time scales. Subspace comparisons performed on the most relevant eigenvectors using several statistical metrics quantify similarity/overlap of high dimensional vector spaces. Free energy landscapes are available for both cPCA and dpPCA.

Conclusion

JEDi is a convenient toolkit that applies best practices in multivariate statistics for comparative studies on the essential dynamics of similar biopolymers. JEDi helps identify functional mechanisms through many integrated tools and visual aids for inspecting and quantifying similarity/differences in mobility and dynamic correlations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04140-5.

Molecular function recognition by supervised projection pursuit machine learning

Identifying mechanisms that control molecular function is a significant challenge in pharmaceutical science and molecular engineering. Here, we present a novel projection pursuit recurrent neural network to identify functional mechanisms in the context of iterative supervised machine learning for discovery-based design optimization. Molecular function recognition is achieved by pairing experiments that categorize systems with digital twin molecular dynamics simulations to generate working hypotheses. Feature extraction decomposes emergent properties of a system into a complete set of basis vectors. Feature selection requires signal-to-noise, statistical significance, and clustering quality to concurrently surpass acceptance levels. Formulated as a multivariate description of differences and similarities between systems, the data-driven working hypothesis is refined by analyzing new systems prioritized by a discovery-likelihood. Utility and generality are demonstrated on several benchmarks, including the elucidation of antibiotic resistance in TEM-52 beta-lactamase. The software is freely available, enabling turnkey analysis of massive data streams found in computational biology and material science.

Machine Learning Classification Model for Functional Binding Modes of TEM-1 β-Lactamase

TEM family of enzymes is one of the most commonly encountered β-lactamases groups with different catalytic capabilities against various antibiotics. Despite the studies investigating the catalytic mechanism of TEM β-lactamases, the binding modes of these enzymes against ligands in different functional catalytic states have been largely overlooked. But the binding modes may play a critical role in the function and even the evolution of these proteins. In this work, a newly developed machine learning analysis approach to the recognition of protein dynamics states was applied to compare the binding modes of TEM-1 β-lactamase with regard to penicillin in different catalytic states. While conventional analysis methods, including principal components analysis (PCA), could not differentiate TEM-1 in different binding modes, the application of a machine learning method led to excellent classification models differentiating these states. It was also revealed that both reactant/product states and apo/product states are more differentiable than the apo/reactant states. The feature importance generated by the training procedure of the machine learning model was utilized to evaluate the contribution from residues at active sites and in different secondary structures. Key active site residues, Ser70 and Ser130, play a critical role in differentiating reactant/product states, while other active site residues are more important for differentiating apo/product states. Overall, this study provides new insights into the different dynamical function states of TEM-1 and may open a new venue for β-lactamases functional and evolutional studies in general.

On-target Resistance to the Mutant-Selective EGFR Inhibitor Osimertinib Can Develop in an Allele-Specific Manner Dependent on the Original EGFR-Activating Mutation

PURPOSE

The third-generation EGFR inhibitor, osimertinib, is the first mutant-selective inhibitor that has received regulatory approval for the treatment of patients with EGFR-mutant lung cancer. Despite the development of highly selective third-generation inhibitors, acquired resistance remains a significant clinical challenge. Recently, we and others have identified a novel osimertinib resistance mutation, G724S, which was not predicted in in vitro screens. Here, we investigate how G724S confers resistance to osimertinib.Experimental Design: We combine structure-based predictive modeling of G724S in combination with the 2 most common EGFR-activating mutations, exon 19 deletion (Ex19Del) and L858R, with in vitro drug-response models and patient genomic profiling.

RESULTS

Our simulations suggest that the G724S mutation selectively reduces osimertinib-binding affinity in the context of Ex19Del. Consistent with our simulations, cell lines transduced with Ex19Del/G724S demonstrate resistance to osimertinib, whereas cells transduced with L858R/G724S are sensitive to osimertinib. Subsequent clinical genomic profiling data further suggest G724S occurs with Ex19Del but not L858R. Furthermore, we demonstrate that Ex19Del/G724S retains sensitivity to afatinib, but not to erlotinib, suggesting a possible therapy for patients at the time of disease relapse.

CONCLUSIONS

Altogether, these data suggest that G724S is an allele-specific resistance mutation emerging in the context of Ex19Del but not L858R. Our results fundamentally reframe the problem of targeted therapy resistance from one focused on the "drug-resistance mutation" pair to one focused on the "activating mutation-drug-resistance mutation" trio. This has broad implications across clinical oncology.

High throughput nonparametric probability density estimation

In high throughput applications, such as those found in bioinformatics and finance, it is important to determine accurate probability distribution functions despite only minimal information about data characteristics, and without using human subjectivity. Such an automated process for univariate data is implemented to achieve this goal by merging the maximum entropy method with single order statistics and maximum likelihood. The only required properties of the random variables are that they are continuous and that they are, or can be approximated as, independent and identically distributed. A quasi-log-likelihood function based on single order statistics for sampled uniform random data is used to empirically construct a sample size invariant universal scoring function. Then a probability density estimate is determined by iteratively improving trial cumulative distribution functions, where better estimates are quantified by the scoring function that identifies atypical fluctuations. This criterion resists under and over fitting data as an alternative to employing the Bayesian or Akaike information criterion. Multiple estimates for the probability density reflect uncertainties due to statistical fluctuations in random samples. Scaled quantile residual plots are also introduced as an effective diagnostic to visualize the quality of the estimated probability densities. Benchmark tests show that estimates for the probability density function (PDF) converge to the true PDF as sample size increases on particularly difficult test probability densities that include cases with discontinuities, multi-resolution scales, heavy tails, and singularities. These results indicate the method has general applicability for high throughput statistical inference.