Accelerated molecular dynamics: a promising and efficient simulation method for biomolecules

Which cryptic sites are feasible drug targets?

Cryptic sites can expand the space of druggable proteins, but the potential usefulness of such sites needs to be investigated before any major effort. Given that the binding pockets are not formed, the druggability of such sites is not well understood. The analysis of proteins and their ligands shows that cryptic sites that are formed primarily by the motion of side chains moving out of the pocket to enable ligand binding generally do not bind drug-sized molecules with sufficient potency. By contrast, sites that are formed by loop or hinge motion are potentially valuable drug targets. Arguments are provided to explain the underlying causes in terms of classical enzyme inhibition theory and the kinetics of side chain motion and ligand binding.

Molecular mechanics studies of factors affecting overall rate in cascade reactions: Multi-enzyme colocalization and environment

Millions of years of evolution have optimized many biosynthetic pathways by use of multi-step catalysis. In addition, multi-step metabolic pathways are commonly found in and on membrane-bound organelles in eukaryotic biochemistry. The fundamental mechanisms that facilitate these reaction processes provide strategies to bioengineer metabolic pathways in synthetic chemistry. Using Brownian dynamics simulations, here we modeled intermediate substrate transportation of colocalized yeast-ester biosynthesis enzymes on the membrane. The substrate acetate ion traveled from the pocket of aldehyde dehydrogenase to its target enzyme acetyl-CoA synthetase, then the substrate acetyl CoA diffused from Acs1 to the active site of the next enzyme, alcohol-O-acetyltransferase. Arranging two enzymes with the smallest inter-enzyme distance of 60 Å had the fastest average substrate association time as compared with anchoring enzymes with larger inter-enzyme distances. When the off-target side reactions were turned on, most substrates were lost, which suggests that native localization is necessary for efficient final product synthesis. We also evaluated the effects of intermolecular interactions, local substrate concentrations, and membrane environment to bring mechanistic insights into the colocalization pathways. The computation work demonstrates that creating spatially organized multi-enzymes on membranes can be an effective strategy to increase final product synthesis in bioengineering systems.

Nuclear Receptor Interdomain Communication is Mediated by the Hinge with Ligand Specificity

Nuclear receptors are ligand-induced transcription factors that bind directly to target genes and regulate their expression. Ligand binding initiates conformational changes that propagate to other domains, allosterically regulating their activity. The nature of this interdomain communication in nuclear receptors is poorly understood, largely owing to the difficulty of experimentally characterizing full-length structures. We have applied computational modeling approaches to describe and study the structure of the full-length farnesoid X receptor (FXR), approximated by the DNA binding domain (DBD) and ligand binding domain (LBD) connected by the flexible hinge region. Using extended molecular dynamics simulations (>10 microseconds) and enhanced sampling simulations, we provide evidence that ligands selectively induce domain rearrangement, leading to interdomain contact. We use protein-protein interaction assays to provide experimental evidence of these interactions, identifying a critical role of the hinge in mediating interdomain contact. Our results illuminate previously unknown aspects of interdomain communication in FXR and provide a framework to enable characterization of other full-length nuclear receptors.

Catalytic mechanism study of ATP-citrate lyase during citryl-CoA synthesis process

ATP-citrate lyase (ACLY) is a critical metabolic enzyme and promising target for drug development. The structure determinations of ACLY have revealed its homotetramer states with various subunit symmetries, but catalytic mechanism of ACLY tetramer and the importance of subunit symmetry have not been clarified. Here, we constructed the free energy landscape of ACLY tetramer with arbitrary subunit symmetries and investigated energetic and conformational coupling of subunits during citryl-CoA synthesis process. The optimal conformational pathway indicates that ACLY tetramer encounters three critical conformational barriers and undergoes a loss of rigid-D2 symmetry to gain an energetic advantage. Energetic coupling of conformational changes and biochemical reactions suggests that these biological events are not independent but rather coupled with each other, showing a comparable energy barrier to the experimental data for the rate-limiting step. These findings could contribute to further research on catalytic mechanism, functional modulation, and inhibitor design of ACLY.

Assessing AF2's ability to predict structural ensembles of proteins

Recent breakthroughs in protein structure prediction have enhanced the precision and speed at which protein configurations can be determined. Additionally, molecular dynamics (MD) simulations serve as a crucial tool for capturing the conformational space of proteins, providing valuable insights into their structural fluctuations. However, the scope of MD simulations is often limited by the accessible timescales and the computational resources available, posing challenges to comprehensively exploring protein behaviors. Recently emerging approaches have focused on expanding the capability of AlphaFold2 (AF2) to predict conformational substates of protein. Here, we benchmark the performance of various workflows that have adapted AF2 for ensemble prediction and compare the obtained structures with ensembles obtained from MD simulations and NMR. We provide an overview of the levels of performance and accessible timescales that can currently be achieved with machine learning (ML) based ensemble generation. Significant minima of the free energy surfaces remain undetected.

ATOMDANCE: Kernel-based denoising and choreographic analysis for protein dynamic comparison

Comparative methods in molecular evolution and structural biology rely heavily upon the site-wise analysis of DNA sequence and protein structure, both static forms of information. However, it is widely accepted that protein function results from nanoscale nonrandom machine-like motions induced by evolutionarily conserved molecular interactions. Comparisons of molecular dynamics (MD) simulations conducted between homologous sites representative of different functional or mutational states can potentially identify local effects on binding interaction and protein evolution. In addition, comparisons of different (i.e., nonhomologous) sites within MD simulations could be employed to identify functional shifts in local time-coordinated dynamics indicative of logic gating within proteins. However, comparative MD analysis is challenged by the large fraction of protein motion caused by random thermal noise in the surrounding solvent. Therefore, properly denoised MD comparisons could reveal functional sites involving these machine-like dynamics with good accuracy. Here, we introduce ATOMDANCE, a user-interfaced suite of comparative machine learning-based denoising tools designed for identifying functional sites and the patterns of coordinated motion they can create within MD simulations. ATOMDANCE-maxDemon4.0 employs Gaussian kernel functions to compute site-wise maximum mean discrepancy between learned features of motion, thereby assessing denoised differences in the nonrandom motions between functional or evolutionary states (e.g., ligand bound versus unbound, wild-type versus mutant). ATOMDANCE-maxDemon4.0 also employs maximum mean discrepancy to analyze potential random amino acid replacements allowing for a site-wise test of neutral versus nonneutral evolution on the divergence of dynamic function in protein homologs. Finally, ATOMDANCE-Choreograph2.0 employs mixed-model analysis of variance and graph network to detect regions where time-synchronized shifts in dynamics occur. Here, we demonstrate ATOMDANCE's utility for identifying key sites involved in dynamic responses during functional binding interactions involving DNA, small-molecule drugs, and virus-host recognition, as well as understanding shifts in global and local site coordination occurring during allosteric activation of a pathogenic protease.

Methods for Classical-Mechanical Molecular Simulation in Chemistry: Achievements, Limitations, Perspectives

More than a half century ago it became feasible to simulate, using classical-mechanical equations of motion, the dynamics of molecular systems on a computer. Since then classical-physical molecular simulation has become an integral part of chemical research. It is widely applied in a variety of branches of chemistry and has significantly contributed to the development of chemical knowledge. It offers understanding and interpretation of experimental results, semiquantitative predictions for measurable and nonmeasurable properties of substances, and allows the calculation of properties of molecular systems under conditions that are experimentally inaccessible. Yet, molecular simulation is built on a number of assumptions, approximations, and simplifications which limit its range of applicability and its accuracy. These concern the potential-energy function used, adequate sampling of the vast statistical-mechanical configurational space of a molecular system and the methods used to compute particular properties of chemical systems from statistical-mechanical ensembles. During the past half century various methodological ideas to improve the efficiency and accuracy of classical-physical molecular simulation have been proposed, investigated, evaluated, implemented in general simulation software or were abandoned. The latter because of fundamental flaws or, while being physically sound, computational inefficiency. Some of these methodological ideas are briefly reviewed and the most effective methods are highlighted. Limitations of classical-physical simulation are discussed and perspectives are sketched.

