Computer simulation studies of self-assembling macromolecules

Coarse-grained (CG) molecular models are now widely used to understand the structure and functionality of macromolecular self-assembling systems. In the last few years, significant efforts have been devoted to construct quantitative CG models based on data from molecular dynamics (MD) simulations with more detailed all-atom (AA) intermolecular force fields as well as experimental thermodynamic data. We review some of the recent progress pertaining to the MD simulation of self-assembling macromolecular systems, using as illustrations the application of CG models to probe surfactant and lipid self-assembly including liposome and dendrimersome formation as well as the interaction of biomembranes with nanoparticles.     

Simulating the vapour–liquid equilibria of 1,4-dioxane

1,4-dioxane, a cyclic ether, is an emerging contaminant which is difficult to remove from water with conventional water treatment methods and resistant to biodegradation. Once a reliable force field is developed for 1,4-dioxane, molecular simulation techniques can be useful to study alternative adsorbents for its removal. For this purpose, we carried out Monte Carlo simulations in a constant volume Gibbs Ensemble to generate a force field which is capable of predicting the vapour–liquid coexistence curve and critical data of 1,4-dioxane. Results are given in comparison with experimental data and results from simulations with other force fields. Liquid densities and critical temperature are predicted in excellent agreement with experimental data using the new force field. At high temperatures, predicted vapour densities are in good agreement with experimental data, however, at lower temperatures the predicted vapour densities deviate about an order of magnitude from the experimental values. The critical density is slightly underestimated with our new force field. However, overall, the results of simulations with the new parameters give much better agreement with experimental data compared to the results obtained using other force fields.     

A force field for dynamic Cu-BTC metal-organic framework

A new force field that can describe the flexibility of Cu-BTC metal-organic framework (MOF) was developed in this work. Part of the parameters were obtained using density functional theory calculations, and the others were taken from other force fields. The new force field could reproduce well the experimental crystal structure, negative thermal expansion, vibrational properties as well as adsorption behavior in Cu-BTC. In addition, the bulk modulus of Cu-BTC was predicted using the new force field. We believe the new force field is useful in understanding the structure-property relationships for MOFs, and the approach can be extended to other MOFs.     

Molecular dynamics methods to predict peptide locations in membranes: LAH4 as a stringent test case

Determining the structure of membrane-active peptides inside lipid bilayers is essential to understand their mechanism of action. Molecular dynamics simulations can easily provide atomistic details, but need experimental validation. We assessed the reliability of self-assembling (or “minimum-bias”) and potential of mean force (PMF) approaches, using all-atom (AA) and coarse-grained (CG) force-fields. The LAH4 peptide was selected as a stringent test case, since it is known to attain different orientations depending on the protonation state of its four histidine residues.In all simulations the histidine side-chains inserted in the membrane when neutral, while they interacted with phospholipid headgroups in their charged state. This led to transmembrane orientations for neutral-His LAH4 in all minimum-bias AA simulations and in most CG trajectories. By contrast, the charged-His peptide stabilized membrane defects in AA simulations, whereas it was located at the membrane surface in some CG trajectories, and interacted with both lipid leaflets in others. This behavior is consistent with the higher antimicrobial activity and membrane-permeabilizing behavior of the charged-His LAH4. In addition, good agreement with solid-state NMR orientational data was observed in AA simulations.PMF calculations correctly predicted a higher membrane affinity for the neutral-His peptide. Interestingly, the structures and relative populations of PMF local free-energy minima corresponded to those determined in the less computationally demanding minimum-bias simulations.These data provide an indication about the possible membrane-perturbation mechanism of the charged-His LAH4 peptide: by interacting with lipid headgroups of both leaflets through its cationic side-chains, it could favor membrane defects and facilitate translocation across the bilayer.     

TOUCHSTONEX: protein structure prediction with sparse NMR data

