SARS-CoV-2-vaccination in primary humoral immunodeficiency - Experience from a German lung clinic.

During the COVID-19 pandemic the effectiveness of vaccines against SARS-CoV-2 in immunodeficient patients did not only affect the individual risk of these vulnerable patients but endangered the selection of new variants of concern due to prolonged virus shedding by these patients. In a tertiary center for pulmonary diseases, we investigated the immune response of 11 patients with primary humoral immunodeficiency and 13 healthy controls on the humoral and cellular level after full vaccination with a mRNA or vector vaccine against SARS-CoV-2. In the majority of patients (73%), we found antibodies against the Spike protein above the thresh-old of positivity.

The Role of General Acid Catalysis in the Mechanism of an Alkyl Transferase Ribozyme.

MTR1 is an in vitro-selected alkyl transferase ribozyme that transfers an alkyl group from -alkylguanine to N1 of the target adenine in the RNA substrate (A63). The structure of the ribozyme suggested a mechanism in which a cytosine (C10) acts as a general acid to protonate -alkylguanine N1. Here, we have analyzed the role of the C10 general acid and the A63 nucleophile by atomic mutagenesis and computation.

Improvements in Precision of Relative Binding Free Energy Calculations Afforded by the Alchemical Enhanced Sampling (ACES) Approach.

Accurate predictions of how strongly small molecules bind to proteins, such as those afforded by relative binding free energy (RBFE) calculations, can greatly increase the efficiency of the hit-to-lead and lead optimization stages of the drug discovery process. The success of such calculations, however, relies heavily on their precision. Here, we show that a recently developed alchemical enhanced sampling (ACES) approach can consistently improve the precision of RBFE calculations on a large and diverse set of proteins and small molecule ligands.

Software Infrastructure for Next-Generation QM/MM-ΔMLP Force Fields.

We present software infrastructure for the design and testing of new quantum mechanical/molecular mechanical and machine-learning potential (QM/MM-ΔMLP) force fields for a wide range of applications. The software integrates Amber's molecular dynamics simulation capabilities with fast, approximate quantum models in the xtb package and machine-learning potential corrections in DeePMD-kit. The xtb package implements the recently developed density-functional tight-binding QM models with multipolar electrostatics and density-dependent dispersion (GFN2-xTB), and the interface with Amber enables their use in periodic boundary QM/MM simulations with linear-scaling QM/MM particle-mesh Ewald electrostatics.

Electronic and Nuclear Quantum Effects on Proton Transfer Reactions of Guanine-Thymine (G-T) Mispairs Using Combined Quantum Mechanical/Molecular Mechanical and Machine Learning Potentials.

Rare tautomeric forms of nucleobases can lead to Watson-Crick-like (WC-like) mispairs in DNA, but the process of proton transfer is fast and difficult to detect experimentally. NMR studies show evidence for the existence of short-time WC-like guanine-thymine (G-T) mispairs; however, the mechanism of proton transfer and the degree to which nuclear quantum effects play a role are unclear. We use a B-DNA helix exhibiting a wGT mispair as a model system to study tautomerization reactions.

Amber free energy tools: Interoperable software for free energy simulations using generalized quantum mechanical/molecular mechanical and machine learning potentials.

We report the development and testing of new integrated cyberinfrastructure for performing free energy simulations with generalized hybrid quantum mechanical/molecular mechanical (QM/MM) and machine learning potentials (MLPs) in Amber. The Sander molecular dynamics program has been extended to leverage fast, density-functional tight-binding models implemented in the DFTB+ and xTB packages, and an interface to the DeePMD-kit software enables the use of MLPs. The software is integrated through application program interfaces that circumvent the need to perform "system calls" and enable the incorporation of long-range Ewald electrostatics into the external software's self-consistent field procedure.

Alchemical Enhanced Sampling with Optimized Phase Space Overlap.

An alchemical enhanced sampling (ACES) method has recently been introduced to facilitate importance sampling in free energy simulations. The method achieves enhanced sampling from Hamiltonian replica exchange within a dual topology framework while utilizing new smoothstep softcore potentials. A common sampling problem encountered in lead optimization is the functionalization of aromatic rings that exhibit distinct conformational preferences when interacting with the protein.

A multi-center interventional study to assess pharmacokinetics, effectiveness, and tolerability of prolonged-release tacrolimus after pediatric kidney transplantation: study protocol for a prospective, open-label, randomized, two-phase, two-sequence, single dose, crossover, phase III b trial.

