Computational models for studying physical instabilities in high concentration biotherapeutic formulations

Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.

Computer Prediction of pK a Values in Small Molecules and Proteins

Accurately determining the acid dissociation constants (K _a or their logarithmic form, pK _a) of small molecules and large biomolecules has proven to be pivotal for the study different biological processes and developing new drugs. This Viewpoint summarizes some of the most common methodologies and recent advances described for pK _a prediction using computational techniques when experimental values are not easily accessible such as in proteins and/or for the screening of large libraries of new compounds.

Computing Protein pKas Using the TABI Poisson–Boltzmann Solver

A common approach to computing protein pKas uses a continuum dielectric model in which the protein is a low dielectric medium with embedded atomic point charges, the solvent is a high dielectric medium with a Boltzmann distribution of ionic charges, and the pKa is related to the electrostatic free energy which is obtained by solving the Poisson–Boltzmann equation. Starting from the model pKa for a titrating residue, the method obtains the intrinsic pKa and then computes the protonation probability for a given pH including site–site interactions. This approach assumes that acid dissociation does not affect protein conformation aside from adding or deleting charges at titratable sites. In this work, we demonstrate our treecode-accelerated boundary integral (TABI) solver for the relevant electrostatic calculations. The pKa computing procedure is enclosed in a convenient Python wrapper which is publicly available at the corresponding author’s website. Predicted results are compared with experimental pKas for several proteins. Among ongoing efforts to improve protein pKa calculations, the advantage of TABI is that it reduces the numerical errors in the electrostatic calculations so that attention can be focused on modeling assumptions.

Enhanced solvation force extrapolation for speeding up molecular dynamics simulations of complex biochemical liquids

We propose an enhanced approach to the extrapolation of mean potential forces acting on atoms of solute macromolecules due to their interactions with solvent atoms in complex biochemical liquids. It improves and extends our previous extrapolation schemes by additionally including new techniques such as an exponential scaling transformation of coordinate space with weights complemented by an automatically adjusted balancing between the least square minimization of force deviations and the norm of expansion coefficients in the approximation. The expensive mean potential forces are treated in terms of the three-dimensional reference interaction site model with Kovalenko-Hirata closure molecular theory of solvation. During the dynamics, they are calculated only after every long (outer) time interval, i.e., quite rarely to reduce the computational costs. At much shorter (inner) time steps, these forces are extrapolated on the basis of their outer values. The equations of motion are then solved using a multiple time step integration within an optimized isokinetic Nosé-Hoover chain thermostat. The new approach is applied to molecular dynamics simulations of various systems consisting of solvated organic and biomolecules of different complexity. For example, we consider hydrated alanine dipeptide, asphaltene in toluene solvent, miniprotein 1L2Y, and protein G in aqueous solution. It is shown that in all these cases, the enhanced extrapolation provides much better accuracy of the solvation force approximation than the existing approaches. As a result, it can be used with much larger outer time steps, leading to a significant speedup of the simulations.

MTS-MD of Biomolecules Steered with 3D-RISM-KH Mean Solvation Forces Accelerated with Generalized Solvation Force Extrapolation

We developed a generalized solvation force extrapolation (GSFE) approach to speed up multiple time step molecular dynamics (MTS-MD) of biomolecules steered with mean solvation forces obtained from the 3D-RISM-KH molecular theory of solvation (three-dimensional reference interaction site model with the Kovalenko-Hirata closure). GSFE is based on a set of techniques including the non-Eckart-like transformation of coordinate space separately for each solute atom, extension of the force-coordinate pair basis set followed by selection of the best subset, balancing the normal equations by modified least-squares minimization of deviations, and incremental increase of outer time step in motion integration. Mean solvation forces acting on the biomolecule atoms in conformations at successive inner time steps are extrapolated using a relatively small number of best (closest) solute atomic coordinates and corresponding mean solvation forces obtained at previous outer time steps by converging the 3D-RISM-KH integral equations. The MTS-MD evolution steered with GSFE of 3D-RISM-KH mean solvation forces is efficiently stabilized with our optimized isokinetic Nosé-Hoover chain (OIN) thermostat. We validated the hybrid MTS-MD/OIN/GSFE/3D-RISM-KH integrator on solvated organic and biomolecules of different stiffness and complexity: asphaltene dimer in toluene solvent, hydrated alanine dipeptide, miniprotein 1L2Y, and protein G. The GSFE accuracy and the OIN efficiency allowed us to enlarge outer time steps up to huge values of 1-4 ps while accurately reproducing conformational properties. Quasidynamics steered with 3D-RISM-KH mean solvation forces achieves time scale compression of conformational changes coupled with solvent exchange, resulting in further significant acceleration of protein conformational sampling with respect to real time dynamics. Overall, this provided a 50- to 1000-fold effective speedup of conformational sampling for these systems, compared to conventional MD with explicit solvent. We have been able to fold the miniprotein from a fully denatured, extended state in about 60 ns of quasidynamics steered with 3D-RISM-KH mean solvation forces, compared to the average physical folding time of 4-9 μs observed in experiment.

