Conformational Search for the Building Block of Proteins Based on the Gradient Gravitational Search Algorithm (ConfGGS) Using Force Fields: CHARMM, AMBER, and OPLS-AA

Proteins are linear polymers built from a repertoire of 20 different amino acids, which are considered building blocks of proteins. The diversity and versatility of these 20 building blocks with regard to their conformations are key to adopting three-dimensional structures that facilitate proteins to undergo important mechanistic biological processes in living systems. The present investigation reports a conformational search of 20 different amino acids, building blocks of proteins, using three different force fields, CHARMM, AMBER, and OPLS-AA, implemented in the gradient gravitational search algorithm. The search technique (ConfGGS) includes the contribution from both bonded and nonbonded terms using Cartesian coordinates. The efficiency of such conformational searches has also been compared with other optimization algorithms: DE/Best, DE/Rand, and PSO algorithms with respect to computational time and accuracy based on the minimum number of iteration steps and computed lowest mean absolute error (MAE) and mean standard deviation (MSD) values for dihedral angles of respective near-optimal structures. Moreover, the ConfGGS technique has also been extended to an ordered protein fragment (PQITL) extracted from HIV-1 protease (PDB ID: 1YTH), an intrinsically disordered protein fragment, i.e., an amyloid-forming segment (AVVTGVTAV), from the NAC domain of Parkinson's disease protein α-synuclein, residues 69-77 (PDB ID: 4RIK), the experimental NMR atomic-resolution structure of α-synuclein fibrils (PDB ID: 2N0A), and a disulfide bond-containing protein fragment sequence (PCYGWPVCY), residues 59-67 (PDB ID: 6Y4F) toward structure prediction as a close homologue compared with experimental accuracy, using the CHARMM force field. The MolProbity validation results for the protein fragment (PQITL) obtained by ConfGGS/CHARMM are in better agreement with the native protein fragment structure of HIV-1 protease (PDB ID: 1YTH). Furthermore, the computed results have also been compared with the coordinates obtained from the AlphaFold network.

Adaptive Synthesis of Functional Amphiphilic Dendrons as a Novel Approach to Artificial Supramolecular Objects

Supramolecular structures, such as micelles, liposomes, polymerosomes or dendrimerosomes, are widely studied and used as drug delivery systems. The behavior of amphiphilic building blocks strongly depends on their spatial distribution and shape of polar and nonpolar component. This report is focused on the development of new versatile synthetic protocols for amphiphilic carbosilane dendrons (amp-CS-DDNs) capable of self-assembly to regular micelles and other supramolecular objects. The presented strategy enables the fine modification of amphiphilic structure in several ways and also enables the facile connection of a desired functionality. DLS experiments demonstrated correlations between structural parameters of amp-CS-DDNs and the size of formed nanoparticles. For detailed information about the organization and spatial distribution of amp-CS-DDNs assemblies, computer simulation models were studied by using molecular dynamics in explicit water.

Carbosilane Glycodendrimers for Anticancer Drug Delivery: Synthetic Route, Characterization, and Biological Effect of Glycodendrimer-Doxorubicin Complexes

The complexity of drug delivery mechanisms calls for the development of new transport system designs. Here, we report a robust synthetic procedure toward stable glycodendrimer (glyco-DDM) series bearing glucose, galactose, and oligo(ethylene glycol)-modified galactose peripheral units. In vitro cytotoxicity assays showed exceptional biocompatibility of the glyco-DDMs. To demonstrate applicability in drug delivery, the anticancer agent doxorubicin (DOX) was encapsulated in the glyco-DDM structure. The anticancer activity of the resulting glyco-DDM/DOX complexes was evaluated on the noncancerous (BJ) and cancerous (MCF-7 and A2780) cell lines, revealing their promising generation- and concentration-dependent effect. The glyco-DDM/DOX complexes show gradual and pH-dependent DOX release profiles. Fluorescence spectra elucidated the encapsulation process. Confocal fluorescence microscopy demonstrated preferential cancer cell internalization of the glyco-DDM/DOX complexes. The conclusions were supported by computer modeling. Overall, our results are consistent with the assumption that novel glyco-DDMs and their drug complexes are very promising in drug delivery and related applications.

Adsorption of Organic Acids and Phosphate to an Iron (Oxyhydr)oxide Mineral: A Combined Experimental and Density Functional Theory Study