TOUCHSTONEX, a new method for folding proteins that uses a small number of long-range contact restraints derived from NMR experimental NOE (nuclear Overhauser enhancement) data, is described. The method employs a new lattice-based, reduced model of proteins that explicitly represents C(alpha), C(beta), and the sidechain centers of mass. The force field consists of knowledge-based terms to produce protein-like behavior, including various short-range interactions, hydrogen bonding, and one-body, pairwise, and multibody long-range interactions. Contact restraints were incorporated into the force field as an NOE-specific pairwise potential. We evaluated the algorithm using a set of 125 proteins of various secondary structure types and lengths up to 174 residues. Using N/8 simulated, long-range sidechain contact restraints, where N is the number of residues, 108 proteins were folded to a C(alpha)-root-mean-square deviation (RMSD) from native below 6.5 A. The average RMSD of the lowest RMSD structures for all 125 proteins (folded and unfolded) was 4.4 A. The algorithm was also applied to limited experimental NOE data generated for three proteins. Using very few experimental sidechain contact restraints, and a small number of sidechain-main chain and main chain-main chain contact restraints, we folded all three proteins to low-to-medium resolution structures. The algorithm can be applied to the NMR structure determination process or other experimental methods that can provide tertiary restraint information, especially in the early stage of structure determination, when only limited data are available.     

Polarizable Water Model for the Coarse-Grained MARTINI Force Field

Coarse-grained (CG) simulations have become an essential tool to study a large variety of biomolecular processes, exploring temporal and spatial scales inaccessible to traditional models of atomistic resolution. One of the major simplifications of CG models is the representation of the solvent, which is either implicit or modeled explicitly as a van der Waals particle. The effect of polarization, and thus a proper screening of interactions depending on the local environment, is absent. Given the important role of water as a ubiquitous solvent in biological systems, its treatment is crucial to the properties derived from simulation studies. Here, we parameterize a polarizable coarse-grained water model to be used in combination with the CG MARTINI force field. Using a three-bead model to represent four water molecules, we show that the orientational polarizability of real water can be effectively accounted for. This has the consequence that the dielectric screening of bulk water is reproduced. At the same time, we parameterized our new water model such that bulk water density and oil/water partitioning data remain at the same level of accuracy as for the standard MARTINI force field. We apply the new model to two cases for which current CG force fields are inadequate. First, we address the transport of ions across a lipid membrane. The computed potential of mean force shows that the ions now naturally feel the change in dielectric medium when moving from the high dielectric aqueous phase toward the low dielectric membrane interior. In the second application we consider the electroporation process of both an oil slab and a lipid bilayer. The electrostatic field drives the formation of water filled pores in both cases, following a similar mechanism as seen with atomistically detailed models.     

Water activity coefficients in aqueous amino acid solutions by molecular dynamics simulation: 1. Force field development

New force fields for molecular dynamics (MD) simulation of aqueous zwitterionic amino acid simulations were developed. These were especially designed to calculate activity coefficient of water in amino acid solutions with high accuracy. For example, aqueous solutions of the following amino acids were considered: glycine, alanine, α-aminobutyric acid, α-aminovalerianic acid, valine and leucine. The force fields were obtained by quantum chemical calculations using B3LYP/6-31G and MP2/6-311(d,p) model theories in combination with the Merz–Kollmann–Singh scheme. To further increase the accuracy of the force field, a polarised continuum was considered in all quantum chemical calculations. Water activity coefficients obtained from MD using different all-purpose literature force fields, namely, OPLS, AMBER ff03 and GROMOS 53A6 as well as experimental data are compared with the results utilising the new force field. The new force field is shown to give better results compared with experimental data than existing force fields.     

Multicomponent mesofluidic system for the detection of veterinary drug residues based on competitive immunoassay

