Implementations of Nosé-Hoover and Nosé-Poincaré thermostats in mesoscopic dynamic simulations with the united-residue model of a polypeptide chain

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Peptide ligands for the affinity purification of adeno-associated viruses from HEK 293 cell lysates

Adeno-associated viruses (AAVs) are the vector of choice for delivering gene therapies that can cure inherited and acquired diseases. Clinical research on various AAV serotypes significantly increased in recent years alongside regulatory approvals of AAV-based therapies. The current AAV purification platform hinges on the capture step, for which several affinity resins are commercially available. These adsorbents rely on protein ligands-typically camelid antibodies-that provide high binding capacity and selectivity, but suffer from low biochemical stability and high cost, and impose harsh elution conditions (pH < 3) that can harm the transduction activity of recovered AAVs. Addressing these challenges, this study introduces peptide ligands that selectively capture AAVs and release them under mild conditions (pH = 6.0). The peptide sequences were identified by screening a focused library and modeled in silico against AAV serotypes 2 and 9 (AAV2 and AAV9) to select candidate ligands that target homologous sites at the interface of the VP1-VP2 and VP2-VP3 virion proteins with mild binding strength (KD  ~ 10-5 -10- 6  M). Selected peptides were conjugated to Toyopearl resin and evaluated via binding studies against AAV2 and AAV9, demonstrating the ability to target both serotypes with values of dynamic binding capacity (DBC10%  > 1013 vp/mL of resin) and product yields (~50%-80%) on par with commercial adsorbents. The peptide-based adsorbents were finally utilized to purify AAV2 from a HEK 293 cell lysate, affording high recovery (50%-80%), 80- to 400-fold reduction of host cell proteins (HCPs), and high transduction activity (up to 80%) of the purified viruses.

Molecular recognition of an aversive odorant by the murine trace amine-associated receptor TAAR7f

There are two main families of G protein-coupled receptors that detect odours in humans, the odorant receptors (ORs) and the trace amine-associated receptors (TAARs). Their amino acid sequences are distinct, with the TAARs being most similar to the aminergic receptors such as those activated by adrenaline, serotonin and histamine. To elucidate the structural determinants of ligand recognition by TAARs, we have determined the cryo-EM structure of a murine receptor, mTAAR7f, coupled to the heterotrimeric G protein Gs and bound to the odorant N,N-dimethylcyclohexylamine (DMCH) to an overall resolution of 2.9 Å. DMCH is bound in a hydrophobic orthosteric binding site primarily through van der Waals interactions and a strong charge-charge interaction between the tertiary amine of the ligand and an aspartic acid residue. This site is distinct and non-overlapping with the binding site for the odorant propionate in the odorant receptor OR51E2. The structure, in combination with mutagenesis data and molecular dynamics simulations suggests that the activation of the receptor follows a similar pathway to that of the β-adrenoceptors, with the significant difference that DMCH interacts directly with one of the main activation microswitch residues.

Molecular engineering of cyclic azobenzene-peptide hybrid ligands for the purification of human blood Factor VIII via photo-affinity chromatography

The use of benign stimuli to control the binding and release of labile biologics for their isolation from complex feedstocks is a key goal of modern biopharmaceutical technology. This study introduces cyclic azobenzene-peptide (CAP) hybrid ligands for the rapid and discrete photo-responsive capture and release of blood coagulation Factor VIII (FVIII). A predictive method - based on amino acid sequence and molecular architecture of CAPs - was developed to correlate the conformation of cis/trans CAP photo-isomers to FVIII binding and release. The combined in silico and in vitro analysis of FVIII:peptide interactions guided the design of a rational approach to optimize isomerization kinetics and biorecognition of CAPs. A photoaffinity adsorbent, prepared by conjugating selected CAP G-cycloAZOB[Lys-YYKHLYN-Lys]-G on translucent chromatographic beads, featured high binding capacity (> 6 mg of FVIII per mL of resin) and rapid photo-isomerization kinetics (τ < 30s) when exposed to 420-450 nm light at the intensity of 0.1 W·cm-2. The adsorbent purified FVIII from a recombinant harvest using a single mobile phase, affording high product yield (>90%), purity (>95%), and blood clotting activity. The CAPs introduced in this report demonstrate a novel route integrating gentle operational conditions in a rapid and efficient bioprocess for the purification of life-saving biotherapeutics.

AmberMDrun: A Scripting Tool for Running Amber MD in an Easy Way