The interaction of soil organic matter with mineral surfaces is a critical reaction involved in many ecosystem services, including stabilization of organic matter in the terrestrial carbon pool and bioavailability of plant nutrients. Using model organic acids typically present in soil solutions, this study couples laboratory adsorption studies with density functional theory (DFT) to provide physical insights into the nature of the chemical bonding between carboxylate functional groups and a model FeOOH cluster. Topological determination of electron density at bond critical points using quantum theory of atoms in molecules (QTAIM) analysis revealed that the presence of multiple bonding paths between the organic acid and the FeOOH cluster is essential in determining the competitive adsorption of organic acids and phosphate for FeOOH surface adsorption sites. The electron density and Laplacian parameter values from QTAIM indicated that the primary carboxylate-FeOOH bond was more ionic than covalent in nature. The experimental and computational results provide molecular-level evidence of the important role of electrostatic forces in the bonding between carboxylic acids and Fe-hydroxides. This knowledge may assist in the formulation of management studies to meet the challenges of maintaining ecosystems services in the face of a changing climate.

A new benchmark illustrates that integration of geometric constraints inferred from enzyme reaction chemistry can increase enzyme active site modeling accuracy

Enzymes play a critical role in a wide array of industrial, medical, and research applications and with the recent explosion of genomic sequencing, we now have sequences for millions of enzymes for which there is no known structure. In order to utilize modern computational design tools for constructing inhibitors or engineering novel catalysts, the ability to accurately model enzymes is critical. A popular approach for modeling enzymes are comparative modeling techniques which can often accurately predict the global structural features. However, achieving atomic accuracy of an active site remains a challenge and is an issue when trying to utilize the molecular details for designing inhibitors or enhanced catalysts. Here we explore integrating knowledge about the required geometric orientation of conserved catalytic residues into the comparative modeling process in order to improve modeling accuracy. In order to investigate the utility of adding this information, we first carefully construct a benchmark set of reference structures to use. Consistent with previous findings, our benchmark demonstrates that the geometry between catalytic residues across an enzyme family is conserved and does not tend to deviate by more than 0.5Å. We then find that by integrating these geometric constraints during modeling, we can double the number of atomic level accuracy models (<1Å RMSD to the crystal structure ligand) within our benchmarking dataset, even for targets with templates as low as 20-30% sequence identity. Catalytic residues within an enzyme family are highly conserved and can often be readily identified through comparative sequence analysis to a known structure within the enzyme family. Therefore utilizing this readily available information has the potential to significantly improve drug design and enzyme engineering efforts for which there is no known structure for the enzyme of interest.

Anharmonic Molecular Mechanics: Ab Initio Based Morse Parametrizations for the Popular MM3 Force Field

Methodologies for creating reactive potential energy surfaces from molecular mechanics force-fields are becoming increasingly popular. To date, molecular mechanics force-fields in biochemistry and small molecule organic chemistry tend to use harmonic expressions to treat bonding stretches, which is a poor approximation in reactive and nonequilibirum molecular dynamics simulations since bonds are often displaced significantly from their equilibrium positions. For such applications there is need for a better treatment of anharmonicity. In this contribution, Morse bonding potentials have been extensively parametrized for the atom types in the MM3 force field of Allinger and co-workers using high level CCSD(T)(F12*) energies. To our knowledge this is among the first instances of a comprehensive parametrization of Morse potentials in a popular organic chemistry force field. In the context of molecular dynamics simulations, these data will: (1) facilitate the fitting of reactive potential energy surfaces using empirical valence bond approaches and (2) enable more accurate treatments of energy transfer.

Anticancer Activity of Dendriplexes against Advanced Prostate Cancer from Protumoral Peptides and Cationic Carbosilane Dendrimers

The interaction of neuropeptides, vasoactive intestinal peptide (VIP), or growth hormone-releasing hormone (GHRH), with a cationic carbosilane dendrimer forms dendriplexes with antitumoral behavior in advanced prostate cancer cells PC3. At the concentrations used for dendriplexes formation, the free peptides were protumoral and prometastatic in advanced prostate cancer, while dendrimer only showed low cytotoxicity, but did not avoid the metastatic behavior of PC3 cells. However, these nanoplexes favored also cell adhesion and avoided cell migration. Also, the dendriplexes were not toxic for no tumoral prostate cells (RPWE-1) or fibroblasts. The use of labeled GHRH peptide (rhodamine labeled) and a dendrimer (fluorescein labeled) allowed us to observe that both systems reach the intracellular milieu after dendriplex formation. The treatment of PC3 cells with the nanoplexes reduced expression of vascular endothelial growth factor (VEGF) and cyclic adenosine monophosphate (cAMP). Molecular modeling analysis highlights the important contribution of the carbosilane framework in the stabilization of the dendriplex, since dendrimer interacts with a peptide region where hydrophobic amino acids are presented.