Technical advances in molecular simulation since the 1980s

This review describes how the theory and practice of molecular simulation have evolved since the beginning of the 1980s when the author started his career in this field. The account is of necessity brief and subjective and highlights the changes that the author considers have had significant impact on his research and mode of working.

Molecular Modeling to Study Dendrimers for Biomedical Applications

Molecular modeling techniques provide a powerful tool to study the properties of molecules and their interactions at the molecular level. The use of computational techniques to predict interaction patterns and molecular properties can inform the design of drug delivery systems and therapeutic agents. Dendrimers are hyperbranched macromolecular structures that comprise repetitive building blocks and have defined architecture and functionality. Their unique structural features can be exploited to design novel carriers for both therapeutic and diagnostic agents. Many studies have been performed to iteratively optimise the properties of dendrimers in solution as well as their interaction with drugs, nucleic acids, proteins and lipid membranes. Key features including dendrimer size and surface have been revealed that can be modified to increase their performance as drug carriers. Computational studies have supported experimental work by providing valuable insights about dendrimer structure and possible molecular interactions at the molecular level. The progress in computational simulation techniques and models provides a basis to improve our ability to better predict and understand the biological activities and interactions of dendrimers. This review will focus on the use of molecular modeling tools for the study and design of dendrimers, with particular emphasis on the efforts that have been made to improve the efficacy of this class of molecules in biomedical applications.

Structure-function relationship of PAMAM dendrimers as robust oil dispersants

PAMAM dendrimers have recently been investigated as efficient and biocompatible oil dispersants utilizing their encapsulation capacity; however, their high cationic charge density has been shown to be cytotoxic. It is therefore imperative to mitigate cationic charge-induced toxicity and understand the effects of such changes. Presented here is a synergistic experimental and computational approach to examine the effects of varying terminal surface charge on the capacity of dendrimers to disperse model liner, polycyclic aromatic, and hybrid hydrocarbons. Uncharged dendrimers collapse by forming intramolecular hydrogen bonds, which reduce the hosting capability. On the other hand, changing the surface charges from positive to negative greatly shifts the pKa of tertiary amines of the PAMAM dendrimer interior. As a result, the negatively charged dendrimers have a significant percentage of tertiary amines protonated, ∼30%. This unexpected change in the interior protonation state causes electrostatic interactions with the anionic terminal groups, leading to contraction and a marked decrease in hydrocarbon hosting capacity. The present work highlights the robust nature of dendrimer oil dispersion and also illuminates potentially unintended or unanticipated effects of varying dendrimer surface chemistry on their encapsulation or hosting efficacy, which is important for their environmental, industrial, and biomedical applications.

Thermodynamics of transport through the ammonium transporter Amt-1 investigated with free energy calculations

Amt-1 from Archaeoglobus fulgidus (AfAmt-1) belongs to the Amt/Rh family of ammonium/ammonia transporting membrane proteins. The transport mode and the precise microscopic permeation mechanism utilized by these proteins are intensely debated. Open questions concern the identity of the transported substrate (ammonia and/or ammonium) and whether the transport is passive or active. To address these questions, we studied the overall thermodynamics of the different transport modes as a function of the environmental conditions. Then, we investigated the thermodynamics of the underlying microscopic transport mechanisms with free energy calculations within a continuum electrostatics model. The formalism developed for this purpose is of general utility in the calculation of binding free energies for ligands with multiple protonation forms or other binding forms. The results of our calculations are compared to the available experimental and theoretical data on Amt/Rh proteins and discussed in light of the current knowledge on the physiological conditions experienced by microorganisms and plants. We found that microscopic models of electroneutral and electrogenic transport modes are in principle thermodynamically viable. However, only the electrogenic variants have a net thermodynamic driving force under the physiological conditions experienced by microorganisms and plants. Thus, the transport mechanism of AfAmt-1 is most likely electrogenic.