MD simulations have been widely applied and become a powerful tool in the field of biomacromolecule simulations and computer-aided drug design, etc., which can estimate binding free energy between receptor and ligand. However, the inputs and force field preparation for performing Amber MD is somewhat complicated, and challenging for beginners. To address this issue, we have developed a script for automatically preparing Amber MD input files, balancing the system, performing Amber MD for production, and predicting receptor-ligand binding free energy. This script is open-source, extensible and can support customization. The core code is written in C++ and has a Python interface, providing both efficient performance and convenience.

Optimization of parallel implementation of UNRES package for coarse-grained simulations to treat large proteins

We report major algorithmic improvements of the UNRES package for physics-based coarse-grained simulations of proteins. These include (i) introduction of interaction lists to optimize computations, (ii) transforming the inertia matrix to a pentadiagonal form to reduce computing and memory requirements, (iii) removing explicit angles and dihedral angles from energy expressions and recoding the most time-consuming energy/force terms to minimize the number of operations and to improve numerical stability, (iv) using OpenMP to parallelize those sections of the code for which distributed-memory parallelization involves unfavorable computing/communication time ratio, and (v) careful memory management to minimize simultaneous access of distant memory sections. The new code enables us to run molecular dynamics simulations of protein systems with size exceeding 100,000 amino-acid residues, reaching over 1 ns/day (1 μs/day in all-atom timescale) with 24 cores for proteins of this size. Parallel performance of the code and comparison of its performance with that of AMBER, GROMACS and MARTINI 3 is presented.

Development of peptide affinity ligands for the purification of polyclonal and monoclonal Fabs from recombinant fluids

Engineered multi-specific monoclonal antibodies (msAbs) and antibody fragments offer valuable therapeutic options against metabolic disorders, aggressive cancers, and viral infections. The advancement in molecular design and recombinant expression of these next-generation drugs, however, is not equaled by the progress in downstream bioprocess technology. The purification of msAbs and fragments requires affinity adsorbents with orthogonal biorecognition of different portions of the antibody structure, namely its Fc (fragment crystallizable) and Fab (fragment antigen-binding) regions or the C_H1-3 and C_L chains. Current adsorbents rely on protein ligands that, while featuring high binding capacity and selectivity, need harsh elution conditions and suffer from high cost, limited biochemical stability, and potential release of immunogenic fragments. Responding to these challenges, we undertook the de novo discovery of peptide ligands that target different regions of human Fab and enable product release under mild conditions. The ligands were discovered by screening a focused library of 12-mer peptides against a feedstock comprising human Fab and Chinese hamster ovary host cell proteins (CHO HCPs). The identified ligands were evaluated via binding studies as well as molecular docking simulations, returning excellent values of binding capacity (Q_max ∼ 20 mg of Fab per mL of resin) and dissociation constant (K_D = 2.16·10^-6 M). Selected ligand FRWNFHRNTFFP and commercial Protein L ligands were further characterized by measuring the dynamic binding capacity (DBC_10%) at different residence times (RT) and performing the purification of polyclonal and monoclonal Fabs from CHO-K1 cell culture fluids. The peptide ligand featured DBC_10% ∼ 6-16 mg/mL (RT of 2 min) and afforded values of yield (93-96%) and purity (89-96%) comparable to those provided by Protein L resins.

Modeling the Structure, Dynamics, and Transformations of Proteins with the UNRES Force Field

The physics-based united-residue (UNRES) model of proteins ( www.unres.pl ) has been designed to carry out large-scale simulations of protein folding. The force field has been derived and parameterized based on the principles of statistical-mechanics, which makes it independent of structural databases and applicable to treat nonstandard situations such as, proteins that contain D-amino-acid residues. Powered by Langevin dynamics and its replica-exchange extensions, UNRES has found a variety of applications, including ab initio and database-assisted protein-structure prediction, simulating protein-folding pathways, exploring protein free-energy landscapes, and solving biological problems. This chapter provides a summary of UNRES and a guide for potential users regarding the application of the UNRES package in a variety of research tasks.

Probing Protein Aggregation Using the Coarse-Grained UNRES Force Field

