15N-Labelling and structure determination of adamantylated azolo-azines in solution.

Determining the accurate chemical structures of synthesized compounds is essential for biomedical studies and computer-assisted drug design. The unequivocal determination of N-adamantylation or N-arylation site(s) in nitrogen-rich heterocycles, characterized by a low density of hydrogen atoms, using NMR methods at natural isotopic abundance is difficult. In these compounds, the heterocyclic moiety is covalently attached to the carbon atom of the substituent group that has no bound hydrogen atoms, and the connection between the two moieties of the compound cannot always be established via conventional 1H-1H and 1H-13C NMR correlation experiments (COSY and HMBC, respectively) or nuclear Overhauser effect spectroscopy (NOESY or ROESY). The selective incorporation of 15N-labelled atoms in different positions of the heterocyclic core allowed for the use of 1H-15N (JHN) and 13C-15N (JCN) coupling constants for the structure determinations of N-alkylated nitrogen-containing heterocycles in solution. This method was tested on the N-adamantylated products in a series of azolo-1,2,4-triazines and 1,2,4-triazolo[1,5-a]pyrimidine. The syntheses of adamantylated azolo-azines were based on the interactions of azolo-azines and 1-adamatanol in TFA solution. For azolo-1,2,4-triazinones, the formation of mixtures of N-adamantyl derivatives was observed. The JHN and JCN values were measured using amplitude-modulated 1D 1H spin-echo experiments with the selective inversion of the 15N nuclei and line-shape analysis in the 1D 13С spectra acquired with selective 15N decoupling, respectively. Additional spin-spin interactions were detected in the 15N-HMBC spectra. NMR data and DFT (density functional theory) calculations permitted to suggest a possible mechanism of isomerization for the adamantylated products of the azolo-1,2,4-triazines. The combined analysis of the JHN and JCN couplings in 15N-labelled compounds provides an efficient method for the structure determination of N-alkylated azolo-azines even in the case of isomer formation. The isomerization of adamantylated tetrazolo[1,5-b][1,2,4]triazin-7-ones in acidic conditions occurs through the formation of the adamantyl cation.

Multiscale modelling of biomolecular corona formation on metallic surfaces.

In the realm of food industry, the choice of non-consumable materials used plays a crucial role in ensuring consumer safety and product quality. Aluminum is widely used in food packaging and food processing applications, including dairy products. However, the interaction between aluminum and milk content requires further investigation to understand its implications. In this work, we present the results of multiscale modelling of the interaction between various surfaces, that is (100), (110), and (111), of fcc aluminum with the most abundant milk proteins and lactose. Our approach combines atomistic molecular dynamics, a coarse-grained model of protein adsorption, and kinetic Monte Carlo simulations to predict the protein corona composition in the deposited milk layer on aluminum surfaces. We consider a simplified model of milk, which is composed of the six most abundant milk proteins found in natural cow milk and lactose, which is the most abundant sugar found in dairy. Through our study, we ranked selected proteins and lactose adsorption affinities based on their corresponding interaction strength with aluminum surfaces and predicted the content of the naturally forming biomolecular corona. Our comprehensive investigation sheds light on the implications of aluminum in food processing and packaging, particularly concerning its interaction with the most abundant milk proteins and lactose. By employing a multiscale modelling approach, we simulated the interaction between metallic aluminum surfaces and the proteins and lactose, considering different crystallographic orientations. The results of our study provide valuable insights into the mechanisms of lactose and protein deposition on aluminum surfaces, which can aid in the general understanding of protein corona formation.

In silico prediction of protein binding affinities onto core-shell PEGylated noble metal nanoparticles for rational design of drug nanocarriers.

