Search

In-silico design of polymeric biomaterials requires molecular dynamics (MD) simulations that retain essential atomistic/molecular details (e.g., explicit water around the biofunctional macromolecule) while simultaneously achieving large length and time scales pertinent to macroscale function. Such large-scale atomistically detailed macromolecular MD simulations with explicit solvent representation are computationally expensive. One way to overcome this limitation is to use an adaptive resolution scheme (AdResS) in which the explicit solvent molecules dynamically adopt either atomistic or coarse-grained resolution depending on their location (e.g., near or far from the macromolecule) in the system. In this study we present the feasibility and the limitations of AdResS methodology for studying polyethylene glycol (PEG) in adaptive resolution water, for varying PEG length and architecture. We first validate the AdResS methodology for such systems, by comparing PEG and solvent structure with that from all-atom simulations. We elucidate the role of the atomistic zone size and the need for calculating thermodynamic force correction within this AdResS approach to correctly reproduce the structure of PEG and water. Lastly, by varying the PEG length and architecture, we study the hydration of PEG, and the effect of PEG architectures on the structural properties of water. Changing the architecture of PEG from linear to multiarm star, we observe reduction in the solvent accessible surface area of the PEG, and an increase in the order of water molecules in the hydration shells.

Stimuli-responsive biomaterials are used to facilitate drug and gene delivery by shielding the drug/gene during circulation times and selectively releasing the cargo at the desired target. Within stimuli-responsive materials, pH-responsive materials are exploited for delivery to specific organs, intracellular compartments, cancer cells, site of inflammation or infection as those sites are characterized by pH that is different from the blood pH. In this paper we use molecular dynamics (MD) simulations to design such pH-responsive biomaterials where the balance between the various intermolecular interactions (e.g., electrostatics, van der Waals) within the biomaterials allow biofunctional molecules to be reversibly shielded and exposed to the environment with change in pH. In our model the shielding aspect is imparted by a polyethylene glycol (PEG) brush and the pH-responsive component is a PEG-tethered oligopeptide that undergoes changes in conformations via protonation of residues upon changes in pH. Starting with a PEG-tethered peptide in a monodisperse short PEG brush, we first vary the composition and sequence of histidine (H), lysine (K), and glutamate (E) along the oligopeptide sequence to find the design parameters that maximize the shielding and exposure of the oligopeptide at pH ∼ 7.0 and pH < 7.0, respectively. Then, we probe the effect of the PEG brush on the conformations of the oligopeptides by simulating PEG-tethered peptide in a bimodal PEG brush containing short PEG and long PEG chains. We characterize the intermolecular interactions involving the PEG, peptide, and solvent that influence the shielded and exposed conformations of the oligopeptides at the two different pHs. In a short monodisperse PEG brush, with a longer PEG-tethered peptide containing large blocks of histidines that undergo change in protonation state as a response to pH change, placed between a protonated lysine and deprotonated glutamate, the PEG brush exhibits maximum shielding and exposure with pH change. This change from shielded to exposed state is driven by electrostatic repulsion upon H protonation. The presence of long PEG chains in a bimodal PEG brush leads to dominating PEG-peptide attractive interactions that reduces the contrast in shielded and exposed conformations of the PEG-tethered peptide upon protonation of histidines.

Molecular dynamics (MD) simulations were used to gain insight on the molecular interactions in a model biological membrane comprised of a bilayer with DPPC (dipalmitoylphosphotidylcholine) and antimicrobial dendritic amphiphile molecules [RCONHC(CH(2)CH(2)COOH)(3), where R is the saturated hydrocarbon tail (R = n-C(n)H(2n+1)), to be abbreviated as 3CAmn]. This study analyzes different biophysical properties of the equilibrated mixed bilayers, at 300 and 325 K, to determine how the presence of the 3CAmn, in varying concentrations and tail lengths, affects the lipid bilayer. Lipid tail order parameter data, bilayer thickness trends, and qualitative lipid tail tilt observations suggest that a molar ratio of 0.2 3CAm19/DPPC is sufficient to induce a phase transition in the bilayer from gel to liquid crystalline at 300 K. These results also imply that the phase transition temperature of the mixed bilayer decreases upon incorporation of higher concentrations of 3CAm19. Hydrogen bonding takes place between the 3CAmn and DPPC at specific sites, as evidenced by the radial distribution function. Increased hydrogen bonding and the smaller headgroup size of the 3CAmn molecule result in a decrease in the total lateral area with higher concentrations of 3CAm19. Diffusion constants of 3CAmn varied with concentration and tail length; diffusion constants of DPPC and 3CAm19 increased with increasing 3CAm19 concentration at 300 K and shorter 3CAmn tails had higher diffusion constants at both temperatures. These computational studies provide a comprehensive understanding of the biophysical changes to model biological membranes by the association of 3CAmn.