Use of caps in the auxiliary basis set formulation of the fragment molecular orbital method

An auxiliary polarization formulation of the fragment molecular orbital (FMO) method is developed, combining a basis set correction computed for capped isolated fragments with a polarization obtained from uncapped fragments. For a set of organic and inorganic test systems, it is shown that the total energy and atomic charges are accurately reproduced with respect to full unfragmented calculations. It is demonstrated that the method is accurate for computing electronic excited states. The developed approach is applied to rank the isomers of chignolin from experimental NMR data (PDB: 1UAO) according to their relative energy. Contributions of polarization and basis set effects to pair interactions between fragments are elucidated.

Hydrophobic substrate binding pocket remodeling of echinocandin B deacylase based on multi-dimensional rational design

Actinoplanes utahensis deacylase (AAC)-catalyzed deacylation of echinocandin B (ECB) is a promising method for the synthesis of anidulafungin, the newest of the echinocandin antifungal agents. However, the low activity of AAC significantly limits its practical application. In this work, we have devised a multi-dimensional rational design strategy for AAC, conducting separate analyses on the substrate-binding pocket's volume, curvature, and length. Furthermore, we quantitatively analyzed substrate properties, particularly on hydrophilic and hydrophobic. Accordingly, we tailored the linoleic acid-binding pocket of AAC to accommodate the extended long lipid chain of ECB. By fine-tuning the key residues, the resulting AAC mutants can accommodate the ECB lipid chain with a lower curvature binding pocket. The D53A/I55F/G57M/F154L/Q661L mutant (MT) displayed 331 % higher catalytic efficiency than the wild-type (WT) enzyme. The MT product conversion was 94.6 %, reaching the highest reported level. Utilizing a multi-dimensional rational design for a customized mutation strategy of the substrate-binding pocket is an effective approach to enhance the catalytic efficiency of enzymes in handling complicated substrates.

Promiscuous recognition of MR1 drives self-reactive mucosal-associated invariant T cell responses

Mucosal-associated invariant T (MAIT) cells use canonical semi-invariant T cell receptors (TCR) to recognize microbial riboflavin precursors displayed by the antigen-presenting molecule MR1. The extent of MAIT TCR crossreactivity toward physiological, microbially unrelated antigens remains underexplored. We describe MAIT TCRs endowed with MR1-dependent reactivity to tumor and healthy cells in the absence of microbial metabolites. MAIT cells bearing TCRs crossreactive toward self are rare but commonly found within healthy donors and display T-helper-like functions in vitro. Experiments with MR1-tetramers loaded with distinct ligands revealed significant crossreactivity among MAIT TCRs both ex vivo and upon in vitro expansion. A canonical MAIT TCR was selected on the basis of extremely promiscuous MR1 recognition. Structural and molecular dynamic analyses associated promiscuity to unique TCRβ-chain features that were enriched within self-reactive MAIT cells of healthy individuals. Thus, self-reactive recognition of MR1 represents a functionally relevant indication of MAIT TCR crossreactivity, suggesting a potentially broader role of MAIT cells in immune homeostasis and diseases, beyond microbial immunosurveillance.

Large scale peptide screening against main protease of SARS CoV-2

The COVID-19 pandemic has been a public health emergency, with deadly forms constantly emerging around the world, highlighting the dire need for highly effective antiviral therapeutics. Peptide therapeutics show significant potential for this viral disease due to their efficiency, safety, and specificity. Here, two thousand seven hundred eight antibacterial peptides were screened computationally targeting the Main protease (Mpro) of SARS CoV-2. Six top-ranked peptides according to their binding scores, binding pose were investigated by molecular dynamics to explore the interaction and binding behavior of peptide-Mpro complexes. The structural and energetic characteristics of Mpro-DRAMP01760 and Mpro-DRAMP01808 complexes fluctuated less during a 250 ns MD simulation. In addition, three peptides (DRAMP01760, DRAMP01808, and DRAMP01342) bind strongly to Mpro protein, according to the free energy landscape and principal component analysis. Peptide helicity and secondary structure analysis are in agreement with our findings. Interaction analysis of protein-peptide complexes demonstrated that Mpro's residue CYS145, HIS41, PRO168, GLU166, GLN189, ASN142, MET49, and THR26 play significant contributions in peptide-protein attachment. Binding free energy analysis (MM-PBSA) demonstrated the energy profile of interacting residues of Mpro in peptide-Mpro complexes. To summarize, the peptides DRAMP01808 and DRAMP01760 may be highly Mpro specific, resulting disruption in a viral replication and transcription. The results of this research are expected to assist future research toward the development of antiviral peptide-based therapeutics for Covid-19 treatment.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

Computational analysis of binding free energies, hotspots and the binding mechanism of Bcl-xL/Bcl-2 binding to Bad/Bax

Hierarchical clustering tree of residues providing contributions to system binding based on the binding free energy of specific residues for (A) Bcl-xL systems (B) Bcl-2 systems.

An accurate free energy estimator: based on MM/PBSA combined with interaction entropy for protein-ligand binding affinity

The molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) method is constantly used to calculate the binding free energy of protein-ligand complexes, and has been shown to effectively balance computational cost against accuracy. The relative binding affinities obtained by the MM/PBSA approach are acceptable, while it usually overestimates the absolute binding free energy. This paper proposes four free energy estimators based on the MM/PBSA for enthalpy change combined with interaction entropy (IE) for entropy change using different weights for individual energy terms. The ΔG_{PBSA_IE} method is determined to be an optimal estimator based on its performance in terms of the correlation between experimental and theoretical values and error estimations. This approach is optimized using high-quality experimental values from a training set containing 84 protein-ligand systems, and the coefficients for the sum of electrostatic energy and polar solvation free energy, van der Waals (vdW) energy, non-polar solvation energy and entropy change are obtained by multivariate linear fitting to the corresponding experimental values. A comparison between the traditional MM/PBSA method and this method shows that the correlation coefficient is improved from 0.46 to 0.72 and the slope of the regression line increases from 0.10 to 1.00. More importantly, the mean absolute error (MAE) is significantly reduced from 22.52 to 1.59 kcal mol^-1. Furthermore, the numerical stability of this method is validated on a test set with a similar correlation coefficient, slope and MAE to those of the training set. Based on the above advantages, the ΔG_{PBSA_IE} method can be a powerful tool for a reliable and accurate estimation of binding free energy and plays a significant role in a detailed energetic investigation of protein-ligand interaction.

Improving the performance of the MM/PBSA and MM/GBSA methods in recognizing the native structure of the Bcl-2 family using the interaction entropy method

Correct discrimination of native structure plays an important role in drug design. IE method significantly improves the performance of MM/PB(GB)SA method in discriminating native and decoy structures in protein–ligand/protein systems of Bcl-2 family.