Regression-Based Active Learning for Accessible Acceleration of Ultra-Large Library Docking.

Structure-based drug discovery is a process for both hit finding and optimization that relies on a validated three-dimensional model of a target biomolecule, used to rationalize the structure-function relationship for this particular target. An ultralarge virtual screening approach has emerged recently for rapid discovery of high-affinity hit compounds, but it requires substantial computational resources. This study shows that active learning with simple linear regression models can accelerate virtual screening, retrieving up to 90% of the top-1% of the docking hit list after docking just 10% of the ligands. The results demonstrate that it is unnecessary to use complex models, such as deep learning approaches, to predict the imprecise results of ligand docking with a low sampling depth. Furthermore, we explore active learning meta-parameters and find that constant batch size models with a simple ensembling method provide the best ligand retrieval rate. Finally, our approach is validated on the ultralarge size virtual screening data set, retrieving 70% of the top-0.05% of ligands after screening only 2% of the library. Altogether, this work provides a computationally accessible approach for accelerated virtual screening that can serve as a blueprint for the future design of low-compute agents for exploration of the chemical space via large-scale accelerated docking. With recent breakthroughs in protein structure prediction, this method can significantly increase accessibility for the academic community and aid in the rapid discovery of high-affinity hit compounds for various targets.

The Role of Tyr-His-Trp Triad and Water Molecule Near the N1-Atom of 2-Hydroperoxycoelenterazine in Bioluminescence of Hydromedusan Photoproteins: Structural and Mutagenesis Study.

Hydromedusan photoproteins responsible for the bioluminescence of a variety of marine jellyfish and hydroids are a unique biochemical system recognized as a stable enzyme-substrate complex consisting of apoprotein and preoxygenated coelenterazine, which is tightly bound in the protein inner cavity. The binding of calcium ions to the photoprotein molecule is only required to initiate the light emission reaction. Although numerous experimental and theoretical studies on the bioluminescence of these photoproteins were performed, many features of their functioning are yet unclear. In particular, which ionic state of dioxetanone intermediate decomposes to yield a coelenteramide in an excited state and the role of the water molecule residing in a proximity to the N1 atom of 2-hydroperoxycoelenterazine in the bioluminescence reaction are still under discussion. With the aim to elucidate the function of this water molecule as well as to pinpoint the amino acid residues presumably involved in the protonation of the primarily formed dioxetanone anion, we constructed a set of single and double obelin and aequorin mutants with substitutions of His, Trp, Tyr, and Ser to residues with different properties of side chains and investigated their bioluminescence properties (specific activity, bioluminescence spectra, stopped-flow kinetics, and fluorescence spectra of Ca2+-discharged photoproteins). Moreover, we determined the spatial structure of the obelin mutant with a substitution of His64, the key residue of the presumable proton transfer, to Phe. On the ground of the bioluminescence properties of the obelin and aequorin mutants as well as the spatial structures of the obelin mutants with the replacements of His64 and Tyr138, the conclusion was made that, in fact, His residue of the Tyr-His-Trp triad and the water molecule perform the "catalytic function" by transferring the proton from solvent to the dioxetanone anion to generate its neutral ionic state in complex with water, as only the decomposition of this form of dioxetanone can provide the highest light output in the light-emitting reaction of the hydromedusan photoproteins.

Crystal structure of semi-synthetic obelin-v after calcium induced bioluminescence implies coelenteramine as the main reaction product.

Coelenterazine-v (CTZ-v), a synthetic vinylene-bridged π-extended derivative, is able to significantly alter bioluminescence spectra of different CTZ-dependent luciferases and photoproteins by shifting them towards longer wavelengths. However, Ca2+-regulated photoproteins activated with CTZ-v display very low bioluminescence activities that hampers its usage as a substrate of photoprotein bioluminescence. Here, we report the crystal structure of semi-synthetic Ca2+-discharged obelin-v bound with the reaction product determined at 2.1 Å resolution. Comparison of the crystal structure of Ca2+-discharged obelin-v with those of other obelins before and after bioluminescence reaction reveals no considerable changes in the overall structure. However, the drastic changes in CTZ-binding cavity are observed owing to the completely different reaction product, coelenteramine-v (CTM-v). Since CTM-v is certainly the main product of obelin-v bioluminescence and is considered to be a product of the "dark" pathway of dioxetanone intermediate decomposition, it explains the low bioluminescence activity of obelin and apparently of other photoproteins with CTZ-v.