Reinforcing Tunnel Network Exploration in Proteins Using Gaussian Accelerated Molecular Dynamics

Tunnels are structural conduits in biomolecules responsible for transporting chemical compounds and solvent molecules from the active site. They have been shown to be present in a wide variety of enzymes across all functional and structural classes. However, the study of such pathways is experimentally challenging, because they are typically transient. Computational methods, such as molecular dynamics (MD) simulations, have been successfully proposed to explore tunnels. Conventional MD (cMD) provides structural details to characterize tunnels but suffers from sampling limitations to capture rare tunnel openings on longer time scales. Therefore, in this study, we explored the potential of Gaussian accelerated MD (GaMD) simulations to improve the exploration of complex tunnel networks in enzymes. We used the haloalkane dehalogenase LinB and its two variants with engineered transport pathways, which are not only well-known for their application potential but have also been extensively studied experimentally and computationally regarding their tunnel networks and their importance in multistep catalytic reactions. Our study demonstrates that GaMD efficiently improves tunnel sampling and allows the identification of all known tunnels for LinB and its two mutants. Furthermore, the improved sampling provided insight into a previously unknown transient side tunnel (ST). The extensive conformational landscape explored by GaMD simulations allowed us to investigate in detail the mechanism of ST opening. We determined variant-specific dynamic properties of ST opening, which were previously inaccessible due to limited sampling of cMD. Our comprehensive analysis supports multiple indicators of the functional relevance of the ST, emphasizing its potential significance beyond structural considerations. In conclusion, our research proves that the GaMD method can overcome the sampling limitations of cMD for the effective study of tunnels in enzymes, providing further means for identifying rare tunnels in enzymes with the potential for drug development, precision medicine, and rational protein engineering.

Stable and accurate atomistic simulations of flexible molecules using conformationally generalisable machine learned potentials

Computational simulation methods based on machine learned potentials (MLPs) promise to revolutionise shape prediction of flexible molecules in solution, but their widespread adoption has been limited by the way in which training data is generated. Here, we present an approach which allows the key conformational degrees of freedom to be properly represented in reference molecular datasets. MLPs trained on these datasets using a global descriptor scheme are generalisable in conformational space, providing quantum chemical accuracy for all conformers. These MLPs are capable of propagating long, stable molecular dynamics trajectories, an attribute that has remained a challenge. We deploy the MLPs in obtaining converged conformational free energy surfaces for flexible molecules via well-tempered metadynamics simulations; this approach provides a hitherto inaccessible route to accurately computing the structural, dynamical and thermodynamical properties of a wide variety of flexible molecular systems. It is further demonstrated that MLPs must be trained on reference datasets with complete coverage of conformational space, including in barrier regions, to achieve stable molecular dynamics trajectories.

Computational insights into the mechanisms underlying structural destabilization and recovery in trafficking-deficient hERG mutants

Cardiovascular diseases are a major global health concern, responsible for a significant number of deaths each year, often linked to cardiac arrhythmias resulting from dysfunction in ion channels. Hereditary Long QT Syndrome (LQTS) is a condition characterized by a prolonged QT interval on ECG, increasing the risk of sudden cardiac death. The most common type of LQTS, LQT2, is caused by mutations in the hERG gene, affecting a potassium ion channel. The majority of these mutations disrupt the channel's trafficking to the cell membrane, leading to intracellular retention. Specific high-affinity hERG blockers (e.g., E-4031) can rescue this mutant phenotype, but the exact mechanism is unknown. This study used accelerated molecular dynamics simulations to investigate how these mutations affect the hERG channel's structure, folding, endoplasmic reticulum (ER) retention, and trafficking. We reveal that these mutations induce structural changes in the channel, narrowing its central pore and altering the conformation of the intracellular domains. These changes expose internalization signals that contribute to ER retention and degradation of the mutant hERG channels. Moreover, the study found that the trafficking rescue drug E-4031 can inhibit these structural changes, potentially rescuing the mutant channels. This research offers valuable insights into the structural issues responsible for the degradation of rescuable transmembrane trafficking mutants. Understanding the defective trafficking structure of the hERG channel could help identify binding sites for small molecules capable of restoring proper folding and facilitating channel trafficking. This knowledge has the potential to lead to mechanism-based therapies that address the condition at the cellular level, which may prove more effective than treating clinical symptoms, ultimately offering hope for individuals with hereditary Long QT Syndrome.

Navigating the landscape of enzyme design: from molecular simulations to machine learning

Global environmental issues and sustainable development call for new technologies for fine chemical synthesis and waste valorization. Biocatalysis has attracted great attention as the alternative to the traditional organic synthesis. However, it is challenging to navigate the vast sequence space to identify those proteins with admirable biocatalytic functions. The recent development of deep-learning based structure prediction methods such as AlphaFold2 reinforced by different computational simulations or multiscale calculations has largely expanded the 3D structure databases and enabled structure-based design. While structure-based approaches shed light on site-specific enzyme engineering, they are not suitable for large-scale screening of potential biocatalysts. Effective utilization of big data using machine learning techniques opens up a new era for accelerated predictions. Here, we review the approaches and applications of structure-based and machine-learning guided enzyme design. We also provide our view on the challenges and perspectives on effectively employing enzyme design approaches integrating traditional molecular simulations and machine learning, and the importance of database construction and algorithm development in attaining predictive ML models to explore the sequence fitness landscape for the design of admirable biocatalysts.

In silico evaluation of the role of Fab glycosylation in cetuximab antibody dynamics

Introduction

N-glycosylation is a post-translational modification that is highly important for the development of monoclonal antibodies (mAbs), as it regulates their biological activity, particularly in terms of immune effector functions. While typically added at the Fc level, approximately 15-25% of circulating antibodies exhibit glycosylation in the Fab domains as well. To the best of our knowledge, cetuximab (Erbitux®) is the only therapeutic antibody presenting Fab glycosylation approved world-wide targeting the epidermal growth factor receptor for the treatment of metastatic-colorectal and head and neck cancers. Additionally, it can trigger antibody-dependent cell cytotoxicity (ADCC), a response that typically is influenced by N-glycosylation at Fc level. However, the role of Fab glycosylation in cetuximab remains poorly understood. Hence, this study aims to investigate the structural role of Fab glycosylation on the conformational behavior of cetuximab.

Methods

The study was performed in silico via accelerated molecular dynamics simulations. The commercial cetuximab was compared to its form without Fab glycosylation and structural descriptors were evaluated to establish conformational differences.

Results

The results clearly show a correlation between the Fab glycosylation and structural descriptors that may modulate the conformational freedom of the antibody, potentially affecting Fc effector functions, and suggesting a negative role of Fab glycosylation on the interaction with FcγRIIIa.

Conclusion

Fab glycosylation of cetuximab is the most critical challenge for biosimilar development, but the differences highlighted in this work with respect to its aglycosylated form can improve the knowledge and represent also a great opportunity to develop novel strategies of biotherapeutics.

Parameter Optimization Method in Multidimensional Umbrella Sampling

