Decline in the Age of Menarche in Istanbul Schoolgirls Over the Last 12 Years

Objective: Menarche is the endpoint of a sequence of maturational events of female puberty. The timing of menarche is a strongly heritable trait. However, secular trends suggest that lifestyle and environmental factors are important. To assess the trend in age at menarche (AAM), and its associated factors in Istanbul over the last 12 years. Methods: A cross-sectional study was carried out between March and April 2022 on schoolgirls aged 9-18 years. A predesigned and self-administered questionnaire was filled out anonymously by the students. The data of AAM was included in the statistical analysis if the time of AAM is remembered in both months and years. A probit model was used to calculate the median AAM. The findings were compared with those from a study performed 12 years ago in the same region of Istanbul. Results: Among 9000 girls to whom the questionnaire was distributed, 1749 (19.5%) responded. The median AAM of 1374 girls whose AAM information was considered valid was 12.04 years (95% confidence interval: 12.01-12.13), 0.7 years lower than was reported 12 years ago (p<0.0001). AAM was correlated positively with maternal AAM, and negatively with body mass index (BMI) standard deviation score and maternal educational status (p<0.0001, p<0.0001 and p=0.002), respectively. There was no correlation between the AAM and birth weight. Girls with BMI percentile ≥85% (n=251) had earlier menarche than the ones with BMI percentile <85% (n=1072) (11.5 vs. 12.1 years, p<0.0001). Among the mother-daughter pairs (n=1162), AAM of girls was 0.91 years (median 0.94 years) earlier than their mothers. Conclusion: The present study demonstrates a significant downward trend in the menarcheal age in Istanbul over the last twelve years. These findings support a strong contribution from genetic factors and BMI on AAM.

Modeling Adsorption, Conformation, and Orientation of the Fis1 Tail Anchor at the Mitochondrial Outer Membrane.

Proteins can be targeted to organellar membranes by using a tail anchor (TA), a stretch of hydrophobic amino acids found at the polypeptide carboxyl-terminus. The Fis1 protein (Fis1p), which promotes mitochondrial and peroxisomal division in the yeast Saccharomyces cerevisiae, is targeted to those organelles by its TA. Substantial evidence suggests that Fis1p insertion into the mitochondrial outer membrane can occur without the need for a translocation machinery. However, recent findings raise the possibility that Fis1p insertion into mitochondria might be promoted by a proteinaceous complex. Here, we have performed atomistic and coarse-grained molecular dynamics simulations to analyze the adsorption, conformation, and orientation of the Fis1(TA). Our results support stable insertion at the mitochondrial outer membrane in a monotopic, rather than a bitopic (transmembrane), configuration. Once inserted in the monotopic orientation, unassisted transition to the bitopic orientation is expected to be blocked by the highly charged nature of the TA carboxyl-terminus and by the Fis1p cytosolic domain. Our results are consistent with a model in which Fis1p does not require a translocation machinery for insertion at mitochondria.

Representation of the conformational ensemble of peptides in coarse grained simulations.

In their native state, many proteins/peptides display an ensemble of conformations, rather than a unique tertiary structure. Novel experimental techniques have enabled a quantitative analysis of this structural heterogeneity. In molecular dynamics simulations, however, capturing this conformational ensemble quantitatively remains a major challenge even with all atom simulations. In coarse grained (CG) simulations, with fewer degrees of freedom, representation of the conformational ensemble becomes more problematic. Here, we revisit a CG model from our group, which was designed to address the conformational transferability problem by using the LKα14 peptide as a model system. The LKα14 peptide transitions from a random/unstructured state in dilute solution to a solely α-helical conformation upon aggregation as evidenced by circular dichroism. Here, we demonstrate that the structure/physics based approach, used in the original parameterization of our CG model, strongly depends on the reference system chosen and excluded volume interactions that are often considered to be of secondary importance. We first tune the excluded volume parameters by using both α-helix and ß-sheet type structures as reference and then update the nonbonded interactions by using a goodness-of-fit metric for representation of the conformational ensemble of LKα14. We demonstrate that the updated model can recover the whole conformational ensemble quantitatively while maintaining the aggregation driven conformational transition. This balanced parametrization with regard to alternative secondary structures opens the door for the generalization of the CG model to other sequences, which we demonstrate on a ß-sheet forming triblock peptide.

A multi-state coarse grained modeling approach for an intrinsically disordered peptide.

