NMR and MD Analysis of the Bonding Interaction of Vancomycin with Muramyl Pentapeptide.

The article describes an NMR spectroscopy study of interactions between vancomycin and a muramyl pentapeptide in two complexes: vancomycin and a native muramyl pentapeptide ended with D-alanine (MPP-D-Ala), and vancomycin and a modified muramyl pentapeptide ended with D-serine (MPP-D-Ser). The measurements were made in a 9:1 mixture of H2O and D2O. The obtained results confirmed the presence of hydrogen bonds previously described in the literature. At the same time, thanks to the pentapeptide model used, we were able to prove the presence of two more hydrogen bonds formed by the side chain amino group of L-lysine and oxygen atoms from the vancomycin carboxyl and amide groups. This type of interaction has not been described before. The existence of these hydrogen bonds was confirmed by the 1H NMR and molecular modeling. The formation of these bonds incurs additional through-space interactions, visible in the NOESY spectrum, between the protons of the L-lysine amino group and a vancomycin-facing hydrogen atom in the benzylic position. The presence of such interactions was also confirmed by molecular dynamics trajectory analysis.

Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions.

In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the ß and the α+ß structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES.

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment.

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

The Quaternization Reaction of 5-O-Sulfonates of Methyl 2,3-o-Isopropylidene-ß-D-Ribofuranoside With Selected Heterocyclic and Aliphatic Amines.

The synthesis of N-((methyl 5-deoxy-2,3-O-isopropylidene-ß-D-ribofuranoside)-5-yl)ammonium salts are presented. To determine the effect of the nucleophile type and outgoing group on the quaternization reaction, selected aliphatic and heterocyclic aromatic amines reacted with: methyl 2,3-O-isopropylidene-5-O-tosyl-ß-D-ribofuranoside or methyl 2,3-O-isopropylidene-5-O-mesyl-ß-D-ribofuranoside or methyl 2,3-O-isopropylidene-5-O-triflyl-ß-D-ribofuranoside were performed on a micro scale. High-resolution 1H- and 13C-NMR spectral data for all new compounds were recorded. Additionally, the single-crystal X-ray diffraction analysis for methyl 2,3-O-isopropylidene-5-O-mesyl-ß-D-ribofuranoside and selected in silico interaction models are reported.

Improved Consensus-Fragment Selection in Template-Assisted Prediction of Protein Structures with the UNRES Force Field in CASP13.

The method for protein-structure prediction, which combines the physics-based coarse-grained UNRES force field with knowledge-based modeling, has been developed further and tested in the 13th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP13). The method implements restraints from the consensus fragments common to server models. In this work, the server models to derive fragments have been chosen on the basis of quality assessment; a fully automatic fragment-selection procedure has been introduced, and Dynamic Fragment Assembly pseudopotentials have been fully implemented. The Global Distance Test Score (GDT_TS), averaged over our "Model 1" predictions, increased by over 10 units with respect to CASP12 for the free-modeling category to reach 40.82. Our "Model 1" predictions ranked 20 and 14 for all and free-modeling targets, respectively (upper 20.2% and 14.3% of all models submitted to CASP13 in these categories, respectively), compared to 27 (upper 21.1%) and 24 (upper 18.9%) in CASP12, respectively. For oligomeric targets, the Interface Patch Similarity (IPS) and Interface Contact Similarity (ICS) averaged over our best oligomer models increased from 0.28 to 0.36 and from 12.4 to 17.8, respectively, from CASP12 to CASP13, and top-ranking models of 2 targets (H0968 and T0997o) were obtained (none in CASP12). The improvement of our method in CASP13 over CASP12 was ascribed to the combined effect of the overall enhancement of server-model quality, our success in selecting server models and fragments to derive restraints, and improvements of the restraint and potential-energy functions.

Extension of the UNRES Coarse-Grained Force Field to Membrane Proteins in the Lipid Bilayer.