Tinker 8: Software Tools for Molecular Design

The Tinker software, currently released as version 8, is a modular molecular mechanics and dynamics package written primarily in a standard, easily portable dialect of Fortran 95 with OpenMP extensions. It supports a wide variety of force fields, including polarizable models such as the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. The package runs on Linux, macOS, and Windows systems. In addition to canonical Tinker, there are branches, Tinker-HP and Tinker-OpenMM, designed for use on message passing interface (MPI) parallel distributed memory supercomputers and state-of-the-art graphical processing units (GPUs), respectively. The Tinker suite also includes a tightly integrated Java-based graphical user interface called Force Field Explorer (FFE), which provides molecular visualization capabilities as well as the ability to launch and control Tinker calculations.

Probing chirality recognition of protonated glutamic acid dimers by gas-phase vibrational spectroscopy and first-principles simulations

We characterize stereospecific aspects of homochiral and heterochiral dimers of glutamic acid by infrared spectroscopy and first-principles molecular dynamics simulations.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Phosphonium carbosilane dendrimers for biomedical applications – synthesis, characterization and cytotoxicity evaluation

Phosphonium carbosilane dendrimers could represent an alternative to ammonium ones in gene therapy applications with high potential of mitochondrial targeting.

Solvation Dynamics of CO₂(g) by Monoethanolamine at the Gas-Liquid Interface: A Molecular Mechanics Approach

A classical force field approach was used to characterize the solvation dynamics of high-density CO₂(g) by monoethanolamine (MEA) at the air–liquid interface. Intra- and intermolecular CO₂ and MEA potentials were parameterized according to the energetics calculated at the MP2 and BLYP-D2 levels of theory. The thermodynamic properties of CO₂ and MEA, such as heat capacity and melting point, were consistently predicted using this classical potential. An approximate interfacial simulation for CO₂(g)/MEA(l) was performed to monitor the depletion of the CO₂(g) phase, which was influenced by amino and hydroxyl groups of MEA. There are more intramolecular hydrogen bond interactions notably identified in the interfacial simulation than the case of bulk MEA(l) simulation. The hydroxyl group of MEA was found to more actively approach CO₂ and overpower the amino group to interact with CO₂ at the air–liquid interface. With artificially reducing the dipole moment of the hydroxyl group, CO₂–amino group interaction was enhanced and suppressed CO₂(g) depletion. The hydroxyl group of MEA was concluded to play dual but contradictory roles for CO₂ capture.

Conformational Interconversions of Amino Acid Derivatives

Exhaustive conformational interconversions including transition structure analyses of N-acetyl-l-glycine-N-methylamide as well as its alanine, serine, and cysteine analogues have been investigated at the MP2/6-31G** level, yielding a total of 142 transition states. Improved estimates of relative energies were obtained by separately extrapolating the Hartree-Fock and MP2 energies to the basis set limit and adding the difference between CCSD(T) and MP2 results with the cc-pVDZ basis set to the extrapolated MP2 results. The performance of eight empirical force fields (AMBER94, AMBER14SB, MM2, MM3, MMFFs, CHARMM22_CMAP, OPLS_2005, and AMOEBAPRO13) in reproducing ab initio energies of transition states was tested. Our results indicate that commonly used class I force fields employing a fixed partial charge model for the electrostatic interaction provide mean errors in the ∼10 kJ/mol range for energies of conformational transition states for amino acid conformers. Modern reparametrized versions, such as CHARMM22_CMAP, and polarizable force fields, such as AMOEBAPRO13, have slightly lower mean errors, but maximal errors are still in the 35 kJ/mol range. There are differences between the force fields in their ability for reproducing conformational transitions classified according to backbone/side-chain or regions in the Ramachandran angles, but the data set is likely too small to draw any general conclusions. Errors in conformational interconversion barriers by ∼10 kJ/mol suggest that the commonly used force field may bias certain types of transitions by several orders of magnitude in rate and thus lead to incorrect dynamics in simulations. It is therefore suggested that information for conformational transition states should be included in parametrizations of new force fields.

Databases of Conformations and NMR Structures of Glycan Determinants

