Integrative modeling of diverse protein-peptide systems using CABS-dock.

The CABS model can be applied to a wide range of protein-protein and protein-peptide molecular modeling tasks, such as simulating folding pathways, predicting structures, docking, and analyzing the structural dynamics of molecular complexes. In this work, we use the CABS-dock tool in two diverse modeling tasks: 1) predicting the structures of amyloid protofilaments and 2) identifying cleavage sites in the peptide substrates of proteolytic enzymes. In the first case, simulations of the simultaneous docking of amyloidogenic peptides indicated that the CABS model can accurately predict the structures of amyloid protofilaments which have an in-register parallel architecture. Scoring based on a combination of symmetry criteria and estimated interaction energy values for bound monomers enables the identification of protofilament models that closely match their experimental structures for 5 out of 6 analyzed systems. For the second task, it has been shown that CABS-dock coarse-grained docking simulations can be used to identify the positions of cleavage sites in the peptide substrates of proteolytic enzymes. The cleavage site position was correctly identified for 12 out of 15 analyzed peptides. When combined with sequence-based methods, these docking simulations may lead to an efficient way of predicting cleavage sites in degraded proteins. The method also provides the atomic structures of enzyme-substrate complexes, which can give insights into enzyme-substrate interactions that are crucial for the design of new potent inhibitors.

Docking of peptides to GPCRs using a combination of CABS-dock with FlexPepDock refinement.

The structural description of peptide ligands bound to G protein-coupled receptors (GPCRs) is important for the discovery of new drugs and deeper understanding of the molecular mechanisms of life. Here we describe a three-stage protocol for the molecular docking of peptides to GPCRs using a set of different programs: (1) CABS-dock for docking fully flexible peptides; (2) PD2 method for the reconstruction of atomistic structures from C-alpha traces provided by CABS-dock and (3) Rosetta FlexPepDock for the refinement of protein-peptide complex structures and model scoring. We evaluated the proposed protocol on the set of seven different GPCR-peptide complexes (including one containing a cyclic peptide), for which crystallographic structures are available. We show that CABS-dock produces high resolution models in the sets of top-scored models. These sets of models, after reconstruction to all-atom representation, can be further improved by Rosetta high-resolution refinement and/or minimization, leading in most of the cases to sub-Angstrom accuracy in terms of interface root-mean-square-deviation measure.

The Effect of Conjugation of Ciprofloxacin and Moxifloxacin with Fatty Acids on Their Antibacterial and Anticancer Activity.

Novel conjugates (CP) of moxifloxacin (MXF) with fatty acids (1m-16m) were synthesized with good yields utilizing amides chemistry. They exhibit a more pronounced cytotoxic potential than the parent drug. They were the most effective for prostate cancer cells with an IC50 below 5 µM for respective conjugates with sorbic (2m), oleic (4m), 6-heptenoic (10m), linoleic (11m), caprylic (15m), and stearic (16m) acids. All derivatives were evaluated against a panel of standard and clinical bacterial strains, as well as towards mycobacteria. The highest activity towards standard isolates was observed for the acetic acid derivative 14m, followed by conjugates of unsaturated crotonic (1m) and sorbic (2m) acids. The activity of conjugates tested against an expanded panel of clinical coagulase-negative staphylococci showed that the compound (14m) was recognized as a leading structure with an MIC of 0.5 µg/mL denoted for all quinolone-susceptible isolates. In the group of CP derivatives, sorbic (2) and geranic (3) acid amides exhibited the highest bactericidal potential against clinical strains. The M. tuberculosis Spec. 210 strain was the most sensitive to sorbic (2m) conjugate and to conjugates with medium- and long-chain polyunsaturated acids. To establish the mechanism of antibacterial action, selected CP and MXF conjugates were examined in both topoisomerase IV decatenation assay and the DNA gyrase supercoiling assay, followed by suitable molecular docking studies.

Design and Synthesis of Menthol and Thymol Derived Ciprofloxacin: Influence of Structural Modifications on the Antibacterial Activity and Anticancer Properties.