Many proteins display a marginally stable tertiary structure, which can be altered via external stimuli. Since a majority of coarse grained (CG) models are aimed at structure prediction, their success for an intrinsically disordered peptide's conformational space with marginal stability and sensitivity to external stimuli cannot be taken for granted. In this study, by using the LKα14 peptide as a test system, we demonstrate a bottom-up approach for constructing a multi-state CG model, which can capture the conformational behavior of this peptide in three distinct environments with a unique set of interaction parameters. LKα14 is disordered in dilute solutions; however, it strictly adopts the α-helix conformation upon aggregation or when in contact with a hydrophobic/hydrophilic interface. Our bottom-up approach combines a generic base model, that is unbiased for any particular secondary structure, with nonbonded interactions which represent hydrogen bonds, electrostatics, and hydrophobic forces. We demonstrate that by using carefully designed all atom potential of mean force calculations from all three states of interest, one can get a balanced representation of the nonbonded interactions. Our CG model behaves intrinsically disordered in bulk water, folds into an α-helix in the presence of an interface or a neighboring peptide, and is stable as a tetrameric unit, successfully reproducing the all atom molecular dynamics simulations and experimental results.

Assembly of Huntingtin headpiece into α-helical bundles.

Protein aggregation is a hallmark of neurodegenerative disorders. In this group of brain-related disorders, a disease-specific "host" protein or fragment misfolds and adopts a metastatic, aggregate-prone conformation. Often, this misfolded conformation is structurally and thermodynamically different from its native state. Intermolecular contacts, which arise in this non-native state, promote aggregation. In this regard, understanding the molecular principles and mechanisms that lead to the formation of such a non-native state and further promote the formation of the critical nucleus for fiber growth is essential. In this study, the authors analyze the aggregation propensity of Huntingtin headpiece (httNT), which is known to facilitate the polyQ aggregation, in relation to the helix mediated aggregation mechanism proposed by the Wetzel group. The authors demonstrate that even though httNT displays a degenerate conformational spectrum on its own, interfaces of macroscopic or molecular origin can promote the α-helix conformation, eliminating all other alternatives in the conformational phase space. Our findings indicate that httNT molecules do not have a strong orientational preference for parallel or antiparallel orientation of the helices within the aggregate. However, a parallel packed bundle of helices would support the idea of increased polyglutamine concentration, to pave the way for cross-ß structures.

Role of Hydrophobic/Aromatic Residues on the Stability of Double-Wall ß-Sheet Structures Formed by a Triblock Peptide.

Bioinspired self-assembling peptides serve as powerful building blocks in the manufacturing of nanomaterials with tailored features. Because of their ease of synthesis, biocompatibility, and tunable activity, this emerging branch of biomolecules has become very popular. The triblock peptide architecture designed by the Hartgerink group is a versatile system that allows control over its assembly and has been shown to demonstrate tunable bioactivity. Three main forces, Coulomb repulsion, hydrogen bonding and hydrophobicity act together to guide the triblock peptides' assembly into one-dimensional objects and hydrogels. It was shown previously that both the nanofiber morphology (e.g., intersheet spacing, formation of antiparallel/parallel ß-sheets) and hydrogel rheology strictly depend on the choice of the core residue where the triblock peptide fibers with aromatic cores in general form shorter fibers and yield poor hydrogels with respect to the ones with aliphatic cores. However, an elaborate understanding of the molecular reasons behind these changes remained unclear. In this study, by using carefully designed computer based free energy calculations, we analyzed the influence of the core residue on the formation of double-wall fibers and single-wall ß-sheets. Our results demonstrate that the aromatic substitution impairs the fiber cores and this impairment is mainly associated with a reduced hydrophobic character of the aromatic side chains. Such weakening is most obvious in tryptophan containing peptides where the fiber core absorbs a significant amount of water. We also show that the ability of tyrosine to form side chain hydrogen bonds plays an indispensable role in the fiber stability. As opposed to the impairment of the fiber cores, single-wall ß-sheets with aromatic faces become more stable compared to the ones with aliphatic faces suggesting that the choice of the core residue can also affect the underlying assembly mechanism. We also provide an in-depth comparison of competing structures (zero-dimensional aggregates, short and long fibers) in the triblock peptides' assembly and show that by adjusting the length of the terminal blocks, the fiber growth can be turned on or off while keeping the nanofiber morphology intact.

Assembly of Triblock Amphiphilic Peptides into One-Dimensional Aggregates and Network Formation.