Umbrella sampling (US) is an effective method for calculating free-energy landscapes (FELs). However, the complexity of controlling the sampling positions complicates multidimensional FEL calculations. In this study, we proposed a method for controlling sampling by optimizing the US parameters. This method comprises the introduction of a target point and the optimization of the parameters to sample a window around this point. We approximated each window to normal distributions using an umbrella integration method and calculated the divergences between the window distributions and the state distributed at the target position by a variationally enhanced sampling method. Thus, the minimization of the divergence facilitated sampling around the target point, after which the parameters could be optimized on the fly while performing equilibration simulation. In practice, our method employs bias potentials with off-diagonal terms, ensuring a more efficient calculation of multidimensional FEL. Additionally, we developed an algorithm for determining the target point for automated FEL search; the algorithm samples in a specified direction while controlling the overlap of distributions. We performed three different FEL calculations as examples: (1) the calculation of the permeation of a water molecule through a lipid bilayer (one-dimensional FEL), (2) the calculation of the internal structural changes in alanine dipeptide in water (two-dimensional FEL), and (3) the calculation of the internal structural changes from a β-strand structure to an α-helix structure in alanine decapeptide (Ala10, 16-dimensional FEL). These results confirmed that our method could control the number of US windows and calculate the high-dimensional FELs that could not be evaluated by the conventional US method.

PepBinding: A Workflow for Predicting Peptide Binding Structures by Combining Peptide Docking and Peptide Gaussian Accelerated Molecular Dynamics Simulations

Predicting protein-peptide interactions is crucial for understanding peptide binding processes and designing peptide drugs. However, traditional computational modeling approaches face challenges in accurately predicting peptide-protein binding structures due to the slow dynamics and high flexibility of the peptides. Here, we introduce a new workflow termed "PepBinding" for predicting peptide binding structures, which combines peptide docking, all-atom enhanced sampling simulations using the Peptide Gaussian accelerated Molecular Dynamics (Pep-GaMD) method, and structural clustering. PepBinding has been demonstrated on seven distinct model peptides. In peptide docking using HPEPDOCK, the peptide backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures ranged from 3.8 to 16.0 Å, corresponding to the medium to inaccurate quality models according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. The Pep-GaMD simulations performed for only 200 ns significantly improved the docking models, resulting in five medium and two acceptable quality models. Therefore, PepBinding is an efficient workflow for predicting peptide binding structures and is publicly available at https://github.com/MiaoLab20/PepBinding.

Multiple Parameter Replica Exchange Gaussian Accelerated Molecular Dynamics for Enhanced Sampling and Free Energy Calculation of Biomolecular Systems

This study introduces a novel method named multiple parameter replica exchange Gaussian accelerated molecular dynamics (MP-Rex-GaMD), building on the Gaussian accelerated molecular dynamics (GaMD) algorithm. GaMD enhances sampling and retrieves free energy information for biomolecular systems by adding a harmonic boost potential to smooth the potential energy surface without the need for predefined reaction coordinates. Our innovative approach advances the acceleration power and energetic reweighting accuracy of GaMD by incorporating a replica exchange algorithm that enables the exchange of multiple parameters, including the GaMD boost parameters of force constant and energy threshold, as well as temperature. Applying MP-Rex-GaMD to the three model systems of dialanine, chignolin, and HIV protease, we demonstrate its superior capability over conventional molecular dynamics and GaMD simulations in exploring protein conformations and effectively navigating various biomolecular states across energy barriers. MP-Rex-GaMD allows users to accurately map free energy landscapes through energetic reweighting, capturing the ensemble of biomolecular states from low-energy conformations to rare high-energy transitions within practical computational time scales.

Exploring Conformational Landscapes Along Anharmonic Low-Frequency Vibrations

We aim to automatize the identification of collective variables to simplify and speed up enhanced sampling simulations of conformational dynamics in biomolecules. We focus on anharmonic low-frequency vibrations that exhibit fluctuations on time scales faster than conformational transitions but describe a path of least resistance toward structural change. A key challenge is that harmonic approximations are ill-suited to characterize these vibrations, which are observed at far-infrared frequencies and are easily excited by thermal collisions at room temperature. Here, we approached this problem with a frequency-selective anharmonic (FRESEAN) mode analysis that does not rely on harmonic approximations and successfully isolates anharmonic low-frequency vibrations from short molecular dynamics simulation trajectories. We applied FRESEAN mode analysis to simulations of alanine dipeptide, a common test system for enhanced sampling simulation protocols, and compared the performance of isolated low-frequency vibrations to conventional user-defined collective variables (here backbone dihedral angles) in enhanced sampling simulations. The comparison shows that enhanced sampling along anharmonic low-frequency vibrations not only reproduces known conformational dynamics but can even further improve the sampling of slow transitions compared to user-defined collective variables. Notably, free energy surfaces spanned by low-frequency anharmonic vibrational modes exhibit lower barriers associated with conformational transitions relative to representations in backbone dihedral space. We thus conclude that anharmonic low-frequency vibrations provide a promising path for highly effective and fully automated enhanced sampling simulations of conformational dynamics in biomolecules.

Ligand Gaussian Accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional molecular dynamics (cMD), due to limited simulation time scales. Based on our previously developed ligand Gaussian accelerated molecular dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3″, in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding, and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as the model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 μs simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 were in agreement with the available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligands and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.

Repurposing of Drug Bank Compounds against Plasmodium falciparum Dihydroorotate Dehydrogenase as novel anti malarial drug candidates by Computational approaches

This study aimed to repurpose Drug Bank Compounds against P. falciparum Dihydroorotate dehydrogenase (Pf-DHODH)a potential molecular target for antimalarial drug development due to its vital role in P. falciparum survival. Initially, the MATGEN server was used to screen drugs against Pf-DHODH (PDB ID 6GJG), followed by revalidating the results through docking by Autodock Vina through PyRx. Based on the docking results, three drugs namely, Talnifumate, Sulfaphenazole, and (3S)-N-[(2S)-1-[2-(1H-indol-3-yl)ethylamino]-1-oxopropan-2-yl]-1-(4-methoxyphenyl)-5-oxopyrrolidine-3-carboxamide-were subjected to molecular dynamics simulation for 100 ns. Molecular dynamics simulation results indicate that (3S)-N-[(2S)-1-[2-(1H-indol-3-yl)ethylamino]-1-oxopropan-2-yl]-1-(4-methoxyphenyl)-5-oxopyrrolidine-3-carboxamide- and Sulfaphenazole may target Pf-DHODH by forming a stable protein-ligand complex as they showed better free binding energy -130.58 kJ/mol, and -79.84 kJ/mol, respectively as compared to the free binding energy 116.255 kJ/mol of the reference compound; 3,6-dimethyl- ~ {N}-[4-(trifluoromethyl)phenyl]-[1,2]oxazolo[5,4-d]pyrimidin-4-amine. Although the studied compounds are drugs, still we applied Lipinski's rules and ADMET analysis that reconfirmed that these drugs have favorable drug-like properties. In conclusion, the results of the study show that Talniflumate and Sulfaphenazole may be potential antimalarial drug candidates.The derivatives of these drugs could be designed and tested to develop better drugs against Plasmodium species.

Graphical Abstract

Structure-Based Identification of Novel Histone Deacetylase 4 (HDAC4) Inhibitors

Histone deacetylases (HDACs) are important cancer drug targets. Existing FDA-approved drugs target the catalytic pocket of HDACs, which is conserved across subfamilies (classes) of HDAC. However, engineering specificity is an important goal. Herein, we use molecular modeling approaches to identify and target potential novel pockets specific to Class IIA HDAC-HDAC4 at the interface between HDAC4 and the transcriptional corepressor component protein NCoR. These pockets were screened using an ensemble docking approach combined with consensus scoring to identify compounds with a different binding mechanism than the currently known HDAC modulators. Binding was compared in experimental assays between HDAC4 and HDAC3, which belong to a different family of HDACs. HDAC4 was significantly inhibited by compound 88402 but not HDAC3. Two other compounds (67436 and 134199) had IC50 values in the low micromolar range for both HDACs, which is comparable to the known inhibitor of HDAC4, SAHA (Vorinostat). However, both of these compounds were significantly weaker inhibitors of HDAC3 than SAHA and thus more selective, albeit to a limited extent. Five compounds exhibited activity on human breast carcinoma and/or urothelial carcinoma cell lines. The present result suggests potential mechanistic and chemical approaches for developing selective HDAC4 modulators.