An automated multicomponent mesofluidic system (MCMS) based on biorecognitions carried out on meso-scale glass beads in polydimethylsiloxane (PDMS) channels was developed. The constructed MCMS consisted of five modules: a bead introduction module, a bioreaction module, a solution handling module, a liquid driving module, and a signal collection module. The integration of these modules enables the assay to be automated and reduces it to a one-step protocol. The MCMS has successfully been applied toward the detection of veterinary drug residues in animal-derived foods. The drug antigen-coated beads (?250 μm) were arrayed in the PDMS channels (?300 μm). The competitive immunoassay was then carried out on the surface of the glass beads. After washing, the Cy3-labeled secondary antibody was introduced to probe the antigen-antibody complex anchored to the beads. The fluorescence intensity of each bead was measured and used to determine the residual drug concentration. The MCMS is highly sensitive, with its detection limits ranging from 0.02 (salbutamol) to 3.5 μg/L (sulfamethazine), and has a short assay time of 45 min or less. The experimental results demonstrate that the MCMS proves to be an economic, efficient, and sensitive platform for multicomponent detection of compound residues for contamination in foods or the environment.     

Label-free Isolation and Enrichment of Cells Through Contactless Dielectrophoresis

Dielectrophoresis (DEP) is the phenomenon by which polarized particles in a non-uniform electric field undergo translational motion, and can be used to direct the motion of microparticles in a surface marker-independent manner. Traditionally, DEP devices include planar metallic electrodes patterned in the sample channel. This approach can be expensive and requires a specialized cleanroom environment. Recently, a contact-free approach called contactless dielectrophoresis (cDEP) has been developed. This method utilizes the classic principle of DEP while avoiding direct contact between electrodes and sample by patterning fluidic electrodes and a sample channel from a single polydimethylsiloxane (PDMS) substrate, and has application as a rapid microfluidic strategy designed to sort and enrich microparticles. Unique to this method is that the electric field is generated via fluidic electrode channels containing a highly conductive fluid, which are separated from the sample channel by a thin insulating barrier. Because metal electrodes do not directly contact the sample, electrolysis, electrode delamination, and sample contamination are avoided. Additionally, this enables an inexpensive and simple fabrication process.cDEP is thus well-suited for manipulating sensitive biological particles. The dielectrophoretic force acting upon the particles depends not only upon spatial gradients of the electric field generated by customizable design of the device geometry, but the intrinsic biophysical properties of the cell. As such, cDEP is a label-free technique that avoids depending upon surface-expressed molecular biomarkers that may be variably expressed within a population, while still allowing characterization, enrichment, and sorting of bioparticles.Here, we demonstrate the basics of fabrication and experimentation using cDEP. We explain the simple preparation of a cDEP chip using soft lithography techniques. We discuss the experimental procedure for characterizing crossover frequency of a particle or cell, the frequency at which the dielectrophoretic force is zero. Finally, we demonstrate the use of this technique for sorting a mixture of ovarian cancer cells and fluorescing microspheres (beads).     

Local measurements of viscoelastic parameters of adherent cell surfaces by magnetic bead microrheometry.

A magnetic bead microrheometer has been designed which allows the generation of forces up to 10(4) pN on 4.5 micron paramagnetic beads. It is applied to measure local viscoelastic properties of the surface of adhering fibroblasts. Creep response and relaxation curves evoked by tangential force pulses of 500-2500 pN (and approximately 1 s duration) on the magnetic beads fixed to the integrin receptors of the cell membrane are recorded by particle tracking. Linear three-phasic creep responses consisting of an elastic deflection, a stress relaxation, and a viscous flow are established. The viscoelastic response curves are analyzed in terms of a series arrangement of a dashpot and a Voigt body, which allows characterization of the viscoelastic behavior of the adhering cell surface in terms of three parameters: an effective elastic constant, a viscosity, and a relaxation time. The displacement field generated by the local tangential forces on the cell surface is visualized by observing the induced motion of assemblies of nonmagnetic colloidal probes fixed to the membrane. It is found that the displacement field decays rapidly with the distance from the magnetic bead. A cutoff radius of Rc approximately 7 micron of the screened elastic field is established. Partial penetration of the shear field into the cytoplasm is established by observing the induced deflection of intracellular compartments. The cell membrane was modeled as a thin elastic plate of shear modulus mu * coupled to a viscoelastic layer, which is fixed to a solid support on the opposite side; the former accounts for the membrane/actin cortex, and the latter for the contribution of the cytoskeleton to the deformation of the cell envelope. It is characterized by the coupling constant chi characterizing the elasticity of the cytoskeleton. The coupling constant chi and the surface shear modulus mu * are obtained from the measured displacements of the magnetic and nonmagnetic beads. By analyzing the experimental data in terms of this model a surface shear modulus of mu * approximately 2 . 10(-3) Pa m to 4 . 10(-3) Pa m is found. By assuming an approximate plate thickness of 0.1 micron one estimates an average bulk shear modulus of mu approximately (2 / 4) . 10(-4) Pa, which is in reasonable agreement with data obtained by atomic force microscopy. The viscosity of the dashpot is related to the apparent viscosity of the cytoplasm, which is obtained by assuming that the top membrane is coupled to the bottom (fixed) membrane by a viscous medium. By application of the theory of diffusion of membrane proteins in supported membranes we find a coefficient of friction of bc approximately 2 . 10(9) Pa s/m corresponding to a cytoplasmic viscosity of 2 . 10(3) Pa s.     

