PathInHydro, a Set of Machine Learning Models to Identify Unbinding Pathways of Gas Molecules in [NiFe] Hydrogenases.

Machine learning (ML) is a powerful tool for the automated data analysis of molecular dynamics (MD) simulations. Recent studies showed that ML models can be used to identify protein-ligand unbinding pathways and understand the underlying mechanism. To expedite the examination of MD simulations, we constructed PathInHydro, a set of supervised ML models capable of automatically assigning unbinding pathways for the dissociation of gas molecules from [NiFe] hydrogenases, using the unbinding trajectories of CO and H from [NiFe] hydrogenase as a training set.

Recent Progress in Modeling and Simulation of Biomolecular Crowding and Condensation Inside Cells.

Macromolecular crowding in the cellular cytoplasm can potentially impact diffusion rates of proteins, their intrinsic structural stability, binding of proteins to their corresponding partners as well as biomolecular organization and phase separation. While such intracellular crowding can have a large impact on biomolecular structure and function, the molecular mechanisms and driving forces that determine the effect of crowding on dynamics and conformations of macromolecules are so far not well understood. At a molecular level, computational methods can provide a unique lens to investigate the effect of macromolecular crowding on biomolecular behavior, providing us with a resolution that is challenging to reach with experimental techniques alone.

Computational screening of the effects of mutations on protein-protein off-rates and dissociation mechanisms by τRAMD.

The dissociation rate, or its reciprocal, the residence time (τ), is a crucial parameter for understanding the duration and biological impact of biomolecular interactions. Accurate prediction of τ is essential for understanding protein-protein interactions (PPIs) and identifying potential drug targets or modulators for tackling diseases. Conventional molecular dynamics simulation techniques are inherently constrained by their limited timescales, making it challenging to estimate residence times, which typically range from minutes to hours.

Characterization of the Bottlenecks and Pathways for Inhibitor Dissociation from [NiFe] Hydrogenase.

[NiFe] hydrogenases can act as efficient catalysts for hydrogen oxidation and biofuel production. However, some [NiFe] hydrogenases are inhibited by gas molecules present in the environment, such as O and CO. One strategy to engineer [NiFe] hydrogenases and achieve O- and CO-tolerant enzymes is by introducing point mutations to block the access of inhibitors to the catalytic site.

A community effort in SARS-CoV-2 drug discovery.

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors.

pH-dependence of the Plasmodium falciparum chloroquine resistance transporter is linked to the transport cycle.

The chloroquine resistance transporter, PfCRT, of the human malaria parasite Plasmodium falciparum is sensitive to acidic pH. Consequently, PfCRT operates at 60% of its maximal drug transport activity at the pH of 5.2 of the digestive vacuole, a proteolytic organelle from which PfCRT expels drugs interfering with heme detoxification.

Molecular characterization of AIFM2/FSP1 inhibition by iFSP1-like molecules.

Ferroptosis is a form of cell death characterized by phospholipid peroxidation, where numerous studies have suggested that the induction of ferroptosis is a therapeutic strategy to target therapy refractory cancer entities. Ferroptosis suppressor protein 1 (FSP1), an NAD(P)H-ubiquinone reductase, is a key determinant of ferroptosis vulnerability, and its pharmacological inhibition was shown to strongly sensitize cancer cells to ferroptosis. A first generation of FSP1 inhibitors, exemplified by the small molecule iFSP1, has been reported; however, the molecular mechanisms underlying inhibition have not been characterized in detail.

Advances in computational methods for ligand binding kinetics.

Binding kinetic parameters can be correlated with drug efficacy, which in recent years led to the development of various computational methods for predicting binding kinetic rates and gaining insight into protein-drug binding paths and mechanisms. In this review, we introduce and compare computational methods recently developed and applied to two systems, trypsin-benzamidine and kinase-inhibitor complexes. Methods involving enhanced sampling in molecular dynamics simulations or machine learning can be used not only to predict kinetic rates, but also to reveal factors modulating the duration of residence times, selectivity, and drug resistance to mutations.

Diffusion of small molecule drugs is affected by surface interactions and crowder proteins.

Crowded environments are known to affect the diffusion of macromolecules, but their effects on the diffusion of small molecules are largely uncharacterized. We investigate how three protein crowders, bovine serum albumin (BSA), hen egg-white lysozyme, and myoglobin, influence the diffusion rates and interactions of four small molecules: fluorescein, and three drugs, doxorubicin, glycogen synthase kinase-3 inhibitor SB216763, and quinacrine. Using Line-FRAP measurements, Brownian dynamics simulations, and molecular docking, we find that the diffusion rates of the small molecules are highly affected by self-aggregation, interactions with the proteins, and surface adsorption.

The Non-phosphorylating Glyceraldehyde-3-Phosphate Dehydrogenase GapN Is a Potential New Drug Target in .

The strict human pathogen causes infections of varying severity, ranging from self-limiting suppurative infections to life-threatening diseases like necrotizing fasciitis or streptococcal toxic shock syndrome. Here, we show that the non-phosphorylating glyceraldehyde-3-phosphate dehydrogenase GapN is an essential enzyme for . GapN converts glyceraldehyde 3-phosphate into 3-phosphoglycerate coupled to the reduction of NADP to NADPH.