Molecular Mechanisms Involved in the B Cell Growth and Clonogenic Activity of HIV-1 Matrix Protein p17 Variants

The human immunodeficiency virus (HIV-1) matrix protein p17 (p17) is released from infected cells as a protein capable of deregulating the biological activity of different cells. P17 variants (vp17s), more frequently detected in the plasma of HIV-1+ patients with rather than without lymphoma and characterized by amino acids insertions in their C-terminal region, were found to trigger B cell growth and clonogenicity. Vp17s endowed with B-cell-growth-promoting activity are drastically destabilized, whereas, in a properly folded state, reference p17 (refp17) does not exert any biological activity on B cell growth and clonogenicity. However, misfolding of refp17 is necessary to expose a masked functional epitope, interacting with the protease-activated receptor 1 (PAR-1), endowed with B cell clonogenicity. Indeed, it is worth noting that changes in the secondary structure can strongly impact the function of a protein. Here, we performed computational studies to show that the gain of function of vp17s is linked to dramatic conformational changes due to structural modification in the secondary-structure elements and in the rearrangement of the hydrogen bond (H-bond) network. In particular, all clonogenic vp17s showed the disengagement of two critical residues, namely Trp16 and Tyr29, from their hydrophobic core. Biological data showed that the mutation of Trp16 and Tyr29 to Ala in the refp17 backbone, alone or in combination, resulted in a protein endowed with B cell clonogenic activity. These data show the pivotal role of the hydrophobic component in maintaining refp17 stability and identify a novel potential therapeutic target to counteract vp17-driven lymphomagenesis in HIV-1+ patients.

Probability Density Reweighting of High-Temperature Molecular Dynamics

Molecular dynamics (MD) simulation is a popular method for elucidating the structures and functions of biomolecules. However, exploring the conformational space, especially for large systems with slow transitions, often requires enhanced sampling methods. Although conducting MD at high temperatures provides a straightforward approach, resulting conformational ensembles diverge significantly from those at low temperatures. To address this discrepancy, we propose a novel probability density-based reweighting (PDR) method. PDR exhibits robust performance across four distinct systems, including a miniprotein, a cyclic peptide, a protein loop, and a protein-peptide complex. It accurately restores the conformational distributions at high temperatures to those at low temperatures. Additionally, we apply PDR to reweight previously studied high-T MD simulations of 12 protein-peptide complexes, enabling a comprehensive investigation of the conformational space of protein-peptide complexes.

Unveiling Conformational States of CDK6 Caused by Binding of Vcyclin Protein and Inhibitor by Combining Gaussian Accelerated Molecular Dynamics and Deep Learning

CDK6 plays a key role in the regulation of the cell cycle and is considered a crucial target for cancer therapy. In this work, conformational transitions of CDK6 were identified by using Gaussian accelerated molecular dynamics (GaMD), deep learning (DL), and free energy landscapes (FELs). DL finds that the binding pocket as well as the T-loop binding to the Vcyclin protein are involved in obvious differences of conformation contacts. This result suggests that the binding pocket of inhibitors (LQQ and AP9) and the binding interface of CDK6 to the Vcyclin protein play a key role in the function of CDK6. The analyses of FELs reveal that the binding pocket and the T-loop of CDK6 have disordered states. The results from principal component analysis (PCA) indicate that the binding of the Vcyclin protein affects the fluctuation behavior of the T-loop in CDK6. Our QM/MM-GBSA calculations suggest that the binding ability of LQQ to CDK6 is stronger than AP9 with or without the binding of the Vcyclin protein. Interaction networks of inhibitors with CDK6 were analyzed and the results reveal that LQQ contributes more hydrogen binding interactions (HBIs) and hot interaction spots with CDK6. In addition, the binding pocket endures flexibility changes from opening to closing states and the Vcyclin protein plays an important role in the stabilizing conformation of the T-loop. We anticipate that this work could provide useful information for further understanding the function of CDK6 and developing new promising inhibitors targeting CDK6.

The mutation of Japanese encephalitis virus envelope protein residue 389 attenuates viral neuroinvasiveness

The envelope (E) protein of the Japanese encephalitis virus (JEV) is a key protein for virus infection and adsorption of host cells, which determines the virulence of the virus and regulates the intensity of inflammatory response. The mutation of multiple aa residues in the E protein plays a critical role in the attenuated strain of JEV. This study demonstrated that the Asp to Gly, Ser, and His mutation of the E389 site, respectively, the replication ability of the viruses in cells was significantly reduced, and the viral neuroinvasiveness was attenuated to different degrees. Among them, the mutation at E389 site enhanced the E protein flexibility contributed to the attenuation of neuroinvasiveness. In contrast, less flexibility of E protein enhanced the neuroinvasiveness of the strain. Our results indicate that the mechanism of attenuation of E389 aa mutation attenuates neuroinvasiveness is related to increased flexibility of the E protein. In addition, the increased flexibility of E protein enhanced the viral sensitivity to heparin inhibition in vitro, which may lead to a decrease in the viral load entering brain. These results suggest that E389 residue is a potential site affecting JEV virulence, and the flexibility of the E protein of aa at this site plays an important role in the determination of neuroinvasiveness.

Recent advances in computational and experimental protein-ligand affinity determination techniques

INTRODUCTION

Modern drug discovery revolves around designing ligands that target the chosen biomolecule, typically proteins. For this, the evaluation of affinities of putative ligands is crucial. This has given rise to a multitude of dedicated computational and experimental methods that are constantly being developed and improved.

AREAS COVERED

In this review, the authors reassess both the industry mainstays and the newest trends among the methods for protein - small-molecule affinity determination. They discuss both computational affinity predictions and experimental techniques, describing their basic principles, main limitations, and advantages. Together, this serves as initial guide to the currently most popular and cutting-edge ligand-binding assays employed in rational drug design.

EXPERT OPINION

The affinity determination methods continue to develop toward miniaturization, high-throughput, and in-cell application. Moreover, the availability of data analysis tools has been constantly increasing. Nevertheless, cross-verification of data using at least two different techniques and careful result interpretation remain of utmost importance.

Advanced computational approaches to understand protein aggregation

Protein aggregation is a widespread phenomenon implicated in debilitating diseases like Alzheimer's, Parkinson's, and cataracts, presenting complex hurdles for the field of molecular biology. In this review, we explore the evolving realm of computational methods and bioinformatics tools that have revolutionized our comprehension of protein aggregation. Beginning with a discussion of the multifaceted challenges associated with understanding this process and emphasizing the critical need for precise predictive tools, we highlight how computational techniques have become indispensable for understanding protein aggregation. We focus on molecular simulations, notably molecular dynamics (MD) simulations, spanning from atomistic to coarse-grained levels, which have emerged as pivotal tools in unraveling the complex dynamics governing protein aggregation in diseases such as cataracts, Alzheimer's, and Parkinson's. MD simulations provide microscopic insights into protein interactions and the subtleties of aggregation pathways, with advanced techniques like replica exchange molecular dynamics, Metadynamics (MetaD), and umbrella sampling enhancing our understanding by probing intricate energy landscapes and transition states. We delve into specific applications of MD simulations, elucidating the chaperone mechanism underlying cataract formation using Markov state modeling and the intricate pathways and interactions driving the toxic aggregate formation in Alzheimer's and Parkinson's disease. Transitioning we highlight how computational techniques, including bioinformatics, sequence analysis, structural data, machine learning algorithms, and artificial intelligence have become indispensable for predicting protein aggregation propensity and locating aggregation-prone regions within protein sequences. Throughout our exploration, we underscore the symbiotic relationship between computational approaches and empirical data, which has paved the way for potential therapeutic strategies against protein aggregation-related diseases. In conclusion, this review offers a comprehensive overview of advanced computational methodologies and bioinformatics tools that have catalyzed breakthroughs in unraveling the molecular basis of protein aggregation, with significant implications for clinical interventions, standing at the intersection of computational biology and experimental research.