Protein aggregation is the cause of many, often lethal, diseases, including the Alzheimer's, Parkinson's, and Huntington's diseases, and familial amyloidosis. Theoretical investigation of the mechanism of this process, including the structures of the oligomeric intermediates which are the most toxic, is difficult because of long time scale of aggregation. Coarse-grained models, which enable us to extend the simulation time scale by three or more orders of magnitude, are, therefore, of great advantage in such studies. In this chapter, we describe the application of the physics-based UNited RESidue (UNRES) force field developed in our laboratory to study protein aggregation, in both free simulations and simulations of aggregation propagation from an existing template (seed), and illustrate it with the examples of Aβ-peptide aggregation and Aβ-peptide-assisted aggregation of the peptides derived from the repeat domains of tau (TauRD).

Characterization of folic acid-functionalized PLA-PEG nanomicelle to deliver Letrozole: A nanoinformatics study

Effective cancer treatment is currently the number one challenge to human health. To date, several treatment methods have been introduced for cancer cell targeting. Among the proposed new methods for attacking cancer cells, nanotechnology has attracted much attention. Hence, various nanocarriers have been developed for targeted delivery of available drugs and improve their effectiveness against malignant cells. The PLA-PEG functionalised with folic acid (PLA-PEG-FA) is one of the nanocarriers with a limited range of applications for targeting cancer cells. In this investigation, different types of in-silico methods such as molecular docking approach, molecular dynamics simulation and free energy calculations are employed to characterise the carriers studied. The effectiveness of PLA-PEG-FA and PLA-PEG in delivering Letrozole as an aromatase inhibitor in cancer cells is examined. It is found that in the presence of folic acid, the stability and cell membrane permeability of nanomicelle are increased. Therefore, PLA-PEG-FA can be considered as a versatile carrier that can increase the effectiveness of aromatase inhibitors (such as Letrozole) and reduce their side effects.

Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems

Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.

Extension of coarse-grained UNRES force field to treat carbon nanotubes

Carbon nanotubes (CNTs) have recently received considerable attention because of their possible applications in various branches of nanotechnology. For their cogent application, knowledge of their interactions with biological macromolecules, especially proteins, is essential and computer simulations are very useful for such studies. Classical all-atom force fields limit simulation time scale and size of the systems significantly. Therefore, in this work, we implemented CNTs into the coarse-grained UNited RESidue (UNRES) force field. A CNT is represented as a rigid infinite-length cylinder which interacts with a protein through the Kihara potential. Energy conservation in microcanonical coarse-grained molecular dynamics simulations and temperature conservation in canonical simulations with UNRES containing the CNT component have been verified. Subsequently, studies of three proteins, bovine serum albumin (BSA), soybean peroxidase (SBP), and α-chymotrypsin (CT), with and without CNTs, were performed to examine the influence of CNTs on the structure and dynamics of these proteins. It was found that nanotubes bind to these proteins and influence their structure. Our results show that the UNRES force field can be used for further studies of CNT-protein systems with 3-4 order of magnitude larger timescale than using regular all-atom force fields. Graphical abstract Bovine serum albumin (BSA), soybean peroxidase (SBP), and α-chymotrypsin (CT), with and without CNTsᅟ.

Charge transfer dynamics at the boron subphthalocyanine chloride/C60 interface: non-adiabatic dynamics study with Libra-X

The Libra-X software for non-adiabatic molecular dynamics is reported. It is used to comprehensively study the charge transfer dynamics at the boron subphtalocyanine chloride (SubPc)/fullerene (C₆₀) interface.

Probing the intermolecular interactions of PPARγ-LBD with polyunsaturated fatty acids and their anti-inflammatory metabolites to infer most potential binding moieties

Background

PPARγ is an isoform of peroxisome proliferator-activated receptor (PPAR) belonging to a super family of nuclear receptors. PPARγ receptor is found to play a crucial role in the modulation of lipid and glucose homeostasis. Its commotion has been reported to play a significant role in a broad spectrum of diseases such as type 2 diabetes mellitus, inflammatory diseases, Alzheimer’s disease, and in some cancers. Hence, PPARγ is an important therapeutic target. Polyunsaturated fatty acids (PUFAs) and their metabolites (henceforth referred to as bioactive lipids) are known to function as agonists of PPARγ. However, agonistic binding modes and affinity of these ligands to PPARγ are yet to be deciphered.

Methods

In this study, we performed a comparative molecular docking, binding free energy calculation and molecular dynamics simulation to infer and rank bioactive lipids based on the binding affinities with the ligand binding domain (LBD) of PPARγ.

Results