The physics-based UNRES coarse-grained force field for the simulations of protein structure and dynamics has been extended to treat membrane proteins. The lipid bilayer has been modeled by introducing a continuous nonpolar phase with the water-interface region of appropriate thickness. The potentials for average electrostatic and correlation interactions of the peptide groups have been rescaled to account for the reduction of the dielectric permittivity compared to the water phase and new potentials for protein side-chain-side-chain interactions inside and across the lipid phase have been introduced. The model was implemented in the UNRES package for coarse-grained simulations of proteins, and the package with the new functionality was tested for total energy conservation and thermostat behavior in microcanonical and canonical molecular dynamics simulations runs, respectively. The method was validated by running unrestricted ab initio blind-prediction tests of 10 short α-helical membrane proteins, all runs started from the extended structures. The modified UNRES force field was able to predict correctly the overall folds of the membrane proteins studied.

Evaluation of the scale-consistent UNRES force field in template-free prediction of protein structures in the CASP13 experiment.

The recent NEWCT-9P version of the coarse-grained UNRES force field for proteins, with scale-consistent formulas for the local and correlation terms, has been tested in the CASP13 experiment of the blind-prediction of protein structure, in the ab initio, contact-assisted, and data-assisted modes. Significant improvement of the performance has been observed with respect to the CASP11 and CASP12 experiments (by over 10 GDT_TS units for the ab initio mode predictions and by over 15 GDT_TS units for the contact-assisted prediction, respectively), which is a result of introducing scale-consistent terms and improved handling of contact-distance restraints. As in previous CASP exercises, UNRES ranked higher in the free modeling category than in the general category that included template based modeling targets. Use of distance restraints from the predicted contacts, albeit many of them were wrong, resulted in the increase of GDT_TS by over 8 units on average and introducing sparse restraints from small-angle X-ray/neutron scattering and chemical cross-link-mass-spectrometry experiments, and ambiguous restraints from nuclear magnetic resonance experiments has also improved the predictions by 8.6, 9.7, and 10.7 GDT_TS units on average, respectively.

An analysis and evaluation of the WeFold collaborative for protein structure prediction and its pipelines in CASP11 and CASP12.

Every two years groups worldwide participate in the Critical Assessment of Protein Structure Prediction (CASP) experiment to blindly test the strengths and weaknesses of their computational methods. CASP has significantly advanced the field but many hurdles still remain, which may require new ideas and collaborations. In 2012 a web-based effort called WeFold, was initiated to promote collaboration within the CASP community and attract researchers from other fields to contribute new ideas to CASP. Members of the WeFold coopetition (cooperation and competition) participated in CASP as individual teams, but also shared components of their methods to create hybrid pipelines and actively contributed to this effort. We assert that the scale and diversity of integrative prediction pipelines could not have been achieved by any individual lab or even by any collaboration among a few partners. The models contributed by the participating groups and generated by the pipelines are publicly available at the WeFold website providing a wealth of data that remains to be tapped. Here, we analyze the results of the 2014 and 2016 pipelines showing improvements according to the CASP assessment as well as areas that require further adjustments and research.

Use of the UNRES force field in template-assisted prediction of protein structures and the refinement of server models: Test with CASP12 targets.

Knowledge-based methods are, at present, the most effective ones for the prediction of protein structures; however, their results heavily depend on the similarity of a target sequence to those of proteins with known structures. On the other hand, the physics-based methods, although still less accurate and more expensive to execute, are independent of databases and give reasonable results where the knowledge-based methods fail because of weak sequence similarity. Therefore, a plausible approach seems to be the use of knowledge-based methods to determine the sections of the structures that correspond to sufficient sequence similarity and physics-based methods to determine the remaining structure. By participating in the 12th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP12) as the KIAS-Gdansk group, we tested our recently developed hybrid approach, in which protein-structure prediction is carried out by using the physics-based UNRES coarse-grained energy function, with restraints derived from the server models. Best predictions among all groups were obtained for 2 targets and 80% of our models were in the upper 50% of the models submitted to CASP. Our method was also able to exclude, with about 70% confidence, the information from the servers that performed poorly on a given target. Moreover, the method resulted in the best models of 2 refinement targets and performed remarkably well on oligomeric targets.

Performance of protein-structure predictions with the physics-based UNRES force field in CASP11.