Illuminating ligand-induced dynamics in nuclear receptors through MD simulations

Nuclear receptors (NRs) regulate gene expression in critical physiological processes, with their functionality finely tuned by ligand-induced conformational changes. While NRs may sometimes undergo significant conformational motions in response to ligand-binding, these effects are more commonly subtle and challenging to study by traditional structural or biophysical methods. Molecular dynamics (MD) simulations are a powerful tool to bridge the gap between static protein-ligand structures and dynamical changes that govern NR function. Here, we summarize a handful of recent studies that apply MD simulations to study NRs. We present diverse methodologies for analyzing simulation data with a detailed examination of the information each method can yield. By delving into the strengths, limitations and unique contributions of these tools, this review provides guidance for extracting meaningful data from MD simulations to advance the goal of understanding the intricate mechanisms by which ligands orchestrate a range of functional outcomes in NRs.

An alternative conformation of the N-terminal loop of human dihydroorotate dehydrogenase drives binding to a potent antiproliferative agent

Over the years, human dihydroorotate dehydrogenase (hDHODH), which is a key player in the de novo pyrimidine-biosynthesis pathway, has been targeted in the treatment of several conditions, including autoimmune disorders and acute myelogenous leukaemia, as well as in host-targeted antiviral therapy. A molecular exploration of its inhibitor-binding behaviours yielded promising candidates for innovative drug design. A detailed description of the enzymatic pharmacophore drove the decoration of well-established inhibitory scaffolds, thus gaining further in vitro and in vivo efficacy. In the present work, using X-ray crystallography, an atypical rearrangement was identified in the binding pose of a potent inhibitor characterized by a polar pyridine-based moiety (compound 18). The crystal structure shows that upon binding compound 18 the dynamics of a protein loop involved in a gating mechanism at the cofactor-binding site is modulated by the presence of three water molecules, thus fine-tuning the polarity/hydrophobicity of the binding pocket. These solvent molecules are engaged in the formation of a hydrogen-bond mesh in which one of them establishes a direct contact with the pyridine moiety of compound 18, thus paving the way for a reappraisal of the inhibition of hDHODH. Using an integrated approach, the thermodynamics of such a modulation is described by means of isothermal titration calorimetry coupled with molecular modelling. These structural insights will guide future drug design to obtain a finer Kd/logD7.4 balance and identify membrane-permeable molecules with a drug-like profile in terms of water solubility.

Protein Oxidative Modifications in Neurodegenerative Diseases: From Advances in Detection and Modelling to Their Use as Disease Biomarkers

Oxidation-reduction post-translational modifications (redox-PTMs) are chemical alterations to amino acids of proteins. Redox-PTMs participate in the regulation of protein conformation, localization and function, acting as signalling effectors that impact many essential biochemical processes in the cells. Crucially, the dysregulation of redox-PTMs of proteins has been implicated in the pathophysiology of numerous human diseases, including neurodegenerative diseases such as Alzheimer's disease and Parkinson's disease. This review aims to highlight the current gaps in knowledge in the field of redox-PTMs biology and to explore new methodological advances in proteomics and computational modelling that will pave the way for a better understanding of the role and therapeutic potential of redox-PTMs of proteins in neurodegenerative diseases. Here, we summarize the main types of redox-PTMs of proteins while providing examples of their occurrence in neurodegenerative diseases and an overview of the state-of-the-art methods used for their detection. We explore the potential of novel computational modelling approaches as essential tools to obtain insights into the precise role of redox-PTMs in regulating protein structure and function. We also discuss the complex crosstalk between various PTMs that occur in living cells. Finally, we argue that redox-PTMs of proteins could be used in the future as diagnosis and prognosis biomarkers for neurodegenerative diseases.

Exploring the activation mechanism of metabotropic glutamate receptor 2

Homomeric dimerization of metabotropic glutamate receptors (mGlus) is essential for the modulation of their functions and represents a promising avenue for the development of novel therapeutic approaches to address central nervous system diseases. Yet, the scarcity of detailed molecular and energetic data on mGlu2 impedes our in-depth comprehension of their activation process. Here, we employ computational simulation methods to elucidate the activation process and key events associated with the mGlu2, including a detailed analysis of its conformational transitions, the binding of agonists, Gi protein coupling, and the guanosine diphosphate (GDP) release. Our results demonstrate that the activation of mGlu2 is a stepwise process and several energy barriers need to be overcome. Moreover, we also identify the rate-determining step of the mGlu2's transition from the agonist-bound state to its active state. From the perspective of free-energy analysis, we find that the conformational dynamics of mGlu2's subunit follow coupled rather than discrete, independent actions. Asymmetric dimerization is critical for receptor activation. Our calculation results are consistent with the observation of cross-linking and fluorescent-labeled blot experiments, thus illustrating the reliability of our calculations. Besides, we also identify potential key residues in the Gi protein binding position on mGlu2, mGlu2 dimer's TM6-TM6 interface, and Gi α5 helix by the change of energy barriers after mutation. The implications of our findings could lead to a more comprehensive grasp of class C G protein-coupled receptor activation.

Ligand Gaussian accelerated Molecular Dynamics 3 (LiGaMD3): Improved Calculations of Binding Thermodynamics and Kinetics of Both Small Molecules and Flexible Peptides

Binding thermodynamics and kinetics play critical roles in drug design. However, it has proven challenging to efficiently predict ligand binding thermodynamics and kinetics of small molecules and flexible peptides using conventional Molecular Dynamics (cMD), due to limited simulation timescales. Based on our previously developed Ligand Gaussian accelerated Molecular Dynamics (LiGaMD) method, we present a new approach, termed "LiGaMD3", in which we introduce triple boosts into three individual energy terms that play important roles in small-molecule/peptide dissociation, rebinding and system conformational changes to improve the sampling efficiency of small-molecule/peptide interactions with target proteins. To validate the performance of LiGaMD3, MDM2 bound by a small molecule (Nutlin 3) and two highly flexible peptides (PMI and P53) were chosen as model systems. LiGaMD3 could efficiently capture repetitive small-molecule/peptide dissociation and binding events within 2 microsecond simulations. The predicted binding kinetic constant rates and free energies from LiGaMD3 agreed with available experimental values and previous simulation results. Therefore, LiGaMD3 provides a more general and efficient approach to capture dissociation and binding of both small-molecule ligand and flexible peptides, allowing for accurate prediction of their binding thermodynamics and kinetics.

The use of collective variables and enhanced sampling in the simulations of existing and emerging microporous materials

Microporous materials, including zeolites, metal-organic frameworks, and cage compounds, offer diverse functionalities due to their unique dynamics and guest confinement properties. These materials play a significant role in separation, catalysis, and sensing, but their complexity hinders exploration using traditional atomistic simulations. This review explores collective variables (CVs) paired with enhanced sampling as a powerful approach to enable efficient investigation of key features in microporous materials. We highlight successful applications of CVs in studying adsorption, diffusion, phase transitions, and mechanical properties, demonstrating their crucial role in guiding material design and optimisation. The future of CVs lies in integration with techniques like machine learning, allowing for enhanced efficiency and accuracy. By tailoring CVs to specific materials and developing multi-scale approaches we can further unlock the intricacies of these fascinating materials. Simulations are a cornerstone in unravelling the complexities of microporous materials and are crucial for our future understanding.