The results inferred affinity in the order of resolvin E1 > neuroprotectin D1 > hydroxy-linoleic acid > docosahexaenoic acid > lipoxin A4 > gamma-linolenic acid, arachidonic acid > alpha-linolenic acid > eicosapentaenoic acid > linoleic acid. Of all the bioactive lipids studied, resolvin E1, neuroprotectin D1 and hydroxy-linoleic acid showed significant affinity comparable to proven PPARγ agonist namely, rosiglitazone, in terms of Glide XP docking score, H-bond formation with the key residues, binding free energy and stable complex formation with LBD favouring co-activator binding, as inferred through Molecular Dynamics trajectory analysis.

Conclusion

Hence, these three bioactive lipids (resolvin E1, neuroprotectin D1 and hydroxy-linoleic acid) may be favourably considered as ideal drug candidates in therapeutic modulation of clinical conditions such as type 2 DM, Alzheimer’s disease and other instances where PPARγ is a key player.

Markov state modeling and dynamical coarse-graining via discrete relaxation path sampling

A method is derived to coarse-grain the dynamics of complex molecular systems to a Markov jump process (MJP) describing how the system jumps between cells that fully partition its state space. The main inputs are relaxation times for each pair of cells, which are shown to be robust with respect to positioning of the cell boundaries. These relaxation times can be calculated via molecular dynamics simulations performed in each cell separately and are used in an efficient estimator for the rate matrix of the MJP. The method is illustrated through applications to Sinai billiards and a cluster of Lennard-Jones discs.

Ligand-based pharmacophore modelling and screening of DNA minor groove binders targeting Staphylococcus aureus

The recognition of DNA by small molecules is of special importance in the design of new drugs. Many natural and synthetic compounds have the ability to interact with the minor groove of DNA. In the present study, identification of minor groove binding compounds was attained by the combined approach of pharmacophore modelling, virtual screening and molecular dynamics approach. Experimentally reported 32 minor groove binding compounds were used to develop the pharmacophore model. Based on the fitness score, best three pharmacophore hypotheses were selected and used as template for screening the compounds from drug bank database. This pharmacophore-based screening provides many compounds with the same pharmacological properties. All these compounds were subjected to four phases of docking protocols with combined Glide-quantum-polarized ligand docking approach. Molecular dynamics results indicated that selected compounds are more active and showed good interaction in the binding site of DNA. Based on the scoring parameters and energy values, the best compounds were selected, and antibacterial activity of these compounds was identified using in vitro antimicrobial techniques.

Validation of potential inhibitors for SrtA against Bacillus anthracis by combined approach of ligand-based and molecular dynamics simulation

The development of SrtA inhibitors targeting the biothreat organism namely Bacillus anthracis was achieved by the combined approach of pharmacophore modeling, binding interactions, electron transferring capacity, ADME, and Molecular dynamics studies. In this study, experimentally reported Ba-SrtA inhibitors (pyridazinone and pyrazolethione derivatives) were considered for the development of enhanced pharmacophoric model. The obtained AAAHR hypothesis was a pure theoretical concept that accounts for common molecular interaction network present in experimentally active pyridazinone and pyrazolethione derivatives. Pharmacophore-based screening of AAAHR hypothesis provides several new compounds, and those compounds were treated with four phases of docking protocols with combined Glide-QPLD docking approach. In this approach, scoring and charge accuracy variations were seen to be dominated by QM/MM approach through the allocation of partial charges. Finally, we reported the best compounds from binding db, Chembridge db, and Toslab based on scoring values, energy parameters, electron transfer reaction, ADME, and cell adhesion inhibition activity. The dynamic state of interaction and binding energy assess that new compounds are more active inside the binding pocket and these compounds on experimental validations will survive as better inhibitors for targeting the cell adhesion mechanism of Ba-SrtA.

Simulation of the opening and closing of Hsp70 chaperones by coarse-grained molecular dynamics