Polymer-coated nanoparticles (NP) are commonly used as drug carriers or theranostic agents. Their uptake rates are modulated by the interactions with essential serum proteins such as transferrin and albumin. Understanding the control parameters of these interactions is crucial for improving the efficiency of these nanoscale devices. In this work, we perform a multiscale computational study of protein adsorption onto polyethylene glycol (PEG) coated gold and silver NPs, producing protein-NP adsorption rankings as a function of PEG grafting density, which are validated against previously reported experimental protein-NP binding constants. Furthermore, the applied nano-docking method provides information on the preferred orientation of proteins immobilised on the surface of NPs. We propose a method of construction of model core-shell NPs in silico. The presented protocol can provide molecular level insights for the experimental development of biosensors, nanocarriers, or other nanoplatforms where information on the preferred orientation of protein at the bio-nano interface is crucial, and enables fast in silico prescreening of assays of various nanocarriers, i.e., combinations of proteins, NPs, and coatings.

Milk Protein Adsorption on Metallic Iron Surfaces.

Food processing and consumption involves multiple contacts between biological fluids and solid materials of processing devices, of which steel is one of the most common. Due to the complexity of these interactions, it is difficult to identify the main control factors in the formation of undesirable deposits on the device surfaces that may affect safety and efficiency of the processes. Mechanistic understanding of biomolecule-metal interactions involving food proteins could improve management of these pertinent industrial processes and consumer safety in the food industry and beyond. In this work, we perform a multiscale study of the formation of protein corona on iron surfaces and nanoparticles in contact with cow milk proteins. By calculating the binding energies of proteins with the substrate, we quantify the adsorption strength and rank proteins by the adsorption affinity. We use a multiscale method involving all-atom and coarse-grained simulations based on generated ab initio three-dimensional structures of milk proteins for this purpose. Finally, using the adsorption energy results, we predict the composition of protein corona on iron curved and flat surfaces via a competitive adsorption model.

Self condensation of enamines mediated by acetylation. A novel approach to 1-(azol-5-yl)-(1E,3Z)-butadiene-4-N,N-dimethylamines.

Novel self-condensation of 3-(azol-5-yl)-1,1-dimethylenamines has been found to form new C-C bonds leading to 2,4-(1,2,3-triazole-1,2,3-thiadiazole-3-phenylisothiazole)-(1E,3Z)-5-yl-butadiene-1-amines. The discovered reaction represents a new example of C-H functionalization in unsaturated systems and can serve an efficient synthetic approach to rational design of new 2,4-(diazole-5-yl)-dieneamines.

Multiscale Modeling of Bio-Nano Interactions of Zero-Valent Silver Nanoparticles.

Understanding the specifics of interaction between the protein and nanomaterial is crucial for designing efficient, safe, and selective nanoplatforms, such as biosensor or nanocarrier systems. Routing experimental screening for the most suitable complementary pair of biomolecule and nanomaterial used in such nanoplatforms might be a resource-intensive task. While a range of computational tools are available for prescreening libraries of proteins for their interactions with small molecular ligands, choices for high-throughput screening of protein libraries for binding affinities to new and existing nanomaterials are very limited. In the current work, we present the results of the systematic computational study of interaction of various biomolecules with pristine zero-valent noble metal nanoparticles, namely, AgNPs, by using the UnitedAtom multiscale approach. A set of blood plasma and dietary proteins for which the interaction with AgNPs was described experimentally were examined computationally to evaluate the performance of the UnitedAtom method. A set of interfacial descriptors (log PNM, adsorption affinities, and adsorption affinity ranking), which can characterize the relative hydrophobicity/hydrophilicity/lipophilicity of the nanosized silver and its ability to form bio(eco)corona, was evaluated for future use in nano-QSAR/QSPR studies.

Representing and describing nanomaterials in predictive nanoinformatics.

Engineered nanomaterials (ENMs) enable new and enhanced products and devices in which matter can be controlled at a near-atomic scale (in the range of 1 to 100 nm). However, the unique nanoscale properties that make ENMs attractive may result in as yet poorly known risks to human health and the environment. Thus, new ENMs should be designed in line with the idea of safe-and-sustainable-by-design (SSbD). The biological activity of ENMs is closely related to their physicochemical characteristics, changes in these characteristics may therefore cause changes in the ENMs activity. In this sense, a set of physicochemical characteristics (for example, chemical composition, crystal structure, size, shape, surface structure) creates a unique 'representation' of a given ENM. The usability of these characteristics or nanomaterial descriptors (nanodescriptors) in nanoinformatics methods such as quantitative structure-activity/property relationship (QSAR/QSPR) models, provides exciting opportunities to optimize ENMs at the design stage by improving their functionality and minimizing unforeseen health/environmental hazards. A computational screening of possible versions of novel ENMs would return optimal nanostructures and manage ('design out') hazardous features at the earliest possible manufacturing step. Safe adoption of ENMs on a vast scale will depend on the successful integration of the entire bulk of nanodescriptors extracted experimentally with data from theoretical and computational models. This Review discusses directions for developing appropriate nanomaterial representations and related nanodescriptors to enhance the reliability of computational modelling utilized in designing safer and more sustainable ENMs.