PaCS-Toolkit: Optimized Software Utilities for Parallel Cascade Selection Molecular Dynamics (PaCS-MD) Simulations and Subsequent Analyses

Parallel cascade selection molecular dynamics (PaCS-MD) is an enhanced conformational sampling method conducted as a "repetition of time leaps in parallel worlds", comprising cycles of multiple molecular dynamics (MD) simulations performed in parallel and selection of the initial structures of MDs for the next cycle. We developed PaCS-Toolkit, an optimized software utility enabling the use of different MD software and trajectory analysis tools to facilitate the execution of the PaCS-MD simulation and analyze the obtained trajectories, including the preparation for the subsequent construction of the Markov state model. PaCS-Toolkit is coded with Python, is compatible with various computing environments, and allows for easy customization by editing the configuration file and specifying the MD software and analysis tools to be used. We present the software design of PaCS-Toolkit and demonstrate applications of PaCS-MD variations: original targeted PaCS-MD to peptide folding; rmsdPaCS-MD to protein domain motion; and dissociation PaCS-MD to ligand dissociation from adenosine A2A receptor.

Nuclear receptor interdomain communication is mediated by the hinge with ligand specificity

Nuclear receptors are ligand-induced transcription factors that bind directly to target genes and regulate their expression. Ligand binding initiates conformational changes that propagate to other domains, allosterically regulating their activity. The nature of this interdomain communication in nuclear receptors is poorly understood, largely owing to the difficulty of experimentally characterizing full-length structures. We have applied computational modeling approaches to describe and study the structure of the full length farnesoid X receptor (FXR), approximated by the DNA binding domain (DBD) and ligand binding domain (LBD) connected by the flexible hinge region. Using extended molecular dynamics simulations (> 10 microseconds) and enhanced sampling simulations, we provide evidence that ligands selectively induce domain rearrangement, leading to interdomain contact. We use protein-protein interaction assays to provide experimental evidence of these interactions, identifying a critical role of the hinge in mediating interdomain contact. Our results illuminate previously unknown aspects of interdomain communication in FXR and provide a framework to enable characterization of other full length nuclear receptors.

Exploring biased activation characteristics by molecular dynamics simulation and machine learning for the μ-opioid receptor

Biased ligands selectively activating specific downstream signaling pathways (termed as biased activation) exhibit significant therapeutic potential. However, the conformational characteristics revealed are very limited for the biased activation, which is not conducive to biased drug development. Motivated by the issue, we combine extensive accelerated molecular dynamics simulations and an interpretable deep learning model to probe the biased activation features for two complex systems constructed by the inactive μOR and two different biased agonists (G-protein-biased agonist TRV130 and β-arrestin-biased agonist endomorphin2). The results indicate that TRV130 binds deeper into the receptor core compared to endomorphin2, located between W2936.48 and D1142.50, and forms hydrogen bonding with D1142.50, while endomorphin2 binds above W2936.48. The G protein-biased agonist induces greater outward movements of the TM6 intracellular end, forming a typical active conformation, while the β-arrestin-biased agonist leads to a smaller extent of outward movements of TM6. Compared with TRV130, endomorphin2 causes more pronounced inward movements of the TM7 intracellular end and more complex conformational changes of H8 and ICL1. In addition, important residues determining the two different biased activation states were further identified by using an interpretable deep learning classification model, including some common biased activation residues across Class A GPCRs like some key residues on the TM2 extracellular end, ECL2, TM5 intracellular end, TM6 intracellular end, and TM7 intracellular end, and some specific important residues of ICL3 for μOR. The observations will provide valuable information for understanding the biased activation mechanism for GPCRs.

Flanking Domains Modulate α-Synuclein Monomer Structure: A Molecular Dynamics Domain Deletion Study

Aggregates of misfolded α-synuclein proteins (asyn) are key markers of Parkinson's disease. Asyn proteins have three domains: an N-terminal domain, a hydrophobic NAC core implicated in aggregation, and a proline-rich C-terminal domain. Proteins with truncated C-terminal domains are known to be prone to aggregation and suggest that studying domain-domain interactions in asyn monomers could help elucidate the role of the flanking domains in modulating protein structure. To this end, we used Gaussian accelerated molecular dynamics (GAMD) to simulate wild-type (WT), N-terminal truncated (DN), C-terminal truncated (ΔC), and isolated NAC domain variants (isoNAC). Using clustering and contact analysis, we found that N- and C-terminal domains interact via electrostatic interactions, while the NAC and N-terminal domains interact through hydrophobic contacts. Our work also suggests that the C-terminal domain does not interact directly with the NAC domain but instead interacts with the N-terminal domain. Removal of the N-terminal domain led to increased contacts between NAC and C-terminal domains and the formation of interdomain β-sheets. Removal of either flanking domain also resulted in increased compactness of every domain. We also found that the contacts between flanking domains results in an electrostatic potential (ESP) that could possibly lead to favorable interactions with anionic lipid membranes. Removal of the C-terminal domain disrupts the ESP in a way that is likely to over-stabilize protein-membrane interactions. All of this suggests that one of the roles of the flanking domains may be to modulate the protein structure in a way that helps maintain elongation, hide hydrophobic residue from the solvent, and maintain an ESP that aids favorable interactions with the membrane.

How nuclear receptors transition between active and inactive forms: An energetic perspective

Nuclear receptors regulate transcriptional programs in response to the binding of natural and synthetic ligands. These ligands modulate the receptor by inducing dynamic changes in the ligand binding domain that shift the C-terminal helix (H12) between active and inactive conformations. Despite decades of study, many questions persist regarding the nature of the inactive state and how ligands shift receptors between different states. Here, we use molecular dynamics (MD) simulations to investigate the timescale and energetic landscape of the conformational transition between inactive and active forms of progesterone receptor (PR) bound to a partial agonist. We observe that the microsecond timescale is insufficient to observe any transitions; only at millisecond timescales achieved via accelerated MD simulations do we find the inactive PR switches to the active state. Energetic analysis reveals that both active and inactive PR states represent energy minima separated by a barrier that can be traversed. In contrast, little or no transition is observed between active and inactive states when an agonist or antagonist is bound, confirming that ligand identity plays a key role in defining the energy landscape of nuclear receptor conformations.

Exploring Chemical Space Using Ab Initio Hyperreactor Dynamics

In recent years, first-principles exploration of chemical reaction space has provided valuable insights into intricate reaction networks. Here, we introduce ab initio hyperreactor dynamics, which enables rapid screening of the accessible chemical space from a given set of initial molecular species, predicting new synthetic routes that can potentially guide subsequent experimental studies. For this purpose, different hyperdynamics derived bias potentials are applied along with pressure-inducing spherical confinement of the molecular system in ab initio molecular dynamics simulations to efficiently enhance reactivity under mild conditions. To showcase the advantages and flexibility of the hyperreactor approach, we present a systematic study of the method's parameters on a HCN toy model and apply it to a recently introduced experimental model for the prebiotic formation of glycinal and acetamide in interstellar ices, which yields results in line with experimental findings. In addition, we show how the developed framework enables the study of complicated transitions like the first step of a nonenzymatic DNA nucleoside synthesis in an aqueous environment, where the molecular fragmentation problem of earlier nanoreactor approaches is avoided.

Fractional Extended Diffusion Theory to capture anomalous relaxation from biased/accelerated molecular simulations