Structural basis for receptor selectivity and inverse agonism in S1P5 receptors.

The bioactive lysophospholipid sphingosine-1-phosphate (S1P) acts via five different subtypes of S1P receptors (S1PRs) - S1P1-5. S1P5 is predominantly expressed in nervous and immune systems, regulating the egress of natural killer cells from lymph nodes and playing a role in immune and neurodegenerative disorders, as well as carcinogenesis. Several S1PR therapeutic drugs have been developed to treat these diseases; however, they lack receptor subtype selectivity, which leads to side effects. In this article, we describe a 2.2 Å resolution room temperature crystal structure of the human S1P5 receptor in complex with a selective inverse agonist determined by serial femtosecond crystallography (SFX) at the Pohang Accelerator Laboratory X-Ray Free Electron Laser (PAL-XFEL) and analyze its structure-activity relationship data. The structure demonstrates a unique ligand-binding mode, involving an allosteric sub-pocket, which clarifies the receptor subtype selectivity and provides a template for structure-based drug design. Together with previously published S1PR structures in complex with antagonists and agonists, our structure with S1P5-inverse agonist sheds light on the activation mechanism and reveals structural determinants of the inverse agonism in the S1PR family.

Nonstationary dynamics of the sine lattice consisting of three pendula (trimer).

The low- and high-amplitude oscillations in the system of three nonlinear coupled pendula (trimer) are analyzed beyond the quasilinear approximation. The considered oscillations are fundamental for many models of the energy exchange processes in physical, mechanical, and biological systems, in particular, for the torsional vibrations of flexible polymers or DNA's double strands. We obtained the conditions of the basic stationary solutions' stability. These solutions correspond to the nonlinear normal modes (NNMs), the instability of which leads to the appearance of localized NNMs (stationary energy localization). Using an asymptotic procedure, we reduce the dimension of the system's phase space that allows us to analyze the energy exchange between pendula in the slow timescale and to reveal periodic interparticle energy exchange and nonstationary energy localization. It has been shown recently that essentially nonstationary resonance processes of this type are adequately described in terms of the limiting phase trajectories (LPTs) corresponding to beatings between the oscillators or coherence domains in the slow timescale. Moreover, it turns out that criteria of the transition to the stationary and nonstationary energy localization can be formulated as the bifurcation conditions for NNMs and LPTs, respectively. The trimer under consideration is a nonintegrable system, and therefore its equations of motion is only after dimensions reduction can be analyzed by the Poincare sections method. Finally, we aim to study the highly nonstationary regimes, which correspond to beatinglike periodic or quasiperiodic recurrent energy exchange between the pendula.

Nonconventional synchronization and energy localization in weakly coupled autogenerators.

The present work follows our previous study dealing with a new type of synchronization in a system of two weakly coupled generalized van der Pol-Duffing autogenerators. The essence of the effect revealed is that the synchronized oscillations are not stationary but accompanied by the most intensive energy exchange between the oscillators. The phase shift between the generators remains constant most of the time, except for vanishingly small transitional intervals. The current analysis deals with a generalized model in order to clarify the frequency detuning effect. We found that varying the frequency detuning, nonlinearity, and dissipation parameters can lead to structural changes in phase diagrams of the energy exchange dynamics, with important transitions from the intensive energy exchange to its localization on one of the two oscillators. The main conclusion is that stationary and nonstationary synchronizations associate with nonlinear normal and local modes, respectively. The analysis uses phase plane diagrams, including the concept of limiting phase trajectories, whose role in nonstationary synchronization appears to be similar to the role of nonlinear normal modes in conventional stationary states.