Sixteen new Ciprofloxacin derivatives were designed and successfully synthesized. In an in silico experiment, lipophilicity was established for obtained compounds. All compounds were screened for antimicrobial activity using standard and clinical strains. As for Gram-positive hospital microorganisms, all tested derivatives were active. Measured MICs were in the range 1-16 µg/mL, confirming high antimicrobial potency. Derivative 12 demonstrated activity against all standard Gram-positive Staphylococci, within the range of 0.8-1.6 µg/mL and was confirmed as the leading structure with MICs 1 µg/mL for S. pasteuri KR 4358 and S. aureus T 5591 (clinical strains). All compounds were screened for their in vitro cytotoxic properties via the MTT method. Three of the examined compounds (3, 11 and 16) showed good activity against cancer cells, and in parallel were found not to be cytotoxic toward normal cells. Doxorubicin SI ranged 0.14-1.11 while the mentioned three ranged 1.9-3.4. Selected Ciprofloxacin derivatives were docked into the crystal structure of topoisomerase II (DNA gyrase) in complex with DNA (PDB ID: 5BTC). In summary, leading structures were established (3, 11, 12 and 16). We have observed poor results in preformed studies for disubstituted derivatives, suggesting that 3-oxo-4-carboxylic acid core is the active DNA-gyrase binding site, and when structural changes were made in this fragment, there was an observed decrease in antibacterial potency.

Multiscale Modeling of Amyloid Fibrils Formed by Aggregating Peptides Derived from the Amyloidogenic Fragment of the A-Chain of Insulin.

Computational prediction of molecular structures of amyloid fibrils remains an exceedingly challenging task. In this work, we propose a multi-scale modeling procedure for the structure prediction of amyloid fibrils formed by the association of ACC1-13 aggregation-prone peptides derived from the N-terminal region of insulin's A-chain. First, a large number of protofilament models composed of five copies of interacting ACC1-13 peptides were predicted by application of CABS-dock coarse-grained (CG) docking simulations. Next, the models were reconstructed to all-atom (AA) representations and refined during molecular dynamics (MD) simulations in explicit solvent. The top-scored protofilament models, selected using symmetry criteria, were used for the assembly of long fibril structures. Finally, the amyloid fibril models resulting from the AA MD simulations were compared with atomic force microscopy (AFM) imaging experimental data. The obtained results indicate that the proposed multi-scale modeling procedure is capable of predicting protofilaments with high accuracy and may be applied for structure prediction and analysis of other amyloid fibrils.

Novel Tetrazole-Based Antimicrobial Agents Targeting Clinical Bacteria Strains: Exploring the Inhibition of Staphylococcus aureus DNA Topoisomerase IV and Gyrase.

Eleven novel imide-tetrazoles were synthesized. In the initial stage of research, in silico structure-based pharmacological prediction was conducted. All compounds were screened for antimicrobial activity using standard and clinical strains. Within the studied group, compounds 1-3 were recognized as leading structures with the most promising results in antimicrobial studies. Minimal inhibitory concentration values for compounds 1, 2, 3 were within the range of 0.8-3.2 µg/mL for standard and clinical Gram-positive and Gram-negative bacterial strains, showing in some cases higher activity than the reference Ciprofloxacin. Additionally, all three inhibited the growth of all clinical Staphylococci panels: Staphylococcus aureus (T5592; T5591) and Staphylococcus epidermidis (5253; 4243) with MIC values of 0.8 µg/mL. Selected compounds were examined in topoisomerase IV decatenation assay and DNA gyrase supercoiling assay, followed by suitable molecular docking studies to explore the possible binding modes. In summary, the presented transition from substrate imide-thioureas to imide-tetrazole derivatives resulted in significant increase of antimicrobial properties. The compounds 1-3 proposed here provide a promising basis for further exploration towards novel antimicrobial drug candidates.

Molecular Dynamics Scoring of Protein-Peptide Models Derived from Coarse-Grained Docking.