Participating as the Cornell-Gdansk group, we have used our physics-based coarse-grained UNited RESidue (UNRES) force field to predict protein structure in the 11th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP11). Our methodology involved extensive multiplexed replica exchange simulations of the target proteins with a recently improved UNRES force field to provide better reproductions of the local structures of polypeptide chains. All simulations were started from fully extended polypeptide chains, and no external information was included in the simulation process except for weak restraints on secondary structure to enable us to finish each prediction within the allowed 3-week time window. Because of simplified UNRES representation of polypeptide chains, use of enhanced sampling methods, code optimization and parallelization and sufficient computational resources, we were able to treat, for the first time, all 55 human prediction targets with sizes from 44 to 595 amino acid residues, the average size being 251 residues. Complete structures of six single-domain proteins were predicted accurately, with the highest accuracy being attained for the T0769, for which the CαRMSD was 3.8 Å for 97 residues of the experimental structure. Correct structures were also predicted for 13 domains of multi-domain proteins with accuracy comparable to that of the best template-based modeling methods. With further improvements of the UNRES force field that are now underway, our physics-based coarse-grained approach to protein-structure prediction will eventually reach global prediction capacity and, consequently, reliability in simulating protein structure and dynamics that are important in biochemical processes. AVAILABILITY AND IMPLEMENTATION: Freely available on the web at http://www.unres.pl/ CONTACT: has5@cornell.edu.

A unified coarse-grained model of biological macromolecules based on mean-field multipole-multipole interactions.

A unified coarse-grained model of three major classes of biological molecules--proteins, nucleic acids, and polysaccharides--has been developed. It is based on the observations that the repeated units of biopolymers (peptide groups, nucleic acid bases, sugar rings) are highly polar and their charge distributions can be represented crudely as point multipoles. The model is an extension of the united residue (UNRES) coarse-grained model of proteins developed previously in our laboratory. The respective force fields are defined as the potentials of mean force of biomacromolecules immersed in water, where all degrees of freedom not considered in the model have been averaged out. Reducing the representation to one center per polar interaction site leads to the representation of average site-site interactions as mean-field dipole-dipole interactions. Further expansion of the potentials of mean force of biopolymer chains into Kubo's cluster-cumulant series leads to the appearance of mean-field dipole-dipole interactions, averaged in the context of local interactions within a biopolymer unit. These mean-field interactions account for the formation of regular structures encountered in biomacromolecules, e.g., α-helices and ß-sheets in proteins, double helices in nucleic acids, and helicoidally packed structures in polysaccharides, which enables us to use a greatly reduced number of interacting sites without sacrificing the ability to reproduce the correct architecture. This reduction results in an extension of the simulation timescale by more than four orders of magnitude compared to the all-atom representation. Examples of the performance of the model are presented.

Molecular modeling of Gram-positive bacteria peptidoglycan layer, selected glycopeptide antibiotics and vancomycin derivatives modified with sugar moieties.

Proper understanding of the mechanisms of binding to Gram-positive bacteria cell wall layers-especially to the peptidoglycan (PG) layer, seems to be crucial for proper development of new drug candidates which are effective against these bacteria. In this work we have constructed two different models of the Gram-positive bacteria PG layer: the layered and the scaffold models. PG conformational changes during geometry optimization, models relaxation, and molecular dynamics were described and discussed. We have found that the border surface of both PG layer models differs from the surface located away from the edge of models and the chains formed by disaccharide units prefer helix-like conformation. This curling of PG chains significantly affects the shape of antibiotic-accessible surface and the process is thus crucial for new drug development. Glycopeptide antibiotics effective against Gram-positive bacteria, such as vancomycin and its semisynthetic derivatives-oritavancin and telavancin, bind to d-alanyl-d-alanine stem termini on the peptidoglycan precursors of the cell wall. This binding inhibits cross-linking between the peptides and subsequently prevents cell wall synthesis. In this study some of the aspects of conformational freedom of vancomycin and restrictions from the modifications of vancomycin structure introduced into oritavancin and telavancin and five other vancomycin derivatives (with addition of 2-acetamido-2-deoxy-ß-d-galactopyranosylamine, 2-acetamido-2-deoxy-ß-d-glucopyranosylamine, 1-amine-1-deoxy-d-glucitol, 2-amino-2-deoxy-d-galactitol, or 2-amino-2-deoxy-d-glucitol to the C-terminal amino acid group in the vancomycin) are presented and discussed. The resulting molecular dynamics trajectories, root mean square deviation changes of aglycon and saccharide moieties as well as a comparative study of possible interactions with cyclic and chain forms of modified groups have been carried out, measured, and analyzed. Energetically advantageous conformations show close similarity to the structures known from the experimental data, but the diversity of others suggest very high conformational freedom of all modeled antibiotics and vancomycin derivatives. Alditol derivatives move closer to the peptidoglycan chain more easily but they also form intramolecular interactions more frequently than their homologous cyclic forms. One of the proposed derivatives seems to be a promising agent which is efficient in treatment of infections caused by Gram-positive bacteria.