Background: Tacrolimus, a calcineurin inhibitor (CNI), is currently the first-line immunosuppressive agent in kidney transplantation. The therapeutic index of tacrolimus is narrow due to due to the substantial impact of minor variations in drug concentration or exposure on clinical outcomes (i.e.

Surface-Accelerated String Method for Locating Minimum Free Energy Paths.

We present a surface-accelerated string method (SASM) to efficiently optimize low-dimensional reaction pathways from the sampling performed with expensive quantum mechanical/molecular mechanical (QM/MM) Hamiltonians. The SASM accelerates the convergence of the path using the aggregate sampling obtained from the current and previous string iterations, whereas approaches like the string method in collective variables (SMCV) or the modified string method in collective variables (MSMCV) update the path only from the sampling obtained from the current iteration. Furthermore, the SASM decouples the number of images used to perform sampling from the number of synthetic images used to represent the path.

AmberTools.

AmberTools is a free and open-source collection of programs used to set up, run, and analyze molecular simulations. The newer features contained within AmberTools23 are briefly described in this Application note.

Multifunctional IgG/IgM antibodies and cellular cytotoxicity are elicited by the full-length MSP1 SumayaVac-1 malaria vaccine.

Radical control of malaria likely requires a vaccine that targets both the asymptomatic liver stages and the disease-causing blood stages of the human malaria parasite Plasmodium falciparum. While substantial progress has been made towards liver stage vaccines, the development of a blood stage vaccine is lagging behind. We have recently conducted a first-in-human clinical trial to evaluate the safety and immunogenicity of the recombinant, full-length merozoite surface protein 1 (MSP1) formulated with GLA-SE as adjuvant.

Estimation of frequency factors for the calculation of kinetic isotope effects from classical and path integral free energy simulations.

We use the modified Bigeleisen-Mayer equation to compute kinetic isotope effect values for non-enzymatic phosphoryl transfer reactions from classical and path integral molecular dynamics umbrella sampling. The modified form of the Bigeleisen-Mayer equation consists of a ratio of imaginary mode vibrational frequencies and a contribution arising from the isotopic substitution's effect on the activation free energy, which can be computed from path integral simulation. In the present study, we describe a practical method for estimating the frequency ratio correction directly from umbrella sampling in a manner that does not require normal mode analysis of many geometry optimized structures.

Catalytic mechanism and pH dependence of a methyltransferase ribozyme (MTR1) from computational enzymology.

A methyltransferase ribozyme (MTR1) was selected in vitro to catalyze alkyl transfer from exogenous O6-methylguanine (O6mG) to a target adenine N1, and recently, high-resolution crystal structures have become available. We use a combination of classical molecular dynamics, ab initio quantum mechanical/molecular mechanical (QM/MM) and alchemical free energy (AFE) simulations to elucidate the atomic-level solution mechanism of MTR1. Simulations identify an active reactant state involving protonation of C10 that hydrogen bonds with O6mG:N1.

Modern semiempirical electronic structure methods and machine learning potentials for drug discovery: Conformers, tautomers, and protonation states.

Modern semiempirical electronic structure methods have considerable promise in drug discovery as universal "force fields" that can reliably model biological and drug-like molecules, including alternative tautomers and protonation states. Herein, we compare the performance of several neglect of diatomic differential overlap-based semiempirical (MNDO/d, AM1, PM6, PM6-D3H4X, PM7, and ODM2), density-functional tight-binding based (DFTB3, DFTB/ChIMES, GFN1-xTB, and GFN2-xTB) models with pure machine learning potentials (ANI-1x and ANI-2x) and hybrid quantum mechanical/machine learning potentials (AIQM1 and QDπ) for a wide range of data computed at a consistent ωB97X/6-31G* level of theory (as in the ANI-1x database). This data includes conformational energies, intermolecular interactions, tautomers, and protonation states.

Dissociative Transition State in Hepatitis Delta Virus Ribozyme Catalysis.

Ribonucleases and small nucleolytic ribozymes are both able to catalyze RNA strand cleavage through 2'--transphosphorylation, provoking the question of whether protein and RNA enzymes facilitate mechanisms that pass through the same or distinct transition states. Here, we report the primary and secondary O kinetic isotope effects for hepatitis delta virus ribozyme catalysis that reveal a dissociative, metaphosphate-like transition state in stark contrast to the late, associative transition states observed for reactions catalyzed by specific base, Zn ions, or ribonuclease A. This new information provides evidence for a discrete ribozyme active site design that modulates the RNA cleavage pathway to pass through an altered transition state.

QDπ: A Quantum Deep Potential Interaction Model for Drug Discovery.

We report QDπ-v1.0 for modeling the internal energy of drug molecules containing H, C, N, and O atoms. The QDπ model is in the form of a quantum mechanical/machine learning potential correction (QM/Δ-MLP) that uses a fast third-order self-consistent density-functional tight-binding (DFTB3/3OB) model that is corrected to a quantitatively high-level of accuracy through a deep-learning potential (DeepPot-SE).

A TREM2-activating antibody with a blood-brain barrier transport vehicle enhances microglial metabolism in Alzheimer's disease models.

Loss-of-function variants of TREM2 are associated with increased risk of Alzheimer's disease (AD), suggesting that activation of this innate immune receptor may be a useful therapeutic strategy. Here we describe a high-affinity human TREM2-activating antibody engineered with a monovalent transferrin receptor (TfR) binding site, termed antibody transport vehicle (ATV), to facilitate blood-brain barrier transcytosis. Upon peripheral delivery in mice, ATV:TREM2 showed improved brain biodistribution and enhanced signaling compared to a standard anti-TREM2 antibody.

AMBER Free Energy Tools: A New Framework for the Design of Optimized Alchemical Transformation Pathways.

We develop a framework for the design of optimized alchemical transformation pathways in free energy simulations using nonlinear mixing and a new functional form for so-called "softcore" potentials. We describe the implementation and testing of this framework in the GPU-accelerated AMBER software suite. The new optimized alchemical transformation pathways integrate a number of important features, including (1) the use of smoothstep functions to stabilize behavior near the transformation end points, (2) consistent power scaling of Coulomb and Lennard-Jones (LJ) interactions with unitless control parameters to maintain balance of electrostatic attractions and exchange repulsions, (3) pairwise form based on the LJ contact radius for the effective interaction distance with separation-shifted scaling, and (4) rigorous smoothing of the potential at the nonbonded cutoff boundary.

The HIF-prolyl hydroxylases have distinct and nonredundant roles in colitis-associated cancer.

Colitis-associated colorectal cancer (CAC) is a severe complication of inflammatory bowel disease (IBD). HIF-prolyl hydroxylases (PHD1, PHD2, and PHD3) control cellular adaptation to hypoxia and are considered promising therapeutic targets in IBD. However, their relevance in the pathogenesis of CAC remains elusive.

AMBER Drug Discovery Boost Tools: Automated Workflow for Production Free-Energy Simulation Setup and Analysis (ProFESSA).

We report an automated workflow for production free-energy simulation setup and analysis (ProFESSA) using the GPU-accelerated AMBER free-energy engine with enhanced sampling features and analysis tools, part of the AMBER Drug Discovery Boost package that has been integrated into the AMBER22 release. The workflow establishes a flexible, end-to-end pipeline for performing alchemical free-energy simulations that brings to bear technologies, including new enhanced sampling features and analysis tools, to practical drug discovery problems. ProFESSA provides the user with top-level control of large sets of free-energy calculations and offers access to the following key functionalities: (1) automated setup of file infrastructure; (2) enhanced conformational and alchemical sampling with the ACES method; and (3) network-wide free-energy analysis with the optional imposition of cycle closure and experimental constraints.

Multireference Generalization of the Weighted Thermodynamic Perturbation Method.

We describe the generalized weighted thermodynamic perturbation (gwTP) method for estimating the free energy surface of an expensive "high-level" potential energy function from the umbrella sampling performed with multiple inexpensive "low-level" reference potentials. The gwTP method is a generalization of the weighted thermodynamic perturbation (wTP) method developed by Li and co-workers [J. Chem.

Pharmacodynamic Monitoring of Ciclosporin and Tacrolimus: Insights From Nuclear Factor of Activated T-Cell-Regulated Gene Expression in Healthy Volunteers.

Background: Although therapeutic drug monitoring of calcineurin inhibitor (CNI) concentrations is performed routinely in clinical practice, an identical concentration may lead to different effects in different patients. Although the quantification of nuclear factor of activated T-cell-regulated gene expression (NFAT-RGE) is a promising method for measuring individual CNI effects, CNI pharmacodynamics are as of yet incompletely understood.

Methods: CNI concentrations and NFAT-RGEs were quantified in 24 healthy volunteers receiving either ciclosporin or tacrolimus in 2 clinical trials.

Incidences of Infectious Events in a Renal Transplant Cohort of the German Center of Infectious Diseases (DZIF).

Background: Infectious complications are a major cause of morbidity and mortality after kidney transplantation.

Methods: In this transplant cohort study at the German Center of Infectious Diseases (DZIF), we evaluated all infections occurring during the first year after renal transplantation. We assessed microbial etiology, incidence rates, and temporal occurrence of these infections.

Combined QM/MM, Machine Learning Path Integral Approach to Compute Free Energy Profiles and Kinetic Isotope Effects in RNA Cleavage Reactions.

We present a fast, accurate, and robust approach for determination of free energy profiles and kinetic isotope effects for RNA 2'-O-transphosphorylation reactions with inclusion of nuclear quantum effects. We apply a deep potential range correction (DPRc) for combined quantum mechanical/molecular mechanical (QM/MM) simulations of reactions in the condensed phase. The method uses the second-order density-functional tight-binding method (DFTB2) as a fast, approximate base QM model.