Peptide assembly plays a key role in both neurological diseases and development of novel biomaterials with well-defined nanostructures. Synthetic model peptides provide a unique platform to explore the role of intermolecular interactions in the assembly process. A triblock peptide architecture designed by the Hartgerink group is a versatile system which relies on Coulomb interactions, hydrogen bonding, and hydrophobicity to guide these peptides' assembly at three different length scales: ß-sheets, double-wall ribbon-like aggregates, and finally a highly porous network structure which can support gels with ≤1% by weight peptide concentration. In this study, by using molecular dynamics simulations of a structure based implicit solvent coarse grained model, we analyzed this hierarchical assembly process. Parametrization of our CG model is based on multiple-state points from atomistic simulations, which enables this model to represent the conformational adaptability of the triblock peptide molecule based on the surrounding medium. Our results indicate that emergence of the double-wall ß-sheet packing mechanism, proposed in light of the experimental evidence, strongly depends on the subtle balance of the intermolecular forces. We demonstrate that, even though backbone hydrogen bonding dominates the early nucleation stages, depending on the strength of the hydrophobic and Coulomb forces, alternative structures such as zero-dimensional aggregates with two ß-sheets oriented orthogonally (which we refer to as a cross-packed structure) and ß-sheets with misoriented hydrophobic side chains are also feasible. We discuss the implications of these competing structures for the three different length scales of assembly by systematically investigating the influence of density, counterion valency, and hydrophobicity.

Conformation and Aggregation of LKα14 Peptide in Bulk Water and at the Air/Water Interface.

Historically, the protein folding problem has mainly been associated with understanding the relationship between amino acid sequence and structure. However, it is known that both the conformation of individual molecules and their aggregation strongly depend on the environmental conditions. Here, we study the aggregation behavior of the model peptide LKα14 (with amino acid sequence LKKLLKLLKKLLKL) in bulk water and at the air/water interface. We start by a quantitative analysis of the conformational space of a single LKα14 in bulk water. Next, in order to analyze the aggregation tendency of LKα14, by using the umbrella sampling technique we calculate the potential of mean force for pulling a single peptide from an n-molecule aggregate. In agreement with the experimental results, our calculations yield the optimal aggregate size as four. This equilibrium state is achieved by two opposing forces: Coulomb repulsion between the lysine side chains and the reduction of solvent accessible hydrophobic surface area upon aggregation. At the vacuum/water interface, however, even dimers of LKα14 become marginally stable, and any larger aggregate falls apart instantaneously. Our results indicate that even though the interface is highly influential in stabilizing the α-helix conformation for a single molecule, it significantly reduces the attraction between two LKα14 peptides, along with their aggregation tendency.

Tipping the Scale from Disorder to Alpha-helix: Folding of Amphiphilic Peptides in the Presence of Macroscopic and Molecular Interfaces.

Secondary amphiphilicity is inherent to the secondary structural elements of proteins. By forming energetically favorable contacts with each other these amphiphilic building blocks give rise to the formation of a tertiary structure. Small proteins and peptides, on the other hand, are usually too short to form multiple structural elements and cannot stabilize them internally. Therefore, these molecules are often found to be structurally ambiguous up to the point of a large degree of intrinsic disorder in solution. Consequently, their conformational preference is particularly susceptible to environmental conditions such as pH, salts, or presence of interfaces. In this study we use molecular dynamics simulations to analyze the conformational behavior of two synthetic peptides, LKKLLKLLKKLLKL (LK) and EAALAEALAEALAE (EALA), with built-in secondary amphiphilicity upon forming an alpha-helix. We use these model peptides to systematically study their aggregation and the influence of macroscopic and molecular interfaces on their conformational preferences. We show that the peptides are neither random coils in bulk water nor fully formed alpha helices, but adopt multiple conformations and secondary structure elements with short lifetimes. These provide a basis for conformation-selection and population-shift upon environmental changes. Differences in these peptides' response to macroscopic and molecular interfaces (presented by an aggregation partner) can be linked to their inherent alpha-helical tendencies in bulk water. We find that the peptides' aggregation behavior is also strongly affected by presence or absence of an interface, and rather subtly depends on their surface charge and hydrophobicity.

A transferable coarse-grained model for diphenylalanine: how to represent an environment driven conformational transition.