Lessons from application of the UNRES force field to predictions of structures of CASP10 targets.

The performance of the physics-based protocol, whose main component is the United Residue (UNRES) physics-based coarse-grained force field, developed in our laboratory for the prediction of protein structure from amino acid sequence, is illustrated. Candidate models are selected, based on probabilities of the conformational families determined by multiplexed replica-exchange simulations, from the 10th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP10). For target T0663, classified as a new fold, which consists of two domains homologous to those of known proteins, UNRES predicted the correct symmetry of packing, in which the domains are rotated with respect to each other by 180° in the experimental structure. By contrast, models obtained by knowledge-based methods, in which each domain is modeled very accurately but not rotated, resulted in incorrect packing. Two UNRES models of this target were featured by the assessors. Correct domain packing was also predicted by UNRES for the homologous target T0644, which has a similar structure to that of T0663, except that the two domains are not rotated. Predictions for two other targets, T0668 and T0684_D2, are among the best ones by global distance test score. These results suggest that our physics-based method has substantial predictive power. In particular, it has the ability to predict domain-domain orientations, which is a significant advance in the state of the art.

Interactions of vasopressin and oxytocin receptors with vasopressin analogues substituted in position 2 with 3,3'-diphenylalanine--a molecular docking study.

Vasopressin and oxytocin receptors belong to the superfamily of G protein-coupled receptors and play an important role in many physiological functions. They are also involved in a number of pathological conditions being important drug targets. In this work, four vasopressin analogues substituted at position 2 with 3,3'-diphenylalanine have been docked into partially flexible vasopressin and oxytocin receptors. The bulky residue at position 2 acts as a structural restraint much stronger in the oxytocin receptor (OTR) than in the vasopressin V2 receptor (V2R), resulting in a different location of the analogues in these receptors. This explains the different, either agonistic or antagonistic, activities of the analogues in V2R and OTR, respectively. In all complexes, the conserved polar residues serve as anchor points for the ligand both in OTR and V2R. Strong interactions of the C-terminus of analogue II ([Mpa(1) ,d-Dpa(2) ,Val(4) ,d-Arg(8) ]VP) with extracellular loop 3 may be responsible for its highest activity at V2R. It also appears that V2R adapts more readily to the docking analogues by conformational changes in the aromatic side chains triggering receptor activation. A weak activity at V1a vasopressin receptor appears to be caused by weak receptor-ligand interactions. Results of this study may facilitate a rational design of new analogues with the highest activity/selectivity at vasopressin and OTRs.

Influence of bulky 3,3'-diphenylalanine enantiomers replacing position 2 of AVP analogues on their conformations: NMR and molecular modeling studies.

In this paper we use NMR spectroscopy and molecular modeling to examine four vasopressin analogues substituted with bulky 3,3'-diphenylalanine (Dpa) enantiomers: [Mpa(1),Dpa(2),Val(4),D-Arg(8)]VP (I), [Mpa(1),D-Dpa(2),Val(4),D-Arg(8)]VP (II), [D-Dpa(2),D-Arg(8)]VP (III) and [Mpa(1),D-Dpa(2)]AVP (IV). All the peptides exhibit a strong and prolonged antidiuretic activity. Additionally, analogues II, III and IV display antiuterotonic activity and analogue II is also a weak V(1a) receptor blocker. The conformational analysis has shown that beta-turns at positions 2,3 and/or 3,4 are characteristic of OT antagonists. In turn, the beta-turn in the Cys(6)-Gly(9) fragment seems to be crucial for enhancement of the antidiuretic activity. The high accessibility of aromatic side chains at positions 2 and 3 plays a crucial role in antagonist-receptor binding. Moreover, orientation of the Phe(3) side chain is claimed to be important for V(1a) receptor affinity.

Conformational stability of the full-atom hexameric model of the ClpB chaperone from Escherichia coli.