One of the major challenges in the computational prediction of protein-peptide complexes is the scoring of predicted models. Usually, it is very difficult to find the most accurate solutions out of the vast number of sometimes very different and potentially plausible predictions. In this work, we tested the protocol for Molecular Dynamics (MD)-based scoring of protein-peptide complex models obtained from coarse-grained (CG) docking simulations. In the first step of the scoring procedure, all models generated by CABS-dock were reconstructed starting from their original C-alpha trace representations to all-atom (AA) structures. The second step included geometry optimization of the reconstructed complexes followed by model scoring based on receptor-ligand interaction energy estimated from short MD simulations in explicit water. We used two well-known AA MD force fields, CHARMM and AMBER, and a CG MARTINI force field. Scoring results for 66 different protein-peptide complexes show that the proposed MD-based scoring approach can be used to identify protein-peptide models of high accuracy. The results also indicate that the scoring accuracy may be significantly affected by the quality of the reconstructed protein receptor structures.

Coarse-Grained Protein Models and Their Applications.

The traditional computational modeling of protein structure, dynamics, and interactions remains difficult for many protein systems. It is mostly due to the size of protein conformational spaces and required simulation time scales that are still too large to be studied in atomistic detail. Lowering the level of protein representation from all-atom to coarse-grained opens up new possibilities for studying protein systems. In this review we provide an overview of coarse-grained models focusing on their design, including choices of representation, models of energy functions, sampling of conformational space, and applications in the modeling of protein structure, dynamics, and interactions. A more detailed description is given for applications of coarse-grained models suitable for efficient combinations with all-atom simulations in multiscale modeling strategies.

Structural Basis of GD2 Ganglioside and Mimetic Peptide Recognition by 14G2a Antibody.

Monoclonal antibodies targeting GD2 ganglioside (GD2) have recently been approved for the treatment of high risk neuroblastoma and are extensively evaluated in clinics in other indications. This study illustrates how a therapeutic antibody distinguishes between different types of gangliosides present on normal and cancer cells and informs how synthetic peptides can imitate ganglioside in its binding to the antibody. Using high resolution crystal structures we demonstrate that the ganglioside recognition by a model antibody (14G2a) is based primarily on an extended network of direct and water molecule mediated hydrogen bonds. Comparison of the GD2-Fab structure with that of a ligand free antibody reveals an induced fit mechanism of ligand binding. These conclusions are validated by directed mutagenesis and allowed structure guided generation of antibody variant with improved affinity toward GD2. Contrary to the carbohydrate, both evaluated mimetic peptides utilize a "key and lock" interaction mechanism complementing the surface of the antibody binding groove exactly as found in the empty structure. The interaction of both peptides with the Fab relies considerably on hydrophobic contacts however, the detailed connections differ significantly between the peptides. As such, the evaluated peptide carbohydrate mimicry is defined primarily in a functional and not in structural manner.

Coarse-Grained Simulations of Membrane Insertion and Folding of Small Helical Proteins Using the CABS Model.

The CABS coarse-grained model is a well-established tool for modeling globular proteins (predicting their structure, dynamics, and interactions). Here we introduce an extension of the CABS representation and force field (CABS-membrane) to the modeling of the effect of the biological membrane environment on the structure of membrane proteins. We validate the CABS-membrane model in folding simulations of 10 short helical membrane proteins not using any knowledge about their structure. The simulations start from random protein conformations placed outside the membrane environment and allow for full flexibility of the modeled proteins during their spontaneous insertion into the membrane. In the resulting trajectories, we have found models close to the experimental membrane structures. We also attempted to select the correctly folded models using simple filtering followed by structural clustering combined with reconstruction to the all-atom representation and all-atom scoring. The CABS-membrane model is a promising approach for further development toward modeling of large protein-membrane systems.

The mechanism of ligand-induced activation or inhibition of µ- and κ-opioid receptors.

G-protein-coupled receptors (GPCRs) are important targets for treating severe diseases. However why certain molecules act as activators whereas others, with similar structures, block GPCR activation, is poorly understood since the same molecule can activate one receptor subtype while blocking another closely related receptor. To shed light on these central questions, we used all-atom, long-time-scale molecular dynamics simulations on the κ-opioid and µ-opioid receptors (κOR and µOR). We found that water molecules penetrating into the receptor interior mediate the activating versus blocking effects of a particular ligand-receptor interaction. Both the size and the flexibility of the bound ligand regulated water influx into the receptor. The solvent-accessible inner surface area was found to be a parameter that can help predict the function of the bound ligand.

