Search

Powerful, facile new ways to create libraries of site-directed mutants are demonstrated. These include: (1) one-pot-PCR, (2) multi-pot-PCR, and (3) split-mix-PCR. One-pot-PCR uses mutant oligonucleotides to generate megaprimers in situ, and it was used to randomly incorporate 28 mutations in a gabT gene in a single reaction. In more difficult cases, multi-pot-PCR can be employed: mutant megaprimers are synthesized individually, then combined in a single mutagenesis PCR. This method was used to incorporate 14 out of 15 mutations in a pabB gene. Split-mix-PCR is a conceptually novel method for creation of site-directed mutant libraries. Separate PCRs for each mutant primer are performed, followed by pooling the products of the individual reactions. The pooled mixture is re-aliquoted into individual mutant oligonucleotide PCRs. These steps are repeated for each cycle. Split-mix-PCR results in a nearly random distribution of mutation sites, and a distribution of number-of-mutations per gene that is computable and narrow. Split-mix-PCR was applied to the directed evolution of aminodeoxychorismate synthase into anthranilate synthase, and easily allowed the determination of the fewest mutations required for introduction of novel activity.

We compute potentials of mean force for bend and twist deformations via force pulling and umbrella sampling experiments for four β-solenoid proteins (BSPs) that show promise in nanotechnology applications. In all cases, we find quasi-Hooke's law behavior until the point of rupture. Bending moduli show modest anisotropy for two-sided and three-sided BSPs, and little anisotropy for a four-sided BSP. There is a slight clockwise/counterclockwise asymmetry in the twist potential of mean force, showing greater stiffness when the applied twist follows the intrinsic twist. When we extrapolate to beam theory appropriate for amyloid fibrils of the BSPs, we find bend/twist moduli which are somewhat smaller than those in the literature for other amyloid fibrils. Twist persistence lengths are on the order of a micron, and bend persistence lengths are several microns. Provided the intrinsic twist can be reversed, these results support the usage of BSPs in biomaterials applications.

Measurements from hydroxyl radical footprinting (HRF) provide rich information about the solvent accessibility of amino acid side chains of a protein. Traditional HRF data analyses focus on comparing the difference in the modification/footprinting rate of a specific site to infer structural changes across two protein states, e.g., between a free and ligand-bound state. However, the rate information itself is not fully used for the purpose of comparing different protein sites within a protein on an absolute scale. To provide such a cross-site comparison, we present a new, to our knowledge, data analysis algorithm to convert the measured footprinting rate constant to a protection factor (PF) by taking into account the known intrinsic reactivity of amino acid side chain. To examine the extent to which PFs can be used for structural interpretation, this PF analysis is applied to three model systems where radiolytic footprinting data are reported in the literature. By visualizing structures colored with the PF values for individual peptides, a rational view of the structural features of various protein sites regarding their solvent accessibility is revealed, where high-PF regions are buried and low-PF regions are more exposed to the solvent. Furthermore, a detailed analysis correlating solvent accessibility and local structural contacts for gelsolin shows a statistically significant agreement between PF values and various structure measures, demonstrating that the PFs derived from this PF analysis readily explain fundamental HRF rate measurements. We also tested this PF analysis on alternative, chemical-based HRF data, showing improved correlations of structural properties of a model protein barstar compared to examining HRF rate data alone. Together, this PF analysis not only permits a novel, to our knowledge, approach of mapping protein structures by using footprinting data, but also elevates the use of HRF measurements from a qualitative, cross-state comparison to a quantitative, cross-site assessment of protein structures in the context of individual conformational states of interest.

We present estimates of ultimate tensile strength (UTS) for two engineered β -solenoid protein mutant fibril structures (spruce budworm and Rhagium inquisitor antifreeze proteins) derived from sonication-based measurements and from force pulling molecular dynamics simulations, both in water. Sonication experiments generate limiting scissioned fibrils with a well-defined length-to-width correlation for the mutant spruce budworm protein and the resultant UTS estimate is 0.66 ± 0.08 GPa. For fibrils formed from engineered R. inquisitor antifreeze protein, depending upon geometry, we estimate UTSs of 3.5 ± 3.2–5.5 ± 5.1 GPa for proteins with interfacial disulfide bonds, and 1.6 ± 1.5–2.5 ± 2.3 GPa for the reduced form. The large error bars for the R. inquisitor structures are intrinsic to the broad distribution of limiting scission lengths. Simulations provide pulling velocity-dependent UTSs increasing from 0.2 to 1 GPa in the available speed range, and 1.5 GPa extrapolated to the speeds expected in the sonication experiments. Simulations yield low-velocity values for the Young's modulus of 6.0 GPa. Without protein optimization, these mechanical parameters are similar to those of spider silk and Kevlar, but in contrast to spider silk, these proteins have a precisely known sequence-structure relationship.