Biased and accelerated molecular simulations (BAMS) are widely used tools to observe relevant molecular phenomena occurring on time scales inaccessible to standard molecular dynamics, but evaluation of the physical time scales involved in the processes is not directly possible from them. For this reason, the problem of recovering dynamics from such kinds of simulations is the object of very active research due to the relevant theoretical and practical implications of dynamics on the properties of both natural and synthetic molecular systems. In a recent paper [A. Rapallo et al., J. Comput. Chem. 42, 586-599 (2021)], it has been shown how the coupling of BAMS (which destroys the dynamics but allows to calculate average properties) with Extended Diffusion Theory (EDT) (which requires input appropriate equilibrium averages calculated over the BAMS trajectories) allows to effectively use the Smoluchowski equation to calculate the orientational time correlation function of the head-tail unit vector defined over a peptide in water solution. Orientational relaxation of this vector is the result of the coupling of internal molecular motions with overall molecular rotation, and it was very well described by correlation functions expressed in terms of weighted sums of suitable time-exponentially decaying functions, in agreement with a Brownian diffusive regime. However, situations occur where exponentially decaying functions are no longer appropriate to capture the actual dynamical behavior, which exhibits persistent long time correlations, compatible with the so called subdiffusive regimes. In this paper, a generalization of EDT will be given, exploiting a fractional Smoluchowski equation (FEDT) to capture the non-exponential character observed in the relaxation of intramolecular distances and molecular radius of gyration, whose dynamics depend on internal molecular motions only. The calculation methods, proper to EDT, are adapted to implement the generalization of the theory, and the resulting algorithm confirms FEDT as a tool of practical value in recovering dynamics from BAMS, to be used in general situations, involving both regular and anomalous diffusion regimes.

Collective Variable-Based Enhanced Sampling: From Human Learning to Machine Learning

Enhanced-sampling algorithms relying on collective variables (CVs) are extensively employed to study complex (bio)chemical processes that are not amenable to brute-force molecular simulations. The selection of appropriate CVs characterizing the slow movement modes is of paramount importance for reliable and efficient enhanced-sampling simulations. In this Perspective, we first review the application and limitations of CVs obtained from chemical and geometrical intuition. We also introduce path-sampling algorithms, which can identify path-like CVs in a high-dimensional free-energy space. Machine-learning algorithms offer a viable approach to finding suitable CVs by analyzing trajectories from preliminary simulations. We discuss both the performance of machine-learning-derived CVs in enhanced-sampling simulations of experimental models and the challenges involved in applying these CVs to realistic, complex molecular assemblies. Moreover, we provide a prospective view of the potential advancements of machine-learning algorithms for the development of CVs in the field of enhanced-sampling simulations.

Exploring the DNA2-PNA heterotriplex formation in targeting the Bcl-2 gene promoter: A structural insight by physico-chemical and microsecond-scale MD investigation

Peptide Nucleic Acids (PNAs) represent a promising tool for gene modulation in anticancer treatment. The uncharged peptidyl backbone and the resistance to chemical and enzymatic degradation make PNAs highly advantageous to form stable hybrid complexes with complementary DNA and RNA strands, providing higher stability than the corresponding natural analogues. Our and other groups' research has successfully shown that tailored PNA sequences can effectively downregulate the expression of human oncogenes using antigene, antisense, or anti-miRNA approaches. Specifically, we identified a seven bases-long PNA sequence, complementary to the longer loop of the main G-quadruplex structure formed by the bcl2midG4 promoter sequence, capable of downregulating the expression of the antiapoptotic Bcl-2 protein and enhancing the anticancer activity of an oncolytic adenovirus. Here, we extended the length of the PNA probe with the aim of including the double-stranded Bcl-2 promoter among the targets of the PNA probe. Our investigation primarily focused on the structural aspects of the resulting DNA2-PNA heterotriplex that were determined by employing conventional and accelerated microsecond-scale molecular dynamics simulations and chemical-physical analysis. Additionally, we conducted preliminary biological experiments using cytotoxicity assays on human A549 and MDA-MB-436 adenocarcinoma cell lines, employing the oncolytic adenovirus delivery strategy.

Molecular mechanism of specific HLA-A mRNA recognition by the RNA-binding-protein hMEX3B to promote tumor immune escape

Immunotherapy, including immune checkpoint inhibitors and adoptive cell transfer, has obtained great progress, but their efficiencies vary among patients due to the genetic and epigenetic differences. Human MEX3B (hMEX3B) protein is an RNA-binding protein that contains two KH domains at the N-terminus and a RING domain at its C-terminus, which has the activity of E3 ubiquitin ligase and is essential for RNA degradation. Current evidence suggests that hMEX3B is involved in many important biological processes, including tumor immune evasion and HLA-A regulation, but the sequence of substrate RNA recognized by hMEX3B and the functional molecular mechanisms are unclear. Here, we first screened the optimized hMEX3B binding sequence on the HLA-A mRNA and reported that the two tandem KH domains can bind with their substrate one hundred times more than the individual KH domains. We systematically investigated the binding characteristics between the two KH domains and their RNA substrates by nuclear magnetic resonance (NMR). Based on this information and the small-angle X-ray scattering (SAXS) data, we used molecular dynamics simulations to obtain structural models of KH domains in complex with their corresponding RNAs. By analyzing the models, we noticed that on the KH domains' variable loops, there were two pairs of threonines and arginines that can disrupt the recognition of the RNA completely, and this influence had also been verified both in vitro and in vivo. Finally, we presented a functional model of the hMEX3B protein, which indicated that hMEX3B regulated the degradation of its substrate mRNAs in many biological processes. Taken together, our research illustrated how the hMEX3B protein played a key role in translation inhibition during the immune response to tumor cells and provided an idea and a lead for the study of the molecular mechanism and function of other MEX3 family proteins.

Structural studies of catalytic peptides using molecular dynamics simulations

Many self-assembling peptides can form amyloid like structures with different sizes and morphologies. Driven by non-covalent interactions, their aggregation can occur through distinct pathways. Additionally, they can bind metal ions to create enzyme like active sites that allow them to catalyze diverse reactions. Due to the non-crystalline nature of amyloids, it is quite challenging to elucidate their structures using experimental spectroscopic techniques. In this aspect, molecular dynamics (MD) simulations provide a useful tool to derive structures of these macromolecules in solution. They can be further validated by comparing with experimentally measured structural parameters. However, these simulations require a multi-step process starting from the selection of the initial structure to the analysis of MD trajectories. There are multiple force fields, parametrization protocols, equilibration processes, software and analysis tools available for this process. Therefore, it is complicated for non-experts to select the most relevant tools and perform these simulations effectively. In this chapter, a systematic methodology that covers all major aspects of modeling of catalytic peptides is provided in a user-friendly manner. It will be helpful for researchers in this critical area of research.

Molecular modeling of ARF6 dysregulation caused by mutations in IQSEC2

IQSEC2 gene mutations are associated with epilepsy, autism, and intellectual disability. The primary function IQSEC2, mediated via its Sec 7 domain, is to act as a guanine nucleotide exchange factor for ARF6. We sought to develop a molecular model, which may explain the aberrant Sec 7 activity on ARF6 of different human IQSEC2 mutations. We integrated experimental data of IQSEC2 mutants with protein structure prediction by the RaptorX server combined with molecular modeling and molecular dynamics simulations. Normally, apocalmodulin (apoCM) binds to IQSEC2 resulting in its N-terminal fragment inhibiting access of its Sec 7 domain to ARF6. An increase in Ca2+ concentration destabilizes the interaction of IQSEC2 with apoCM and removes steric hindrance of Sec 7 binding with ARF6. Mutations at amino acid residue 350 of IQSEC2 result in loss of steric hindrance of Sec 7 binding with ARF6 leading to constitutive activation of ARF6 by Sec 7. On the other hand, a mutation at amino acid residue 359 of IQSEC2 results in constitutive hindrance of Sec 7 binding to ARF6 leading to the loss of the ability of IQSEC2 to activate ARF6. These studies provide a model for dysregulation of IQSEC2 Sec 7 activity by mutant IQSEC2 proteins.Communicated by Ramaswamy H. Sarma.