Structure prediction of the second extracellular loop in G-protein-coupled receptors.

G-protein-coupled receptors (GPCRs) play key roles in living organisms. Therefore, it is important to determine their functional structures. The second extracellular loop (ECL2) is a functionally important region of GPCRs, which poses significant challenge for computational structure prediction methods. In this work, we evaluated CABS, a well-established protein modeling tool for predicting ECL2 structure in 13 GPCRs. The ECL2s (with between 13 and 34 residues) are predicted in an environment of other extracellular loops being fully flexible and the transmembrane domain fixed in its x-ray conformation. The modeling procedure used theoretical predictions of ECL2 secondary structure and experimental constraints on disulfide bridges. Our approach yielded ensembles of low-energy conformers and the most populated conformers that contained models close to the available x-ray structures. The level of similarity between the predicted models and x-ray structures is comparable to that of other state-of-the-art computational methods. Our results extend other studies by including newly crystallized GPCRs.

Structural elucidation of specific noncovalent association of folic acid with native cyclodextrins using an ion mobility mass spectrometry and theoretical approach.

The combination of ion mobility mass spectrometry studies and theoretical calculations including docking studies permitted a detailed structural description of noncovalent complexes of folic acid (FA) and native cyclodextrins (α-CD, ß-CD, and γ-CD). The mode of noncovalent association depended on the cavity size of the cyclodextrin. The structure of FA/α-CD represented the exclusion complex in which the aminobenzoic moiety and the aromatic pteridine ring of folic acid remain outside the cyclodextrin cavity, while the glutamate residue is anchored in the interior of the α-cyclodextrin. A rotaxane-type structure was proposed for the FA/ß-CD complex with the aminobenzoic part of FA being trapped in the central cavity of ß-CD. The glutamate residue and the aromatic pteridine ring interact with the primary and secondary rim hydroxyl residues, respectively, enhancing complex stability. Two possible structures of FA/γ-CD were suggested, the first one being analogous to the FA/ß-CD complex and the second one being more stable-in which the aromatic pteridine ring penetrates into the CD cavity while the glutamate residue with the aminobenzoic part of FA is exposed to the cone exterior of CD at its wider edge. Further insight into the association behavior of the folic acid toward cyclodextrins evaluated by thermodynamic calculations indicates that the process is highly exothermic. The complex stability increased in the order FA/α-CD < FA/ß-CD < FA/γ-CD. This order is consistent with the previously determined relative gas-phase stability established based on the dissociation efficiency curves of the FA/CD complexes.

Structural investigation of the C-terminal catalytic fragment of presenilin 1.

The gamma-secretase complex has a decisive role in the development of Alzheimer's disease, in that it cleaves a precursor to create the amyloid beta peptide whose aggregates form the senile plaques encountered in the brains of patients. Gamma-secretase is a member of the intramembrane-cleaving proteases which process their transmembrane substrates within the bilayer. Many of the mutations encountered in early onset familial Alzheimer's disease are linked to presenilin 1, the catalytic component of gamma-secretase, whose active form requires its endoproteolytic cleavage into N-terminal and C-terminal fragments. Although there is general agreement regarding the topology of the N-terminal fragment, studies of the C-terminal fragment have yielded ambiguous and contradictory results that may be difficult to reconcile in the absence of structural information. Here we present the first structure of the C-terminal fragment of human presenilin 1, as obtained from NMR studies in SDS micelles. The structure reveals a topology where the membrane is likely traversed three times in accordance with the more generally accepted nine transmembrane domain model of presenilin 1, but contains unique structural features adapted to accommodate the unusual intramembrane catalysis. These include a putative half-membrane-spanning helix N-terminally harboring the catalytic aspartate, a severely kinked helical structure toward the C terminus as well as a soluble helix in the assumed-to-be unstructured N-terminal loop.

A comprehensive pharmacological analysis of fenoterol and its derivatives to unravel the role of ß2-adrenergic receptor in zebrafish.