Structural refinement of the hERG1 pore and voltage-sensing domains with ROSETTA-membrane and molecular dynamics simulations.

The hERG1 gene (Kv11.1) encodes a voltage-gated potassium channel. Mutations in this gene lead to one form of the Long QT Syndrome (LQTS) in humans. Promiscuous binding of drugs to hERG1 is known to alter the structure/function of the channel leading to an acquired form of the LQTS. Expectably, creation and validation of reliable 3D model of the channel have been a key target in molecular cardiology and pharmacology for the last decade. Although many models were built, they all were limited to pore domain. In this work, a full model of the hERG1 channel is developed which includes all transmembrane segments. We tested a template-driven de-novo design with ROSETTA-membrane modeling using side-chain placements optimized by subsequent molecular dynamics (MD) simulations. Although backbone templates for the homology modeled parts of the pore and voltage sensors were based on the available structures of KvAP, Kv1.2 and Kv1.2-Kv2.1 chimera channels, the missing parts are modeled de-novo. The impact of several alignments on the structure of the S4 helix in the voltage-sensing domain was also tested. Herein, final models are evaluated for consistency to the reported structural elements discovered mainly on the basis of mutagenesis and electrophysiology. These structural elements include salt bridges and close contacts in the voltage-sensor domain; and the topology of the extracellular S5-pore linker compared with that established by toxin foot-printing and nuclear magnetic resonance studies. Implications of the refined hERG1 model to binding of blockers and channels activators (potent new ligands for channel activations) are discussed.

Interactions of H562 in the S5 helix with T618 and S621 in the pore helix are important determinants of hERG1 potassium channel structure and function.

hERG1 is a member of the cyclic nucleotide binding domain family of K(+) channels. Alignment of cyclic nucleotide binding domain channels revealed an evolutionary conserved sequence HwX(A/G)C in the S5 domain. We reasoned that histidine 562 in hERG1 could play an important structure-function role. To explore this role, we created in silica models of the hERG1 pore domain based on the KvAP crystal structure with Rosetta-membrane modeling and molecular-dynamics simulations. Simulations indicate that the H562 residue in the S5 helix spans the gap between the S5 helix and the pore helix, stabilizing the pore domain, and that mutation at the H562 residue leads to a disruption of the hydrogen bonding to T618 and S621, resulting in distortion of the selectivity filter. Analysis of the simulated point mutations at positions 562/618/621 showed that the reciprocal double mutations H562W/T618I would partially restore the orientation of the 562 residue. Matching hydrophobic interactions between mutated W562 residue and I618 partially compensate for the disrupted hydrogen bonding. Complementary in vitro electrophysiological studies confirmed the results of the molecular-dynamics simulations on single mutations at positions 562, 618, and 621. Experimentally, mutations of the H562 to tryptophan produced a functional channel, but with slowed deactivation and shifted V(1/2) of activation. Furthermore, the double mutation T618I/H562W rescued the defects seen in activation, deactivation, and potassium selectivity seen with the H562W mutation. In conclusion, interactions between H562 in the S5 helix and amino acids in the pore helix are important determinants of hERG1 potassium channel function, as confirmed by theory and experiment.

Molecular mechanism of ion-ion and ion-substrate coupling in the Na+-dependent leucine transporter LeuT.