Accelerated Molecular Dynamics and AlphaFold Uncover a Missing Conformational State of Transporter Protein OxlT

Transporter proteins change their conformations to carry their substrate across the cell membrane. The conformational dynamics is vital to understanding the transport function. We have studied the oxalate transporter (OxlT), an oxalate:formate antiporter from Oxalobacter formigenes, significant in avoiding kidney stone formation. The atomic structure of OxlT has been recently solved in the outward-open and occluded states. However, the inward-open conformation is still missing, hindering a complete understanding of the transporter. Here, we performed a Gaussian accelerated molecular dynamics simulation to sample the extensive conformational space of OxlT and successfully predicted the inward-open conformation where cytoplasmic substrate formate binding was preferred over oxalate binding. We also identified critical interactions for the inward-open conformation. The results were complemented by an AlphaFold2 structure prediction. Although AlphaFold2 solely predicted OxlT in the outward-open conformation, mutation of the identified critical residues made it partly predict the inward-open conformation, identifying possible state-shifting mutations.

Molecular Dynamics Simulations for Rationalizing Polymer Bioconjugation Strategies: Challenges, Recent Developments, and Future Opportunities

The covalent modification of proteins with polymers is a well-established method for improving the pharmacokinetic properties of therapeutically valuable biologics. The conjugated polymer chains of the resulting hybrid represent highly flexible macromolecular structures. As the dynamics of such systems remain rather elusive for established experimental techniques from the field of protein structure elucidation, molecular dynamics simulations have proven as a valuable tool for studying such conjugates at an atomistic level, thereby complementing experimental studies. With a focus on new developments, this review aims to provide researchers from the polymer bioconjugation field with a concise and up to date overview of such approaches. After introducing basic principles of molecular dynamics simulations, as well as methods for and potential pitfalls in modeling bioconjugates, the review illustrates how these computational techniques have contributed to the understanding of bioconjugates and bioconjugation strategies in the recent past and how they may lead to a more rational design of novel bioconjugates in the future.

Accelerating Molecular Dynamics Simulations for Drug Discovery

Accurate prediction of ligand binding thermodynamics and kinetics is crucial in drug design. However, it remains challenging for conventional molecular dynamics (MD) simulations due to sampling issues. Gaussian accelerated MD (GaMD) is an enhanced sampling method that adds a harmonic boost to overcome energy barriers, which has demonstrated significant benefits in exploring protein-ligand interactions. Especially, the ligand GaMD (LiGaMD) applies a selective boost potential to the ligand nonbonded potential energy, significantly improving sampling for ligand binding and dissociation. Furthermore, a selective boost potential is applied to the potential of both ligand and protein residues around binding pocket in LiGaMD2 to further increase the sampling of protein-ligand interaction. LiGaMD and LiGaMD2 simulations could capture repetitive ligand binding and unbinding events within microsecond simulations, allowing to simultaneously characterize ligand binding thermodynamics and kinetics, which is expected to greatly facilitate drug design. In this chapter, we provide a brief review of the status of LiGaMD in drug discovery and outline its usage.

Roles of Accelerated Molecular Dynamics Simulations in Predictions of Binding Kinetic Parameters

Rational predictions on binding kinetics parameters of drugs to targets play significant roles in future drug designs. Full conformational samplings of targets are requisite for accurate predictions of binding kinetic parameters. In this review, we mainly focus on the applications of enhanced sampling technologies in calculations of binding kinetics parameters and residence time of drugs. The methods involved in molecular dynamics simulations are applied to not only probe conformational changes of targets but also reveal calculations of residence time that is significant for drug efficiency. For this review, special attention are paid to accelerated molecular dynamics (aMD) and Gaussian aMD (GaMD) simulations that have been adopted to predict the association or disassociation rate constant. We also expect that this review can provide useful information for future drug design.

Addressing Challenges of Macrocyclic Conformational Sampling in Polar and Apolar Solvents: Lessons for Chameleonicity

We evaluated a workflow to reliably sample the conformational space of a set of 47 peptidic macrocycles. Starting from SMILES strings, we use accelerated molecular dynamics simulations to overcome high energy barriers, in particular, the cis-trans isomerization of peptide bonds. We find that our approach performs very well in polar solvents like water and dimethyl sulfoxide. Interestingly, the protonation state of a secondary amine in the ring only slightly influences the conformational ensembles of our test systems. For several of the macrocycles, determining the conformational distribution in chloroform turns out to be considerably more challenging. Especially, the choice of partial charges crucially influences the ensembles in chloroform. We address these challenges by modifying initial structures and the choice of partial charges. Our results suggest that special care has to be taken to understand the configurational distribution in apolar solvents, which is a key step toward a reliable prediction of membrane permeation of macrocycles and their chameleonic properties.

Deep-Learning-Assisted Enhanced Sampling for Exploring Molecular Conformational Changes

We present a novel strategy to explore conformational changes and identify stable states of molecular objects, eliminating the need for a priori knowledge. The approach applies a deep learning method to extract information about the movement modes of the molecular object from a short, high-dimensional, and parameter-free preliminary enhanced-sampling simulation. The gathered information is described by a small set of deep-learning-based collective variables (dCVs), which steer the production-enhanced-sampling simulation. Considering the challenge of adequately exploring the configurational space using the low-dimensional, suboptimal dCVs, we incorporate a method designed for ergodic sampling, namely, Gaussian-accelerated molecular dynamics (MD), into the framework of CV-based enhanced sampling. MD simulations on both toy models and nontrivial examples demonstrate the remarkable computational efficiency of the strategy in capturing the conformational changes of molecular objects without a priori knowledge. Specifically, we achieved the blind folding of two fast folders, chignolin and villin, within a time scale of hundreds of nanoseconds and successfully reconstructed the free-energy landscapes that characterize their reversible folding. All in all, the presented strategy holds significant promise for investigating conformational changes in macromolecules, and it is anticipated to find extensive applications in the fields of chemistry and biology.

QMrebind: incorporating quantum mechanical force field reparameterization at the ligand binding site for improved drug-target kinetics through milestoning simulations

Understanding the interaction of ligands with biomolecules is an integral component of drug discovery and development. Challenges for computing thermodynamic and kinetic quantities for pharmaceutically relevant receptor-ligand complexes include the size and flexibility of the ligands, large-scale conformational rearrangements of the receptor, accurate force field parameters, simulation efficiency, and sufficient sampling associated with rare events. Our recently developed multiscale milestoning simulation approach, SEEKR2 (Simulation Enabled Estimation of Kinetic Rates v.2), has demonstrated success in predicting unbinding (koff) kinetics by employing molecular dynamics (MD) simulations in regions closer to the binding site. The MD region is further subdivided into smaller Voronoi tessellations to improve the simulation efficiency and parallelization. To date, all MD simulations are run using general molecular mechanics (MM) force fields. The accuracy of calculations can be further improved by incorporating quantum mechanical (QM) methods into generating system-specific force fields through reparameterizing ligand partial charges in the bound state. The force field reparameterization process modifies the potential energy landscape of the bimolecular complex, enabling a more accurate representation of the intermolecular interactions and polarization effects at the bound state. We present QMrebind (Quantum Mechanical force field reparameterization at the receptor-ligand binding site), an ORCA-based software that facilitates reparameterizing the potential energy function within the phase space representing the bound state in a receptor-ligand complex. With SEEKR2 koff estimates and experimentally determined kinetic rates, we compare and interpret the receptor-ligand unbinding kinetics obtained using the newly reparameterized force fields for model host-guest systems and HSP90-inhibitor complexes. This method provides an opportunity to achieve higher accuracy in predicting receptor-ligand koff rate constants.