The present study reports a comprehensive nuclear magnetic resonance (NMR) characterization and a systematic conformational sampling of the conformational preferences of 170 glycan moieties of glycosphingolipids as produced in large-scale quantities by bacterial fermentation. These glycans span across a variety of families including the blood group antigens (A, B and O), core structures (Types 1, 2 and 4), fucosylated oligosaccharides (core and lacto-series), sialylated oligosaccharides (Types 1 and 2), Lewis antigens, GPI-anchors and globosides. A complementary set of about 100 glycan determinants occurring in glycoproteins and glycosaminoglycans has also been structurally characterized using molecular mechanics-based computation. The experimental and computational data generated are organized in two relational databases that can be queried by the user through a user-friendly search engine. The NMR ((1)H and (13)C, COSY, TOCSY, HMQC, HMBC correlation) spectra and 3D structures are available for visualization and download in commonly used structure formats. Emphasis has been given to the use of a common nomenclature for the structural encoding of the carbohydrates and each glycan molecule is described by four different types of representations in order to cope with the different usages in chemistry and biology. These web-based databases were developed with non-proprietary software and are open access for the scientific community available at http://glyco3d.cermav.cnrs.fr.

Carbohydrate-protein interactions: molecular modeling insights

The article reviews the significant contributions to, and the present status of, applications of computational methods for the characterization and prediction of protein-carbohydrate interactions. After a presentation of the specific features of carbohydrate modeling, along with a brief description of the experimental data and general features of carbohydrate-protein interactions, the survey provides a thorough coverage of the available computational methods and tools. At the quantum-mechanical level, the use of both molecular orbitals and density-functional theory is critically assessed. These are followed by a presentation and critical evaluation of the applications of semiempirical and empirical methods: QM/MM, molecular dynamics, free-energy calculations, metadynamics, molecular robotics, and others. The usefulness of molecular docking in structural glycobiology is evaluated by considering recent docking- validation studies on a range of protein targets. The range of applications of these theoretical methods provides insights into the structural, energetic, and mechanistic facets that occur in the course of the recognition processes. Selected examples are provided to exemplify the usefulness and the present limitations of these computational methods in their ability to assist in elucidation of the structural basis underlying the diverse function and biological roles of carbohydrates in their dialogue with proteins. These test cases cover the field of both carbohydrate biosynthesis and glycosyltransferases, as well as glycoside hydrolases. The phenomenon of (macro)molecular recognition is illustrated for the interactions of carbohydrates with such proteins as lectins, monoclonal antibodies, GAG-binding proteins, porins, and viruses.

Cationic carbosilane dendrimers and oligonucleotide binding: an energetic affair

Generation 2 cationic carbosilane dendrimers hold great promise as internalizing agents for gene therapy as they present low toxicity and retain and internalize the genetic material as an oligonucleotide or siRNA. In this work we carried out complete in silico structural and energetical characterization of the interactions of a set of G2 carbosilane dendrimers, showing different affinity towards two single strand oligonucleotide (ODN) sequences in vitro. Our simulations predict that these four dendrimers and the relevant ODN complexes are characterized by similar size and shape, and that the molecule-specific ODN binding ability can be rationalized only by considering a critical molecular design parameter: the normalized effective binding energy ΔG(bind,eff)/N(eff), i.e. the performance of each active individual dendrimer branch directly involved in a binding interaction.

All-atom empirical potential for molecular modeling and dynamics studies of proteins

New protein parameters are reported for the all-atom empirical energy function in the CHARMM program. The parameter evaluation was based on a self-consistent approach designed to achieve a balance between the internal (bonding) and interaction (nonbonding) terms of the force field and among the solvent-solvent, solvent-solute, and solute-solute interactions. Optimization of the internal parameters used experimental gas-phase geometries, vibrational spectra, and torsional energy surfaces supplemented with ab initio results. The peptide backbone bonding parameters were optimized with respect to data for N-methylacetamide and the alanine dipeptide. The interaction parameters, particularly the atomic charges, were determined by fitting ab initio interaction energies and geometries of complexes between water and model compounds that represented the backbone and the various side chains. In addition, dipole moments, experimental heats and free energies of vaporization, solvation and sublimation, molecular volumes, and crystal pressures and structures were used in the optimization. The resulting protein parameters were tested by applying them to noncyclic tripeptide crystals, cyclic peptide crystals, and the proteins crambin, bovine pancreatic trypsin inhibitor, and carbonmonoxy myoglobin in vacuo and in crystals. A detailed analysis of the relationship between the alanine dipeptide potential energy surface and calculated protein φ, χ angles was made and used in optimizing the peptide group torsional parameters. The results demonstrate that use of ab initio structural and energetic data by themselves are not sufficient to obtain an adequate backbone representation for peptides and proteins in solution and in crystals. Extensive comparisons between molecular dynamics simulations and experimental data for polypeptides and proteins were performed for both structural and dynamic properties. Energy minimization and dynamics simulations for crystals demonstrate that the latter are needed to obtain meaningful comparisons with experimental crystal structures. The presented parameters, in combination with the previously published CHARMM all-atom parameters for nucleic acids and lipids, provide a consistent set for condensed-phase simulations of a wide variety of molecules of biological interest.