The Escherichia coli heat shock protein ClpB, a member of the Hsp100 family, plays a crucial role in cellular thermotolerance. In co-operation with the Hsp70 chaperone system, it is able to solubilize proteins aggregated by heat shock conditions and refold them into the native state in an ATP-dependent way. It was established that the mechanism of ClpB action depends on the formation of a ring-shaped hexameric structure and the translocation of a protein substrate through an axial channel. The structural aspects of this process are not fully known. By means of homology modeling and protein-protein docking, we obtained a model of the hexameric arrangement of the full-length ClpB protein complexed with ATP. A molecular dynamics simulation of this model was performed to assess its flexibility and conformational stability. The high mobility of the "linker" M-domain, essential for the renaturing activity of ClpB, was demonstrated, and the size and shape of central channel were analyzed. In this model, we propose the coordinates for a loop between b4 and B6 structural elements, not defined in previous structural research, which faces the inside of the channel and may therefore play a role in substrate translocation.

Molecular dynamics study of the internal water molecules in vasopressin and oxytocin receptors.

The role of the internal water molecules in vasopressin and oxytocin receptors has been investigated via molecular dynamics simulations in hydrated membrane model. Several water molecules have been identified within the binding pockets of receptors, where they interact with the conserved residues. In all unliganded receptors, the water molecules bound to the highly conserved D2.50 cluster have been observed. It has been proposed which water molecules may significantly contribute to the stability of overall receptor structure. In receptor-ligand complexes the water molecules are involved in the receptor-ligand interactions by forming water-mediated hydrogen bonds at their contact surface.

Molecular docking-based study of vasopressin analogues modified at positions 2 and 3 with N-methylphenylalanine: influence on receptor-bound conformations and interactions with vasopressin and oxytocin receptors.

In this study, four cyclic vasopressin (CYFQNCPRG-NH(2), AVP) analogues substituted at positions 2 and 3 with four combinations of enantiomers of N-methylphenylalanine have been investigated. Three-dimensional structures of analogues have been formerly determined using NMR spectroscopy in dimethyl sulfoxide. Three-dimensional models of the vasopressin and oxytocin receptors were constructed by combining the multiple sequence alignment and the RD crystal structure as a template. The analogues have been docked into the receptor using the AutoDock program. The relaxation of the receptor-ligand complexes using energy minimization, followed by the constrained simulated annealing protocols (CSA), has been performed. The receptor-bound conformations of the investigated analogues have been proposed. We concluded that the N-methylated residues at positions 2 and 3 act as a structural restraint, determining the conformation of analogues, their location inside the receptor cavity, and mutual arrangement of the aromatic side chains. The conserved polar residues constitute the handles keeping the biologically active analogues inside the binding cavity. The Arg(8)-D(2.65) salt bridge might be responsible for analogue-selective binding in OTR and V1aR versus V2R, where the positively charged K(2.65) 100 is present at the equivalent position.

Investigation of cis/trans ratios of peptide bonds in AVP analogues containing N-methylphenylalanine enantiomers.

The solution conformation of vasopressin analogues, modified at positions 2 and 3 with N-methylphenylalanine or its enantiomer, [D-MePhe2,MePhe3]AVP and [MePhe2,D-MePhe3]AVP, were studied by 2D NMR spectroscopy in H2O/D2O and theoretical calculations (EDMC/ANALYZE). In the case of [MePhe2,D-MePhe3]AVP, the synthesis afforded two products, A and B, which had identical molecular ions and similar retention times on HPLC. This finding was explained by racemization of Cys1, which gave an additional analogue, [D-Cys1,MePhe2,D-MePhe3]AVP (B). The possibility is not excluded of racemization of Cys1 in the remaining analogues of this series. However, only in the case of [MePhe2,D-MePhe3]AVP was this process so distinct that two strong peaks in the HPLC chromatogram were observed. The NMR spectra of all the analogues showed several distinct sets of residual proton resonances. This suggests that the peptides adopt more than two groups of conformations in H2O/D2O. This fact is due to cis/trans isomerization. Two more populated isomers arise from the cis/trans isomerization across the 2-3 peptide bond in [D-MePhe2,MePhe3]AVP and [MePhe2,D-MePhe3]AVP (A) and across the 1-2 peptide bond in [D-Cys1,MePhe2,D-MePhe3]AVP (B).