The use of beta solenoid proteins as functionalizable, nanoscale, self-assembling molecular building blocks may have many applications, including templating the growth of wires or higher-dimensional structures. By understanding their mechanical strengths, we can efficiently design the proteins for specific functions. We present a study of the mechanical properties of seven beta solenoid proteins using GROMACS molecular dynamics software to produce force/torque-displacement data, implement umbrella sampling of bending/twisting trajectories, produce Potentials of Mean Force (PMFs), extract effective spring constants, and calculate rigidities for two bending and two twisting directions for each protein. We examine the differences between computing the strength values from force/torque-displacement data alone and PMF data, and show how higher precision estimates can be obtained from the former. In addition to the analysis of the methods, we report estimates for the bend/twist persistence lengths for each protein, which range from 0.5–3.4 μm. We note that beta solenoid proteins with internal disulfide bridges do not enjoy enhanced bending or twisting strength, and that the strongest correlate with bend/twist rigidity is the number of hydrogen bonds per turn. In addition, we compute estimates of the Young's modulus (Y) for each protein, which range from Y = 3.5 to 7.2 GPa.

The ligand-binding domain (LBD) of an estrogen receptor undergoes a large conformational switching from an inactive to active state in response to hormone stimuli. Very recently, a novel D538G mutant has been identified to be active in advanced breast cancer tumors. Here, we ask if molecular simulations can provide insight on its mechanistic impact on the receptor's activation status. It has been challenging for ab initio modeling to identify two distinct conformations of a single amino acid sequence as large as that of the LBD. Using a coarse-grained (CG) model, we are able to correctly reproduce this LBD conformational switching. Furthermore, we found that the D538G mutation reshapes the energy landscape by stabilizing both active and inactive conformations, but preferring the active by 1.5 kcal/mol. This observation is consistent with the concept of a mutation-shifting landscape and provides a structural explanation for the oncogenic D538G mutation at the detailed conformational level.

Estrogen receptor alpha (ERa) is a hormone-responsive transcription factor that contains several discrete functional domains, including a ligand-binding domain (LBD) and a DNA-binding domain (DBD). Despite a wealth of knowledge about the behaviors of individual domains, the molecular mechanisms of cross-talk between LBD and DBD during signal transduction from hormone to DNA-binding of ERa remain elusive. Here, we apply a multiscale approach combining coarse-grained (CG) and atomistically detailed simulations to characterize this cross-talk mechanism via an investigation of the ERa conformational landscape. First, a CG model of ERa is built based on crystal structures of individual LBDs and DBDs, with more emphasis on their interdomain interactions. Second, molecular dynamics simulations are implemented and enhanced sampling is achieved via the “push-pull-release” strategy in the search for different LBD–DBD orientations. Third, multiple energetically stable ERa conformations are identified on the landscape. A key finding is that estradiol-bound LBDs utilize the well-described activation helix H12 to pack and stabilize LBD–DBD interactions. Our results suggest that the estradiol-bound LBDs can serve as a scaffold to position and stabilize the DBD–DNA complex, consistent with experimental observations of enhanced DNA binding with the LBD. Final assessment using atomic-level simulations shows that these CG-predicted models are significantly stable within a 15-ns simulation window and that specific pairs of lysine residues in close proximity at the domain interfaces could serve as candidate sites for chemical cross-linking studies. Together, these simulation results provide a molecular view of the role of ERa domain interactions in response to hormone binding.