One of the major challenges in the development of coarse grained (CG) simulation models that aim at biomolecular structure formation processes is the correct representation of an environment-driven conformational change, for example, a folding/unfolding event upon interaction with an interface or upon aggregation. In the present study, we investigate this transferability challenge for a CG model using the example of diphenylalanine. This dipeptide displays a transition from a trans-like to a cis-like conformation upon aggregation as well as upon transfer from bulk water to the cyclohexane/water interface. Here, we show that one can construct a single CG model that can reproduce both the bulk and interface conformational behavior and the segregation between hydrophobic/hydrophilic medium. While the general strategy to obtain nonbonded interactions in the present CG model is to reproduce solvation free energies of small molecules representing the CG beads in the respective solvents, the success of the model strongly depends on nontrivial decisions one has to make to capture the delicate balance between the bonded and nonbonded interactions. In particular, we found that the peptide's conformational behavior is qualitatively affected by the cyclohexane/water interaction potential, an interaction that does not directly involve the peptide at all but merely influences the properties of the hydrophobic/hydrophilic interface. Furthermore, we show that a small modification to improve the structural/conformational properties of the CG model could dramatically alter the thermodynamic properties.

Self-assembling multidomain peptide fibers with aromatic cores.

Self-assembling multidomain peptides have been shown to have desirable properties, such as the ability to form hydrogels that rapidly recover following shear-thinning and the potential to be tailored by amino acid selection to vary their elasticity and encapsulate and deliver proteins and cells. Here we describe the effects of substitution of aliphatic hydrophobic amino acids in the central domain of the peptide for the aromatic amino acids phenylalanine, tyrosine, and tryptophan. While the basic nanofibrous morphology is retained in all cases, selection of the particular core residues results in switching from antiparallel hydrogen bonding to parallel hydrogen bonding in addition to changes in nanofiber morphology and in hydrogel rheological properties. Peptide nanofiber assemblies are investigated by circular dichroism polarimetry, infrared spectroscopy, atomic force microscopy, transmission and scanning electron microscopy, oscillatory rheology, and molecular dynamics simulations. Results from this study will aid in designing next generation cell scaffolding materials.

Adsorption, folding, and packing of an amphiphilic peptide at the air/water interface.

Peptide oligomers play an essential role as model compounds for identifying key motifs in protein structure formation and protein aggregation. Here, we present our results, based on extensive molecular dynamics simulations, on adsorption, folding, and packing within a surface monolayer of an amphiphilic peptide at the air/water interface. Experimental results suggest that these molecules spontaneously form ordered monolayers at the interface, adopting a ß-hairpin-like structure within the surface layer. Our results reveal that the ß-hairpin structure can be observed both in bulk and at the air/water interface. However, the presence of an interface leads to ideal partitioning of the hydrophobic and hydrophilic residues, and therefore reduces the conformational space for the molecule and increases the stability of the hairpin structure. We obtained the adsorption free energy of a single ß-hairpin at the air/water interface, and analyzed the enthalpic and entropic contributions. The adsorption process is favored by two main factors: (1) Free-energy reduction due to desolvation of the hydrophobic side chains of the peptide and release of the water molecules which form a cage around these hydrophobic groups in bulk water. (2) Reduction of the total air/water contact area at the interface upon adsorption of the peptide amphiphile. By performing mutations on the original molecule, we demonstrated the relative role of key design features of the peptide. Finally, by analyzing the potential of mean force among two peptides at the interface, we investigated possible packing mechanisms for these molecules within the surface monolayer.

Driving forces for adsorption of amphiphilic peptides to the air-water interface.

We have studied the partitioning of amphiphilic peptides at the air-water interface. The free energy of adsorption from bulk to interface was calculated by determining the potential of mean force via atomistic molecular dynamics simulations. To this end a method is introduced to restrain or constrain the center of mass of a group of molecules in a periodic system. The model amphiphilic peptides are composed of alternating valine and asparagine residues. The decomposition of the free energy difference between the bulk and interface is studied for different peptide block lengths. Our analysis revealed that for short amphiphilic peptides the surface driving force dominantly stems from the dehydration of hydrophobic side chains. The only opposing force is associated with the loss of orientational freedom of the peptide at the interface. For the peptides studied, the free energy difference scales linearly with the size of the molecule, since the peptides mainly adopt extended conformations both in bulk and at the interface. The free energy difference depends strongly on the water model, which can be rationalized through the hydration thermodynamics of hydrophobic solutes. Finally, we measured the reduction of the surface tension associated with complete coverage of the interface with peptides.

Twist-writhe partitioning in a coarse-grained DNA minicircle model.