A less‐biased analysis of metalloproteins reveals novel zinc coordination geometries

Zinc metalloproteins are involved in many biological processes and play crucial biochemical roles across all domains of life. Local structure around the zinc ion, especially the coordination geometry (CG), is dictated by the protein sequence and is often directly related to the function of the protein. Current methodologies in characterizing zinc metalloproteins' CG consider only previously reported CG models based mainly on nonbiological chemical context. Exceptions to these canonical CG models are either misclassified or discarded as “outliers.” Thus, we developed a less‐biased method that directly handles potential exceptions without pre‐assuming any CG model. Our study shows that numerous exceptions could actually be further classified and that new CG models are needed to characterize them. Also, these new CG models are cross‐validated by strong correlation between independent structural and functional annotation distance metrics, which is partially lost if these new CGs models are ignored. Furthermore, these new CG models exhibit functional propensities distinct from the canonical CG models. Proteins 2015; 83:1470–1487. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc.     

Packing of transmembrane helices in bacteriorhodopsin folding: structure and thermodynamics

We propose a coarse-grained (CG) model to study the native structure and physical properties of helical membrane proteins (HMPs) using off-lattice computer simulations. Instead of considering sequence heterogeneity explicitly, we model its effect on the packing of helices by employing a mean packing parameter r(0), which is calculated from an all-atom (AA) model. Specifically, this CG model is applied to investigate the packing of helices in bacteriorhodopsin (BR), and predicts the seven helix bundle structure of BR with a root mean square deviation (RMSD) in coordinates of helix backbone atoms (N, C, C(alpha)) of 3.99 A from its crystal structure. This predicted structure is further refined in an AA model by Amber and the refined structure has a RMSD (in coordinates of helix backbone atoms) of 2.64 A. The predicted packing position, tilting angle, and orientation angle of each helix in the refined structure are consistent with experimental data and their physical origins can be well understood in our model. Our results show that a reasonably good structure of BR can be predicted by using such a dual-scale approach, provided that its secondary structure is known. Starting from a random initial configuration, the folded structure can be obtained in days using a regular desktop computer. Various thermodynamic properties of helix packing of BR are also investigated in this CG model.     

Disposable polydimethylsiloxane/silicon hybrid chips for protein detection

This paper presents disposable protein analysis chips with single- or four-chamber-constructed from poly(dimethylsiloxane) (PDMS) and silicon. The chips are composed of a multilayer stack of PDMS layers that sandwich a silicon microchip. This inner silicon chip features an etched array of micro-cavities hosting polymeric beads. The sample is introduced into the fluid network through the top PDMS layer, where it is directed to the bead chamber. After reaction of the analyte with the probe beads, the signal generated on the beads is captured with a CCD camera, digitally processed, and analyzed. An established bead-based fluorescent assay for C-reactive protein (CRP) was used here to characterize these hybrid chips. The detection limit of the single-chamber protein chip was found to be 1 ng/ml. Additionally, using a back pressure compensation method, the signals from each chamber of the four-chamber chip were found to fall within 10% of each other.     