ß-adrenergic receptors (ßARs) belong to a key molecular targets that regulate the most important processes occurring in the human organism. Although over the last decades a zebrafish model has been developed as a model complementary to rodents in biomedical research, the role of ß2AR in regulation of pathological and toxicological effects remains to elucidate. Therefore, the study aimed to clarify the role of ß2AR with a particular emphasis on the distinct role of subtypes A and B of zebrafish ß2AR. As model compounds selective ß2AR agonists - (R,R)-fenoterol ((R,R)-Fen) and its new derivatives: (R,R)-4'-methoxyfenoterol ((R,R)-MFen) and (R,R)-4'-methoxy-1-naphtylfenoterol ((R,R)-MNFen) - were tested. We described dose-dependent changes observed after fenoterols exposure in terms of general toxicity, cardiotoxicity and neurobehavioural responses. Subsequently, to better characterise the role of ß2-adrenergic stimulation in zebrafish, we have performed a series of molecular docking simulations. Our results indicate that (R,R)-Fen displays the highest affinity for subtype A of zebrafish ß2AR and ß2AAR might be involved in pigment depletion. (R,R)-MFen shows the lowest affinity for zebrafish ß2ARs out of the tested fenoterols and this might be associated with its cardiotoxic and anxiogenic effects. (R,R)-MNFen displays the highest affinity for subtype B of zebrafish ß2AR and modulation of this receptor might be associated with the development of malformations, increases locomotor activity and induces a negative chronotropic effect. Taken together, the presented data offer insights into the functional responses of the zebrafish ß2ARs confirming their intraspecies conservation, and support the translation of the zebrafish model in pharmacological and toxicological research.

Menthol- and thymol-based ciprofloxacin derivatives against Mycobacterium tuberculosis: in vitro activity, lipophilicity, and computational studies.

In this work, we investigated the antitubercular properties of Ciprofloxacin derivatives conjugated with menthol and thymol moieties. For the sixteen derivatives, we established minimal inhibitory concentrations (MIC) using isolates of Mycobacterium tuberculosis that were resistant or susceptible to other antibiotics. For the most potent compound 1-cyclopropyl-6-fluoro-7-{4-[6-((1R,2S,5R)-2-isopropyl-5-methylcyclohexyloxy)-6-oxohexyl]piperazin-1-yl}-4-oxo-1,4-dihydroquinoline-3-carboxylic acid (6), we determined fractional inhibitory concentration index (FICI) values to confirm antibacterial susceptibility and synergistic effects with other reference drugs. In addition, chromatographic studies of all the derivatives demonstrated a significant three to four-fold increase in lipophilicity and affinity to phospholipids compared to Ciprofloxacin. Finally, we conducted structure-based studies of the investigated compounds using molecular docking and taking into account protein target mutations associated with fluoroquinolone resistance. In summary, our findings indicate that the investigated compounds possess tuberculostatic properties, with some showing similar or even better activity against resistant strains compared to reference drugs. Increased lipophilicity and affinity to phospholipids of the new derivatives can offer several advantages for new drug candidates, beyond just improved cell membrane penetration. However, further studies are needed to fully understand their safety, efficacy, and mechanism of action.

N-Acylated Ciprofloxacin Derivatives: Synthesis and In Vitro Biological Evaluation as Antibacterial and Anticancer Agents.

A novel series of N-acylated ciprofloxacin (CP) conjugates 1-21 were synthesized and screened as potential antimicrobial agents. Conjugates 1 and 2 were 1.25-10-fold more potent than CP toward all Staphylococci (minimal inhibitory concentration 0.05-0.4 µg/mL). Most of the chloro- (3-7), bromo- (8-11), and CF3-alkanoyl (14-16) derivatives expressed higher or comparable activity to CP against selected Gram-positive strains. A few CP analogues (5, 10, and 11) were also more effective toward the chosen clinical Gram-negative rods. Conjugates 5, 10, and 11 considerably influenced the phases of the bacterial growth cycle over 18 h. Additionally, compounds 2, 4-7, 9-12, and 21 exerted stronger tuberculostatic action against three Mycobacterium tuberculosis isolates than the first-line antitubercular drugs. Amides 1, 2, 5, 6, 10, and 11 targeted gyrase and topoisomerase IV at 2.7-10.0 µg/mL, which suggests a mechanism of antibacterial action related to CP. These findings were confirmed by molecular docking studies. In addition, compounds 3 and 15 showed high antiproliferative activities against prostate PC3 cells (IC50 2.02-4.8 µM), up to 6.5-2.75 stronger than cisplatin. They almost completely reduced the growth and proliferation rates in these cells, without a cytotoxic action against normal HaCaT cell lines. Furthermore, derivatives 3 and 21 induced apoptosis/necrosis in PC3 cells, probably by increasing the intracellular ROS amount, as well as they diminished the IL-6 level in tumor cells.