Here we present a systematic study of supercoil formation in DNA minicircles under varying linking number by using molecular-dynamics simulations of a two-bead coarse-grained model. Our model is designed with the purpose of simulating long chains without sacrificing the characteristic structural properties of the DNA molecule, such as its helicity, backbone directionality, and the presence of major and minor grooves. The model parameters are extracted directly from full-atomistic simulations of DNA oligomers via Boltzmann inversion; therefore, our results can be interpreted as an extrapolation of those simulations to presently inaccessible chain lengths and simulation times. Using this model, we measure the twist/writhe partitioning in DNA minicircles, in particular its dependence on the chain length and excess linking number. We observe an asymmetric supercoiling transition consistent with experiments. Our results suggest that the fraction of the linking number absorbed as twist and writhe is nontrivially dependent on chain length and excess linking number. Beyond the supercoiling transition, chains of the order of one persistence length carry equal amounts of twist and writhe. For longer chains, an increasing fraction of the linking number is absorbed by the writhe.

Equilibrium polyelectrolyte bundles with different multivalent counterion concentrations.

We present the results of molecular-dynamics simulations on the salt concentration dependence of the formation of polyelectrolyte bundles in thermodynamic equilibrium. Extending our results on salt-free systems we investigate here deficiency or excess of trivalent counterions in solution. Our results reveal that the trivalent counterion concentration significantly alters the bundle size and size distribution. The onset of bundle formation takes place at earlier Bjerrum length values with increasing trivalent counterion concentration. For the cases of 80%, 95%, and 100% charge compensation via trivalent counterions, the net charge of the bundles decreases with increasing size. We suggest that competition among two different mechanisms, counterion condensation and merger of bundles, leads to a nonmonotonic change in line-charge density with increasing Bjerrum length. The investigated case of having an abundance of trivalent counterions by 200% prohibits such a behavior. In this case, we find that the difference in effective line-charge density of different size bundles diminishes. In fact, the system displays an isoelectric point, where all bundles become charge neutral.

The introduction of hydrogen bond and hydrophobicity effects into the rotational isomeric states model for conformational analysis of unfolded peptides.

Relative contributions of local and non-local interactions to the unfolded conformations of peptides are examined by using the rotational isomeric states model which is a Markov model based on pairwise interactions of torsion angles. The isomeric states of a residue are well described by the Ramachandran map of backbone torsion angles. The statistical weight matrices for the states are determined by molecular dynamics simulations applied to monopeptides and dipeptides. Conformational properties of tripeptides formed from combinations of alanine, valine, tyrosine and tryptophan are investigated based on the Markov model. Comparison with molecular dynamics simulation results on these tripeptides identifies the sequence-distant long-range interactions that are missing in the Markov model. These are essentially the hydrogen bond and hydrophobic interactions that are obtained between the first and the third residue of a tripeptide. A systematic correction is proposed for incorporating these long-range interactions into the rotational isomeric states model. Preliminary results suggest that the Markov assumption can be improved significantly by renormalizing the statistical weight matrices to include the effects of the long-range correlations.

Buckling of stiff polymers: influence of thermal fluctuations.

The buckling of biopolymers is a frequently studied phenomenon The influence of thermal fluctuations on the buckling transition is, however, often ignored and not completely understood. A quantitative theory of the buckling of a wormlike chain based on a semiclassical approximation of the partition function is presented. The contribution of thermal fluctuations to the force-extension relation that allows one to go beyond the classical Euler buckling is derived in the linear and nonlinear regimes as well. It is shown that the thermal fluctuations in the nonlinear buckling regime increase the end-to-end distance of the semiflexible rod if it is confined to two dimensions as opposed to the three-dimensional case. The transition to a buckled state softens at finite temperature. We derive the scaling behavior of the transition shift with increasing ratio of contour length versus persistence length.

Assembly of one-dimensional supramolecular objects: from monomers to networks.

One-dimensional supramolecular aggregates can form networks at exceedingly low concentrations. Recent experiments in several laboratories, including our own, have demonstrated the formation of gels by these systems at concentrations well under 1% by weight. The systems of interest in our laboratory form either cylindrical nanofibers or ribbons as a result of strong noncovalent interactions among monomers. The stiffness and interaction energies among these thread-like objects can vary significantly depending on the chemical structure of the monomers used. We have used Monte Carlo simulations to study the structure of the threads and their ability to form networks through bundle formation. The persistence length of the threads was found to be strongly affected not only by stiffness, but also by the strength of attractive two-body interactions among thread segments. The relative values of stiffness and attractive two-body interaction strength determine if threads collapse or create bundles. Only in the presence of sufficiently long threads and bundle formation can these systems assemble into networks of high connectivity.