Ligand unbinding mechanisms and kinetics for T4 lysozyme mutants from τRAMD simulations.

The protein-ligand residence time, τ, influences molecular function in biological networks and has been recognized as an important determinant of drug efficacy. To predict τ, computational methods must overcome the problem that τ often exceeds the timescales accessible to conventional molecular dynamics (MD) simulation. Here, we apply the τ-Random Acceleration Molecular Dynamics (τRAMD) method to a set of kinetically characterized complexes of T4 lysozyme mutants with small, engineered binding cavities.

Prediction of the Drug-Target Binding Kinetics for Flexible Proteins by Comparative Binding Energy Analysis.

There is growing consensus that the optimization of the kinetic parameters for drug-protein binding leads to improved drug efficacy. Therefore, computational methods have been developed to predict kinetic rates and to derive quantitative structure-kinetic relationships (QSKRs). Many of these methods are based on crystal structures of ligand-protein complexes.

A Protocol to Use Comparative Binding Energy Analysis to Estimate Drug-Target Residence Time.

Comparative Binding Energy (COMBINE) analysis is an approach for deriving a target-specific scoring function to compute binding free energy, drug-binding kinetics, or a related property by exploiting the information contained in the three-dimensional structures of receptor-ligand complexes. Here, we describe the process of setting up and running COMBINE analysis to derive a Quantitative Structure-Kinetics Relationship (QSKR) for the dissociation rate constants (k) of inhibitors of a drug target. The derived QSKR model can be used to estimate residence times (τ, τ=1/k) for similar inhibitors binding to the same target, and it can also help to identify key receptor-ligand interactions that distinguish inhibitors with short and long residence times.

Line-FRAP, A Versatile Method to Measure Diffusion Rates In Vitro and In Vivo.

The crowded cellular milieu affect molecular diffusion through hard (occluded space) and soft (weak, non-specific) interactions. Multiple methods have been developed to measure diffusion coefficients at physiological protein concentrations within cells, each with its limitations. Here, we show that Line-FRAP, combined with rigours data analysis, is able to determine diffusion coefficients in a variety of environments, from in vitro to in vivo.

Structure-kinetic relationship reveals the mechanism of selectivity of FAK inhibitors over PYK2.

There is increasing evidence of a significant correlation between prolonged drug-target residence time and increased drug efficacy. Here, we report a structural rationale for kinetic selectivity between two closely related kinases: focal adhesion kinase (FAK) and proline-rich tyrosine kinase 2 (PYK2). We found that slowly dissociating FAK inhibitors induce helical structure at the DFG motif of FAK but not PYK2.

Recent progress in molecular simulation methods for drug binding kinetics.

Due to the contribution of drug-target binding kinetics to drug efficacy, there is a high level of interest in developing methods to predict drug-target binding kinetic parameters. During the review period, a wide range of enhanced sampling molecular dynamics simulation-based methods has been developed for computing drug-target binding kinetics and studying binding and unbinding mechanisms. Here, we assess the performance of these methods considering two benchmark systems in detail: mutant T4 lysozyme-ligand complexes and a large set of N-HSP90-inhibitor complexes.

Myotubularin-related protein 7 activates peroxisome proliferator-activated receptor-gamma.

Peroxisome proliferator-activated receptor-gamma (PPARγ) is a transcription factor drugable by agonists approved for treatment of type 2 diabetes, but also inhibits carcinogenesis and cell proliferation in vivo. Activating mutations in the Kirsten rat sarcoma viral oncogene homologue (KRAS) gene mitigate these beneficial effects by promoting a negative feedback-loop comprising extracellular signal-regulated kinase 1/2 (ERK1/2) and mitogen-activated kinase kinase 1/2 (MEK1/2)-dependent inactivation of PPARγ. To overcome this inhibitory mechanism, we searched for novel post-translational regulators of PPARγ.

Structural characterization of an Arf dimer interface: molecular mechanism of Arf-dependent membrane scission.

Dimerization of the small GTPase Arf is prerequisite for the scission of COPI-coated transport vesicles. Here, we quantify the monomer/dimer equilibrium of Arf within the membrane and show that after membrane scission, Arf dimers are restricted to donor membranes. By hydrogen exchange mass spectrometry, we define the interface of activated dimeric Arf within its switch II region.

From Brazil to Germany: Challenges and Advantages.

Moving to a new country, with a different culture and a new environment, is not an easy decision. In this perspective, I present some reasons that made me, a Brazilian computational biochemist, move abroad to do postdoctoral research and some of the challenges I faced before and after moving.

Escape of a Small Molecule from Inside T4 Lysozyme by Multiple Pathways.

The T4 lysozyme L99A mutant is often used as a model system to study small-molecule binding to proteins, but pathways for ligand entry and exit from the buried binding site and the associated protein conformational changes have not been fully resolved. Here, molecular dynamics simulations were employed to model benzene exit from its binding cavity using the weighted ensemble (WE) approach to enhance sampling of low-probability unbinding trajectories. Independent WE simulations revealed four pathways for benzene exit, which correspond to transient tunnels spontaneously formed in previous simulations of apo T4 lysozyme.