Multiscale Coarse-Grained Approach to Investigate Self-Association of Antibodies

Self-association of therapeutic monoclonal antibodies (mabs) are thought to modulate the undesirably high viscosity observed in their concentrated solutions. Computational prediction of such a self-association behavior is advantageous early during mab drug candidate selection when material availability is limited. Here, we present a coarse-grained (CG) simulation method that enables microsecond molecular dynamics simulations of full-length antibodies at high concentrations. The proposed approach differs from others in two ways: first, charges are assigned to CG beads in an effort to reproduce molecular multipole moments and charge asymmetry of full-length antibodies instead of only localized charges. This leads to great improvements in the agreement between CG and all-atom electrostatic fields. Second, the distinctive hydrophobic character of each antibody is incorporated through empirical adjustments to the short-range van der Waals terms dictated by cosolvent all-atom molecular dynamics simulations of antibody variable regions. CG simulations performed on a set of 15 different mabs reveal that diffusion coefficients in crowded environments are markedly impacted by intermolecular interactions. Diffusion coefficients computed from the simulations are in correlation with experimentally measured observables, including viscosities at a high concentration. Further, we show that the evaluation of electrostatic and hydrophobic characters of the mabs is useful in predicting the nonuniform effect of salt on the viscosity of mab solutions. This CG modeling approach is particularly applicable as a material-free screening tool for selecting antibody candidates with desirable viscosity properties.     

Theoretical study of the surface properties of poly(dimethylsiloxane) and poly(tetrafluoroethylene)

Molecular dynamics (MD) simulations of poly(dimethylsiloxane) (PDMS) and poly(tetrafluoroethylene) (PTFE) were carried out to determine their surface properties and energies. This study helps to gain better insight into the molecular modeling of PDMS and PTFE, in particular how different approaches affect calculations of surface energy. Current experimental and theoretical data were used to further understand the surface properties of PDMS and PTFE as well as to validate and verify results obtained from the combination of density functional theory (DFT) calculations (including periodic boundary conditions) and MD simulations. Detailed analysis of the structure and electronic properties (by calculation of the projected density of states) of the bulk and surface models of PDMS and PTFE was performed. The sensitivity of the surface energy calculation of these two polymers to the chemistry and model preparation was indicated. The balance between the molecular density, weight (which also reflects bond orientation in the surface region), bond flexibility, and intramolecular interactions including bond stretching was revealed to govern the results obtained. In modeling, the structural organization of polymer near a given surface (types and number of end groups and broken bonds due to application of different cut offs of the periodic structure) also significantly affects the final results. Besides the structural differences, certain simulation parameters, such the DFT functionals and simulation boxes utilized, play an important role in determining surface energy. The models used here were shown to be sufficient due to their good agreement with experimental and other theoretical data related to surface properties and surface energies.     

An improved generalized AMBER force field (GAFF) for urea

We describe an improved force field parameter set for the generalized AMBER force field (GAFF) for urea. Quantum chemical computations were used to obtain geometrical and energetic parameters of urea dimers and larger oligomers using AM1 semiempirical MO theory, density functional theory at the B3LYP/6-31G(d,p) level, MP2 and CCSD ab initio calculations with the 6-311++G(d,p), aug-cc-pVDZ, aug-cc-pVTZ, and aug-cc-pVQZ basis sets, and with the CBS-QB3 and CBS-APNO complete basis set methods. Seven different urea dimer structures were optimized at the MP2/aug-cc-pVDZ level to obtain accurate interaction energies. Atomic partial charges were calculated at the MP2/aug-cc-pVDZ level with the restrained electrostatic potential (RESP) fitting approach. The interaction energies computed with these new RESP charges in the force field are consistent with those obtained from CCSD and MP2 calculations. The linear dimer structure calculated using the force field with modified geometrical parameters and the new RESP charge set agrees well with available experimental data.