Ion-coupled transport of neurotransmitter molecules by neurotransmitter:sodium symporters (NSS) play an important role in the regulation of neuronal signaling. One of the major events in the transport cycle is ion-substrate coupling and formation of the high-affinity occluded state with bound ions and substrate. Molecular mechanisms of ion-substrate coupling and the corresponding ion-substrate stoichiometry in NSS transporters has yet to be understood. The recent determination of a high-resolution structure for a bacterial homolog of Na(+)/Cl(-)-dependent neurotransmitter transporters, LeuT, offers a unique opportunity to analyze the functional roles of the multi-ion binding sites within the binding pocket. The binding pocket of LeuT contains two metal binding sites. The first ion in site NA1 is directly coupled to the bound substrate (Leu) with the second ion in the neighboring site (NA2) only approximately 7 A away. Extensive, fully atomistic, molecular dynamics, and free energy simulations of LeuT in an explicit lipid bilayer are performed to evaluate substrate-binding affinity as a function of the ion load (single versus double occupancy) and occupancy by specific monovalent cations. It was shown that double ion occupancy of the binding pocket is required to ensure substrate coupling to Na(+) and not to Li(+) or K(+) cations. Furthermore, it was found that presence of the ion in site NA2 is required for structural stability of the binding pocket as well as amplified selectivity for Na(+) in the case of double ion occupancy.

Reactivity of polychlorinated biphenyls in nucleophilic and electrophilic substitutions.

To explain the chemical reactivity of polychlorinated biphenyls in nucleophilic (S(N)) and electrophilic (S(E)) substitutions, quantum chemical calculations were carried out at the B3LYP/6-31G(d) level of the Density Functional Theory in gas phase. Carbon atomic charges in biphenyl structure were calculated by the Atoms-in-Molecules method. Chemical hardness and global electrophilicity index parameters were determined for congeners. A comparison of calculated descriptors and experimental data for congener reactivity in the S(N) and S(E) reactions was made. It is shown that interactions in the S(N) mechanism are reactions of the hard acid-hard base type, these are the most effective in case of highly chlorinated substrates. To explain the congener reactivity in the SE reactions, correct descriptors were not established. The obtained results can be used to carry out chemical transformations of the polychlorinated biphenyls in order to prepare them for microbiological destruction or preservation.

Halothane solvation in water and organic solvents from molecular simulations with new polarizable potential function.

The partitioning of a substrate from one phase into another is a complex process with widespread applications: from chemical technology to the pharmaceutical industry. One particularly well-known and well-studied example is 2-bromo-2-chloro-1,1,1-trifluoroethane (halothane) trafficking through the lipid bilayer. Halothane is a model volatile anesthetic known to impact functions of model lipid bilayers, altering the structure and thickness upon its partitioning from the bulk phase. A number of theoretical and experimental investigations suggest the importance of electronic polarizability, determining a preference for halothane to partition in the interfacial systems as in lipid bilayers or binary solvents. The recently published protocol for the development of polarizable force fields based on the classical Drude model has provided fresh impetus to efforts directed at understanding the molecular principles governing complex thermodynamics of the hydrophobic hydration. Here, molecular simulations were combined with free energy simulations to study solvation of halothane in polarizable water and methanol. The absolute free energy of halothane solvation in different solvents (water, methanol, and n-hexane) has been evaluated for additive and polarizable models. It was found that both additive and polarizable models provide an adequate description of the halothane solvation in high-dielectric (polar) solvents such as water, but explicit accounting for electronic polarization is imperative for a correct description of the solvation thermodynamics in nonpolar systems. To study halothane dynamics in binary mixtures, all-atom molecular dynamics (MD) simulations for halothane-methanol mixtures in a wide range of concentrations were performed alongside an analysis of structural organization, dynamics, and thermodynamic properties to dissect the molecular determinants of the halothane solvation in polar and amphiphilic liquids such as methanol. Additionally, a theoretical test of the hypothesis on the weak hydrogen bonding of halothane and methanol in the condensed phase is provided, which was presented on the basis of spectroscopic analysis of the C-H vibrations in different gas-phase complexes. The simulations performed in the condensed phase suggest that hydrophobic interactions between halothane and methanol play a dominant role in preferential solvation.