Heat-shock proteins 70 (Hsp70s) are key molecular chaperones which assist in the folding and refolding/disaggregation of proteins. Hsp70s, which consist of a nucleotide-binding domain (NBD, consisting of NBD-I and NBD-II subdomains) and a substrate-binding domain [SBD, further split into the β-sheet (SBD-β) and α-helical (SBD-α) subdomains], occur in two major conformations having (a) a closed SBD, in which the SBD and NBD domains do not interact, (b) an open SBD, in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II. In the SBD-closed conformation, SBD is bound to a substrate protein, with release occurring after transition to the open conformation. While the transition from the closed to the open conformation is triggered efficiently by binding of adenosine triphosphate (ATP) to the NBD, it also occurs, although less frequently, in the absence of ATP. The reverse transition occurs after ATP hydrolysis. Here, we report canonical and multiplexed replica exchange simulations of the conformational dynamics of Hsp70s using a coarse-grained molecular dynamics approach with the UNRES force field. The simulations were run in the following three modes: (i) with the two halves of the NBD unrestrained relative to each other, (ii) with the two halves of the NBD restrained in an "open" geometry as in the SBD-closed form of DnaK (2KHO), and (iii) the two halves of NBD restrained in a "closed" geometry as in known experimental structures of ATP-bound NBD forms of Hsp70. Open conformations, in which the SBD interacted strongly with the NBD, formed spontaneously during all simulations; the number of transitions was largest in simulations carried out with the "closed" NBD domain, and smallest in those carried out with the "open" NBD domain; this observation is in agreement with the experimentally-observed influence of ATP-binding on the transition of Hsp70's from the SBD-closed to the SBD-open form. Two kinds of open conformations were observed: one in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II (as observed in the structures of nucleotide exchange factors), and another one in which this interaction pattern is swapped. A third type of motion, in which SBD-α binds to NBD without dissociating from SBD-β was also observed. It was found that the first stage of interdomain communication (approach of SBD-β, to NBD) is coupled with the rotation of the long axes of NBD-I and NBD-II towards each other. To the best of our knowledge, this is the first successful simulation of the full transition of an Hsp70 from the SBD-closed to the SBD-open conformation.

Coarse-grained force field: general folding theory

We review the coarse-grained UNited RESidue (UNRES) force field for the simulations of protein structure and dynamics, which is being developed in our laboratory over the last several years. UNRES is a physics-based force field, the prototype of which is defined as a potential of mean force of polypeptide chains in water, where all the degrees of freedom except the coordinates of α-carbon atoms and side-chain centers have been integrated out. We describe the initial implementation of UNRES to protein-structure prediction formulated as a search for the global minimum of the potential-energy function and its subsequent molecular dynamics and extensions of molecular-dynamics implementation, which enabled us to study protein-folding pathways and thermodynamics, as well as to reformulate the protein-structure prediction problem as a search for the conformational ensemble with the lowest free energy at temperatures below the folding-transition temperature. Applications of UNRES to study biological problems are also described.

Investigation of protein folding by coarse-grained molecular dynamics with the UNRES force field

Coarse-grained molecular dynamics simulations offer a dramatic extension of the time-scale of simulations compared to all-atom approaches. In this article, we describe the use of the physics-based united-residue (UNRES) force field, developed in our laboratory, in protein-structure simulations. We demonstrate that this force field offers about a 4000-times extension of the simulation time scale; this feature arises both from averaging out the fast-moving degrees of freedom and reduction of the cost of energy and force calculations compared to all-atom approaches with explicit solvent. With massively parallel computers, microsecond folding simulation times of proteins containing about 1000 residues can be obtained in days. A straightforward application of canonical UNRES/MD simulations, demonstrated with the example of the N-terminal part of the B-domain of staphylococcal protein A (PDB code: 1BDD, a three-alpha-helix bundle), discerns the folding mechanism and determines kinetic parameters by parallel simulations of several hundred or more trajectories. Use of generalized-ensemble techniques, of which the multiplexed replica exchange method proved to be the most effective, enables us to compute thermodynamics of folding and carry out fully physics-based prediction of protein structure, in which the predicted structure is determined as a mean over the most populated ensemble below the folding-transition temperature. By using principal component analysis of the UNRES folding trajectories of the formin-binding protein WW domain (PDB code: 1E0L; a three-stranded antiparallel beta-sheet) and 1BDD, we identified representative structures along the folding pathways and demonstrated that only a few (low-indexed) principal components can capture the main structural features of a protein-folding trajectory; the potentials of mean force calculated along these essential modes exhibit multiple minima, as opposed to those along the remaining modes that are unimodal. In addition, a comparison between the structures that are representative of the minima in the free-energy profile along the essential collective coordinates of protein folding (computed by principal component analysis) and the free-energy profile projected along the virtual-bond dihedral angles gamma of the backbone revealed the key residues involved in the transitions between the different basins of the folding free-energy profile, in agreement with existing experimental data for 1E0L .