Current state-of-the-art molecular dynamics methods and applications

Molecular dynamics simulations are used to describe the patterns, strength, and properties of protein behavior, drug-receptor interactions, the solvation of molecules, the conformational changes that a protein or molecule may undergo under various conditions, and other events that require the systematic evaluation of molecular properties in dynamic molecular systems. Only few years ago proteins were considered to be rigid body structures with very limited conformational flexibility. However, it is now clear that proteins are highly dynamic structures, the internal organization of which is the key to their 3D spatial arrangement and hence biological function. The study of protein dynamics in the lab is a very complicated, expensive, and time-consuming process. Therefore, a lot of effort and hope lies with the computers and the in silico study of protein structure and molecular dynamics. Herein, an effort has been made to describe the ever-evolving field of molecular dynamics, the different algorithms, and force fields that are being used as well as to provide some insight on what the near future holds for this auspicious field of computational structural biology.

Synthesis, structure and molecular modelling of anionic carbosilane dendrimers

Anionic carbosilane dendrimers of generations 1-3 have been synthesized containing carboxylate G(n)X(C(2)H(4)CO(2)Na)(m) and sulfonate G(n)X(C(2)H(4)SO(3)Na)(m) peripheral groups and derived from two different cores, 1,3,5-(HO)(3)C(6)H(3) (X = O(3)) and Si(C(3)H(5))(4) (X = Si). The peripheral anionic groups make these dendrimers water soluble, despite their highly hydrophobic framework. These dendrimers present a net negative charge in water, which was influenced by the pH of the medium. This characteristic was studied by pH titration. Also molecular modeling calculations have been performed to study differences in an aqueous medium between carboxylate and sulfonate dendrimers of different cores. The results obtained were also compared with those obtained from DOSY NMR experiments and zeta-potential measurements.

Nature of protein-CO2 interactions as elucidated via molecular dynamics

Rising global temperatures require innovative measures to reduce atmospheric concentrations of CO(2). The most successful carbon capture technology on Earth is the enzymatic capture of CO(2) and its sequestration in the form of glucose. Efforts to improve upon or mimic this naturally occurring process will require a rich understanding of protein-CO(2) interactions. Toward that end, extensive all-atom molecular dynamics (MD) simulations were performed on the CO(2)-utilizing enzyme phosphoenolpyruvate carboxykinase (PEPCK). Preliminary simulations were performed using implicit and explicit solvent models, which yielded similar results: arginine, lysine, tyrosine, and asparagine enhance the ability of a protein to bind carbon dioxide. Extensive explicit solvent simulations were performed for both wild-type PEPCK and five single-point PEPCK mutants, revealing three prevalent channels by which CO(2) enters (or exits) the active site cleft, as well as a fourth channel (observed only once), the existence of which can be rationalized in terms of the position of a key Arg residue. The strongest CO(2) binding sites in these simulations consist of appropriately positioned hydrogen bond donors and acceptors. Interactions between CO(2) and both Mn(2+) and Mg(2+) present in PEPCK are minimal due to the stable protein- and solvent-based coordination environments of these cations. His 232, suggested by X-ray crystallography as being a potential important CO(2) binding site, is indeed found to be particularly "CO(2)-philic" in these simulations. Finally, a recent mechanism, proposed on the basis of X-ray crystallography, for PEPCK active site lid closure is discussed in light of the MD trajectories. Overall, the results of this work will prove useful not only to scientists investigating PEPCK, but also to groups seeking to develop an environmentally benign, protein-based carbon capture, sequestration, and utilization system.

Complexes of sulfuric acid with N,N-dimethylformamide: An ab initio investigation

The (H2SO4)2, H2SO4-DMF, and (H2SO4)2-DMF complexes have been investigated, using the B3LYP functional with cc-pVQZ basis set. The characteristics of structure and energetics for binary complexes of sulfuric acid with dimethylformamide (DMF) have been obtained for the first time. The H-bond formation both between molecules of sulfuric acid as well as sulfuric acid-DMF were studied, on the basis of Weinhold’s natural bond orbital (NBO) analysis. It was shown that the H-bond formation between sulfuric acid and DMF molecules is stronger than ones for the acids dimer. The value of charge transfer from lone pair (LP) orbitals of DMF oxygen to the antibonding orbital of acid OH-bond significantly exceeds the criterion of H-bond existance (0.01 e). As follows from energy, among the complexes under investigation the most preferable one was found to be (H2SO4)2-DMF in which sulfuric acid molecules are linked with each other by three H-bonds.