Understanding protein-protein association is crucial in revealing the molecular basis of many biological processes. Here, we describe a theoretical simulation pipeline to study protein-protein association from an energy landscape perspective. First, a coarse-grained model is implemented and its applications are demonstrated via molecular dynamics simulations for several protein complexes. Second, an enhanced search method is used to efficiently sample a broad range of protein conformations. Third, multiple conformations are identified and clustered from simulation data and further projected on a three-dimensional globe specifying protein orientations and interacting energies. Results from several complexes indicate that the crystal-like conformation is favorable on the energy landscape even if the landscape is relatively rugged with metastable conformations. A closer examination on molecular forces shows that the formation of associated protein complexes can be primarily electrostatics-driven, hydrophobics-driven, or a combination of both in stabilizing specific binding interfaces. Taken together, these results suggest that the coarse-grained simulations and analyses provide an alternative toolset to study protein-protein association occurring in functional biomolecular complexes.

In protein self-assembly, types of surfaces determine the force between them. Yet the extent to which the surrounding water contributes to this force remains as a fundamental question. Here we study three self-assembling filament systems that respectively have hydrated (collagen), dry nonpolar, and dry polar (amyloid) interfaces. Using molecular dynamics simulations, we calculate and compare local hydration maps and hydration forces. We find that the primary hydration shells are formed all over the surface, regardless of the types of the underlying amino acids. The weakly oscillating hydration force arises from coalescence and depletion of hydration shells as two filaments approach, whereas local water diffusion, orientation, or hydrogen-bonding events have no direct effect. Hydration forces between hydrated, polar, and nonpolar interfaces differ in the amplitude and phase of the oscillation relative to the equilibrium surface separation. Therefore, water-mediated interactions between these protein surfaces, ranging in character from "hydrophobic" to "hydrophilic", have a common molecular origin based on the robustly formed hydration shells, which is likely applicable to a broad range of biomolecular assemblies whose interfacial geometry is similar in length scale to those of the present study.

Conformational stability of the collagen triple helix affects its turnover and determines tissue homeostasis. Although it is known that the presence of imino acids (prolines or hydroxyprolines) confer stability to the molecule, little is known regarding the stability of the imino-poor region lacking imino acids, which plays a key role in collagen cleavage. In particular, there have been continuing debates about the role of water in collagen stability. We addressed these issues using molecular dynamics simulations on 30-residue long collagen triple helices, including a structure that has a biologically relevant 9-residue imino-poor region from type III collagen (PDB ID: 1BKV). A torsional map approach was used to characterize the conformational motion of the molecule that differ between imino-rich and imino-poor regions. At temperatures 300 K and above, unwinding initiates at a common cleavage site, the glycine-isoleucine bond in the imino-poor region. This provides a linkage between previous observations that unwinding of the imino-poor region is a requirement for collagenase cleavage, and that isolated collagen molecules are unstable at body temperature. We found that unwinding of the imino-poor region is controlled by dynamic water bridges between backbone atoms with average lifetimes on the order of a few picoseconds, as the degree of unwinding strongly correlated with the loss of water bridges, and unwinding could be either prevented or enhanced, respectively by enforcing or forbidding water bridge formation. While individual water bridges were short-lived in the imino-poor region, the hydration shell surrounding the entire molecule was stable even at 330 K. The diameter of the hydrated collagen including the first hydration shell was about 14 A, in good agreement with the experimentally measured inter-collagen distances. These results elucidate the general role of water in collagen turnover: water not only affects collagen cleavage by controlling its torsional motion, but it also forms a larger-scale lubrication layer mediating collagen self-assembly.

We examine the energetics of β-strand packing in a fibril-forming SH3 domain using a simple sequence-based energy model. First, we describe this packing energy function and then apply it to three model systems: Aβ, HET-s prion, and SH3 domain. The packing results of Aβ and HET-s are compared to and are consistent with available experimental and computational results. Moreover, our results show that a native β-strand in SH3 is strongly disfavored to pack with any other strand, in accord with recent NMR data. Finally, based on packing energy calculations, several SH3 models of β-strand packing are proposed that fit well with known electron microscopy maps.