Interactions of benzo[a]pyrene and diesel exhaust particulate matter with the lung surfactant system.

One of the reasons for cellular changes in the lung tissue exposed to the diesel exhaust composed of soot particles with adsorbed volatile organic molecules is the reduction of the clearance rate in the pulmonary region of the respiratory system. The interaction of the fractal-like particles and organic substances with a surfactant monolayer limits its dynamic activity. The surface properties of Survanta, a purified extract of bovine lung surfactant (LS), which interacted with carbon particles (200 nm aggregates) and benzo[a]pyrene (BaP), molecules were measured with the oscillating bubble technique. The results showed a significant lowering of the dynamics of the surfactant monolayer compared to the control case (no exposure). Additional measurements of surface pressure during the monotonic compression of the air-water interface containing the major LS phospholipid, dipalmitoylphosphatidylcholine (DPPC), showed that the presence of BaP molecules in the system influenced its stability. The experimental results were supplemented with a theoretical molecular dynamics model of the interaction between BaP and DPPC molecules. The simulation results indicated the insertion of BaP molecules into the lipid layer, which explained the measured effects.

Generation and characterization of a novel, permanently active S100P mutant.

S100 proteins function as Ca2+ signal transducers by regulating cellular targets in their Ca2+ bound conformation. S100P is a member of the S100 protein family that can activate the membrane and F-actin binding protein ezrin in a Ca2+ dependent manner at least in vitro. Here we generated a novel tool to elucidate directly the S100P-ezrin interaction in vivo. This was achieved by constructing a S100P derivative that contained mutations in the two EF hand loops predicted to lock the protein in a permanently active state. The resulting S100P mutant, termed here S100P pa, could be purified as a soluble protein and showed biochemical properties displayed by wild-type S100P only in the presence of Ca2+. Importantly, S100P pa bound to the N-terminal domain of ezrin in the absence of Ca2+ showing an affinity only slightly reduced as compared to that of Ca2+-bound WT S100P. In line with this permanent complex formation, S100P pa colocalized with ezrin to plasma membrane protrusions of epithelial cells even in the absence of intracellular Ca2+ transients. Thus, S100P pa is a novel type of S100 protein mutant locked in a permanently active state that shows an unregulated complex formation with its cellular target ezrin.

Rapid self-association of highly amyloidogenic H-fragments of insulin: Experiment and molecular dynamics simulations.

The so-called 'H-fragment' of insulin is an extremely amyloidogenic double chain peptide consisting of the N-terminal parts of A-chain and B-chain linked by a disulfide bond between Cys-7A and Cys-7B. Here, we conduct a detailed investigation of the self-association behavior of H-fragment monomers into amyloid-like fibrils using kinetic assays, infrared spectroscopy, circular dichroism (CD), atomic force microscopy (AFM) and molecular dynamics (MD) simulations. Unlike the intact predominantly α-helical insulin, H-fragment remains in a disordered state in aqueous solutions. Its aggregation accelerates with acidification of the environment leading, at pH 1.9, to the formation of thin and structurally homogenous fibrils with the infrared features typical for parallel ß-sheet conformation. According to time-lapse AFM morphological analysis both secondary nucleation and fragmentation are involved in later stages of H-fibrils' self-assembly. Based on the low nucleation order (two) obtained from the global fitting of kinetic data, realistic all-atom MD simulations of pairs of interacting H-fragment monomers were subsequently carried out. The molecular self-association scenario emerging from these simulations implicates the intrinsic conformational instability of H-monomer in its tendency to aggregate and form intermolecular ß-sheet structure. Our findings provide the new mechanistic context for studies of insulin misfolding and aggregation.