Estrogen receptor alpha (ERα) is a hormone-responsive transcription factor that contains several discrete functional domains, including a ligand-binding domain (LBD) and a DNA-binding domain (DBD). Despite a wealth of knowledge about the behaviors of individual domains, the molecular mechanisms of cross-talk between LBD and DBD during signal transduction from hormone to DNA-binding of ERα remain elusive. Here, we apply a multiscale approach combining coarse-grained (CG) and atomistically detailed simulations to characterize this cross-talk mechanism via an investigation of the ERα conformational landscape. First, a CG model of ERα is built based on crystal structures of individual LBDs and DBDs, with more emphasis on their interdomain interactions. Second, molecular dynamics simulations are implemented and enhanced sampling is achieved via the "push-pull-release" strategy in the search for different LBD-DBD orientations. Third, multiple energetically stable ERα conformations are identified on the landscape. A key finding is that estradiol-bound LBDs utilize the well-described activation helix H12 to pack and stabilize LBD-DBD interactions. Our results suggest that the estradiol-bound LBDs can serve as a scaffold to position and stabilize the DBD-DNA complex, consistent with experimental observations of enhanced DNA binding with the LBD. Final assessment using atomic-level simulations shows that these CG-predicted models are significantly stable within a 15-ns simulation window and that specific pairs of lysine residues in close proximity at the domain interfaces could serve as candidate sites for chemical cross-linking studies. Together, these simulation results provide a molecular view of the role of ERα domain interactions in response to hormone binding.

Molecular interactions of optical clearing agents were investigated using a combination of molecular dynamics (MD) simulations and optical spectroscopy. For a series of sugar alcohols with low to high optical clearing potential, Raman spectroscopy and integrating sphere measurements were used to quantitatively characterize tissue water loss and reduction in light scattering following agent exposures. The rate of tissue water loss was found to correlate with agent optical clearing potential, but equivalent tissue optical clearing was measured in native and fixed tissue in vitro, given long-enough exposure times to the polyol series. MD simulations showed that the rate of tissue optical clearing correlated with the preferential formation of hydrogen bond bridges between agent and collagen. Hydrogen bond bridge formation disrupts the collagen hydration layer and facilitates replacement by a chemical agent to homogenize tissue refractive index. However, the reduction in tissue light scattering did not correlate with the agent index of refraction. Our results suggest that a necessary property of optical clearing agents is hyperosmolarity to tissue, but that the most effective agents with the highest rates of optical clearing are a subset with the highest collagen solubilities.

We investigate the hydration structure of collagen using molecular dynamics simulations. The water density maximum, and orientation near the surface of collagen and their specificity to the proximal collagen sequence are clearly visible in our analysis. These results provide insights into the effect of hydration in collagen self-assembly, and more generally, the anomalous behavior of water near protein surfaces.Copyright © 2010 by ASME

A generalized method, termed Fast-SAXS-pro, for computing small angle x-ray scattering (SAXS) profiles of proteins, nucleic acids, and their complexes is presented. First, effective coarse-grained structure factors of DNA nucleotides are derived using a simplified two-particle-per-nucleotide representation. Second, SAXS data of a 18-bp double-stranded DNA are measured and used for the calibration of the scattering contribution from excess electron density in the DNA solvation layer. Additional test on a 25-bp DNA duplex validates this SAXS computational method and suggests that DNA has a different contribution from its hydration surface to the total scattering compared to RNA and protein. To account for such a difference, a sigmoidal function is implemented for the treatment of non-uniform electron density across the surface of a protein/nucleic-acid complex. This treatment allows differential scattering from the solvation layer surrounding protein/nucleic-acid complexes. Finally, the applications of this Fast-SAXS-pro method are demonstrated for protein/DNA and protein/RNA complexes.

Collagens are triple-helical molecules that self-assemble into higher order fibers forming the major component of extracellular matrix. We had previously reported the role of the first hydration layer in controlling the conformational behavior of the collagen triple helix. Here we perform explicit-water molecular dynamics simulations to elucidate the structural features of water in mediating the interaction between collagen triple helices (PDB ID: 1A3I, 2D3F). By dividing the simulation box into cells, we quantified local water density, diffusion coefficient, and water orientation at atomistic resolution. Around a single collagen triple helix the reduction in diffusion coefficient and density fluctuation extend up to 11 Angstroms from the collagen backbone, and the circumferential and radial orientation of water near hydrophobic and hydrophilic groups, respectively, were clearly distinguishable. When two triple helices were held at a radial separation that is a few Angstroms larger than their crystalline packing distance, water in between them had reduced diffusion coefficient and constrained angular orientation. This indicates that the experimentally observed attractive force between collagens at small distances may have an entropic origin. We also tested three-collagen systems where one collagen is radially translated from its original position in the crystal packing by 4 or 7 Angstroms, and found that it moves towards the other two within 3-ns of simulation, nearly restoring the crystal packing. These results illustrate the microscopic origins of water mediated attraction between collagen molecules.